/** * Returns a new array of random values. * @param n1 1st array dimension. */ public static double[] randdouble(int n1) { return randdouble(_random,n1); }
/** * Returns a new array of random values. * @param n1 1st array dimension. * @param n2 2nd array dimension. */ public static double[][] randdouble(int n1, int n2) { return randdouble(_random,n1,n2); }
/** * Returns a new array of random values. * @param n1 1st array dimension. * @param n2 2nd array dimension. * @param n3 3rd array dimension. */ public static double[][][] randdouble(int n1, int n2, int n3) { return randdouble(_random,n1,n2,n3); }
private static void benchDouble( ArrayFile af, int n) throws IOException { double[] a = randdouble(n); double[] b = zerodouble(n); int nio; Stopwatch sw = new Stopwatch(); sw.start(); for (nio=0; sw.time()<5.0; ++nio) { af.seek(0); af.writeDoubles(a); af.seek(0); af.readDoubles(b); } sw.stop(); for (int i=0; i<n; ++i) if (a[i]!=b[i]) throw new RuntimeException("double: i/o failure"); double time = sw.time(); double rate = 2.0*8.0e-6*nio*n/time; System.out.println("double: rate="+rate+" MB/s"); } }
public void testSymmetric22() { int nrand = 10000; double[][] v = new double[2][2]; double[] d = new double[2]; for (int irand=0; irand<nrand; ++irand) { double[][] a = randdouble(2,2); a = add(a, transpose(a)); Eigen.solveSymmetric22(a,v,d); check(a,v,d); } }
double[][][] a = new double[nrand][][]; for (int irand=0; irand<nrand; ++irand) { a[irand] = randdouble(3,3); a[irand] = add(a[irand],transpose(a[irand]));
for (int is=0; is<ns; ++is) { for (int it=0; it<nt; ++it) { DMatrix matrix = new DMatrix(sub(randdouble(3,3),0.5));
private static void testDouble(ArrayFile af, int n) throws IOException double[] a = randdouble(n); double[] b = zerodouble(n);