Advertisement
FSZ

Mathematica LLL

FSZ
Feb 21st, 2020
117
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.92 KB | None | 0 0
  1. emptwob = {};
  2. For[j = 5, j < 7, j++,
  3. For[hone = 1, hone < j, hone++,
  4. For[htwo = 1, htwo < j, htwo++,
  5. If[hone == (j - 1)/2, Continue,
  6. If[htwo == (j - 1)/2, Continue,
  7. If[hone != htwo, R = 6.0000000000000005*(10)^(15);
  8. Q = 1.01*(10)^(32); SuM = Q + R^2; P = (10)^(50);
  9. B = LatticeReduce[{ {1, 0, Round[P*Log[j]]}, {0, 1,
  10. Round[P*Log[2*hone/(j - 1)]]}, {0, 0,
  11. Round[P*Log[2*htwo/(j - 1)]]} }]; D1 = B[[1]];
  12. D2 = B[[2]]; D3 = B[[3]]; F2 = D1;
  13. F3 = D2 - ((D2.F2)/(F2.F2))*F2;
  14. F4 = D3 - ((D3.F2)/(F2.F2))*F2 - ((D3.F3)/(F3.F3))*F3;
  15. epsloc = {Norm[D1]/Norm[F2], Norm[D1]/Norm[F3],
  16. Norm[D1]/Norm[F4] }; ceoone = Max[epsloc];
  17. teeta = Norm[D1]/ceoone;
  18. If[teeta^2 > SuM, LB = (Sqrt[teeta^2 - Q] - R)/P,
  19. Print[{j, hone, htwo}]]; AppendTo[emptwob, LB],
  20. Continue] ] ] ]]]; Print[Min[emptwob]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement