(*.f64 #s(literal 2 binary64) (/.f64 (/.f64 (-.f64 (/.f64 #s(literal 1 binary64) (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal -4 binary64))) (*.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) (*.f64 (PI.f64) (PI.f64))) #s(literal 1/256 binary64))) (-.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))) (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 (PI.f64) #s(literal -1/4 binary64)))))
(/.f64 (/.f64 #s(literal -1 binary64) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (atan.f64 (exp.f64 (/.f64 y kR))))) (/.f64 #s(literal 1/2 binary64) (-.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)))))
(/.f64 #s(literal 1 binary64) (/.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.f64 (-.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64))) #s(literal 2 binary64))))
(+.f64 (*.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal -1 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (*.f64 y #s(approx (+ (* 1/2 (/ y (pow kR 2))) (/ 1 kR)) (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) (/.f64 y kR))) kR)))))) (/.f64 (PI.f64) #s(literal 4 binary64)))))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal -1 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (*.f64 (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 (*.f64 (/.f64 y (*.f64 kR kR)) #s(literal 1/2 binary64)) y)) (/.f64 y kR)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal -1 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64))))
(if.f64 (<=.f64 (/.f64 y kR) #s(literal -5394370239545291/8299031137761985917024815727382322302024892464484873799991314659381305622825816292799414097894207588576395773222601578364790302150823550615773749668227927374122363606803019047370752 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64))))
(+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))
(-.f64 phip delp)
(asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 (*.f64 (cos.f64 phipp) sinp0) (cos.f64 lampp))))
(asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 (*.f64 (cos.f64 phipp) (cos.f64 lampp)))))
(if.f64 (<=.f64 lampp #s(literal 8920298079412249/44601490397061246283071436545296723011960832 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) (*.f64 sinp0 (cos.f64 phipp))))) (if.f64 (<=.f64 lampp #s(literal 44999999999999997231817817289718218865550436393229298544109048042861511490682889182116529782467557879139928255947822356773107810893791361384810651376521734876411352979114488472976249497974830135699722668048647712322016346347867631619545561838951632079272894529536 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 #s(approx (* (cos phipp) (cos lampp)) (cos.f64 lampp))))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (* sinp0 (+ (* (cos lampp) (cos phipp)) (/ 1 (/ sinp0 (* (sin phipp) cosp0))))) (*.f64 sinp0 (*.f64 (cos.f64 lampp) (cos.f64 phipp))))))))
(if.f64 (<=.f64 phipp #s(literal -4980620899901579/36893488147419103232 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) (*.f64 sinp0 (cos.f64 phipp))))) (if.f64 (<=.f64 phipp #s(literal 1491238505072343/12855504354071922204335696738729300820177623950262342682411008 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (+.f64 (*.f64 cosp0 phipp) (*.f64 sinp0 (cos.f64 lampp))))) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) (*.f64 sinp0 (cos.f64 phipp)))))))
(if.f64 (<=.f64 cosp0 #s(literal -6298479312360475/449891379454319638281053850768598185886969711830191663310075557261183758067148787031904068610389085714992091063352089512320826605549429968900851518086516385513472 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 (+.f64 #s(approx (* (* (cos phipp) sinp0) (/ (cos lampp) cosp0)) (/.f64 #s(approx (* sinp0 (cos phipp)) sinp0) cosp0)) (sin.f64 phipp)) cosp0))) (if.f64 (<=.f64 cosp0 #s(literal 3679531545830007/109836762562089755439710412785302291476310964802292886550311415346968690934362496833960954250583272879636740982263693728593951807995466301001184452657840914432 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (* sinp0 (+ (* (cos lampp) (cos phipp)) (/ 1 (/ sinp0 (* (sin phipp) cosp0))))) (*.f64 sinp0 (*.f64 (cos.f64 lampp) (cos.f64 phipp)))))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 (+.f64 #s(approx (* (* (cos phipp) sinp0) (/ (cos lampp) cosp0)) (/.f64 #s(approx (* sinp0 (cos phipp)) sinp0) cosp0)) (sin.f64 phipp)) cosp0)))))
(if.f64 (<=.f64 phipp #s(literal -760928193040519/2305843009213693952 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 (+.f64 #s(approx (* (* (cos phipp) sinp0) (/ (cos lampp) cosp0)) (/.f64 #s(approx (* sinp0 (cos phipp)) sinp0) cosp0)) (sin.f64 phipp)) cosp0))) (if.f64 (<=.f64 phipp #s(literal 7673845534663173/295147905179352825856 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (+.f64 (*.f64 cosp0 phipp) (*.f64 sinp0 (cos.f64 lampp))))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 (+.f64 #s(approx (* (* (cos phipp) sinp0) (/ (cos lampp) cosp0)) (/.f64 #s(approx (* sinp0 (cos phipp)) sinp0) cosp0)) (sin.f64 phipp)) cosp0)))))
(asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (+.f64 (*.f64 cosp0 phipp) (*.f64 sinp0 (cos.f64 lampp)))))
(if.f64 (<=.f64 cosp0 #s(literal -5224175567749775/18014398509481984 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (+.f64 sinp0 (*.f64 cosp0 phipp))))) (if.f64 (<=.f64 cosp0 #s(literal 5363734393030611/1532495540865888858358347027150309183618739122183602176 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (*.f64 sinp0 (cos.f64 lampp))))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (+.f64 sinp0 (*.f64 cosp0 phipp)))))))
(asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (+.f64 sinp0 (*.f64 cosp0 phipp)))))
(*.f64 ecc (sin.f64 phip))
#s(approx (* ecc (sin phip)) (*.f64 ecc phip))
(asin.f64 sinp0)
#s(approx (/ k0 (/ (- 1 (* (* ecc (sin phi0)) (* ecc (sin phi0)))) (sqrt one_es))) (*.f64 k0 (sqrt.f64 one_es)))
(sqrt.f64 (-.f64 (*.f64 (*.f64 rone_es es) #s(approx (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) (-.f64 #s(literal 1 binary64) (*.f64 (+.f64 phi0 phi0) phi0)))) #s(literal -1 binary64)))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 es #s(approx (* (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) rone_es) rone_es))))
#s(approx (sqrt (+ 1 (* es (* (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) rone_es)))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) #s(approx (* es (* rone_es (pow (cos phi0) 4))) (*.f64 es rone_es)))))
#s(approx (- (log (tan (+ (/ (PI) 4) (* 1/2 phip0)))) (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi0)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))))))) (*.f64 c #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 (*.f64 ecc hlf_e) (*.f64 #s(literal 2 binary64) (sin.f64 phi0))))))
#s(approx (- (log (tan (+ (/ (PI) 4) (* 1/2 phip0)))) (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi0)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))))))) (*.f64 c #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 ecc (*.f64 hlf_e #s(approx (- (sin phi0) (* -1 (sin phi0))) (*.f64 #s(literal 2 binary64) phi0)))))))
#s(approx (- (log (tan (+ (/ (PI) 4) (* 1/2 phip0)))) (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi0)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))))))) (*.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 (+.f64 ecc ecc) phi0)) hlf_e) c))
#s(approx (- (log (tan (+ (/ (PI) 4) (* 1/2 phip0)))) (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi0)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))))))) (*.f64 (*.f64 c #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 (+.f64 ecc ecc) phi0))) hlf_e))
#s(approx (- (log (tan (+ (/ (PI) 4) (* 1/2 phip0)))) (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi0)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))))))) (*.f64 (*.f64 c hlf_e) #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 (+.f64 ecc ecc) phi0))))
(if.f64 (<=.f64 phipp #s(literal -6402987333448277/3369993333393829974333376885877453834204643052817571560137951281152 binary64)) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (if.f64 (<=.f64 phipp #s(literal 4276979833124165/1989292945639146568621528992587283360401824603189390869761855907572637988050133502132224 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 (+.f64 (+.f64 (*.f64 phipp #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (PI.f64))) (PI.f64)))))))
(if.f64 (<=.f64 phipp #s(literal -6402987333448277/3369993333393829974333376885877453834204643052817571560137951281152 binary64)) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (if.f64 (<=.f64 phipp #s(literal 1906546502743257/10141204801825835211973625643008 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(*.f64 kR (/.f64 (-.f64 (pow.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))))) #s(literal 3 binary64)) (pow.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))) #s(literal 3 binary64))) (+.f64 (*.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))) (log.f64 (*.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp)))) (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))))) (pow.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))))) #s(literal 2 binary64)))))
(*.f64 kR (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))))) (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))))))
(*.f64 kR (neg.f64 (log.f64 (/.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))) (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))))))))
(*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(*.f64 kR (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(*.f64 kR lampp)
#s(approx (log (- (/ (* (cos (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))) (/ (* (sin (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))))) (+.f64 (log.f64 (*.f64 #s(literal 1/2 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (*.f64 phi (+.f64 (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (sqrt.f64 #s(literal 2 binary64)))) #s(approx (* phi (+ (* 1/6 (* phi (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (* (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (+ (* 1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4)))))) (* 1/2 (- (* 1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* 1/16 (/ (sqrt 2) (cos (* -1/4 (PI))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3)))))))) (sqrt 2))) (* 16 (/ (* (pow (cos (* -1/4 (PI))) 3) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 3)) (pow (sqrt 2) 3))))))) (* 1/2 (+ (* -4 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 2)) (pow (sqrt 2) 2))) (* 4 (/ (* (cos (* -1/4 (PI))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))))) (sqrt 2))))))) (*.f64 #s(literal 1/6 binary64) (*.f64 (pow.f64 phi #s(literal 2 binary64)) (+.f64 (*.f64 #s(literal -24 binary64) (/.f64 (*.f64 (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))) (*.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (-.f64 (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (pow.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64)))) (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64)))))))))) (pow.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)))) (+.f64 (*.f64 #s(literal 12 binary64) (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal -1/2 binary64) (+.f64 (*.f64 #s(literal -1/48 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (+.f64 (*.f64 #s(literal 1/16 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (pow.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64)) (sqrt.f64 #s(literal 2 binary64))) (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 4 binary64))))))) (*.f64 #s(literal 1/2 binary64) (-.f64 (*.f64 #s(literal 1/48 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (+.f64 (*.f64 #s(literal 1/16 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (pow.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64))))))))) (sqrt.f64 #s(literal 2 binary64)))) (*.f64 #s(literal 16 binary64) (/.f64 (*.f64 (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64)) (pow.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) #s(literal 3 binary64))) (pow.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 3 binary64)))))))))))))
#s(approx (log (- (/ (* (cos (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))) (/ (* (sin (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))))) (+.f64 (log.f64 (*.f64 #s(literal 1/2 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (*.f64 phi #s(approx (+ (* 2 (/ (* (cos (* -1/4 (PI))) (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI))))))) (sqrt 2))) (* phi (+ (* 1/6 (* phi (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (* (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (+ (* 1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4)))))) (* 1/2 (- (* 1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* 1/16 (/ (sqrt 2) (cos (* -1/4 (PI))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3)))))))) (sqrt 2))) (* 16 (/ (* (pow (cos (* -1/4 (PI))) 3) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 3)) (pow (sqrt 2) 3))))))) (* 1/2 (+ (* -4 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 2)) (pow (sqrt 2) 2))) (* 4 (/ (* (cos (* -1/4 (PI))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))))) (sqrt 2)))))))) (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (sqrt.f64 #s(literal 2 binary64))))))))
(log.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 phi #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi)))) (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))
(neg.f64 (log.f64 (/.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (fabs.f64 (sin.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi)))) #s(literal -2028240960365167/20282409603651670423947251286016 binary64)) (log.f64 (tan.f64 (+.f64 #s(approx (+ (* phi 1/2) (+ (* 17/4 (PI)) (PI))) (+.f64 (PI.f64) (*.f64 #s(literal 17/4 binary64) (PI.f64)))) (PI.f64)))) (log.f64 (tan.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi)))) #s(literal -2028240960365167/20282409603651670423947251286016 binary64)) (log.f64 (tan.f64 (+.f64 #s(approx (+ (* phi 1/2) (+ (* 17/4 (PI)) (PI))) (+.f64 (PI.f64) (*.f64 #s(literal 17/4 binary64) (PI.f64)))) (PI.f64)))) (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi)))))
(log.f64 (tan.f64 (+.f64 (+.f64 (*.f64 phi #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 21/4 binary64) (PI.f64)) (PI.f64))) (PI.f64))))
(log.f64 (tan.f64 (+.f64 (*.f64 phi #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (PI.f64)))))
(log.f64 (tan.f64 (+.f64 (PI.f64) #s(approx (+ (* phi 1/2) (+ (* 5/4 (PI)) (PI))) (*.f64 #s(literal 9/4 binary64) (PI.f64))))))
(log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phi)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))
(log.f64 (tan.f64 #s(approx (+ (* phi 1/2) (+ (* 5/4 (PI)) (PI))) (*.f64 #s(literal 9/4 binary64) (PI.f64)))))
(asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 (*.f64 (cos.f64 phip) sinp0) (cos.f64 lamp))))
(asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 (*.f64 (cos.f64 phip) (cos.f64 lamp)))))
(if.f64 (<=.f64 (cos.f64 lamp) #s(literal 4458563631096791/4503599627370496 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 #s(approx (* (cos phip) (cos lamp)) (cos.f64 lamp))))) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) (*.f64 sinp0 (cos.f64 phip))))))
(if.f64 (<=.f64 phip #s(literal -1300000 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) (*.f64 sinp0 (cos.f64 phip))))) (if.f64 (<=.f64 phip #s(literal 7673845534663173/295147905179352825856 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp))))) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) (*.f64 sinp0 (cos.f64 phip)))))))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp)))))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 (neg.f64 (cos.f64 lamp)) sinp0)))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) #s(approx (* -1 (* sinp0 (cos lamp))) (neg.f64 sinp0))))
(if.f64 (<=.f64 phi #s(literal -6687265054221059/361473786714651839609485931802192366508973300717001923159475447150424810286233407987951861887389439612274926783780351561999781998832434041296198795326329101623141899709787663433296905279066051548640942013290819886814068736 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 #s(approx (- (* 1/4 (PI)) (* -1/2 phi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal 2 binary64) (atan.f64 (exp.f64 (+.f64 (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (+.f64 (PI.f64) (PI.f64)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64)))) (+.f64 (PI.f64) (PI.f64))) (PI.f64)))) (*.f64 hlf_e #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 phi (-.f64 ecc (*.f64 #s(literal -1 binary64) ecc))))))) K)))) (/.f64 (PI.f64) #s(literal 2 binary64))))
(if.f64 (<=.f64 phi #s(literal -6687265054221059/361473786714651839609485931802192366508973300717001923159475447150424810286233407987951861887389439612274926783780351561999781998832434041296198795326329101623141899709787663433296905279066051548640942013290819886814068736 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 #s(approx (- (* 1/4 (PI)) (* -1/2 phi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal 2 binary64) (atan.f64 (exp.f64 (+.f64 (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 #s(approx (+ (+ (- (* 1/4 (PI)) (- (* -1/2 phi) (PI))) (+ (PI) (PI))) (PI)) (+.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal 4 binary64) (PI.f64)))) (PI.f64)))) (*.f64 hlf_e #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 phi (-.f64 ecc (*.f64 #s(literal -1 binary64) ecc))))))) K)))) (/.f64 (PI.f64) #s(literal 2 binary64))))
(if.f64 (<=.f64 phi #s(literal -6687265054221059/361473786714651839609485931802192366508973300717001923159475447150424810286233407987951861887389439612274926783780351561999781998832434041296198795326329101623141899709787663433296905279066051548640942013290819886814068736 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 #s(approx (- (* 1/4 (PI)) (* -1/2 phi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 (+.f64 (+.f64 (PI.f64) (PI.f64)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64)))))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(if.f64 (<=.f64 K #s(literal -8820118950590725/21778071482940061661655974875633165533184 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))) (if.f64 (<=.f64 K #s(literal 6749139785887759/1606938044258990275541962092341162602522202993782792835301376 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi)))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))
(if.f64 (<=.f64 K #s(literal -8820118950590725/21778071482940061661655974875633165533184 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))) (if.f64 (<=.f64 K #s(literal 6749139785887759/1606938044258990275541962092341162602522202993782792835301376 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 #s(approx (- (* 1/4 (PI)) (* -1/2 phi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))
(if.f64 (<=.f64 K #s(literal 2331202670670875/1942668892225729070919461906823518906642406839052139521251812409738904285205208498176 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))
(-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (+.f64 (*.f64 #s(approx (* ecc (+ (* 1/6 (* ecc (+ (* -6 (* (sin phi) (pow (- (sin phi) (* -1 (sin phi))) 2))) (+ (* 2 (pow (- (sin phi) (* -1 (sin phi))) 3)) (* 6 (* (pow (sin phi) 2) (- (sin phi) (* -1 (sin phi))))))))) (* 1/2 (+ (* -1 (pow (- (sin phi) (* -1 (sin phi))) 2)) (* 2 (* (sin phi) (- (sin phi) (* -1 (sin phi))))))))) (*.f64 (-.f64 (*.f64 (*.f64 #s(literal 4 binary64) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (*.f64 #s(literal -6 binary64) (sin.f64 phi))) (-.f64 (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64))) (*.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi)) #s(literal -6 binary64))) (*.f64 (pow.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi)) #s(literal 3 binary64)) #s(literal 2 binary64)))) (*.f64 (*.f64 ecc ecc) #s(literal 1/6 binary64)))) ecc) (*.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi)) ecc)))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (-.f64 #s(approx (* ecc (+ (* 1/6 (* ecc (+ (* -6 (* (sin phi) (pow (- (sin phi) (* -1 (sin phi))) 2))) (+ (* 2 (pow (- (sin phi) (* -1 (sin phi))) 3)) (* 6 (* (pow (sin phi) 2) (- (sin phi) (* -1 (sin phi))))))))) (* 1/2 (+ (* -1 (pow (- (sin phi) (* -1 (sin phi))) 2)) (* 2 (* (sin phi) (- (sin phi) (* -1 (sin phi))))))))) (*.f64 (-.f64 (*.f64 (*.f64 #s(literal 4 binary64) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (*.f64 #s(literal -6 binary64) (sin.f64 phi))) (-.f64 (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64))) (*.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi)) #s(literal -6 binary64))) (*.f64 (pow.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi)) #s(literal 3 binary64)) #s(literal 2 binary64)))) (*.f64 (*.f64 ecc ecc) #s(literal 1/6 binary64)))) (*.f64 #s(literal -2 binary64) (sin.f64 phi))) ecc))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 ecc (-.f64 (sin.f64 phi) (*.f64 #s(literal -1 binary64) (sin.f64 phi)))))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) #s(approx (* (+ (+ (* ecc (+ (* 1/6 (* ecc (+ (* -6 (* (sin phi) (pow (- (sin phi) (* -1 (sin phi))) 2))) (+ (* 2 (pow (- (sin phi) (* -1 (sin phi))) 3)) (* 6 (* (pow (sin phi) 2) (- (sin phi) (* -1 (sin phi))))))))) (* 1/2 (+ (* -1 (pow (- (sin phi) (* -1 (sin phi))) 2)) (* 2 (* (sin phi) (- (sin phi) (* -1 (sin phi))))))))) (sin phi)) (sin phi)) ecc) (*.f64 #s(literal 2 binary64) (*.f64 ecc (sin.f64 phi)))))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi))