ge_p2_dbl.q 625 B

1234567891011121314151617181920212223242526272829303132333435363738394041
  1. :name:fe:r->X:r->Y:r->Z:r->T:t0:t1:t2:t3:t4:t5:p->X:p->Y:p->Z:
  2. fe r:var/r=fe:
  3. enter f:enter/f:>X1=fe#11:>Y1=fe#12:>Z1=fe#13:
  4. return:nofallthrough:<X3=fe#1:<Y3=fe#2:<Z3=fe#3:<T3=fe#4:leave:
  5. h=f+g:<f=fe:<g=fe:>h=fe:asm/fe_add(>h,<f,<g);:
  6. h=f-g:<f=fe:<g=fe:>h=fe:asm/fe_sub(>h,<f,<g);:
  7. h=f*g:<f=fe:<g=fe:>h=fe:asm/fe_mul(>h,<f,<g);:
  8. h=f^2:<f=fe:>h=fe:asm/fe_sq(>h,<f);:
  9. h=2*f^2:<f=fe:>h=fe:asm/fe_sq2(>h,<f);:
  10. h=2*g:<g=fe:>h=fe:asm/fe_add(>h,<g,<g);:
  11. :
  12. enter ge_p2_dbl
  13. fe X1
  14. fe Y1
  15. fe Z1
  16. fe A
  17. fe AA
  18. fe XX
  19. fe YY
  20. fe B
  21. fe X3
  22. fe Y3
  23. fe Z3
  24. fe T3
  25. XX=X1^2
  26. YY=Y1^2
  27. B=2*Z1^2
  28. A=X1+Y1
  29. AA=A^2
  30. Y3=YY+XX
  31. Z3=YY-XX
  32. X3=AA-Y3
  33. T3=B-Z3
  34. return