D1 geodesic shooting, 200 arcs before cut time: max rel |d0(end) - length| = 9.92e-15
D1b past cut time: length 8.1681, d0(end) 5.2768 (must be smaller)
D2 max |d(p,q) - d(q,p)| over 500 pairs = 0.00e+00
S0 sphere sampler: max |d0(q) - r|/r = 1.70e-14
HF tau at p0=(0.7, -0.4, 0.9): r=0.2: 1.00000; r=0.05: 1.00000; r=0.0125: 1.00000; r=0.003: 1.00000
HF F (Siegel inversion, CR) at p0=(0.7, -0.4, 0.9): r=0.2: 1.34052; r=0.05: 1.07535; r=0.0125: 1.01832; r=0.003: 1.00437
HF tau o F (anti-CR) at p0=(0.7, -0.4, 0.9): r=0.2: 1.34052; r=0.05: 1.07535; r=0.0125: 1.01832; r=0.003: 1.00437
HF A=(2x,y,2t) control at p0=(0.7, -0.4, 0.9): r=0.2: 2.00000; r=0.05: 2.00000; r=0.0125: 2.00000; r=0.003: 2.00000
