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 17b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm// THIS TYPICALLY TAKES > 4 MINUTES TO RUN! 18b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm// It should generate no output during that time. 19a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 20a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmpackage com.hp.creals; 21a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 22a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport junit.framework.AssertionFailedError; 23a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport junit.framework.TestCase; 24a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 25a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport java.math.BigInteger; 26a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmimport java.util.Random; 27a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 28a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehmpublic class SlowCRTest extends TestCase { 29a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void check(boolean x, String s) { 30a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (!x) throw new AssertionFailedError(s); 31a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 32a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static int TEST_PREC = -200; // 200 bits to the right of 33a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // binary point. 34a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static int NRANDOM = 100; // Number of random values to 35a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // test. Bigger ==> slower 36a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkEq(CR x, CR y, String s) { 37e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm check(x.compareTo(y, TEST_PREC) == 0, s); 38a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 39a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkApprEq(double x, double y, String s) { 40e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm check(Math.abs(x - y) <= 0.000001, s); 41a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 42a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static BigInteger MASK = 43a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm BigInteger.ONE.shiftLeft(-TEST_PREC).subtract(BigInteger.ONE); 44a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static boolean isApprInt(CR x) { 45a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm BigInteger appr = x.get_appr(TEST_PREC); 46a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm return appr.and(MASK).signum() == 0; 47a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 48a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 49a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR ZERO = CR.valueOf(0); 50a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR ONE = CR.valueOf(1); 51a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR TWO = CR.valueOf(2); 52a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR BIG = CR.valueOf(200).exp(); 53a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final static CR SMALL = BIG.inverse(); 54b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm final static CR HALF_PI = CR.PI.divide(CR.valueOf(2)); 55a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 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); 60b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm final static UnaryCRFunction ARCSINE = 61b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm UnaryCRFunction.sinFunction.inverseMonotone(HALF_PI.negate(), 62b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm HALF_PI); 63a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 64a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Perform some consistency checks on trig functions at x 65a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // We assume that x is within floating point range. 66a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkTrig(CR x) { 67a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double xAsDouble = x.doubleValue(); 68a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) < 1000000.0) { 69a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.sin().doubleValue(), Math.sin(xAsDouble), 70a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "sin float compare:" + xAsDouble); 71a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.cos().doubleValue(), Math.cos(xAsDouble), 72a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "cos float compare:" + xAsDouble); 73a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(TAN.execute(x).doubleValue(), Math.tan(xAsDouble), 74a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "tan float compare:" + xAsDouble); 75a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(ATAN.execute(x).doubleValue(), Math.atan(xAsDouble), 76a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "atan float compare:" + xAsDouble); 77a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 78a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) < 1.0) { 79b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkApprEq(x.asin().doubleValue(), Math.asin(xAsDouble), 80a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "asin float compare:" + xAsDouble); 81b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkApprEq(x.acos().doubleValue(), Math.acos(xAsDouble), 82b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm "acos float compare:" + xAsDouble); 83b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ARCSINE.execute(x), x.asin(), 84b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm "inverse(sin) compare:" + xAsDouble); 85a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 86a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (xAsDouble < 3.1415926535 && xAsDouble > 0.0) { 87a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(COSINE.execute(x).doubleValue(), Math.cos(xAsDouble), 88a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "deriv(sin) float compare:" + xAsDouble); 89e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(COSINE.execute(x), x.cos(), 90a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "deriv(sin) float compare:" + xAsDouble); 91a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 92a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Check that sin(x+v) = sin(x)cos(v) + cos(x)sin(v) 93a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // for a couple of different values of v. 94a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 1; i <= 5; ++i) { 95a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR v = CR.valueOf(i); 96e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq( 97a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.add(v).sin(), 98a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.sin().multiply(v.cos()).add(x.cos().multiply(v.sin())), 99a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "Angle sum formula failed for " + xAsDouble + " + " + i); 100a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 101e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.cos().multiply(x.cos()).add(x.sin().multiply(x.sin())), 102a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR.valueOf(1), 103a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "sin(x)^2 + cos(x)^2 != 1:" + xAsDouble); 104a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Check that inverses are consistent 105e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x, TAN.execute(ATAN.execute(x)), 106b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm "tan(atan(" + xAsDouble + ")"); 107b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm CR xcos = x.cos(); 108b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm CR tmp = xcos.acos(); 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); 114b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm CR xsin = x.sin(); 115b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm tmp = ARCSINE.execute(xsin); 116b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm CR tmp2 = xsin.asin(); 117b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(tmp, tmp2, "Asin(sin) computations differ:" + xAsDouble); 118a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Result or its inverse should differ from x by an 119a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // exact multiple of pi. 120a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm check(isApprInt(tmp.subtract(x).divide(CR.PI)) 121a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm || isApprInt(tmp.add(x).divide(CR.PI)), 122a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "acos(cos):" + xAsDouble); 123a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 124a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 125a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkExpLn(CR x) { 126a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double xAsDouble = x.doubleValue(); 127a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) < 10.0) { 128a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.exp().doubleValue(), Math.exp(xAsDouble), 129a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "exp float compare:" + xAsDouble); 130a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 131a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (Math.abs(xAsDouble) <= 1000.0) { 132e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x, x.exp().ln(), "ln(exp) failed:" + xAsDouble); 133e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.multiply(CR.valueOf(2)).exp(), 134a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 135a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.exp().multiply(x.exp()), 136a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "exp^2 failed:" + xAsDouble); 137a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 138a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (xAsDouble > 0.000000001) { 139a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkApprEq(x.ln().doubleValue(), Math.log(xAsDouble), 140a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "exp float compare:" + xAsDouble); 141e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x, x.ln().exp(), "exp(ln) failed:" + xAsDouble); 142e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.ln().divide(CR.valueOf(2)), x.sqrt().ln(), 143a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "ln(sqrt) failed:" + xAsDouble); 144a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // Check that ln(xv) = ln(x) + ln(v) for various v 145a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 1; i <= 5; ++i) { 146a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR v = CR.valueOf(i); 147e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq( 148a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.ln().add(v.ln()), 149a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm x.multiply(v).ln(), 150a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "ln(product) formula failed for:" + xAsDouble + "," + i); 151a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 152a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 153a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 154a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 155a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm private static void checkBasic(CR x) { 156e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.abs().sqrt().multiply(x.abs().sqrt()), x.abs(), 157a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "sqrt*sqrt:" + x.doubleValue()); 158a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (!x.get_appr(TEST_PREC).equals(BigInteger.ZERO)) { 159e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(x.inverse().inverse(), x, 160a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm "inverse(inverse):" + x.doubleValue()); 161a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 162a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 163a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 164a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm public void testSlowTrig() { 165b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ZERO.acos(), CR.PI.divide(TWO), "acos(0)"); 166b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ONE.acos(), ZERO, "acos(1)"); 167b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ONE.negate().acos(), CR.PI, "acos(-1)"); 168b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ZERO.asin(), ZERO, "asin(0)"); 169b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ONE.asin(), CR.PI.divide(TWO), "asin(1)"); 170b849a8eb275670d2bd9b2dc535231cf333414f5eHans Boehm checkEq(ONE.negate().asin(), CR.PI.divide(TWO).negate(), "asin(-1)"); 171a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(ZERO); 172a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR BIG = CR.valueOf(200).exp(); 173a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(BIG); 174a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(BIG.negate()); 175a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(SMALL); 176a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(SMALL.negate()); 177a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI); 178a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.subtract(SMALL)); 179a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.add(SMALL)); 180a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.negate()); 181a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.negate().subtract(SMALL)); 182a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(CR.PI.negate().add(SMALL)); 183a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm Random r = new Random(); // Random seed! 184a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < NRANDOM; ++i) { 185a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(2.0 * r.nextDouble() - 1.0); 186a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 187a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 188a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(x); 189a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 190a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // And a few big ones 191a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < 10; ++i) { 192a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(200.0 * r.nextDouble()); 193a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 194a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 195a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkTrig(x); 196a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 197a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 198a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 199a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm public void testSlowExpLn() { 200e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(CR.valueOf(1).ln(), CR.valueOf(0), "ln(1) != 0"); 201a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(CR.valueOf(0)); 202a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm CR BIG = CR.valueOf(200).exp(); 203a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(BIG); 204a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(BIG.negate()); 205a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(SMALL); 206a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(SMALL.negate()); 207a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(CR.PI); 208a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(ONE); 209a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(ONE.subtract(SMALL)); 210a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(ONE.negate().subtract(SMALL)); 211a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm Random r = new Random(); // Random seed! 212a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < NRANDOM; ++i) { 213a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(10.0 * r.nextDouble() - 1.0); 214a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 215a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 216a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(x); 217a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 218a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // And a few big ones 219a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < 10; ++i) { 220a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(200.0 * r.nextDouble()); 221a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 222a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 223a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkExpLn(x); 224a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 225a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 226a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm 227a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm public void testSlowBasic() { 228e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ZERO.sqrt(), ZERO, "sqrt(0)"); 229e96bec4f09552f6f85fc6d763a7fec938a008feaHans Boehm checkEq(ZERO.abs(), ZERO, "abs(0)"); 230a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm Random r = new Random(); // Random seed! 231a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < NRANDOM; ++i) { 232a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(10.0 * r.nextDouble() - 1.0); 233a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 234a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 235a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkBasic(x); 236a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 237a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm // And a few very big ones, but within IEEE double range 238a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm for (int i = 0; i < 10; ++i) { 239a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm double d = Math.exp(600.0 * r.nextDouble()); 240a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm if (r.nextBoolean()) d = -d; 241a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm final CR x = CR.valueOf(d); 242a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm checkBasic(x); 243a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 244a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm } 245a8e45a3255089266ccf6500188a2d31af8634e4aHans Boehm} 246