Examples of CholeskyDecompositionInner

org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

This implementation of a Cholesky decomposition using the inner-product form. For large matrices a block implementation is better. On larger matrices the lower triangular decomposition is significantly faster. This is faster on smaller matrices than {@link CholeskyDecompositionBlock}but much slower on larger matrices.
@author Peter Abeles

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

     *
     * @return A new solver for symmetric positive definite matrices.
     */
    public static LinearSolver<DenseMatrix64F> symmPosDef( int matrixWidth ) {
        if(matrixWidth < EjmlParameters.SWITCH_BLOCK64_CHOLESKY )  {
            CholeskyDecompositionCommon decomp = new CholeskyDecompositionInner(true);
            return new LinearSolverChol(decomp);
        } else {
            if( EjmlParameters.MEMORY == EjmlParameters.MemoryUsage.FASTER )
                return new LinearSolverCholBlock64();
            else {
                CholeskyDecompositionCommon decomp = new CholeskyDecompositionInner(true);
                return new LinearSolverChol(decomp);
            }
        }
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

public class TestLinearSolverChol {


    @Test
    public void standardTests() {


        CholeskyDecompositionInner alg = new CholeskyDecompositionInner(true);
        LinearSolverChol solver = new LinearSolverChol(alg);


        BaseCholeskySolveTests tests = new BaseCholeskySolveTests();
        tests.standardTests(solver);
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

     * @return A new CholeskyDecomposition.
     */
    public static CholeskyDecomposition<DenseMatrix64F> chol( int matrixSize , boolean lower )
    {
        if( matrixSize < EjmlParameters.SWITCH_BLOCK64_CHOLESKY ) {
            return new CholeskyDecompositionInner(lower);
        } else if( EjmlParameters.MEMORY == EjmlParameters.MemoryUsage.FASTER ){
            return new CholeskyDecompositionBlock64(lower);
        } else {
            return new CholeskyDecompositionBlock(EjmlParameters.BLOCK_WIDTH_CHOL);
        }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

     * @param cov The covariance of the stribution.  Not modified.
     */
    public CovarianceRandomDraw( Random rand , DenseMatrix64F cov )
    {
        r = new DenseMatrix64F(cov.numRows,1);
        CholeskyDecompositionInner choleky = new CholeskyDecompositionInner( true);


        choleky.decompose(cov);


        A = choleky.getT();
        this.rand = rand;
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

     */
    public static boolean isPositiveDefinite( DenseMatrix64F A ) {
        if( !isSquare(A))
           return false;


        CholeskyDecompositionInner chol = new CholeskyDecompositionInner(true);
        if( chol.inputModified() )
            A = A.copy();


        return chol.decompose(A);
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

    }


    private static void runAlgorithms( DenseMatrix64F A , DenseMatrix64F b ,int numTrials )
    {
        System.out.println("Solve Cholesky         = "+solve(
                new LinearSolverChol(new CholeskyDecompositionInner(true)),
                A,b,numTrials));
        System.out.println("Solve Cholesky LDL     = "+solve(
                new LinearSolverCholLDL(new CholeskyDecompositionLDL()),
                A,b,numTrials));
    }

View Full Code Here

Examples of org.ejml.alg.dense.decomposition.chol.CholeskyDecompositionInner

//        System.out.println("invert GJ No Pivot     = "+ invertGJ_NoPivot(mat,numTrials));
//        System.out.println("invert GJ               = "+ invertGJ(mat,numTrials));
//        System.out.println("invert LU-NR            = "+ invertLU_nr(mat,numTrials));
//        System.out.println("invert LU-Alt           = "+ invertLU_alt(mat,numTrials));
        System.out.println("invert Cholesky Inner       = "+ invertCholesky(
                new LinearSolverChol(new CholeskyDecompositionInner( true)),
                mat,numTrials));
        System.out.println("invert Cholesky Block Dense = "+ invertCholesky(
                new LinearSolverChol(new CholeskyDecompositionBlock( EjmlParameters.BLOCK_WIDTH_CHOL)),
                mat,numTrials));
//        System.out.println("invert default              = "+ invertCholesky(

View Full Code Here

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.