Package org.ejml.data

Examples of org.ejml.data.RowD1Matrix64F


        DenseMatrix64F B = new DenseMatrix64F(A);

        DenseMatrix64F C = new DenseMatrix64F(m,n);
        DenseMatrix64F u = new DenseMatrix64F(m,1);

        RowD1Matrix64F UBV = alg.getUBV();

        for( int i = 0; i < n; i++ ) {
            alg.computeU(i);

            SpecializedOps.subvector(UBV,i+1,i,m-i-1,false,i+1,u);
            u.data[i] = 1;

            DenseMatrix64F Q = SpecializedOps.createReflector(u,alg.getGammasU()[i]);

            CommonOps.mult(Q,B,C);

//            u.print();
//            B.print();
//            UBV.print();
//            C.print();

            B.set(C);

            // make sure everything is as expected
            for( int j = i+1; j < m; j++ ) {
                assertEquals(0,C.get(j,i),1e-8);
            }

            for( int j = i+1; j < n; j++ ) {
                assertEquals(UBV.get(i,j),C.get(i,j),1e-8);
            }
            u.data[i] = 0;
        }
    }
View Full Code Here


        DenseMatrix64F B = new DenseMatrix64F(A);

        DenseMatrix64F C = new DenseMatrix64F(m,n);
        DenseMatrix64F u = new DenseMatrix64F(n,1);

        RowD1Matrix64F UBV = alg.getUBV();

//        A.print();

        for( int i = 0; i < n-2; i++ ) {
            alg.computeV(i);

            u.zero();
            SpecializedOps.subvector(UBV,i,i+2,n-i-2,true,i+2,u);
            u.data[i+1] = 1;

            DenseMatrix64F Q = SpecializedOps.createReflector(u,alg.getGammasV()[i]);

//            Q.print();

            CommonOps.mult(B,Q,C);

//            u.print();
//            B.print();
//            UBV.print();
//            C.print();

            B.set(C);

            // make sure everything is as expected
            for( int j = i+2; j < n; j++ ) {
                assertEquals(0,C.get(i,j),1e-8);
            }

            for( int j = i+2; j < m; j++ ) {
                assertEquals(UBV.get(j,i),C.get(j,i),1e-8);
            }
            u.data[i] = 0;
        }

    }
View Full Code Here

TOP

Related Classes of org.ejml.data.RowD1Matrix64F

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.