public void test3Random() { java.util.Random r = new java.util.Random(); int ntest = 1000; for (int itest = 0; itest < ntest; ++itest) { SymmetricTridiagonalFilter stf = makeRandomFilter(); boolean inplace = r.nextBoolean(); // apply in-place? int n1 = 2 + r.nextInt(11); int n2 = 2 + r.nextInt(12); int n3 = 2 + r.nextInt(13); float[][][] t = randfloat(r, n1, n2, n3); float[][][] x = copy(t); float[][][] y = inplace ? x : zerofloat(n1, n2, n3); float[][][] z = inplace ? x : zerofloat(n1, n2, n3); stf.apply1(x, y); stf.applyInverse1(y, z); assertEqual(t, x); assertEqual(t, z); stf.apply2(x, y); stf.applyInverse2(y, z); assertEqual(t, x); assertEqual(t, z); stf.apply3(x, y); stf.applyInverse3(y, z); assertEqual(t, x); assertEqual(t, z); } }
public void test3Simple() { int n1 = 11; int n2 = 12; int n3 = 13; float[][][] r = randfloat(n1, n2, n3); float[][][] x = copy(r); float[][][] y = copy(r); SymmetricTridiagonalFilter stf = new SymmetricTridiagonalFilter(2.6, 2.5, 2.7, 1.2); stf.apply1(x, x); stf.apply2(x, x); stf.apply3(x, x); stf.apply1(y, y); y = transpose12(y); stf.apply1(y, y); y = transpose12(y); y = transpose23(y); stf.apply2(y, y); y = transpose23(y); assertEqual(x, y); }