Package org.ejml.alg.dense.decomposition.lu

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

49
50
51
52
53
54
55
56
57
58
59
            DenseMatrix64F expected = new DenseMatrix64F(N,N);
            DenseMatrix64F found = new DenseMatrix64F(N,N);

            // first compute inverse by LU
            LUDecompositionAlt alg = new LUDecompositionAlt();
            LinearSolverLu solver = new LinearSolverLu(alg);

            assertTrue( solver.setA(A));
            solver.invert(expected);
View Full Code Here
41
42
43
44
45
46
47
48
49
50
51
52
    public void testAll() {
        for( int N = 2; N <= UnrolledDeterminantFromMinor.MAX; N++ ) {
            DenseMatrix64F A = RandomMatrices.createRandom(N,N,rand);

            double unrolled = UnrolledDeterminantFromMinor.det(A);
            LUDecompositionAlt alg = new LUDecompositionAlt();
            assertTrue( alg.decompose(A) );
            double expected = alg.computeDeterminant();

            assertEquals(expected,unrolled,1e-8);
        }
    }
View Full Code Here
87
88
89
90
91
92
93
94
95
96
97
98
99
100
        Random rand = new Random(0xfff);

        for( int width = 5; width < 12; width++ ) {
            DenseMatrix64F A = RandomMatrices.createRandom(width,width,rand);

            LUDecompositionAlt lu = new LUDecompositionAlt();
            lu.decompose(A);

            double luVal = lu.computeDeterminant();

            DeterminantFromMinor minor = new DeterminantFromMinor(width);
            double minorVal = minor.compute(new DenseMatrix64F(width,width, true, A.data));

            assertEquals(luVal,minorVal,1e-6);
View Full Code Here
46
47
48
49
50
51
52
53
54
55
56
    @Test
    public void invert() {
        DenseMatrix64F A = new DenseMatrix64F(3,3, true, 0, 1, 2, -2, 4, 9, 0.5, 0, 5);
        DenseMatrix64F A_inv = RandomMatrices.createRandom(3,3,rand);

        LUDecompositionAlt decomp = new LUDecompositionAlt();
        LinearSolver solver = new LinearSolverLu(decomp);

        solver.setA(A);
        InvertUsingSolve.invert(solver,A,A_inv);
View Full Code Here
57
58
59
60
61
62
63
     * Creates a solver for linear systems.  The A matrix will have dimensions (m,m).
     *
     * @return A new linear solver.
     */
    public static LinearSolver<DenseMatrix64F> linear( int matrixSize ) {
        return new LinearSolverLu(new LUDecompositionAlt());
    }
View Full Code Here
100
101
102
103
104
105
106
     *
     * @parm matrixWidth The matrix size that the decomposition should be optimized for.
     * @return LUDecomposition
     */
    public static LUDecomposition<DenseMatrix64F> lu( int numRows , int numCol ) {
        return new LUDecompositionAlt();
    }
View Full Code Here
628
629
630
631
632
633
634
635
636
637
638
                return UnrolledDeterminantFromMinor.det(mat);
            } else {
                return mat.get(0);
            }
        } else {
            LUDecompositionAlt alg = new LUDecompositionAlt();
            if( !alg.decompose(mat) )
                return 0.0;
            return alg.computeDeterminant();
        }
    }
View Full Code Here
664
665
666
667
668
669
670
671
672
673
674
                UnrolledInverseFromMinor.inv(mat,mat);
            } else {
                mat.set(0, 1.0/mat.get(0));
            }
        } else {
            LUDecompositionAlt alg = new LUDecompositionAlt();
            LinearSolverLu solver = new LinearSolverLu(alg);
            if( solver.setA(mat) ) {
                solver.invert(mat);
            } else {
                return false;
View Full Code Here
710
711
712
713
714
715
716
717
718
719
720
                UnrolledInverseFromMinor.inv(mat,result);
            } else {
                result.set(0,  1.0/mat.get(0));
            }
        } else {
            LUDecompositionAlt alg = new LUDecompositionAlt();
            LinearSolverLu solver = new LinearSolverLu(alg);

            if( !solver.setA(mat))
                return false;
            solver.invert(result);
View Full Code Here
46
47
48
49
50
51
52
53
54
55
56
        DenseMatrix64F A = new DenseMatrix64F(3,3, true, 0, 1, 2, -2, 4, 9, 0.5, 0, 5);
        DenseMatrix64F b = new DenseMatrix64F(3,1, true, 8, 33, 15.5);
        DenseMatrix64F x = RandomMatrices.createRandom(3,1,rand);
        DenseMatrix64F x_improved = new DenseMatrix64F(3,1);

        LUDecompositionAlt alg = new LUDecompositionAlt();

        x_improved.set(x);

        LinearSolverLu solver = new LinearSolverLu(alg);
        assertTrue(solver.setA(A));
View Full Code Here
TOP

Related Classes of org.ejml.alg.dense.decomposition.lu.LUDecompositionAlt

  • org.ejml.alg.dense.linsol.BenchmarkInvertSquare
  • org.ejml.alg.dense.linsol.BenchmarkSolveEq
  • org.ejml.alg.dense.linsol.lu.TestLinearSolverLu
  • org.ejml.alg.dense.linsol.lu.TestLinearSolverLuBase
  • org.ejml.alg.dense.linsol.lu.TestLinearSolverLuKJI
  • org.ejml.alg.dense.linsol.TestInvertUsingSolve
  • org.ejml.alg.dense.misc.BenchmarkDeterminant
  • org.ejml.alg.dense.misc.TestDeterminantFromMinor
  • org.ejml.alg.dense.misc.TestUnrolledDeterminantFromMinor
  • org.ejml.alg.dense.misc.TestUnrolledInverseFromMinor
    • Copyright © 2018 www.massapicom. All rights reserved.
    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.