SlowCRTest.java revision e96bec4f09552f6f85fc6d763a7fec938a008fea
1a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm/* 2a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * Copyright (C) 2015 The Android Open Source Project 3a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * 4a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * Licensed under the Apache License, Version 2.0 (the "License"); 5a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * you may not use this file except in compliance with the License. 6a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * You may obtain a copy of the License at 7a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * 8a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * http://www.apache.org/licenses/LICENSE-2.0 9a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * 10a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * Unless required by applicable law or agreed to in writing, software 11a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * distributed under the License is distributed on an "AS IS" BASIS, 12a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 13a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * See the License for the specific language governing permissions and 14a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm * limitations under the License. 15a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm */ 16a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 17a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm/* THIS TYPICALLY TAKES > 10 MINUTES TO RUN! It shold generate no output during that time. */ 18a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 19a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmpackage com.hp.creals; 20a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 21a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport junit.framework.AssertionFailedError; 22a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport junit.framework.TestCase; 23a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 24a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport java.math.BigInteger; 25a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport java.util.Random; 26a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 27a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmpublic class SlowCRTest extends TestCase { 28a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void check(boolean x, String s) { 29a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (!x) throw new AssertionFailedError(s); 30a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 31a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static int TEST_PREC = -200; // 200 bits to the right of 32a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // binary point. 33a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static int NRANDOM = 100; // Number of random values to 34a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // test. Bigger ==> slower 35a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkEq(CR x, CR y, String s) { 36e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm check(x.compareTo(y, TEST_PREC) == 0, s); 37a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 38a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkApprEq(double x, double y, String s) { 39e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm check(Math.abs(x - y) <= 0.000001, s); 40a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 41a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static BigInteger MASK = 42a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm BigInteger.ONE.shiftLeft(-TEST_PREC).subtract(BigInteger.ONE); 43a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static boolean isApprInt(CR x) { 44a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm BigInteger appr = x.get_appr(TEST_PREC); 45a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm return appr.and(MASK).signum() == 0; 46a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 47a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 48a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR ZERO = CR.valueOf(0); 49a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR ONE = CR.valueOf(1); 50a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR TWO = CR.valueOf(2); 51a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR BIG = CR.valueOf(200).exp(); 52a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR SMALL = BIG.inverse(); 53a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 54a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static UnaryCRFunction ASIN = UnaryCRFunction.asinFunction; 55a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static UnaryCRFunction ACOS = UnaryCRFunction.acosFunction; 56a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static UnaryCRFunction ATAN = UnaryCRFunction.atanFunction; 57a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static UnaryCRFunction TAN = UnaryCRFunction.tanFunction; 58a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static UnaryCRFunction COSINE = UnaryCRFunction.sinFunction 59a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm .monotoneDerivative(ZERO, CR.PI); 60a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 61a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Perform some consistency checks on trig functions at x 62a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // We assume that x is within floating point range. 63a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkTrig(CR x) { 64a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double xAsDouble = x.doubleValue(); 65a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) < 1000000.0) { 66a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.sin().doubleValue(), Math.sin(xAsDouble), 67a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "sin float compare:" + xAsDouble); 68a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.cos().doubleValue(), Math.cos(xAsDouble), 69a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "cos float compare:" + xAsDouble); 70a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(TAN.execute(x).doubleValue(), Math.tan(xAsDouble), 71a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "tan float compare:" + xAsDouble); 72a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(ATAN.execute(x).doubleValue(), Math.atan(xAsDouble), 73a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "atan float compare:" + xAsDouble); 74a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 75a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) < 1.0) { 76a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(ASIN.execute(x).doubleValue(), Math.asin(xAsDouble), 77a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "asin float compare:" + xAsDouble); 78a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(ACOS.execute(x).doubleValue(), Math.acos(xAsDouble), 79a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "acos float compare:" + xAsDouble ); 80a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 81a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (xAsDouble < 3.1415926535 && xAsDouble > 0.0) { 82a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(COSINE.execute(x).doubleValue(), Math.cos(xAsDouble), 83a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "deriv(sin) float compare:" + xAsDouble); 84e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(COSINE.execute(x), x.cos(), 85a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "deriv(sin) float compare:" + xAsDouble); 86a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 87a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Check that sin(x+v) = sin(x)cos(v) + cos(x)sin(v) 88a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // for a couple of different values of v. 89a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 1; i <= 5; ++i) { 90a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR v = CR.valueOf(i); 91e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq( 92a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.add(v).sin(), 93a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.sin().multiply(v.cos()).add(x.cos().multiply(v.sin())), 94a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "Angle sum formula failed for " + xAsDouble + " + " + i); 95a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 96e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.cos().multiply(x.cos()).add(x.sin().multiply(x.sin())), 97a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR.valueOf(1), 98a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "sin(x)^2 + cos(x)^2 != 1:" + xAsDouble); 99a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Check that inverses are consistent 100e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x, TAN.execute(ATAN.execute(x)), 101a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "tan(atan(" + xAsDouble + ")" ); 102a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR tmp = ACOS.execute(x.cos()); 103a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Result or its inverse should differ from x by an 104a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // exact multiple of pi. 105a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm check(isApprInt(tmp.subtract(x).divide(CR.PI)) 106a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm || isApprInt(tmp.add(x).divide(CR.PI)), 107a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "acos(cos):" + xAsDouble); 108a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm tmp = ASIN.execute(x.sin()); 109a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Result or its inverse should differ from x by an 110a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // exact multiple of pi. 111a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm check(isApprInt(tmp.subtract(x).divide(CR.PI)) 112a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm || isApprInt(tmp.add(x).divide(CR.PI)), 113a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "acos(cos):" + xAsDouble); 114a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 115a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 116a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkExpLn(CR x) { 117a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double xAsDouble = x.doubleValue(); 118a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) < 10.0) { 119a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.exp().doubleValue(), Math.exp(xAsDouble), 120a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "exp float compare:" + xAsDouble); 121a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 122a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) <= 1000.0) { 123e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x, x.exp().ln(), "ln(exp) failed:" + xAsDouble); 124e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.multiply(CR.valueOf(2)).exp(), 125a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 126a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.exp().multiply(x.exp()), 127a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "exp^2 failed:" + xAsDouble); 128a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 129a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (xAsDouble > 0.000000001) { 130a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.ln().doubleValue(), Math.log(xAsDouble), 131a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "exp float compare:" + xAsDouble); 132e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x, x.ln().exp(), "exp(ln) failed:" + xAsDouble); 133e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.ln().divide(CR.valueOf(2)), x.sqrt().ln(), 134a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "ln(sqrt) failed:" + xAsDouble); 135a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Check that ln(xv) = ln(x) + ln(v) for various v 136a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 1; i <= 5; ++i) { 137a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR v = CR.valueOf(i); 138e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq( 139a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.ln().add(v.ln()), 140a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.multiply(v).ln(), 141a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "ln(product) formula failed for:" + xAsDouble + "," + i); 142a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 143a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 144a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 145a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 146a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkBasic(CR x) { 147e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.abs().sqrt().multiply(x.abs().sqrt()), x.abs(), 148a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "sqrt*sqrt:" + x.doubleValue()); 149a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (!x.get_appr(TEST_PREC).equals(BigInteger.ZERO)) { 150e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.inverse().inverse(), x, 151a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "inverse(inverse):" + x.doubleValue()); 152a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 153a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 154a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 155a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm public void testSlowTrig() { 156e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ACOS.execute(ZERO), CR.PI.divide(TWO), "acos(0)"); 157e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ACOS.execute(ONE), ZERO, "acos(1)"); 158e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ACOS.execute(ONE.negate()), CR.PI, "acos(-1)"); 159e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ASIN.execute(ZERO), ZERO, "asin(0)"); 160e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ASIN.execute(ONE), CR.PI.divide(TWO), "asin(1)"); 161e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ASIN.execute(ONE.negate()), 162a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR.PI.divide(TWO).negate(), "asin(-1)"); 163a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(ZERO); 164a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR BIG = CR.valueOf(200).exp(); 165a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(BIG); 166a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(BIG.negate()); 167a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(SMALL); 168a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(SMALL.negate()); 169a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI); 170a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.subtract(SMALL)); 171a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.add(SMALL)); 172a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.negate()); 173a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.negate().subtract(SMALL)); 174a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.negate().add(SMALL)); 175a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm Random r = new Random(); // Random seed! 176a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < NRANDOM; ++i) { 177a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(2.0 * r.nextDouble() - 1.0); 178a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 179a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 180a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(x); 181a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 182a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // And a few big ones 183a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < 10; ++i) { 184a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(200.0 * r.nextDouble()); 185a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 186a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 187a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(x); 188a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 189a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 190a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 191a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm public void testSlowExpLn() { 192e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(CR.valueOf(1).ln(), CR.valueOf(0), "ln(1) != 0"); 193a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(CR.valueOf(0)); 194a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR BIG = CR.valueOf(200).exp(); 195a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(BIG); 196a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(BIG.negate()); 197a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(SMALL); 198a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(SMALL.negate()); 199a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(CR.PI); 200a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(ONE); 201a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(ONE.subtract(SMALL)); 202a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(ONE.negate().subtract(SMALL)); 203a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm Random r = new Random(); // Random seed! 204a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < NRANDOM; ++i) { 205a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(10.0 * r.nextDouble() - 1.0); 206a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 207a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 208a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(x); 209a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 210a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // And a few big ones 211a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < 10; ++i) { 212a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(200.0 * r.nextDouble()); 213a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 214a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 215a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(x); 216a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 217a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 218a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 219a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm public void testSlowBasic() { 220e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ZERO.sqrt(), ZERO, "sqrt(0)"); 221e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ZERO.abs(), ZERO, "abs(0)"); 222a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm Random r = new Random(); // Random seed! 223a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < NRANDOM; ++i) { 224a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(10.0 * r.nextDouble() - 1.0); 225a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 226a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 227a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkBasic(x); 228a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 229a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // And a few very big ones, but within IEEE double range 230a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < 10; ++i) { 231a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(600.0 * r.nextDouble()); 232a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 233a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 234a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkBasic(x); 235a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 236a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 237a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm} 238