double[][] u = svd.getU(); double[] lambda = svd.getSingularValues(); lambda[0] = FastMath.sqrt(Math.abs(lambda[0]));
Matrix u = svd.getU(); double[] lambda = svd.getSingularValues(); for(int i = 0; i < tdim; i++) {
@Test public void testJamaExample() { double[][] M = transpose(new double[][] { { 1., 2., 3., 4. }, { 5., 6., 7., 8. }, { 9., 10., 11., 12. } }); SingularValueDecomposition svd = new SingularValueDecomposition(M); // M = U * S * V^T double[][] M2 = times(svd.getU(), timesTranspose(svd.getS(), svd.getV())); assertTrue("Not a proper decomposition.", almostEquals(M, M2, 1e-14)); }
double[][] u = svd.getU(); double[] lambda = svd.getSingularValues();
double[][] u = svd.getU(); double[] lambda = svd.getSingularValues();