Examples of LUDecompositionAlt

org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

An LU decomposition algorithm that originally came from Jama. In general this is faster than what is in NR since it creates a cache of a column, which makes a big difference in larger matrices.
@author Peter Abeles

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

        shouldFailSingular = false;
    }


    @Override
    protected LinearSolver createSolver( DenseMatrix64F A ) {
        LUDecompositionAlt decomp = new LUDecompositionAlt();


        return new LinearSolverLu(decomp);
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

        shouldFailSingular = false;
    }


    @Override
    protected LinearSolver<DenseMatrix64F> createSolver( DenseMatrix64F A ) {
        LUDecompositionAlt decomp = new LUDecompositionAlt();


        return new LinearSolverLuKJI(decomp);
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

        DenseMatrix64F b = RandomMatrices.createRandom(2,5,rand);
        DenseMatrix64F c = RandomMatrices.createRandom(2,5,rand);
        DenseMatrix64F c_exp = RandomMatrices.createRandom(2,5,rand);


        assertTrue(CommonOps.solve(a,b,c));
        LUDecompositionAlt alg = new LUDecompositionAlt();
        LinearSolverLu solver = new LinearSolverLu(alg);
        assertTrue(solver.setA(a));


        solver.solve(b,c_exp);

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

    @Test
    public void invert() {
        for( int i = 1; i <= 10; i++ ) {
            DenseMatrix64F a = RandomMatrices.createRandom(i,i,rand);


            LUDecompositionAlt lu = new LUDecompositionAlt();
            LinearSolverLu solver = new LinearSolverLu(lu);
            assertTrue(solver.setA(a));


            DenseMatrix64F a_inv = new DenseMatrix64F(i,i);
            DenseMatrix64F a_lu = new DenseMatrix64F(i,i);

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

//        System.out.println("invert GJ No Pivot     = "+ invertBenchmark(
//                new GaussJordanNoPivot(),mat,numTrials));
//        System.out.println("invert GJ              = "+ invertBenchmark(
//                new GaussJordan(mat.numRows),mat,numTrials));
        System.out.println("invert LU              = "+ invertBenchmark(
                new LinearSolverLu(new LUDecompositionAlt()),mat,numTrials));
        System.out.println("invert LU  NR          = "+ invertBenchmark(
                new LinearSolverLu(new LUDecompositionNR()),mat,numTrials));
        System.out.println("invert Ops             = "+
                invertOpsBenchmark(mat,numTrials));
//        System.out.println("unrolled               = "+

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

    }


    private static void runAlgorithms( int numTrials )
    {
        System.out.println("solve LU A            = "+ solveBenchmark(
                new LinearSolverLu(new LUDecompositionAlt()),numTrials));
        System.out.println("solve LU B            = "+ solveBenchmark(
                new LinearSolverLuKJI(new LUDecompositionAlt()),numTrials));
        System.out.println("solve QR house        = "+ solveBenchmark(
                new LinearSolverQrHouse(),numTrials));
        System.out.println("solve QR house Col    = "+ solveBenchmark(
                new LinearSolverQrHouseCol(),numTrials));
        System.out.println("solve PInv            = "+ solveBenchmark(

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt


    public static long computeLU( DenseMatrix64F mat , int numTrials )
    {
        long before = System.currentTimeMillis();


        LUDecompositionAlt alg = new LUDecompositionAlt();


        double total = 0;


        for( int i = 0; i < numTrials; i++ ) {
            alg.decompose(mat);
            total += alg.computeDeterminant();
        }


//        System.out.println("   total = "+total);
        long after = System.currentTimeMillis();

View Full Code Here

0 1

TOP

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.