public void testSymmetric33Special() { double[][] v = new double[3][3]; double[] d = new double[3]; double[][][] as = {ASMALL,A100,A110,A111,ATEST1}; for (double[][] a:as) { Eigen.solveSymmetric33(a,v,d); check(a,v,d); } }
for (nloop=0; s.time()<maxtime; ++nloop) { for (int irand=0; irand<nrand; ++irand) { Eigen.solveSymmetric33(a[irand],v,d);
a[2][1] = g23[i3][i2][i1]; a[2][2] = g33[i3][i2][i1]; Eigen.solveSymmetric33(a,z,e); float u1i = (float)z[0][0]; float u2i = (float)z[0][1];
public void testSymmetric33() { double[][] v = new double[3][3]; double[] d = new double[3]; int nrand = 10000; for (int irand=0; irand<nrand; ++irand) { //double[][] a = randdouble(3,3); //a = add(a,transpose(a)); double[][] a = makeRandomSymmetric33(); Eigen.solveSymmetric33(a,v,d); check(a,v,d); } }
double[] ev = new double[3]; Eigen.solveSymmetric33(aa,vv,ev); // slow but accurate! double[] u = vv[0]; double[] w = vv[2];