blob: 33e49ea2d60960bd7f2052d1433e44442ef9d0e5 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
source 2007 Bernstein--Lange
parameter c2
assume c2 = 2*c
compute A = Z1 Z2
compute B = A^2
compute C = X1 X2
compute D = Y1 Y2
compute E = d C D
compute BB = B^2
compute EE = E^2
compute H = (A+B)^2-BB
compute I = (A+E)^2-EE
compute X3 = (H-I)((X1+Y1)(X2+Y2)-C-D)
compute Y3 = (H+I-2 B)(D-C)
compute Z3 = c2(BB-EE)
|