[PASS] rem:bernoulli-general: the three points lie in the upper half disk D^+  (exact rational arithmetic)
[PASS] rem:bernoulli-general: pair terms a_ij = rho1 rho1 - rho2 are positive at the three points  (exact rational arithmetic)  -- {(0, 1): 0.046106106266290345, (0, 2): 0.12675905469487603, (1, 2): 0.018762056328987687}
[PASS] rem:bernoulli-general: T^2/(4 a12 a13 a23) > 1, so no Hermitian determinantal kernel  (exact rational arithmetic)  -- ratio = 9517717790593398661531421149299339000162892188214447234593931084345745631210802384823607/2430837625185067815869656635768991603454329006925359469903536830670798654926728852823607 = 3.9154066450
[PASS] rem:bernoulli-general: the ratio is 3.9154... as printed  (numerical (floating point))  -- 3.9154066450
[PASS] rem:bernoulli-general: the four points lie in D^+ and the configuration is symmetric under z -> -conj z  (exact rational arithmetic)
[PASS] rem:bernoulli-general: rho_2(z,z) = 0 (correlation functions vanish on the diagonal)  (exact rational arithmetic)
[PASS] rem:bernoulli-general: P has degree 4, constant term 1, rational coefficients  (exact rational arithmetic)  -- tau^0: 1; tau^1: -4700693125000/123602168041; tau^2: 10820185818448421389171672037573001707261081651677011994738137938888246429303249856667207922697424893172321491810879687500000000000000/20231885199339142839241952040952947540597893918462506791516297010849354523143059180486650317100961353055451574680866345962330862777; tau^3: -66532559179492487326137906520943706905425506307161450969312245689146291239043310038784189821840667724609375000000000000000000000000000/20231885199339142839241952040952947540597893918462506791516297010849354523143059180486650317100961353055451574680866345962330862777; tau^4: 79660451379755257838986212118509190067717726596122533854029546862707140407566810173493254964209074584319883188815026293396949768066406250000000000000000000000000000000000000000000000000/10705972508193725305582367148243079497497393826458078353455991693219374301588182796030580204208806023224673663204031494082654571683603194065413790649563458912862357089034957447356529
[PASS] rem:bernoulli-general: coefficients match the printed P = 1 - 38.031 tau + 534.81 tau^2 - 3288.5 tau^3 + 7440.7 tau^4  (numerical (floating point))  -- [1.0, -38.0308, 534.8086, -3288.5002, 7440.7487]
[PASS] rem:bernoulli-general: discriminant of P is negative (exact)  (exact rational arithmetic)  -- -1.628590e+04
[PASS] rem:bernoulli-general: P is squarefree: gcd(P, P') = 1 (exact)  (exact rational arithmetic)
[PASS] rem:bernoulli-general: P has exactly two real roots (Sturm), hence one pair of simple nonreal roots  (exact rational arithmetic)  -- real roots: 2
[PASS] rem:bernoulli-general: the nonreal roots are 0.08879 +- 0.00176 i as printed  (numerical (floating point))  -- (0.088790337 - 0.0017573251j), (0.088790337 + 0.0017573251j)
[PASS] rem:bernoulli-general: every Pfaffian used (27) is of an antisymmetric matrix and satisfies Pf^2 = det  (exact rational arithmetic)
== part4_exact_checks: 13 passed, 0 failed, 99.4 s