Package org.ejml.data

Examples of org.ejml.data.D1Submatrix64F

53
54
55
56
57
58
59
60
61
62
            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
77
78
79
80
81
82
83
84
85
86
            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
108
109
110
111
112
113
114
115
116
117
118
119
            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);

            checkVector_row(rowC, end, offset, A, Cb, c);
        }
    }
View Full Code Here
137
138
139
140
141
142
143
144
145
146
            BlockMatrix64F Ab = BlockMatrixOps.convert(A.getMatrix(),r);
            BlockMatrix64F Bb = BlockMatrixOps.convert(B.getMatrix(),r);

            double expected = a.dot(b);

            double found = BlockVectorOps.dot_row(r,new D1Submatrix64F(Ab),rowA, new D1Submatrix64F(Bb),rowB,offset,end);

            assertEquals(expected,found,1e-8);
        }
    }
View Full Code Here
166
167
168
169
170
171
172
173
174
175
176
177
            BlockMatrix64F Bb = BlockMatrixOps.convert(B.getMatrix(),r);

            double expected = a.dot(b);

            double found = BlockVectorOps.dot_row_col(r,
                    new D1Submatrix64F(Ab),rowA,
                    new D1Submatrix64F(Bb),colB,
                    offset,end);

            assertEquals(expected,found,1e-8);
        }
    }
View Full Code Here
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
            BlockMatrix64F T = BlockMatrixOps.createRandom(size,size,-1,1,rand,r);
            BlockMatrixOps.zeroTriangle(true,T);

            BlockMatrix64F T_inv = T.copy();

            BlockTriangularSolver.invert(r,false,new D1Submatrix64F(T),new D1Submatrix64F(T_inv),temp);

            BlockMatrix64F C = new BlockMatrix64F(size,size,r);

            BlockMatrixOps.mult(T,T_inv,C);

            assertTrue(GenericMatrixOps.isIdentity(C,1e-8));

            // see if passing in the same matrix instance twice messes it up or not
            BlockTriangularSolver.invert(r,false,new D1Submatrix64F(T),new D1Submatrix64F(T),temp);

            assertTrue(BlockMatrixOps.isEquals(T,T_inv,1e-8));
        }
    }
View Full Code Here
78
79
80
81
82
83
84
85
86
87
88
            BlockMatrix64F T = BlockMatrixOps.createRandom(size,size,-1,1,rand,r);
            BlockMatrixOps.zeroTriangle(true,T);

            BlockMatrix64F T_inv = T.copy();

            BlockTriangularSolver.invert(r,false,new D1Submatrix64F(T_inv),temp);

            BlockMatrix64F C = new BlockMatrix64F(size,size,r);

            BlockMatrixOps.mult(T,T_inv,C);
View Full Code Here
143
144
145
146
147
148
149
150
151
            // Compute Y directly from the expected result B
            BlockMatrixOps.mult(T,B,Y);
        }

        // Y is overwritten with the solution
        BlockTriangularSolver.solve(r,upper,new D1Submatrix64F(T),new D1Submatrix64F(Y),transT);

        assertTrue( BlockMatrixOps.isEquals(B,Y,1e-8));
    }
View Full Code Here
177
178
179
180
181
182
183
184
185
186
        }

        int size = T.numRows;

        // Y is overwritten with the solution
        BlockTriangularSolver.solve(r,upper,new D1Submatrix64F(T2,0,size,0,size),new D1Submatrix64F(Y),transT);

        assertTrue( "Failed upper = "+upper+" transT = "+transT+" T.length "+T.numRows+" B.cols "+B.numCols,
                BlockMatrixOps.isEquals(B,Y,1e-8));
    }
View Full Code Here
222
223
224
225
226
227
228
229
230
231
232
233
        // do it again using block matrices
        BlockMatrix64F b_L = BlockMatrixOps.convert(L,3);
        BlockMatrix64F b_B = BlockMatrixOps.convert(B,3);

        D1Submatrix64F sub_L = new D1Submatrix64F(b_L,0, 3, 0, 3);
        D1Submatrix64F sub_B = new D1Submatrix64F(b_B,0, 3, 0, 5);

        if( transT ) {
            sub_L.original = BlockMatrixOps.transpose((BlockMatrix64F)sub_L.original,null);
            TestBlockMultiplication.transposeSub(sub_L);
        }
View Full Code Here
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
    • 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.