com.joptimizer.util.UtilsTest.java Source code

Introduction

Here is the source code for com.joptimizer.util.UtilsTest.java

Source

/*
 * Copyright 2011-2016 joptimizer.com
 *
 * This work is licensed under the Creative Commons Attribution-NoDerivatives 4.0 
 * International License. To view a copy of this license, visit 
 *
 *        http://creativecommons.org/licenses/by-nd/4.0/ 
 *
 * or send a letter to Creative Commons, PO Box 1866, Mountain View, CA 94042, USA.
 */
package com.joptimizer.util;

import junit.framework.TestCase;

import org.apache.commons.lang3.ArrayUtils;
import org.apache.commons.logging.Log;
import org.apache.commons.logging.LogFactory;
import org.apache.commons.math3.linear.MatrixUtils;
import org.apache.commons.math3.linear.QRDecomposition;

import cern.colt.matrix.DoubleFactory1D;
import cern.colt.matrix.DoubleFactory2D;
import cern.colt.matrix.DoubleMatrix1D;
import cern.colt.matrix.DoubleMatrix2D;
import cern.colt.matrix.linalg.Algebra;

/**
 * @author alberto trivellato (alberto.trivellato@gmail.com)
 */
public class UtilsTest extends TestCase {

    private Algebra ALG = Algebra.DEFAULT;
    private DoubleFactory1D F1 = DoubleFactory1D.dense;
    private DoubleFactory2D F2 = DoubleFactory2D.dense;
    private Log log = LogFactory.getLog(this.getClass().getName());

    public void testDummy() throws Exception {
        log.debug("testDummy");
    }

    public void testNorm() throws Exception {
        log.debug("testNorm");
        double[] d = new double[] { 1.1, 2.2, 3.3, 4.4, 5.5, 6.6, 7.7, 8.8, 9.9, 10.1 };
        DoubleMatrix1D v = F1.make(d);
        double n0 = v.zDotProduct(v);
        //double n1 = ALG.norm1(v);
        double n2 = ALG.norm2(v);
        log.debug("n0: " + n0);
        //log.debug("n1: " + n1);
        log.debug("n2: " + n2);
        assertEquals(n0, n2);
    }

    public void testCalculateDeterminant() throws Exception {
        log.debug("testCalculateDeterminant");
        double[][] A = new double[][] { { 1, 0, 0 }, { 2, 2, 0 }, { 3, 3, 3 } };
        double det = Utils.calculateDeterminant(A, 3);
        assertEquals(6.0, det);
    }

    public void testCalculateDeterminantBig() throws Exception {
        log.debug("testCalculateDeterminantBig");
        DoubleMatrix2D A = Utils.randomValuesSparseMatrix(40, 40, -1, 1, 0.9, 12345L);
        log.debug("A: " + ArrayUtils.toString(A.toArray()));
        double det = Utils.calculateDeterminant(A.toArray(), A.rows());
        assertEquals(0., det, 1.e-6);
    }

    public void testQRdecomposition() throws Exception {
        log.debug("testQRdecomposition");
        double[][] A = new double[][] { { 1, 0, 0 }, { 2, 2, 0 }, { 3, 3, 3 } };
        QRDecomposition qrFact = new QRDecomposition(MatrixUtils.createRealMatrix(A));
        qrFact.getQ();
    }

    public void testGetExpAndMantissa() throws Exception {
        log.debug("testGetExpAndMantissa");
        float myFloat = -0.22f;
        Utils.getExpAndMantissa(myFloat);

        double myDouble = -0.22e-3;
        Utils.getExpAndMantissa(myDouble);

        myDouble = Math.pow(2, 17);
        Utils.getExpAndMantissa(myDouble);
    }

    public void testRound() throws Exception {
        log.debug("testRound");
        double d = 0.1000000000000009;
        double precision = 1.e10;
        double d2 = Utils.round(d, precision);
        log.debug("d2. " + d2);
        assertEquals(0.1, d2);
    }
}