// Copyright (c) 1999, Silicon Graphics, Inc. -- ALL RIGHTS RESERVED
//
// Permission is granted free of charge to copy, modify, use and distribute
// this software provided you include the entirety of this notice in all
// copies made.
//
// THIS SOFTWARE IS PROVIDED ON AN AS IS BASIS, WITHOUT WARRANTY OF ANY
// KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, WITHOUT LIMITATION,
// WARRANTIES THAT THE SUBJECT SOFTWARE IS FREE OF DEFECTS, MERCHANTABLE, FIT
// FOR A PARTICULAR PURPOSE OR NON-INFRINGING. SGI ASSUMES NO RISK AS TO THE
// QUALITY AND PERFORMANCE OF THE SOFTWARE. SHOULD THE SOFTWARE PROVE
// DEFECTIVE IN ANY RESPECT, SGI ASSUMES NO COST OR LIABILITY FOR ANY
// SERVICING, REPAIR OR CORRECTION. THIS DISCLAIMER OF WARRANTY CONSTITUTES
// AN ESSENTIAL PART OF THIS LICENSE. NO USE OF ANY SUBJECT SOFTWARE IS
// AUTHORIZED HEREUNDER EXCEPT UNDER THIS DISCLAIMER.
//
// UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, WHETHER TORT (INCLUDING,
// WITHOUT LIMITATION, NEGLIGENCE OR STRICT LIABILITY), CONTRACT, OR
// OTHERWISE, SHALL SGI BE LIABLE FOR ANY DIRECT, INDIRECT, SPECIAL,
// INCIDENTAL, OR CONSEQUENTIAL DAMAGES OF ANY CHARACTER WITH RESPECT TO THE
// SOFTWARE INCLUDING, WITHOUT LIMITATION, DAMAGES FOR LOSS OF GOODWILL, WORK
// STOPPAGE, LOSS OF DATA, COMPUTER FAILURE OR MALFUNCTION, OR ANY AND ALL
// OTHER COMMERCIAL DAMAGES OR LOSSES, EVEN IF SGI SHALL HAVE BEEN INFORMED OF
// THE POSSIBILITY OF SUCH DAMAGES. THIS LIMITATION OF LIABILITY SHALL NOT
// APPLY TO LIABILITY RESULTING FROM SGI's NEGLIGENCE TO THE EXTENT APPLICABLE
// LAW PROHIBITS SUCH LIMITATION. SOME JURISDICTIONS DO NOT ALLOW THE
// EXCLUSION OR LIMITATION OF INCIDENTAL OR CONSEQUENTIAL DAMAGES, SO THAT
// EXCLUSION AND LIMITATION MAY NOT APPLY TO YOU.
//
// These license terms shall be governed by and construed in accordance with
// the laws of the United States and the State of California as applied to
// agreements entered into and to be performed entirely within California
// between California residents. Any litigation relating to these license
// terms shall be subject to the exclusive jurisdiction of the Federal Courts
// of the Northern District of California (or, absent subject matter
// jurisdiction in such courts, the courts of the State of California), with
// venue lying exclusively in Santa Clara County, California.
// Superficial sanity test for the constructive reals package.
// Added test for division by negative number. Hans_Boehm@hp.com, 8/13/01
// Modified to use AssertionFailedError. hboehm@google.com, 6/6/14
// Modified further for Android/JUnit testing framework, 12/15/14
// Added basic asin and acos tests, improved messages,
// hboehm@google.com, 5/22/15
package com.hp.creals;
import java.math.BigInteger;
import junit.framework.AssertionFailedError;
import junit.framework.TestCase;
public class CRTest extends TestCase {
private static void check(boolean x, String s) {
if (!x) throw new AssertionFailedError(s);
}
private static void check_eq(CR x, CR y, String s) {
if (x.compareTo(y, -50) != 0) {
throw new AssertionFailedError(s + "(" + x + " vs. " + y + ")");
}
}
private static void check_appr_eq(double x, double y, String s) {
if (Math.abs(x - y) > 0.000001) {
throw new AssertionFailedError(s + "(" + x + " vs. " + y + ")");
}
}
// TODO: Break this up into meaningful smaller tests.
public void testCR() {
CR zero = CR.valueOf(0);
CR one = CR.valueOf(1);
CR two = CR.valueOf(2);
check(one.signum() == 1, "signum(1) failed");
check(one.negate().signum() == -1, "signum(-1) failed");
check(zero.signum(-100) == 0, "signum(0) failed");
check(one.compareTo(two, -10) == -1, "comparison failed");
check(two.toString(4).compareTo("2.0000") == 0, "toString failed");
check_eq(one.shiftLeft(1),two, "shiftLeft failed");
check_eq(two.shiftRight(1),one, "shiftRight failed");
check_eq(one.add(one),two, "add failed 1");
check_eq(one.max(two),two, "max failed");
check_eq(one.min(two),one, "min failed");
check_eq(one.abs(),one, "abs failed 1");
check_eq(one.negate().abs(),one, "abs failed 2");
CR three = two.add(one);
CR four = two.add(two);
check_eq(CR.valueOf(4), four, "2 + 2 failed");
check_eq(CR.valueOf(3), three, "2 + 1 failed");
check_eq(one.negate().add(two), one, "negate failed");
check(one.negate().signum() == -1, "signum(-1) failed");
check_eq(two.multiply(two), four, "multiply failed");
check_eq(one.divide(four).shiftLeft(4), four, "divide failed 1");
check_eq(two.divide(one.negate()), two.negate(), "divide(neg) failed");
CR thirteen = CR.valueOf(13);
check_eq(one.divide(thirteen).multiply(thirteen), one,
"divide failed 2");
check(thirteen.floatValue() == 13.0, "floatValue failed");
check(thirteen.intValue() == 13, "intValue failed");
check(thirteen.longValue() == 13, "longValue failed");
check(thirteen.doubleValue() == 13.0, "doubleValue failed");
check_eq(zero.exp(), one, "exp(0) failed");
CR e = one.exp();
check(e.toString(20).substring(0,17)
.compareTo("2.718281828459045") == 0,
"exp(1) failed");
check_eq(e.ln(), one, "ln(e) failed");
CR half_pi = CR.PI.divide(two);
CR half = one.divide(two);
BigInteger million = BigInteger.valueOf(1000*1000);
BigInteger thousand = BigInteger.valueOf(1000);
CR huge = CR.valueOf(million.multiply(million).multiply(thousand));
UnaryCRFunction asin = UnaryCRFunction.asinFunction;
UnaryCRFunction acos = UnaryCRFunction.acosFunction;
UnaryCRFunction atan = UnaryCRFunction.atanFunction;
UnaryCRFunction tan = UnaryCRFunction.tanFunction;
check_eq(half_pi.sin(), one, "sin(pi/2) failed");
check_eq(asin.execute(one),half_pi, "asin(1) failed");
check_eq(asin.execute(one.negate()),
half_pi.negate(), "asin(-1) failed");
check_eq(asin.execute(zero), zero, "asin(0) failed");
check_eq(asin.execute(half.sin()), half, "asin(sin(0.5)) failed");
UnaryCRFunction cosine = UnaryCRFunction.sinFunction
.monotoneDerivative(zero, CR.PI);
check_eq(cosine.execute(one), one.cos(), "monotoneDerivative failed");
check_eq(cosine.execute(three), three.cos(),
"monotoneDerivative failed 2");
check_eq(asin.execute(one.sin()), one, "asin(sin(1) failed");
check_eq(acos.execute(one.cos()), one, "acos(cos(1) failed");
check_eq(atan.execute(tan.execute(one)), one, "atan(tan(1) failed");
check_eq(atan.execute(tan.execute(one.negate())), one.negate(),
"atan(tan(-1) failed");
check_eq(tan.execute(atan.execute(huge)), huge,
"tan(atan(10**15)) failed");
CR sqrt13 = thirteen.sqrt();
check_eq(sqrt13.multiply(sqrt13), thirteen, "sqrt(13)*sqrt(13) failed");
CR tmp = CR.PI.add(CR.valueOf(-123).exp());
CR tmp2 = tmp.subtract(CR.PI);
check(tmp2.ln().intValue() == -123, "intValue(...) failed");
check(tmp2.ln().longValue() == -123, "longValue(...) failed");
check(tmp2.ln().floatValue() == -123.0, "floatValue(...) failed");
check(tmp2.ln().doubleValue() == -123.0, "doubleValue(...) failed");
for (double n = -10.0; n < 10.0; n += 2.0) {
check_appr_eq(Math.sin(n), CR.valueOf(n).sin().doubleValue(),
"sin failed at " + n);
check_appr_eq(Math.cos(n), CR.valueOf(n).cos().doubleValue(),
"cos failed at " + n);
check_appr_eq(Math.exp(n), CR.valueOf(n).exp().doubleValue(),
"exp failed at " + n);
check_appr_eq(Math.asin(0.1*n),
CR.valueOf(0.1*n).asin().doubleValue(),
"asin failed at " + 0.1*n);
check_appr_eq(Math.acos(0.1*n),
CR.valueOf(0.1*n).acos().doubleValue(),
"acos failed at " + 0.1*n);
if (n > 0.0) {
check_appr_eq(Math.log(n), CR.valueOf(n).ln().doubleValue(),
"ln failed at " + n);
}
}
check_appr_eq(Math.cos(12345678.0),
CR.valueOf(12345678).cos().doubleValue(),
"cos failed at " + 12345678);
}
}