assertEquals("angular Distance=" + res2, 1.0, res1 / res2, 1e-10);
/* Check known problem with vecmath for small vector differences */
Vector3d v3 = new Vector3d(1.5, 1.6, 0.);
Vector3d v4 = new Vector3d(1.5, 1.601, 0.);
double res3 = pt.AngDist(v3, v4);
double res4 = v3.angle(v4);
double expected = res4 - Math.PI / 2.;
System.out.println("res3 = " + res3 + " res4=" + res4);
assertEquals("angular Distance=" + res4, 1., (Math.PI / 2. + res3)
/ res4, 1e-3);