Examples of SimpleMatrix

com.alibaba.security.simpleimage.analyze.sift.matrix.SimpleMatrix
类SimpleMatrix.java的实现描述：TODO 类实现描述 @author axman 2013-3-25 下午4:41:17
org.ejml.simple.SimpleMatrix
nist.gov/javanumerics/jama/
@author Peter Abeles

Examples of org.ejml.simple.SimpleMatrix

            }
        }
    }


    private void checkNaive(int m, int n) {
        SimpleMatrix A = SimpleMatrix.wrap(RandomMatrices.createRandom(m,n,rand));


        BidiagonalDecompositionRow decomp = new BidiagonalDecompositionRow();
        BidiagonalDecompositionNaive naive = new BidiagonalDecompositionNaive();


        assertTrue(decomp.decompose(A.getMatrix().copy()));
        assertTrue(naive.decompose(A.getMatrix()));


        SimpleMatrix U = SimpleMatrix.wrap(decomp.getU(null,false,false));
        SimpleMatrix B = SimpleMatrix.wrap(decomp.getB(null,false));
        SimpleMatrix V = SimpleMatrix.wrap(decomp.getV(null,false,false));


//        U.print();
//        B.print();
//        naive.getB().print();
//        V.print();
//        naive.getV().print();


//        naive.getVTran().print();


        assertTrue(naive.getB().isIdentical(B,1e-8));
        assertTrue(naive.getU().isIdentical(U,1e-8));
        assertTrue(naive.getV().isIdentical(V,1e-8));


        // check the decomposition
        DenseMatrix64F foundA = U.mult(B).mult(V.transpose()).getMatrix();


//        A.print();
//        foundA.print();


        assertTrue(MatrixFeatures.isIdentical(A.getMatrix(),foundA,1e-8));

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        for( int width = 1; width <= 3*r; width++ ) {
//            System.out.println("width "+width);
            int end = width;
            int offset = width > 1 ? 1 : 0;


            SimpleMatrix A = SimpleMatrix.random(r,width,-1,1,rand);
            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Bb = Ab.copy();


            SimpleMatrix b = A.extractVector(true,rowA).scale(alpha);


            BlockVectorOps.scale_row(r,new D1Submatrix64F(Ab),rowA, alpha, new D1Submatrix64F(Bb),rowB,offset,end);


            checkVector_row(rowB, end, offset, A, Bb, b);
        }

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        for( int width = 1; width <= 3*r; width++ ) {
//            System.out.println("width "+width);
            int end = width;
            int offset = width > 1 ? 1 : 0;


            SimpleMatrix A = SimpleMatrix.random(r,width,-1,1,rand);
            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Bb = Ab.copy();


            SimpleMatrix b = A.extractVector(true,rowA).divide(alpha);


            BlockVectorOps.div_row(r,new D1Submatrix64F(Ab),rowA, alpha, new D1Submatrix64F(Bb),rowB,offset,end);


            checkVector_row(rowB, end, offset, A, Bb, b);
        }

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        checkAgainstRandom(5, 7);
        checkAgainstRandom(2, 3);
    }


    private void checkAgainstRandom(int m, int n) {
        SimpleMatrix A = SimpleMatrix.wrap(RandomMatrices.createRandom(m,n,rand));


        BidiagonalDecompositionNaive decomp = new BidiagonalDecompositionNaive();


        assertTrue(decomp.decompose(A.getMatrix()));


        SimpleMatrix U = decomp.getU();
        SimpleMatrix B = decomp.getB();
        SimpleMatrix V = decomp.getV();


//        U.print();
//        B.print();
//        V.print();


//        U.mult(A).mult(V).print();


        // check the decomposition
        DenseMatrix64F foundA = U.mult(B).mult(V.transpose()).getMatrix();


//        A.print();
//        foundA.print();


        assertTrue(MatrixFeatures.isIdentical(A.getMatrix(),foundA,1e-8));

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        for( int width = 1; width <= 3*r; width++ ) {
//            System.out.println("width "+width);
            int end = width;
            int offset = width > 1 ? 1 : 0;


            SimpleMatrix A = SimpleMatrix.random(r,width,-1,1,rand);
            SimpleMatrix B = SimpleMatrix.random(r,width,-1,1,rand);
            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Bb = BlockMatrixOps.convert(B.getMatrix(),r);
            BlockMatrix64F Cb = Ab.copy();


            SimpleMatrix a = A.extractVector(true,rowA).scale(alpha);
            SimpleMatrix b = B.extractVector(true,rowB).scale(beta);
            SimpleMatrix c = a.plus(b);


            BlockVectorOps.add_row(r,
                    new D1Submatrix64F(Ab),rowA, alpha,
                    new D1Submatrix64F(Bb),rowB, beta ,
                    new D1Submatrix64F(Cb),rowC, offset,end);

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        for( int width = 1; width <= 3*r; width++ ) {
//            System.out.println("width "+width);
            int end = width;
            int offset = width > 1 ? 1 : 0;


            SimpleMatrix A = SimpleMatrix.random(r,width,-1,1,rand);
            SimpleMatrix a = A.extractMatrix(rowA,rowA+1,offset,SimpleMatrix.END);
            SimpleMatrix B = SimpleMatrix.random(r,width,-1,1,rand);
            SimpleMatrix b = B.extractMatrix(rowB,rowB+1,offset,SimpleMatrix.END);


            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Bb = BlockMatrixOps.convert(B.getMatrix(),r);


            double expected = a.dot(b);

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

//            System.out.println("width "+width);
            int end = width;
            int offset = width > 1 ? 1 : 0;
            if( colB >= width) colB = 0;


            SimpleMatrix A = SimpleMatrix.random(width,width,-1,1,rand);
            SimpleMatrix a = A.extractMatrix(rowA,rowA+1,offset,SimpleMatrix.END);
            SimpleMatrix B = SimpleMatrix.random(width,width,-1,1,rand);
            SimpleMatrix b = B.extractMatrix(offset,SimpleMatrix.END,colB,colB+1);


            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Bb = BlockMatrixOps.convert(B.getMatrix(),r);


            double expected = a.dot(b);

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        return num;
    }


    private void checkComponents( SingularValueDecomposition<DenseMatrix64F> svd , DenseMatrix64F expected )
    {
        SimpleMatrix U = SimpleMatrix.wrap(svd.getU(null,false));
        SimpleMatrix Vt = SimpleMatrix.wrap(svd.getV(null,true));
        SimpleMatrix W = SimpleMatrix.wrap(svd.getW(null));


        assertTrue( !U.hasUncountable() );
        assertTrue( !Vt.hasUncountable() );
        assertTrue( !W.hasUncountable() );


        if( svd.isCompact() ) {
            assertEquals(W.numCols(),W.numRows());
            assertEquals(U.numCols(),W.numRows());
            assertEquals(Vt.numRows(),W.numCols());
        } else {
            assertEquals(U.numCols(),W.numRows());
            assertEquals(W.numCols(),Vt.numRows());
            assertEquals(U.numCols(),U.numRows());
            assertEquals(Vt.numCols(),Vt.numRows());
        }


        DenseMatrix64F found = U.mult(W).mult(Vt).getMatrix();

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        }
    }


    @Test
    public void testIdentity() {
        SimpleMatrix A = SimpleMatrix.identity(5);


        BidiagonalDecomposition<DenseMatrix64F> decomp = createQRDecomposition();


        assertTrue(decomp.decompose(A.getMatrix().copy()));


        checkGeneric(A.getMatrix(), decomp);
    }

View Full Code Here

Examples of org.ejml.simple.SimpleMatrix

        checkGeneric(A.getMatrix(), decomp);
    }


    @Test
    public void testZero() {
        SimpleMatrix A = new SimpleMatrix(5,5);


        BidiagonalDecomposition<DenseMatrix64F> decomp = createQRDecomposition();


        assertTrue(decomp.decompose(A.getMatrix().copy()));


        checkGeneric(A.getMatrix(), decomp);
    }

View Full Code Here

0 1 2 3 4 5

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.