constraint: d1 + d2 + d3 = 0
e2 = d1*d2 + d2*d3 + d3*d1
delta = d1*d2*d3
s_1 = 0
s_3 = 3*delta
s_5 = -5*e2*delta
s_7 = 7*e2^2*delta
s_9 = (-9*e2^3 + 3*delta^2)*delta
coefficients of ell^(k) in H/delta:
  k=3: -1
  k=5: e2/12
  k=7: -e2^2/360
  k=9: e2^3/20160 - delta^2/60480
normalized root block H/delta in Weierstrass derivatives:
  wp^(1): 1
  wp^(3): -e2/12
  wp^(5): e2^2/360
  wp^(7): -e2^3/20160 + delta^2/60480
all exact identities passed