Examples of org.ejml.data.D1Submatrix64F

org.ejml.data.D1Submatrix64F

Describes a rectangular submatrix inside of a {@link D1Matrix64F}.

Rows are row0 <= i < row1 and Columns are col0 <= j < col1
@author Peter Abeles

            SimpleMatrix v = y.plus(u.scale(-(gamma/2.0)*u.dot(y)));


            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Vb = BlockMatrixOps.convert(V.getMatrix(),r);


            TridiagonalBlockHelper.computeRowOfV(r,new D1Submatrix64F(Ab),new D1Submatrix64F(Vb),
                    row,gamma);


            for( int i = row+1; i < A.numCols(); i++ ) {
                assertEquals(Vb.get(row,i),v.get(i),1e-8);
            }

View Full Code Here

            double val = A.get(i,col)/max;
            expected += val*val;
        }
        expected = Math.sqrt(expected);


        double found = BlockHouseHolder.computeTauAndDivideCol(r,new D1Submatrix64F(A),col,max);


        assertEquals(expected,found,1e-8);


        for( int i = col; i < A.numRows; i++ ) {
            assertEquals(A_orig.get(i,col)/max , A.get(i,col),1e-8);

View Full Code Here

            double val = A.get(row,j)/max;
            expected += val*val;
        }
        expected = Math.sqrt(expected);


        double found = BlockHouseHolder.computeTauAndDivideRow(r,new D1Submatrix64F(A),row,colStart,max);


        assertEquals(expected,found,1e-8);


        for( int j = colStart; j < A.numCols; j++ ) {
            assertEquals(A_orig.get(row,j)/max , A.get(row,j),1e-8);

View Full Code Here


        // make sure it ignores the first element
        A.set(0,1,100000);
        A.set(5,1,-2346);


        double max = BlockHouseHolder.findMaxCol(r,new D1Submatrix64F(A),1);


        assertEquals(2346,max,1e-8);
    }

View Full Code Here


        // make sure it ignores the first element
        A.set(1,1,100000);
        A.set(1,4,-2346);


        double max = BlockHouseHolder.findMaxRow(r,new D1Submatrix64F(A),1,2);


        assertEquals(2346,max,1e-8);
    }

View Full Code Here

        double temp[] = new double[ r ];


        BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
        BlockMatrix64F Wb = new BlockMatrix64F(Ab.numRows,Ab.numCols,r);


        D1Submatrix64F Ab_sub = new D1Submatrix64F(Ab);
        D1Submatrix64F Wb_sub = new D1Submatrix64F(Wb);


        BlockHouseHolder.computeW_Column(r,Ab_sub,Wb_sub,temp,betas,0);


        // see if the result is the same
        assertTrue(GenericMatrixOps.isEquivalent(Wb,W.getMatrix(),1e-8));

View Full Code Here

        double beta = 1.5;


        BlockMatrix64F Wb = BlockMatrixOps.convert(W.getMatrix(),r);
        BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);


        D1Submatrix64F Wb_sub = new D1Submatrix64F(Wb,0, W.numRows(), 0, r);
        D1Submatrix64F Yb_sub = new D1Submatrix64F(Ab,0, A.numRows(), 0, r);


        BlockHouseHolder.initializeW(r,Wb_sub,Yb_sub,r,beta);


        assertEquals(-beta,Wb.get(0,0),1e-8);

View Full Code Here

        SimpleMatrix T = SimpleMatrix.random(M,1,-1,1,rand);


        // -beta * (V + W*T)
        SimpleMatrix expected = V.plus(W.mult(T)).scale(-beta);


        BlockHouseHolder.computeZ(r,new D1Submatrix64F(Ab,0, A.numRows(), 0, r),new D1Submatrix64F(Aw,0, A.numRows(), 0, r),
                M,T.getMatrix().data,beta);


        for( int i = 0; i < A.numRows(); i++ ) {
            assertEquals(expected.get(i),Aw.get(i,M),1e-8);
        }

View Full Code Here

        SimpleMatrix expected = Y.transpose().mult(V);


        BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
        double found[] = new double[ M ];


        BlockHouseHolder.computeY_t_V(r,new D1Submatrix64F(Ab,0, A.numRows(), 0, r),M,found);


        for( int i = 0; i < M; i++ ) {
            assertEquals(expected.get(i),found[i],1e-8);
        }
    }

View Full Code Here

        if( A.blockLength != B.blockLength || A.blockLength != C.blockLength )
            throw new IllegalArgumentException("Block lengths are not all the same.");


        final int blockLength = A.blockLength;


        D1Submatrix64F Asub = new D1Submatrix64F(A,0, A.numRows, 0, A.numCols);
        D1Submatrix64F Bsub = new D1Submatrix64F(B,0, B.numRows, 0, B.numCols);
        D1Submatrix64F Csub = new D1Submatrix64F(C,0, C.numRows, 0, C.numCols);


        BlockMultiplication.mult(blockLength,Asub,Bsub,Csub);
    }

View Full Code Here

0 1 2 3 4 5

TOP

Related Classes of org.ejml.data.D1Submatrix64F

org.ejml.alg.block.BlockMatrixOps

org.ejml.alg.block.BlockTriangularSolver

org.ejml.alg.block.decomposition.chol.BlockCholeskyOuterForm

org.ejml.alg.block.decomposition.hessenberg.TestTridiagonalBlockHelper

org.ejml.alg.block.decomposition.hessenberg.TestTridiagonalDecompositionBlockHouseholder

org.ejml.alg.block.decomposition.hessenberg.TridiagonalDecompositionBlockHouseholder

org.ejml.alg.block.decomposition.qr.BlockMatrix64HouseholderQR

org.ejml.alg.block.decomposition.qr.TestBlockHouseHolder

org.ejml.alg.block.linsol.chol.BlockCholeskyOuterSolver

org.ejml.alg.block.linsol.qr.BlockQrHouseHolderSolver

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.