4546474849505152535455
fail("*** SVD: Inconsistent solution vector"); } // System.out.println("svd: rank = " << svd.rank() << " nullity = " << svd.nullity() << endl; localflag = (svd.rank() != 50 || svd.nullity() != 0); globalflag = globalflag || localflag; if (localflag) { fail("*** SVD: Unexpected rank or nullity in original matrix"); }
5859606162636465666768
a[3][i]=a[2][i]-a[1][i]; a[4][i]=0.0; } SVD svd2 = new SVD(a); // System.out.println("svd: rank = " << svd2.rank() << " nullity = " << svd2.nullity() << endl; localflag = (svd2.rank() != 48 || svd2.nullity() != 2); globalflag = globalflag || localflag; if (localflag) { fail("*** SVD: Unexpected rank or nullity in modified matrix"); }
1213141516171819202122
SVD svd = new SVD(m); double b[] = { 3, 1, 2 }; double[] x = new double[3]; System.out.println(svd.inv_condition()); System.out.println(svd.rank()); System.out.println(svd.nullity()); svd.solve(b, x); System.out.println(NRUtil.toString(x)); } }