(*.f64 esp (-.f64 angle phi0))
#s(approx (* esp (- angle phi0)) (*.f64 angle esp))
(if.f64 (<=.f64 (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))) #s(literal 5 binary64)) (atan2.f64 g (+.f64 (/.f64 (-.f64 (cos.f64 (-.f64 (/.f64 y esp) (-.f64 (-.f64 phi0 (*.f64 (PI.f64) #s(literal -1/2 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))))) (cos.f64 (+.f64 (/.f64 y esp) (-.f64 (-.f64 phi0 (*.f64 (PI.f64) #s(literal -1/2 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64)))))) #s(literal 2 binary64)) (*.f64 (cos.f64 phi0) (cos.f64 (/.f64 y esp))))) (atan2.f64 g (sin.f64 #s(approx (- (/ (neg y) esp) (- phi0 (* 1/2 (PI)))) (*.f64 #s(literal -1 binary64) phi0)))))
(if.f64 (<=.f64 (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))) #s(literal 5 binary64)) (atan2.f64 g (/.f64 (-.f64 (*.f64 (*.f64 (cos.f64 (/.f64 y esp)) (cos.f64 phi0)) #s(literal 2 binary64)) #s(approx (- (cos (- (/ y esp) phi0)) (cos (+ (/ y esp) phi0))) (-.f64 (cos.f64 (neg.f64 phi0)) (cos.f64 phi0)))) #s(literal 2 binary64))) (atan2.f64 g (sin.f64 #s(approx (- (/ (neg y) esp) (- phi0 (* 1/2 (PI)))) (*.f64 #s(literal -1 binary64) phi0)))))
(if.f64 (<=.f64 (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))) #s(literal 5 binary64)) (atan2.f64 g (-.f64 (*.f64 (cos.f64 (/.f64 y esp)) (cos.f64 phi0)) (*.f64 (sin.f64 (/.f64 y esp)) (sin.f64 phi0)))) (atan2.f64 g (sin.f64 #s(approx (- (/ (neg y) esp) (- phi0 (* 1/2 (PI)))) (*.f64 #s(literal -1 binary64) phi0)))))
(if.f64 (<=.f64 (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))) #s(literal 5 binary64)) (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))) (atan2.f64 g (sin.f64 #s(approx (- (/ (neg y) esp) (- phi0 (* 1/2 (PI)))) (*.f64 #s(literal -1 binary64) phi0)))))
(if.f64 (<=.f64 (/.f64 y esp) #s(literal 19999999999999999781741223642818392253569612520802717890360030929450604798220516297708225612915260122593317856641907797168065523046908674225209344 binary64)) (atan2.f64 g (cos.f64 #s(approx (+ phi0 (/ y esp)) phi0))) (atan2.f64 g (sin.f64 #s(approx (- (/ (neg y) esp) (- phi0 (* 1/2 (PI)))) (*.f64 #s(literal -1 binary64) phi0)))))
(atan2.f64 g (cos.f64 #s(approx (+ phi0 (/ y esp)) phi0)))
(atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64))))
(if.f64 (<=.f64 g #s(literal 4000000 binary64)) (asin.f64 (/.f64 (fabs.f64 (+.f64 (*.f64 (sin.f64 (/.f64 y esp)) (cos.f64 phi0)) (*.f64 (cos.f64 (/.f64 y esp)) (sin.f64 phi0)))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g)))))) (asin.f64 #s(approx (/ (fabs (+ (* (sin (/ y esp)) (cos phi0)) (* (cos (/ y esp)) (sin phi0)))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (+.f64 (*.f64 (cos.f64 phi0) (sin.f64 (/.f64 y esp))) (*.f64 (cos.f64 (/.f64 y esp)) (sin.f64 phi0)))) g))))
(if.f64 (<=.f64 g #s(literal 499999999999999990417798086218687295286560007015159396545582405077050056101839291488149134308110575981351030133088002720283516165604201974116686757888 binary64)) (asin.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (/.f64 #s(literal 1 binary64) (/.f64 esp y)) phi0))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g)))))) (asin.f64 #s(approx (/ (fabs (+ (* (sin (/ y esp)) (cos phi0)) (* (cos (/ y esp)) (sin phi0)))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (+.f64 (*.f64 (cos.f64 phi0) (sin.f64 (/.f64 y esp))) (*.f64 (cos.f64 (/.f64 y esp)) (sin.f64 phi0)))) g))))
(if.f64 (<=.f64 g #s(literal 9000000 binary64)) (asin.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (/.f64 #s(literal 1 binary64) (/.f64 esp y)) phi0))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g)))))) (asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)))) g))))
(if.f64 (<=.f64 g #s(literal 4000000 binary64)) (asin.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (/.f64 y esp) phi0))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g)))))) (asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)))) g))))
(if.f64 (<=.f64 g #s(literal 13 binary64)) (asin.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (/.f64 y esp) phi0))) (sqrt.f64 (fabs.f64 #s(approx (- -1 (* g g)) #s(literal -1 binary64)))))) (asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)))) g))))
(if.f64 (<=.f64 g #s(literal 13 binary64)) (asin.f64 (/.f64 #s(literal 1 binary64) #s(approx (/ (sqrt (fabs (- (* g g) -1))) (fabs (sin (+ (/ y esp) phi0)))) (/.f64 #s(literal 1 binary64) (fabs.f64 (sin.f64 (+.f64 phi0 (/.f64 y esp)))))))) (asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)))) g))))
(if.f64 (<=.f64 (+.f64 #s(literal 1 binary64) (*.f64 g g)) #s(literal 5000000000000 binary64)) (asin.f64 (/.f64 (fabs.f64 #s(approx (sin (+ (/ y esp) phi0)) (sin.f64 phi0))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g)))))) (asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)))) g))))
(asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)))) g)))
(asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 (+.f64 phi0 (/.f64 y esp)))) g)))
(asin.f64 #s(approx (/ (fabs (sin (+ (/ y esp) phi0))) (sqrt (fabs (- -1 (* g g))))) (/.f64 (fabs.f64 (sin.f64 #s(approx (+ phi0 (/ y esp)) phi0))) g)))
(sinh.f64 (/.f64 x esp))
(if.f64 (<=.f64 (/.f64 x esp) #s(literal 50000000000000000 binary64)) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 #s(approx (- 1 (/ (* -1/6 (* x x)) (* esp esp))) #s(literal 1 binary64)) (/.f64 x esp))) (if.f64 (<=.f64 (/.f64 x esp) #s(literal 200000000000000008675459394923837214658058664990387862358354757867223362577936222188264751104 binary64)) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 x #s(approx (+ (* 1/6 (/ (pow x 2) (pow esp 3))) (/ 1 esp)) (*.f64 esp (/.f64 esp (*.f64 (*.f64 esp esp) esp)))))) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 (-.f64 #s(literal 1 binary64) (/.f64 (*.f64 #s(literal -1/6 binary64) (*.f64 x x)) (*.f64 esp esp))) (/.f64 x esp)))))
#s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 x (+.f64 (*.f64 #s(literal 1/6 binary64) (/.f64 (/.f64 (*.f64 (/.f64 x esp) x) (neg.f64 esp)) (neg.f64 esp))) (/.f64 #s(literal 1 binary64) esp))))
(if.f64 (<=.f64 (*.f64 #s(literal 1/2 binary64) (-.f64 (exp.f64 (/.f64 x esp)) (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 x esp))))) #s(literal 6189700196426901/309485009821345068724781056 binary64)) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 #s(approx (- 1 (/ (* -1/6 (* x x)) (* esp esp))) #s(literal 1 binary64)) (/.f64 x esp))) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 x (+.f64 (*.f64 #s(literal 1/6 binary64) (*.f64 x (/.f64 x (*.f64 (*.f64 esp esp) esp)))) (/.f64 #s(literal 1 binary64) esp)))))
(if.f64 (<=.f64 (*.f64 #s(literal 1/2 binary64) (-.f64 (exp.f64 (/.f64 x esp)) (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 x esp))))) #s(literal 0 binary64)) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 #s(approx (- 1 (/ (* -1/6 (* x x)) (* esp esp))) #s(literal 1 binary64)) (/.f64 x esp))) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 (/.f64 x esp) (-.f64 (*.f64 #s(literal 1/6 binary64) (*.f64 x (/.f64 x (*.f64 esp esp)))) #s(literal -1 binary64)))))
#s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 x (+.f64 (*.f64 #s(literal 1/6 binary64) (*.f64 (/.f64 (/.f64 x esp) esp) (/.f64 x esp))) (/.f64 #s(literal 1 binary64) esp))))
#s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 (-.f64 #s(literal 1 binary64) (/.f64 (*.f64 (*.f64 (/.f64 x esp) x) #s(literal -1/6 binary64)) esp)) (/.f64 x esp)))
(if.f64 (<=.f64 (*.f64 #s(literal 1/2 binary64) (-.f64 (exp.f64 (/.f64 x esp)) (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 x esp))))) #s(literal 7737125245533627/77371252455336267181195264 binary64)) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 #s(approx (- 1 (/ (* -1/6 (* x x)) (* esp esp))) #s(literal 1 binary64)) (/.f64 x esp))) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 (-.f64 (*.f64 esp esp) (*.f64 #s(literal -1/6 binary64) (*.f64 x x))) (/.f64 x (*.f64 (*.f64 esp esp) esp)))))
(if.f64 (<=.f64 (/.f64 x esp) #s(literal 50000000000000000 binary64)) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 #s(approx (- 1 (/ (* -1/6 (* x x)) (* esp esp))) #s(literal 1 binary64)) (/.f64 x esp))) #s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 x #s(approx (+ (* 1/6 (/ (pow x 2) (pow esp 3))) (/ 1 esp)) (*.f64 esp (/.f64 esp (*.f64 (*.f64 esp esp) esp)))))))
#s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 #s(approx (- 1 (/ (* -1/6 (* x x)) (* esp esp))) #s(literal 1 binary64)) (/.f64 x esp)))
#s(approx (* 1/2 (- (exp (/ x esp)) (/ 1 (exp (/ x esp))))) (*.f64 x #s(approx (+ (* 1/6 (/ (pow x 2) (pow esp 3))) (/ 1 esp)) (/.f64 #s(literal 1 binary64) esp))))
(*.f64 #s(literal 2 binary64) (+.f64 (/.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 3 binary64)) (-.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)) (*.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (atan.f64 (exp.f64 (/.f64 y kR)))) (atan.f64 (exp.f64 (/.f64 y kR)))))) (/.f64 (*.f64 #s(literal -1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (-.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)) (*.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (atan.f64 (exp.f64 (/.f64 y kR)))) (atan.f64 (exp.f64 (/.f64 y kR))))))))
(*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 y) kR)))) (*.f64 (PI.f64) #s(literal 1/4 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 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 #s(literal 1 binary64) #s(approx (exp (/ (neg y) kR)) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1 binary64) (/.f64 y kR)))))) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) (*.f64 (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (-.f64 #s(literal 1 binary64) (/.f64 (-.f64 (*.f64 (*.f64 y (/.f64 y kR)) #s(literal -1/2 binary64)) 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 -1152921504606847/576460752303423488 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 #s(literal 1 binary64) #s(approx (exp (/ (neg y) kR)) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1 binary64) (/.f64 y kR)))))) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(approx (+ 1 (/ y kR)) (*.f64 y (+.f64 (/.f64 #s(literal 1 binary64) kR) (/.f64 #s(literal 1 binary64) y)))))) (/.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 -1152921504606847/576460752303423488 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 #s(literal 1 binary64) #s(approx (exp (/ (neg y) kR)) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1 binary64) (/.f64 y kR)))))) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) (*.f64 (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (-.f64 (/.f64 y kR) #s(literal -1 binary64)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)))
(*.f64 (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (-.f64 (/.f64 y kR) #s(literal -1 binary64)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(*.f64 (-.f64 (atan.f64 #s(approx (exp (/ y kR)) #s(literal 1 binary64))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(-.f64 phip delp)
(asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 (*.f64 (cos.f64 phipp) (cos.f64 lampp)))))
(if.f64 (<=.f64 phipp #s(literal -1770887431076117/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 132000000 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 (cos.f64 lampp) #s(literal 9007199254740979/9007199254740992 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 #s(approx (* (cos phipp) (cos lampp)) (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 sinp0 #s(literal -4943758783237121/4052261297735344686047273304385899561535592023674254785152009111026028136145418111718463914987406049109568248643848426935932764722081811824108276205189417663145685354884286644224 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (* sinp0 (+ (* (cos lampp) (cos phipp)) (* (/ (sin phipp) sinp0) cosp0))) (*.f64 sinp0 (*.f64 (cos.f64 lampp) (cos.f64 phipp)))))) (if.f64 (<=.f64 sinp0 #s(literal 4370447970442389/390218568789499028922699653724145788218574767833121393857394619953171467352470702515038262882936496394978366390175827861930996959911035663286272 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 sinp0 #s(approx (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0)) (+.f64 #s(approx (cos phipp) #s(literal 1 binary64)) (/.f64 (*.f64 cosp0 (sin.f64 phipp)) sinp0)))))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (* sinp0 (+ (* (cos lampp) (cos phipp)) (* (/ (sin phipp) sinp0) cosp0))) (*.f64 sinp0 (*.f64 (cos.f64 lampp) (cos.f64 phipp))))))))
(if.f64 (<=.f64 phipp #s(literal -1770887431076117/36893488147419103232 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 sinp0 #s(approx (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0)) (+.f64 #s(approx (cos phipp) #s(literal 1 binary64)) (/.f64 (*.f64 cosp0 (sin.f64 phipp)) sinp0)))))) (if.f64 (<=.f64 phipp #s(literal 132000000 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 sinp0 #s(approx (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0)) (+.f64 #s(approx (cos phipp) #s(literal 1 binary64)) (/.f64 (*.f64 cosp0 (sin.f64 phipp)) sinp0))))))))
(if.f64 (<=.f64 phipp #s(literal -1770887431076117/36893488147419103232 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 sinp0 #s(approx (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0)) (+.f64 #s(approx (cos phipp) #s(literal 1 binary64)) (*.f64 (/.f64 (sin.f64 phipp) sinp0) cosp0)))))) (if.f64 (<=.f64 phipp #s(literal 132000000 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 sinp0 #s(approx (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0)) (+.f64 #s(approx (cos phipp) #s(literal 1 binary64)) (*.f64 (/.f64 (sin.f64 phipp) sinp0) cosp0))))))))
(if.f64 (<=.f64 phipp #s(literal -280 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 cosp0 (sin.f64 phipp)))) (if.f64 (<=.f64 phipp #s(literal 132000000 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 cosp0 (sin.f64 phipp))))))
(if.f64 (<=.f64 sinp0 #s(literal -4943758783237121/4052261297735344686047273304385899561535592023674254785152009111026028136145418111718463914987406049109568248643848426935932764722081811824108276205189417663145685354884286644224 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))))) (if.f64 (<=.f64 sinp0 #s(literal 4370447970442389/390218568789499028922699653724145788218574767833121393857394619953171467352470702515038262882936496394978366390175827861930996959911035663286272 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) (*.f64 cosp0 (sin.f64 phipp)))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (*.f64 sinp0 (cos.f64 lampp)))))))
(if.f64 (<=.f64 lampp #s(literal 149999999999999998557159424 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))))) (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)))))
(if.f64 (<=.f64 cosp0 #s(literal -569999999999999993189233259358191616 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (*.f64 phipp #s(approx (+ cosp0 (/ (* sinp0 (cos lampp)) phipp)) cosp0))))) (if.f64 (<=.f64 cosp0 #s(literal 2455322750715917/11692013098647223345629478661730264157247460343808 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) #s(approx (+ sinp0 (* cosp0 phipp)) sinp0)))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (*.f64 phipp #s(approx (+ cosp0 (/ (* sinp0 (cos lampp)) phipp)) cosp0)))))))
(asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) #s(approx (+ sinp0 (* cosp0 phipp)) sinp0))))
(*.f64 ecc (sin.f64 phip))
#s(approx (* ecc (sin phip)) (*.f64 (+.f64 (*.f64 (-.f64 (*.f64 (-.f64 (*.f64 #s(literal 1/120 binary64) ecc) (*.f64 #s(literal 1/5040 binary64) (*.f64 (*.f64 phip ecc) phip))) (*.f64 phip phip)) (*.f64 #s(literal 1/6 binary64) ecc)) (*.f64 phip phip)) ecc) phip))
#s(approx (* ecc (sin phip)) (*.f64 (+.f64 (*.f64 (-.f64 (*.f64 #s(approx (- (* 1/120 ecc) (* 1/5040 (* (* phip ecc) phip))) (*.f64 #s(literal 1/120 binary64) ecc)) (*.f64 phip phip)) (*.f64 #s(literal 1/6 binary64) ecc)) (*.f64 phip phip)) ecc) phip))
(*.f64 #s(approx (sin phip) (*.f64 (-.f64 (*.f64 (*.f64 phip phip) #s(literal -1/6 binary64)) #s(literal -1 binary64)) phip)) ecc)
#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 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64))) es) #s(literal -1 binary64)))
#s(approx (sqrt (+ 1 (* es (* (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) rone_es)))) (-.f64 (sqrt.f64 (-.f64 (*.f64 rone_es es) #s(literal -1 binary64))) (/.f64 (*.f64 (*.f64 (*.f64 phi0 phi0) rone_es) es) (sqrt.f64 (-.f64 (*.f64 rone_es es) #s(literal -1 binary64))))))
(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))))
(if.f64 (<=.f64 hlf_e #s(literal -3399999999999999970641029280647652253504989982395107886759936 binary64)) #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))))))) (if.f64 (<=.f64 hlf_e #s(literal 3500000000000000040423454640501799690528344422718440438020279986078058282437574656 binary64)) (-.f64 (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (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 (*.f64 ecc hlf_e) (*.f64 #s(literal 2 binary64) (sin.f64 phi0))))))))
(if.f64 (<=.f64 ecc #s(literal 2152829579878715/299004108316488185848867822351174292268202623451041623610437960853046897358187881202051502002495252055893102400215083077513294952620060704653486412249315029041282722978874323768631243063073425063936 binary64)) #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 (-.f64 (+.f64 (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))) (/.f64 (log.f64 (/.f64 (sin.f64 (+.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal 1/2 binary64) phip0))) (cos.f64 (+.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal 1/2 binary64) phip0))))) c)) (log.f64 (/.f64 (sin.f64 (+.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal 1/2 binary64) phi0))) (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))))))) #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 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))
(if.f64 (<=.f64 ecc #s(literal 7415301886248907/1196016433265952743395471289404697169072810493804166494441751843412187589432751524808206008009981008223572409600860332310053179810480242818613945648997260116165130891915497295074524972252293700255744 binary64)) (-.f64 (log.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (PI.f64))) (PI.f64)) (PI.f64)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (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 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))
(if.f64 (<=.f64 (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))))))))) #s(literal -4835703278458517/9671406556917033397649408 binary64)) (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (-.f64 (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 (-.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi0)) ecc)) hlf_e) (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0))))) 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 c #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 ecc (*.f64 hlf_e (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))
(if.f64 (<=.f64 (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))))))))) #s(literal -4835703278458517/9671406556917033397649408 binary64)) (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 ecc (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 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 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))
(if.f64 (<=.f64 (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))))))))) #s(literal -4835703278458517/9671406556917033397649408 binary64)) (-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (PI.f64))))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 ecc (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 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 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))
(if.f64 (<=.f64 (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))))))))) #s(literal -4835703278458517/9671406556917033397649408 binary64)) (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 ecc #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 c #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 ecc (*.f64 hlf_e (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))
(if.f64 (<=.f64 (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phi0)))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))))))))) #s(literal -4835703278458517/9671406556917033397649408 binary64)) (-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64))) (PI.f64)))) (*.f64 hlf_e #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 ecc #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 c #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 ecc (*.f64 hlf_e (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 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 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 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 (*.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))
(*.f64 kR #s(approx (log (- (/ (* (cos (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 1/2)))) (/ (* (sin (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 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 phipp (+.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 (* phipp (+ (* 1/6 (* phipp (+ (* -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 phipp #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))))))))))))))
(*.f64 kR #s(approx (log (- (/ (* (cos (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 1/2)))) (/ (* (sin (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 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 phipp #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))) (* phipp (+ (* 1/6 (* phipp (+ (* -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)))))))))
(*.f64 kR #s(approx (log (- (/ (* (cos (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 1/2)))) (/ (* (sin (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 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 (*.f64 #s(literal 2 binary64) (*.f64 (*.f64 phipp (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 #s(literal -1/4 binary64) (-.f64 (/.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 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))) (sqrt.f64 #s(literal 2 binary64))))))
(*.f64 #s(approx (log (- (/ (* (cos (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 1/2)))) (/ (* (sin (* -1/2 phipp)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phipp 1/2)))))) (-.f64 (log.f64 (/.f64 (*.f64 #s(literal 1/2 binary64) (sqrt.f64 #s(literal 2 binary64))) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (*.f64 #s(literal -2 binary64) (*.f64 (*.f64 (*.f64 phipp (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 #s(literal -1/4 binary64) (-.f64 (/.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 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) kR)
(*.f64 (+.f64 (log.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))))) (log.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) phipp)) (cos.f64 (*.f64 phipp #s(literal 1/2 binary64))))))) kR)
(*.f64 kR (-.f64 (log.f64 (fabs.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 phipp #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal -1/2 binary64) phipp)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(*.f64 (log.f64 (*.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (-.f64 (cos.f64 (*.f64 phipp #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal -1/2 binary64) phipp))))) kR)
(*.f64 (log.f64 (*.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) #s(approx (- (cos (* phipp 1/2)) (sin (* -1/2 phipp))) (+.f64 #s(literal 1 binary64) (*.f64 phipp (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/8 binary64) phipp))))))) kR)
(if.f64 (<=.f64 phipp #s(literal -29811744551187/709803441694928604052074031140629428079727891296209043243642772637343054798240159498233447962659731992932150006119314388217384402944 binary64)) (*.f64 kR (log.f64 (tan.f64 (/.f64 (-.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) phipp) phipp) (*.f64 #s(literal 1/16 binary64) (*.f64 (PI.f64) (PI.f64)))) (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.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)))))))
(if.f64 (<=.f64 phipp #s(literal -29811744551187/709803441694928604052074031140629428079727891296209043243642772637343054798240159498233447962659731992932150006119314388217384402944 binary64)) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (*.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 (log.f64 (*.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) #s(approx (- (cos (* phipp 1/2)) (sin (* -1/2 phipp))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) phipp))))) kR)
(*.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)
(if.f64 (<=.f64 phi #s(literal 2102928824402889/18446744073709551616 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)))))))) (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 phi #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (*.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) phi)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))) (log.f64 (fabs.f64 (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))))))
(if.f64 (<=.f64 phi #s(literal 8500259669165361/18889465931478580854784 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 (sin.f64 (*.f64 (-.f64 (*.f64 phi #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))) (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 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 (+.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)))))))))))))
(if.f64 (<=.f64 phi #s(literal 7131979884014207/129672361527531029953512745740348785969138944757576153124864291552832900356653379574990845279596993571506183956603149661949848471106617978371464838566061365220661931356297172615168 binary64)) (log.f64 (tan.f64 (+.f64 (+.f64 #s(approx (+ (* phi 1/2) (+ (* 17/4 (PI)) (PI))) (+.f64 (PI.f64) (*.f64 #s(literal 17/4 binary64) (PI.f64)))) (PI.f64)) (PI.f64)))) (log.f64 (tan.f64 (+.f64 (*.f64 #s(literal 17/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 binary64))))))
(-.f64 (+.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (log.f64 (fabs.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) phi)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi))))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))
(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))))))
(log.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (*.f64 phi #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 17/4 binary64) (PI.f64)) (PI.f64))) (PI.f64)) (PI.f64))))
(log.f64 (tan.f64 (+.f64 (*.f64 #s(literal 17/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 binary64)))))
(log.f64 (tan.f64 (+.f64 #s(approx (+ (* phi 1/2) (+ (* 13/4 (PI)) (PI))) (+.f64 (PI.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 sinp0 (*.f64 (cos.f64 phip) (cos.f64 lamp)))))
(if.f64 (<=.f64 phip #s(literal -731966804844795/4722366482869645213696 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 230000000000000 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)))))))
(if.f64 (<=.f64 (cos.f64 lamp) #s(literal 9007192838293625/9007199254740992 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 #s(approx (* (cos phip) (cos lamp)) (cos.f64 lamp))))) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 sinp0 #s(approx (- (/ (* cosp0 (sin phip)) sinp0) (* (cos lamp) (cos phip))) (-.f64 (/.f64 (*.f64 cosp0 (sin.f64 phip)) sinp0) (cos.f64 phip)))))))
(if.f64 (<=.f64 (cos.f64 lamp) #s(literal 9007192838293625/9007199254740992 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 cosp0 #s(literal -4455508415646675/278469275977917188637766821636980671685377716506870112950390168328813181592326500641230029987632531724770648777685205524334027387779262573135538410627576549575927409323746640019706676772864 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (if.f64 (<=.f64 cosp0 #s(literal 6992716752707823/24973988402527937851052777838345330445988785141319769206873255677002973910558124960962448824507935769278615448971252983163583805434306282450321408 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 (*.f64 (neg.f64 (cos.f64 phip)) (cos.f64 lamp)) sinp0))) (if.f64 (<=.f64 cosp0 #s(literal 4183246980753157/2535301200456458802993406410752 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp))))))))
(if.f64 (<=.f64 cosp0 #s(literal -4455508415646675/278469275977917188637766821636980671685377716506870112950390168328813181592326500641230029987632531724770648777685205524334027387779262573135538410627576549575927409323746640019706676772864 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (if.f64 (<=.f64 cosp0 #s(literal 6992716752707823/24973988402527937851052777838345330445988785141319769206873255677002973910558124960962448824507935769278615448971252983163583805434306282450321408 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 (*.f64 (neg.f64 sinp0) (cos.f64 lamp)) (cos.f64 phip)))) (if.f64 (<=.f64 cosp0 #s(literal 4183246980753157/2535301200456458802993406410752 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp))))))))
(if.f64 (<=.f64 phip #s(literal -731966804844795/4722366482869645213696 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (if.f64 (<=.f64 phip #s(literal 16500000000 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))) #s(approx (* sinp0 (cos phip)) sinp0))))))
(if.f64 (<=.f64 cosp0 #s(literal -3331240651246993/9797766621314684873895700802803279209044463565243731922466831101232640732633100491228823617617764419367505179450247842283955649007454149170085442756585554871624752266571753841250508572690789992495054848 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (if.f64 (<=.f64 cosp0 #s(literal 4926807918903085/117304950450073441093299338992332138457996243649210992760592177980666118165925495436678284352817677825758656549761022853424729541660972550965022826666248518125620733165482852770884817967017897067499683873717228533661411547573197142488509591988118713532416 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) #s(approx (* -1 (* sinp0 (* (cos lamp) (cos phip)))) (*.f64 #s(literal -1 binary64) (*.f64 sinp0 (cos.f64 phip)))))) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0))))))
(if.f64 (<=.f64 sinp0 #s(literal -317043563211239/31082702275611665134711390509176302506278509424834232340028998555822468563283335970816 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 (neg.f64 (cos.f64 lamp)) sinp0))) (if.f64 (<=.f64 sinp0 #s(literal 3794275180128377/3794275180128377091639574036764685364535950857523710002444946112771297432041422848 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) #s(approx (* sinp0 (cos phip)) sinp0)))) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 (neg.f64 (cos.f64 lamp)) sinp0)))))
(if.f64 (<=.f64 phip #s(literal -1558380939346983/9444732965739290427392 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 cosp0 (sin.f64 phip)))) (if.f64 (<=.f64 phip #s(literal 220 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) #s(approx (- (* cosp0 phip) (* sinp0 (cos lamp))) (-.f64 (*.f64 cosp0 phip) sinp0)))) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 cosp0 (sin.f64 phip))))))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) #s(approx (- (* cosp0 phip) (* sinp0 (cos lamp))) (-.f64 (*.f64 cosp0 phip) 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 5596146821100995/2945340432158418383223693624588738123559693482299075088767878449688292160397327779966295692450325070170031945807812908771881611572255401942922812303597144053805349165872996110766935565946816006053119311086960734516644260779498911850068592403100913453684334767056261910363295677456051671938422478104563288264146944 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (-.f64 K (*.f64 (-.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e) (log.f64 (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/4 binary64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))) (/.f64 (-.f64 (pow.f64 (*.f64 (atan.f64 (exp.f64 (-.f64 K (*.f64 (-.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e) (log.f64 (tan.f64 (+.f64 (+.f64 (+.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)) (PI.f64))))) c)))) #s(literal 2 binary64)) #s(literal 2 binary64)) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (*.f64 (atan.f64 (exp.f64 (-.f64 K (*.f64 (-.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e) (log.f64 (tan.f64 (+.f64 (+.f64 (+.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)) (PI.f64))))) c)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))
(if.f64 (<=.f64 phi #s(literal 5596146821100995/2945340432158418383223693624588738123559693482299075088767878449688292160397327779966295692450325070170031945807812908771881611572255401942922812303597144053805349165872996110766935565946816006053119311086960734516644260779498911850068592403100913453684334767056261910363295677456051671938422478104563288264146944 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (-.f64 K (*.f64 (-.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e) (log.f64 (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/4 binary64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))) (-.f64 (*.f64 #s(literal 2 binary64) (atan.f64 (exp.f64 (+.f64 (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (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 K #s(literal -3332663724254167/4503599627370496 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 (-.f64 K (*.f64 (-.f64 (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e) (log.f64 (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/4 binary64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))
(if.f64 (<=.f64 K #s(literal 50000000000000 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 (fabs (- -1 (* (sin phi) ecc)))) (log (fabs (- (* (sin phi) ecc) 1)))) (+.f64 (*.f64 (*.f64 (*.f64 (*.f64 ecc ecc) #s(literal 2/3 binary64)) (pow.f64 (sin.f64 phi) #s(literal 3 binary64))) ecc) (*.f64 (*.f64 #s(literal 2 binary64) (sin.f64 phi)) ecc)))
#s(approx (- (log (fabs (- -1 (* (sin phi) ecc)))) (log (fabs (- (* (sin phi) ecc) 1)))) (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 ecc ecc) #s(literal 2/3 binary64)) (pow.f64 (sin.f64 phi) #s(literal 3 binary64))) (*.f64 #s(literal 2 binary64) (sin.f64 phi))) ecc))
#s(approx (- (log (fabs (- -1 (* (sin phi) ecc)))) (log (fabs (- (* (sin phi) ecc) 1)))) (*.f64 #s(approx (+ (+ (* (* (pow (sin phi) 3) 2/3) (* ecc ecc)) (sin phi)) (sin phi)) (*.f64 #s(literal 2 binary64) (sin.f64 phi))) ecc))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi))
(*.f64 (sqrt.f64 (fabs.f64 (/.f64 (-.f64 #s(literal 1 binary64) Up) (-.f64 #s(literal -1 binary64) Up)))) Esc)
(*.f64 #s(approx (sqrt (fabs (/ (- 1 Up) (- -1 Up)))) (+.f64 #s(literal 1 binary64) (*.f64 Up (-.f64 (*.f64 Up (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) Up))) #s(literal 1 binary64))))) Esc)
#s(approx (/ Esc (sqrt (/ (+ 1 Up) (- 1 Up)))) (-.f64 Esc (*.f64 (+.f64 Esc (*.f64 (*.f64 #s(literal -1/2 binary64) Esc) Up)) Up)))
#s(approx (/ Esc (sqrt (/ (+ 1 Up) (- 1 Up)))) (*.f64 (-.f64 #s(literal 1 binary64) Up) Esc))
#s(approx (/ Esc (sqrt (/ (+ 1 Up) (- 1 Up)))) Esc)
(if.f64 (<=.f64 (*.f64 BrA u) #s(literal 2000000000000000137172103703564102993414188346625929973381646791516038639747754425505775838752679231688970493666459275394749789597812172229456459932366192699143082941239010020801269538891555886778514937042106442934926263917068257100320412740354048 binary64)) (neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (/.f64 #s(approx (- (pow (* 0 (sin (* u BrA))) 2) (pow (* 1 (cos (* u BrA))) 2)) (*.f64 #s(literal -1 binary64) (pow.f64 (cos.f64 (*.f64 BrA u)) #s(literal 2 binary64)))) (sin.f64 (+.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 u BrA)) (PI.f64))))))) (neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Vp (-.f64 (*.f64 cosgam (/.f64 Sp Vp)) singam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 (*.f64 BrA u) #s(literal 2000000000000000137172103703564102993414188346625929973381646791516038639747754425505775838752679231688970493666459275394749789597812172229456459932366192699143082941239010020801269538891555886778514937042106442934926263917068257100320412740354048 binary64)) (neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (sin.f64 (-.f64 (*.f64 u BrA) (*.f64 (PI.f64) #s(literal -1/2 binary64))))))) (neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Vp (-.f64 (*.f64 cosgam (/.f64 Sp Vp)) singam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (cos.f64 (*.f64 BrA u)))) #s(literal 49999999999999995601277750478615906956426432484762865091230684279338790788450641385479969549606017377053487170299935055586674081792 binary64)) (neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (cos.f64 (*.f64 BrA u))))) (neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Vp #s(approx (- (/ (* Sp cosgam) Vp) singam) (*.f64 singam (-.f64 (/.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) #s(literal 1 binary64)))))) #s(approx (cos (* BrA u)) #s(literal 1 binary64))))))
(neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) #s(approx (cos (* BrA u)) #s(literal 1 binary64)))))
(if.f64 (<=.f64 (*.f64 Sp cosgam) #s(literal -4742843975160471/118571099379011784113736688648896417641748464297615937576404566024103044751294464 binary64)) (neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Sp #s(approx (+ cosgam (* -1 (/ (* Vp singam) Sp))) cosgam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64))))) (if.f64 (<=.f64 (*.f64 Sp cosgam) #s(literal 200000000000 binary64)) (*.f64 (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 (neg.f64 singam) Vp)) #s(approx (cos (* BrA u)) #s(literal 1 binary64))) (neg.f64 rB)) (neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Sp #s(approx (+ cosgam (* -1 (/ (* Vp singam) Sp))) cosgam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64)))))))
(neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Sp #s(approx (+ cosgam (* -1 (/ (* Vp singam) Sp))) cosgam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64)))))
(pow.f64 (/.f64 #s(literal 1 binary64) (*.f64 (E.f64) (E.f64))) (*.f64 (*.f64 v BrA) #s(literal 1/2 binary64)))
(pow.f64 (/.f64 #s(literal 1 binary64) (E.f64)) (*.f64 v BrA))
(pow.f64 (E.f64) (*.f64 (neg.f64 BrA) v))
(exp.f64 (neg.f64 (*.f64 BrA v)))
(if.f64 (<=.f64 (neg.f64 (*.f64 BrA v)) #s(literal -200 binary64)) #s(approx (exp (neg (* BrA v))) #s(approx (+ 1 (* -1 (* BrA v))) (*.f64 v #s(approx (+ (* -1 BrA) (/ 1 v)) (*.f64 (/.f64 #s(literal 1 binary64) (*.f64 v v)) v))))) (if.f64 (<=.f64 (neg.f64 (*.f64 BrA v)) #s(literal 3602879701896397/72057594037927936 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) (if.f64 (<=.f64 (neg.f64 (*.f64 BrA v)) #s(literal 1000000000000000007630473539575035660514778335511710750780086664439969510636494954611131549135839186513983455555395220895687860544809584999829725260594873271087399626486606146442550988840016917394626449536395208620267012778077787723395914064607119962069483324573977857832138825282954985472 binary64)) #s(approx (exp (neg (* BrA v))) (/.f64 (-.f64 #s(literal 1 binary64) (*.f64 BrA (*.f64 BrA (*.f64 v v)))) (+.f64 (*.f64 v BrA) #s(literal 1 binary64)))) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))))))
(if.f64 (<=.f64 (exp.f64 (neg.f64 (*.f64 BrA v))) #s(literal 6156563468186637/1231312693637327475383720003129487931408741852202045208373384168882678805359287831606695820465153613775207124697088 binary64)) #s(approx (exp (neg (* BrA v))) #s(approx (+ 1 (* -1 (* BrA v))) (*.f64 v #s(approx (+ (* -1 BrA) (/ 1 v)) (*.f64 (/.f64 #s(literal 1 binary64) (*.f64 v v)) v))))) (if.f64 (<=.f64 (exp.f64 (neg.f64 (*.f64 BrA v))) #s(literal 2 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) #s(approx (exp (neg (* BrA v))) (/.f64 (/.f64 (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 v BrA)) (*.f64 v v)) (neg.f64 v)) (neg.f64 v)))))
(if.f64 (<=.f64 (neg.f64 (*.f64 BrA v)) #s(literal -200 binary64)) #s(approx (exp (neg (* BrA v))) #s(approx (+ 1 (* -1 (* BrA v))) (*.f64 v #s(approx (+ (* -1 BrA) (/ 1 v)) (*.f64 (/.f64 #s(literal 1 binary64) (*.f64 v v)) v))))) #s(approx (exp (neg (* BrA v))) (/.f64 (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 v BrA)) v) v)))
#s(approx (exp (neg (* BrA v))) (/.f64 (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 v BrA)) v) v))
#s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA)))
#s(approx (exp (neg (* BrA v))) #s(literal 1 binary64))
(if.f64 (<=.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4503577109372359/4503599627370496 binary64)) (*.f64 ArB (log.f64 (-.f64 (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64)) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64)) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))))) (if.f64 (<=.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4548635623644201/4503599627370496 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (/ (* (sin (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 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 gamma0 #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))) (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -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))))))))) (*.f64 ArB (log.f64 (-.f64 (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64)) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64)) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))))))))
(if.f64 (<=.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4278419646001971/4503599627370496 binary64)) (*.f64 ArB (log.f64 (-.f64 (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64)) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64)) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))))) (*.f64 ArB #s(approx (log (- (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (/ (* (sin (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 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 gamma0 (+.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 (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -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 gamma0 #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)))))))))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5902958103587057/1180591620717411303424 binary64)) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 4722366482869645/9444732965739290427392 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (/ (* (sin (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 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 gamma0 #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))) (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -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))))))))) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 gamma0 #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) gamma0)))) #s(literal -5902958103587057/1180591620717411303424 binary64)) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 4722366482869645/9444732965739290427392 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (/ (* (sin (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 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)))))) #s(approx (* gamma0 (+ (* 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))) (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -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 gamma0 (*.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)))))))) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))))
(*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(*.f64 ArB (log.f64 (*.f64 (-.f64 #s(approx (* (* (/ (* (sqrt 2) 1/2) (cos (- (* gamma0 1/2) (* -1/4 (PI))))) (cos (* gamma0 1/2))) (cos (- (* gamma0 1/2) (* -1/4 (PI))))) (*.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (*.f64 #s(literal 1/2 binary64) gamma0)) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sqrt.f64 #s(literal 2 binary64))) #s(literal 1/2 binary64))) (/.f64 #s(literal 1 binary64) (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(*.f64 (log.f64 (/.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) ArB)
(*.f64 ArB (log.f64 (*.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 7737125245533627/154742504910672534362390528 binary64)) (*.f64 ArB (log.f64 (fabs.f64 #s(approx (tan (+ (/ (PI) 4) (* 1/2 gamma0))) (-.f64 (*.f64 #s(literal 1 binary64) gamma0) #s(literal -1 binary64)))))) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -2028240960365167/20282409603651670423947251286016 binary64)) (*.f64 ArB (log.f64 (fabs.f64 #s(approx (tan (+ (/ (PI) 4) (* 1/2 gamma0))) (-.f64 (*.f64 #s(literal 1 binary64) gamma0) #s(literal -1 binary64)))))) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 gamma0 #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) gamma0)))) #s(literal 944473296573929/18889465931478580854784 binary64)) (*.f64 ArB (log.f64 (fabs.f64 #s(approx (tan (+ (/ (PI) 4) (* 1/2 gamma0))) (-.f64 (*.f64 #s(literal 1 binary64) gamma0) #s(literal -1 binary64)))))) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0)))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 944473296573929/18889465931478580854784 binary64)) (*.f64 ArB (log.f64 (fabs.f64 #s(approx (tan (+ (/ (PI) 4) (* 1/2 gamma0))) (-.f64 (*.f64 #s(literal 1 binary64) gamma0) #s(literal -1 binary64)))))) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64))))))))
(*.f64 ArB (log.f64 (fabs.f64 #s(approx (tan (+ (/ (PI) 4) (* 1/2 gamma0))) (-.f64 (*.f64 #s(literal 1 binary64) gamma0) #s(literal -1 binary64))))))
(*.f64 (log.f64 #s(approx (tan (+ (/ (PI) 4) (* 1/2 gamma0))) (-.f64 gamma0 #s(literal -1 binary64)))) ArB)
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5764607523034235/576460752303423488 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 5902958103587057/295147905179352825856 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (/ (sqrt 2) 2) (cos (* -1/2 gamma0))) (cos (* (+ (* gamma0 -2) (PI)) -1/4))) (/ (* (sin (* gamma0 1/2)) (/ (sqrt 2) 2)) (cos (* (+ (* gamma0 -2) (PI)) -1/4))))) (+.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 gamma0 (+.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 (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -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))))))))) (- (* 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)))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4))) (+ (* -1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (* 1/2 (- (* -1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* -1/16 (/ (sqrt 2) (cos (* -1/4 (PI)))))))))) (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 gamma0 #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/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)))))))) (-.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)))))))) (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/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 (*.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/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 #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/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/16 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))))) (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)))))))))))))) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))) (*.f64 (*.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))))) (log.f64 (fabs.f64 (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5764607523034235/576460752303423488 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 6189700196426901/618970019642690137449562112 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (/ (sqrt 2) 2) (cos (* -1/2 gamma0))) (cos (* (+ (* gamma0 -2) (PI)) -1/4))) (/ (* (sin (* gamma0 1/2)) (/ (sqrt 2) 2)) (cos (* (+ (* gamma0 -2) (PI)) -1/4))))) (+.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 gamma0 #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))) (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -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))))))))) (- (* 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)))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4))) (+ (* -1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (* 1/2 (- (* -1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* -1/16 (/ (sqrt 2) (cos (* -1/4 (PI)))))))))) (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))))))))) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5764607523034235/576460752303423488 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 6189700196426901/618970019642690137449562112 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (/ (sqrt 2) 2) (cos (* -1/2 gamma0))) (cos (* (+ (* gamma0 -2) (PI)) -1/4))) (/ (* (sin (* gamma0 1/2)) (/ (sqrt 2) 2)) (cos (* (+ (* gamma0 -2) (PI)) -1/4))))) (+.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 gamma0 #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))) (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -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))))))))) (- (* 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)))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4))) (+ (* -1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (* 1/2 (- (* -1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* -1/16 (/ (sqrt 2) (cos (* -1/4 (PI)))))))))) (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))))))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (fabs.f64 (/.f64 (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) ArB)))
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5764607523034235/576460752303423488 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 6189700196426901/618970019642690137449562112 binary64)) (*.f64 ArB #s(approx (log (- (/ (* (/ (sqrt 2) 2) (cos (* -1/2 gamma0))) (cos (* (+ (* gamma0 -2) (PI)) -1/4))) (/ (* (sin (* gamma0 1/2)) (/ (sqrt 2) 2)) (cos (* (+ (* gamma0 -2) (PI)) -1/4))))) (+.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)))))) #s(approx (* gamma0 (+ (* 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))) (* gamma0 (+ (* 1/6 (* gamma0 (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -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))))))))) (- (* 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)))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4))) (+ (* -1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (* 1/2 (- (* -1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* -1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* -1/16 (/ (sqrt 2) (cos (* -1/4 (PI)))))))))) (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 gamma0 (*.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)))))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (fabs.f64 (/.f64 (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) ArB)))
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5764607523034235/576460752303423488 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (fabs.f64 (/.f64 (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) ArB))
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -5764607523034235/288230376151711744 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))) (sin.f64 (+.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)) (/.f64 (PI.f64) #s(literal 2 binary64)))))))) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 2206763817411543/2251799813685248 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))) #s(approx (cos (* (+ (* -2 gamma0) (PI)) 1/4)) (cos.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)))))))) (*.f64 ArB (log.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (sin.f64 (+.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64))))))))
(*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))))))
(*.f64 (log.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))) ArB)
(if.f64 (<=.f64 (*.f64 ArB (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -2024022533073/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* gamma0 1/2))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) gamma0)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))))) (log.f64 (fabs.f64 (cos.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal 6189700196426901/618970019642690137449562112 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* gamma0 1/2))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) gamma0)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (neg.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64))))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))))
(if.f64 (<=.f64 (*.f64 ArB (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -2024022533073/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* gamma0 1/2))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) gamma0)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))))
(if.f64 (<=.f64 gamma0 #s(literal -4951760157141521/39614081257132168796771975168 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) (if.f64 (<=.f64 gamma0 #s(literal 8741742960768907/12855504354071922204335696738729300820177623950262342682411008 binary64)) (*.f64 #s(approx ArB #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (- (/ (PI) 4) (* 1/2 gamma0)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))))))))
(if.f64 (<=.f64 gamma0 #s(literal 8741742960768907/12855504354071922204335696738729300820177623950262342682411008 binary64)) (*.f64 #s(approx ArB #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (- (/ (PI) 4) (* 1/2 gamma0)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) ArB))
(*.f64 #s(approx ArB #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (- (/ (PI) 4) (* 1/2 gamma0)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(*.f64 ArB (log.f64 (tan.f64 #s(approx (- (/ (PI) 4) (* 1/2 gamma0)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(-.f64 (*.f64 (*.f64 (-.f64 (cos.f64 (+.f64 phi0 phi0)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) es) #s(literal -1 binary64))
(-.f64 #s(literal 1 binary64) (*.f64 #s(approx (* (sin phi0) (sin phi0)) (*.f64 phi0 phi0)) es))
#s(approx (- 1 (* es (* (sin phi0) (sin phi0)))) #s(literal 1 binary64))
(/.f64 (*.f64 (sqrt.f64 (/.f64 #s(literal 1 binary64) (/.f64 one_es (+.f64 one_es (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es))))) com) (*.f64 (cos.f64 phi0) (sqrt.f64 (-.f64 #s(literal 1 binary64) (*.f64 es (*.f64 (sin.f64 phi0) (sin.f64 phi0)))))))
(*.f64 (/.f64 com (cos.f64 phi0)) (sqrt.f64 #s(approx (/ (- (/ (* (pow (cos phi0) 4) es) one_es) -1) (- 1 (* (* (- 1 (cos (+ phi0 phi0))) 1/2) es))) (+.f64 #s(literal 1 binary64) (/.f64 #s(literal 1 binary64) (/.f64 one_es (+.f64 (*.f64 (+.f64 (*.f64 (-.f64 (/.f64 es one_es) #s(literal -1 binary64)) es) (*.f64 (/.f64 es one_es) #s(literal -2 binary64))) (*.f64 (*.f64 phi0 phi0) one_es)) es)))))))
(/.f64 (*.f64 (sqrt.f64 #s(approx (/ (- (/ (* (pow (cos phi0) 4) es) one_es) -1) (- 1 (* (* (- 1 (cos (+ phi0 phi0))) 1/2) es))) (+.f64 (*.f64 (+.f64 (*.f64 (-.f64 (/.f64 es one_es) #s(literal -1 binary64)) es) (*.f64 (/.f64 es one_es) #s(literal -2 binary64))) (*.f64 phi0 phi0)) (-.f64 (/.f64 es one_es) #s(literal -1 binary64))))) com) (cos.f64 phi0))
#s(approx (/ (* (sqrt (+ 1 (/ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es))) com) (* (cos phi0) (sqrt (- 1 (* es (* (sin phi0) (sin phi0))))))) (*.f64 com (sqrt.f64 (/.f64 #s(literal 1 binary64) (/.f64 one_es (+.f64 (*.f64 #s(literal 1 binary64) one_es) es))))))
#s(approx (/ (* (sqrt (+ 1 (/ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es))) com) (* (cos phi0) (sqrt (- 1 (* es (* (sin phi0) (sin phi0))))))) (*.f64 com (sqrt.f64 (+.f64 #s(literal 1 binary64) (/.f64 es one_es)))))
#s(approx (/ (* (sqrt (+ 1 (/ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es))) com) (* (cos phi0) (sqrt (- 1 (* es (* (sin phi0) (sin phi0))))))) (*.f64 com #s(approx (sqrt (+ 1 (/ es one_es))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) (/.f64 es one_es))))))
#s(approx (/ (* (sqrt (+ 1 (/ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es))) com) (* (cos phi0) (sqrt (- 1 (* es (* (sin phi0) (sin phi0))))))) #s(approx (* com (sqrt (+ 1 (/ es one_es)))) (+.f64 com (*.f64 #s(literal 1/2 binary64) (/.f64 (*.f64 com es) one_es)))))
#s(approx (/ (* (sqrt (+ 1 (/ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es))) com) (* (cos phi0) (sqrt (- 1 (* es (* (sin phi0) (sin phi0))))))) #s(approx (/ com (cos phi0)) com))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (/.f64 (*.f64 es #s(approx (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -2 binary64) (pow.f64 phi0 #s(literal 2 binary64)))))) one_es)))
(sqrt.f64 (/.f64 (+.f64 (*.f64 #s(approx (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) (-.f64 #s(literal 1 binary64) (*.f64 (+.f64 phi0 phi0) phi0))) es) one_es) one_es))
(sqrt.f64 (/.f64 #s(literal 1 binary64) (/.f64 one_es (+.f64 #s(approx (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) es) one_es))))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (/.f64 #s(approx (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) es) one_es)))
(/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 #s(approx (/ one_es (+ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es)) #s(literal 1 binary64)))))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (*.f64 (-.f64 (*.f64 (-.f64 (*.f64 (-.f64 (*.f64 #s(literal -1/7 binary64) (*.f64 (*.f64 U ArB) U)) (*.f64 #s(literal 1/5 binary64) ArB)) (*.f64 U U)) (*.f64 #s(literal 1/3 binary64) ArB)) (*.f64 U U)) ArB) U))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (-.f64 (*.f64 (*.f64 (*.f64 U U) U) (-.f64 (*.f64 (*.f64 (*.f64 U U) ArB) #s(literal -1/5 binary64)) (*.f64 #s(literal 1/3 binary64) ArB))) (*.f64 U ArB)))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (*.f64 (-.f64 (*.f64 (-.f64 (*.f64 #s(literal -1/5 binary64) (*.f64 (*.f64 U ArB) U)) (*.f64 #s(literal 1/3 binary64) ArB)) (*.f64 U U)) ArB) U))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (*.f64 (-.f64 (*.f64 #s(approx (- (* -1/5 (* (* U ArB) U)) (* 1/3 ArB)) (*.f64 #s(literal -1/3 binary64) ArB)) (*.f64 U U)) ArB) U))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (*.f64 (neg.f64 U) ArB))
(*.f64 A lam)
(if.f64 (<=.f64 (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam))) #s(literal 5 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (cos.f64 #s(approx (- (* (neg B) lam) (* (PI) -1/2)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))))
(if.f64 (<=.f64 (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam))) #s(literal 5 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (cos.f64 #s(approx (- (* (neg B) lam) (* (PI) -1/2)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))))
(if.f64 (<=.f64 (*.f64 B lam) #s(literal -1000000000000000000000 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (cos.f64 #s(approx (- (* (neg B) lam) (* (PI) -1/2)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))) (if.f64 (<=.f64 (*.f64 B lam) #s(literal 499999999999999978409886320820907879202552238629189140897698107811441303810555744076971465473716161022374445056 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) #s(approx (* (sin (* B lam)) singam) (*.f64 (*.f64 lam B) singam))) #s(approx (cos (* B lam)) #s(literal 1 binary64)))) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (cos.f64 #s(approx (- (* (neg B) lam) (* (PI) -1/2)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64))))))
(*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) #s(approx (* (sin (* B lam)) singam) (*.f64 B (*.f64 lam singam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64))))
(-.f64 (*.f64 (-.f64 u u0) cosrot) (*.f64 v sinrot))
(if.f64 (<=.f64 (*.f64 v sinrot) #s(literal -1000000000000000 binary64)) (-.f64 #s(approx (* (- u u0) cosrot) (*.f64 cosrot u)) (*.f64 v sinrot)) (if.f64 (<=.f64 (*.f64 v sinrot) #s(literal 213598703592091/4271974071841820164790043412339104229205409044713305539894083215644439451561281100045924173873152 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot (-.f64 u u0))) (-.f64 #s(approx (* (- u u0) cosrot) (*.f64 cosrot u)) (*.f64 v sinrot))))
(if.f64 (<=.f64 (*.f64 v sinrot) #s(literal -19999999999999999098582133569958947190600450174767048236959251965035770900582349244308780304596114601737544754773898621832134656 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 (neg.f64 v) sinrot)) (if.f64 (<=.f64 (*.f64 v sinrot) #s(literal 9999999999999999932209486743616279764617084419440640 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot (-.f64 u u0))) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 (neg.f64 v) sinrot))))
(if.f64 (<=.f64 (*.f64 v sinrot) #s(literal -2894802230932905/57896044618658097711785492504343953926634992332820282019728792003956564819968 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 (neg.f64 v) sinrot)) (if.f64 (<=.f64 (*.f64 v sinrot) #s(literal 213598703592091/4271974071841820164790043412339104229205409044713305539894083215644439451561281100045924173873152 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot u)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 (neg.f64 v) sinrot))))
#s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot u))
(+.f64 (*.f64 v cosrot) (*.f64 (-.f64 u u0) sinrot))
(if.f64 (<=.f64 u #s(literal -19000000000000000820569080407792879310496361049729294918786304442556826584315303900423442216236123079491081822021877760 binary64)) (+.f64 (*.f64 v cosrot) #s(approx (* (- u u0) sinrot) (*.f64 sinrot u))) (if.f64 (<=.f64 u #s(literal 1000000000000000050555427725995033814228237030803003279020481474722232763977085405824233377105062219252417113236701184 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (-.f64 (*.f64 cosrot v) (*.f64 sinrot u0))) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0)))))
(if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal -99999999999999997704951326524533662844684271992415000612999597473199345218078991130326129448151154688 binary64)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0))) (if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal 4000000000000000070142166264077602166149767460708800344583192419745366289222052773513135086095057460819601312121498137127444044423471504910343963196864201302268135998643045972226530032988244005617664 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (-.f64 (*.f64 cosrot v) (*.f64 sinrot u0))) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0)))))
(if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal -1000000000000000000 binary64)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0))) (if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal 4999999999999999966104743371808139882308542209720320 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 cosrot v)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0)))))
(if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal -99999999999999997704951326524533662844684271992415000612999597473199345218078991130326129448151154688 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 #s(approx (- u u0) u) sinrot)) (if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal 200000000000000003111883225893368604853640279384212286673954116086166756232951140652997077983009489535241256173568000 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 cosrot v)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 #s(approx (- u u0) (neg.f64 u0)) sinrot))))
(if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal -99999999999999997704951326524533662844684271992415000612999597473199345218078991130326129448151154688 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 #s(approx (- u u0) u) sinrot)) (if.f64 (<=.f64 (*.f64 (-.f64 u u0) sinrot) #s(literal 4999999999999999772223133475743019061733712700409539130496607211544840259226135691611880105565210303017104153796972357853870064153456670293082673807209411155434429495479368482882719667688996710696271289138913738752 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 cosrot v)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 #s(approx (- u u0) u) sinrot))))
#s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (*.f64 cosrot v))
(sqrt.f64 (+.f64 (pow.f64 (sinh.f64 Ce) #s(literal 2 binary64)) #s(literal 1 binary64)))
(if.f64 (<=.f64 Ce #s(literal 164999999999999990088603704975509451249150407477010779866365775715099928952832 binary64)) #s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) (/.f64 (-.f64 (sqrt.f64 (*.f64 (*.f64 (*.f64 (*.f64 (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce) Ce) Ce) (*.f64 (*.f64 (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce) Ce) Ce)) #s(literal 1/64 binary64))) #s(literal -1 binary64)) (-.f64 (*.f64 (*.f64 Ce Ce) (+.f64 (*.f64 #s(literal 1/4 binary64) (*.f64 Ce Ce)) #s(literal -1/2 binary64))) #s(literal -1 binary64)))) #s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) (+.f64 #s(literal 1 binary64) (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) (*.f64 Ce Ce)) Ce) Ce)))))
#s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) (+.f64 #s(literal 1 binary64) (sqrt.f64 (sqrt.f64 (*.f64 (*.f64 #s(literal 1/16 binary64) (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce)) (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce))))))
#s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) (+.f64 #s(literal 1 binary64) (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) (*.f64 Ce Ce)) Ce) Ce))))
#s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) (-.f64 (*.f64 (*.f64 Ce Ce) #s(literal 1/2 binary64)) #s(literal -1 binary64)))
#s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) #s(approx (- (* (* Ce Ce) 1/2) -1) #s(literal 1 binary64)))
(exp.f64 (+.f64 Ce Ce))
(if.f64 (<=.f64 Ce #s(literal -48999163945791/144115188075855872 binary64)) #s(approx (exp (* 2 Ce)) (/.f64 #s(approx (- (* (* (* (- (+ Ce Ce) -2) Ce) Ce) (- (+ Ce Ce) -2)) 1) #s(literal -1 binary64)) (-.f64 (*.f64 (-.f64 (+.f64 Ce Ce) #s(literal -2 binary64)) Ce) #s(literal 1 binary64)))) #s(approx (exp (* 2 Ce)) (/.f64 (-.f64 (*.f64 (*.f64 (*.f64 (-.f64 (+.f64 Ce Ce) #s(literal -2 binary64)) Ce) Ce) (-.f64 (+.f64 Ce Ce) #s(literal -2 binary64))) #s(literal 1 binary64)) #s(approx (- (* (- (+ Ce Ce) -2) Ce) 1) (-.f64 (*.f64 #s(literal 2 binary64) Ce) #s(literal 1 binary64))))))
(if.f64 (<=.f64 Ce #s(literal -48999163945791/144115188075855872 binary64)) #s(approx (exp (* 2 Ce)) (/.f64 #s(approx (- (* (* (* (- (+ Ce Ce) -2) Ce) Ce) (- (+ Ce Ce) -2)) 1) #s(literal -1 binary64)) (-.f64 (*.f64 (-.f64 (+.f64 Ce Ce) #s(literal -2 binary64)) Ce) #s(literal 1 binary64)))) #s(approx (exp (* 2 Ce)) (+.f64 #s(literal 1 binary64) (*.f64 Ce (+.f64 #s(literal 2 binary64) (*.f64 Ce (+.f64 #s(literal 2 binary64) (*.f64 #s(literal 4/3 binary64) Ce))))))))
#s(approx (exp (* 2 Ce)) (+.f64 #s(literal 1 binary64) (*.f64 Ce (+.f64 #s(literal 2 binary64) (*.f64 Ce (+.f64 #s(literal 2 binary64) (*.f64 #s(literal 4/3 binary64) Ce)))))))
#s(approx (exp (* 2 Ce)) (-.f64 (*.f64 (-.f64 (+.f64 Ce Ce) #s(literal -2 binary64)) Ce) #s(literal -1 binary64)))
#s(approx (exp (* 2 Ce)) (-.f64 (+.f64 Ce Ce) #s(literal -1 binary64)))
#s(approx (exp (* 2 Ce)) (+.f64 (+.f64 #s(literal 1 binary64) Ce) Ce))
#s(approx (exp (* 2 Ce)) (+.f64 #s(approx (+ 1 Ce) #s(literal 1 binary64)) Ce))
#s(approx (exp (* 2 Ce)) #s(literal 1 binary64))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (*.f64 (sin.f64 Cn) (sin.f64 Cn)) (*.f64 (*.f64 (cos.f64 Cn) (cos.f64 Cn)) (*.f64 (cos.f64 lam) (cos.f64 lam))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 lam lam)))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn))))))))
(/.f64 #s(literal 2 binary64) #s(approx (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam))))) (fabs.f64 (cos.f64 lam))))
(/.f64 #s(literal 2 binary64) #s(approx (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam))))) (fabs.f64 #s(approx (cos Cn) #s(literal 1 binary64)))))
(/.f64 (*.f64 (sin.f64 lam) (cos.f64 Cn)) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64))) (pow.f64 (*.f64 (cos.f64 lam) (cos.f64 Cn)) #s(literal 2 binary64)))))
(/.f64 (*.f64 (sin.f64 lam) (cos.f64 Cn)) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn))) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 lam lam)))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn)))))))))
(*.f64 (cos.f64 Cn) (/.f64 (sin.f64 lam) (sqrt.f64 (+.f64 (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 lam lam)))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn))))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn))))))))
(/.f64 (*.f64 (sin.f64 lam) (cos.f64 Cn)) #s(approx (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam))))) (sqrt.f64 (pow.f64 (cos.f64 lam) #s(literal 2 binary64)))))
(*.f64 (cos.f64 Cn) (/.f64 (sin.f64 lam) #s(approx (sqrt (+ (* (- 1/2 (* -1/2 (cos (+ lam lam)))) (- 1/2 (* -1/2 (cos (+ Cn Cn))))) (- 1/2 (* 1/2 (cos (+ Cn Cn)))))) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) lam))))))))
#s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (/.f64 (sin.f64 lam) (fabs.f64 (cos.f64 lam))))
#s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (/.f64 (sin.f64 lam) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) #s(approx (* 1/2 (cos (* 2 (- Cn (* -1/2 (PI)))))) #s(literal -1/2 binary64))))))
#s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) #s(approx (/ (sin lam) (sqrt (pow (cos lam) 2))) (-.f64 lam (*.f64 (*.f64 #s(literal -1/3 binary64) (*.f64 lam lam)) lam))))
#s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) #s(approx (/ (sin lam) (sqrt (pow (cos lam) 2))) (*.f64 lam (+.f64 #s(literal 1 binary64) (*.f64 (*.f64 #s(literal 1/3 binary64) lam) lam)))))
(atan2.f64 (sin.f64 Cn) (*.f64 (cos.f64 Cn) (cos.f64 lam)))
(atan2.f64 (sin.f64 Cn) #s(approx (* (cos Cn) (cos lam)) (cos.f64 lam)))
(atan2.f64 #s(approx (sin Cn) (-.f64 Cn (*.f64 (*.f64 #s(literal 1/6 binary64) (*.f64 Cn Cn)) Cn))) #s(approx (* (cos Cn) (cos lam)) (cos.f64 lam)))
(*.f64 (*.f64 #s(literal 2 binary64) cos_arg_r) cosh_arg_i)
(*.f64 (+.f64 cosh_arg_i cosh_arg_i) cos_arg_r)
(if.f64 (<=.f64 (*.f64 (*.f64 #s(literal 2 binary64) cos_arg_r) cosh_arg_i) #s(literal 8379879956214123/1675975991242824637446753124775730765934920727574049172215445180465220503759193372100234287270862928461253982273310756356719235351493321243304206125760512 binary64)) (+.f64 #s(approx cos_arg_r #s(literal 0 binary64)) #s(approx cos_arg_r #s(literal 0 binary64))) (+.f64 cosh_arg_i cosh_arg_i))
(+.f64 cos_arg_r cos_arg_r)
(*.f64 (*.f64 #s(literal -2 binary64) sin_arg_r) sinh_arg_i)
(if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 60000000000000001896764932096 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 sin_arg_r (+.f64 (*.f64 #s(literal -1 binary64) (/.f64 (*.f64 cos_arg_r (*.f64 hi sinh_arg_i)) sin_arg_r)) (*.f64 cosh_arg_i hr)))))
(if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 200000000000000002881189517449054771166223724485662526027424629870997854138252263253726515251452912161010874366592467075072 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 (-.f64 (*.f64 hr cosh_arg_i) (*.f64 (/.f64 (*.f64 hi sinh_arg_i) sin_arg_r) cos_arg_r)) sin_arg_r)))
(if.f64 (<=.f64 (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) #s(literal 2000000000000000132872935482496206237094312341172584908970922214753713493255768101167089780693375139608812241567134921336075488584322101781755750774742240399521541760156078250259598945212267909879768657149226586411367871939134697181462712041438531269934236247503274787037183937480902858990682112 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 (-.f64 (*.f64 hr cosh_arg_i) (*.f64 (/.f64 (*.f64 hi cos_arg_r) sin_arg_r) sinh_arg_i)) sin_arg_r)))
(if.f64 (<=.f64 hi #s(literal 290000000000000006817283078234158399161900034285931781579881244597883946948870697187951259765402912698399484979538317518878193459495448315479821738931480425208217072926982144 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 hi (-.f64 (*.f64 sin_arg_r (/.f64 (*.f64 hr cosh_arg_i) hi)) (*.f64 cos_arg_r sinh_arg_i)))))
(if.f64 (<=.f64 hi #s(literal 4290498537581631/214524926879081553593184399971293538039669853129478294043576983099548224481176751628829988770670454843040509730983776813660062124991145119142938384097345536 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)))
(if.f64 (<=.f64 sin_arg_r #s(literal -8183476519740355/818347651974035467503297424206899788054160511510766197370822842024033449101168638720817523081476039287721671031890017752304314136471348263332131897344 binary64)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) (-.f64 (*.f64 (*.f64 hr sin_arg_r) cosh_arg_i) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)))
(-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi))
(if.f64 (<=.f64 cosh_arg_i #s(literal 30000000000000000000000 binary64)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 cosh_arg_i (*.f64 hr sin_arg_r)))))
(if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal -5521397077432451/55213970774324510299478046898216203619608871777363092441300193790394368 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr))) (if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal 1000000000000000008493621433689702976148869924598760615894999102702796905906176 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 (*.f64 (neg.f64 sinh_arg_i) cos_arg_r) hi)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr)))))
(if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal -5521397077432451/55213970774324510299478046898216203619608871777363092441300193790394368 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr))) (if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal 1000000000000000008493621433689702976148869924598760615894999102702796905906176 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 (*.f64 (neg.f64 cos_arg_r) hi) sinh_arg_i)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr)))))
(if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal -5521397077432451/55213970774324510299478046898216203619608871777363092441300193790394368 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr))) (if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal 1000000000000000008493621433689702976148869924598760615894999102702796905906176 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 (*.f64 hi sinh_arg_i) (neg.f64 cos_arg_r))) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr)))))
#s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 (*.f64 cosh_arg_i sin_arg_r) hr)))
#s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) #s(approx (* sin_arg_r (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr))) (*.f64 cosh_arg_i (*.f64 hr sin_arg_r))))
(if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal -200000000000000018380567016286756476168068919431369064448 binary64)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) (if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal 5902958103587057/295147905179352825856 binary64)) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (*.f64 sinh_arg_i (+.f64 (*.f64 cos_arg_r hr) (/.f64 (*.f64 cosh_arg_i (*.f64 hi sin_arg_r)) sinh_arg_i)))) (if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal +inf.0 binary64)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (/.f64 #s(approx (* (+ (* hr (* sinh_arg_i cos_arg_r)) (* hi (* cosh_arg_i sin_arg_r))) hr) (*.f64 (*.f64 (*.f64 hr sin_arg_r) hi) cosh_arg_i)) hr)))))
(if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal -3999999999999999869344675216466765095398132743226221891671847117885183383686911450434959473821876224 binary64)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) (if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal 5902958103587057/295147905179352825856 binary64)) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (*.f64 sinh_arg_i (+.f64 (*.f64 cos_arg_r hr) (*.f64 (*.f64 hi sin_arg_r) (/.f64 cosh_arg_i sinh_arg_i))))) (if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal +inf.0 binary64)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (/.f64 #s(approx (* (+ (* hr (* sinh_arg_i cos_arg_r)) (* hi (* cosh_arg_i sin_arg_r))) hr) (*.f64 (*.f64 (*.f64 hr sin_arg_r) hi) cosh_arg_i)) hr)))))
(if.f64 (<=.f64 (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(literal +inf.0 binary64)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (/.f64 #s(approx (* (+ (* hr (* sinh_arg_i cos_arg_r)) (* hi (* cosh_arg_i sin_arg_r))) hr) (*.f64 (*.f64 (*.f64 hr sin_arg_r) hi) cosh_arg_i)) hr)))
(if.f64 (<=.f64 cos_arg_r #s(literal -199999999999999983222784 binary64)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 hr cos_arg_r) sinh_arg_i)))
(if.f64 (<=.f64 sin_arg_r #s(literal -2000000000000000115715919885453939654786757378350080876345294848 binary64)) (+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) (+.f64 (*.f64 (*.f64 hi sin_arg_r) cosh_arg_i) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)))
(+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr))
(+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr))
#s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (*.f64 (*.f64 hr cos_arg_r) sinh_arg_i))
#s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (*.f64 cos_arg_r (*.f64 hr sinh_arg_i)))
(if.f64 (>.f32 #s(literal 1 binary32) #s(literal 1/10000000000 binary32)) (tan.f64 phi) #s(literal 0 binary64))
(*.f64 (/.f64 lam (sqrt.f64 (-.f64 (*.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi phi))) #s(literal 1/2 binary64)) es) #s(literal -1 binary64)))) (cos.f64 phi))
(*.f64 (/.f64 lam #s(approx (sqrt (- (* (- (* 1/2 (cos (+ phi phi))) 1/2) es) -1)) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) (*.f64 es (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) phi))) #s(literal 1/2 binary64))))))) (cos.f64 phi))
(*.f64 (cos.f64 phi) #s(approx (/ lam (sqrt (- (* (- (* 1/2 (cos (+ phi phi))) 1/2) es) -1))) (-.f64 lam (*.f64 (*.f64 (*.f64 es lam) #s(literal 1/2 binary64)) (-.f64 (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)) #s(literal 1/2 binary64))))))
(*.f64 #s(approx (/ lam (sqrt (- (* (- (* 1/2 (cos (+ phi phi))) 1/2) es) -1))) #s(approx (+ lam (* -1/2 (* es (* lam (- (* 1/2 (cos (* 2 phi))) 1/2))))) lam)) (cos.f64 phi))
#s(approx (/ (* (cos phi) lam) (sqrt (- 1 (* es (* (sin phi) (sin phi)))))) (+.f64 lam (*.f64 (*.f64 #s(approx (- (* -1/2 lam) (* -1/2 (* es lam))) (*.f64 (*.f64 es lam) #s(literal 1/2 binary64))) #s(approx phi #s(literal 0 binary64))) #s(approx phi #s(literal 0 binary64)))))
(/.f64 #s(literal 1 binary64) #s(approx (/ (sqrt (- (* (- (* 1/2 (cos (+ phi phi))) 1/2) es) -1)) (* lam (cos phi))) (/.f64 #s(literal 1 binary64) lam)))
#s(approx (/ (* (cos phi) lam) (sqrt (- 1 (* es (* (sin phi) (sin phi)))))) (+.f64 lam (*.f64 #s(approx (* (- (* -1/2 lam) (* -1/2 (* es lam))) phi) (*.f64 #s(literal -1/2 binary64) (*.f64 lam phi))) phi)))
(-.f64 #s(literal 1 binary64) (*.f64 es (*.f64 (sin.f64 phi) (sin.f64 phi))))
(-.f64 (*.f64 (*.f64 (-.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) es) #s(literal -1 binary64))
#s(approx (- 1 (* es (* (sin phi) (sin phi)))) #s(literal 1 binary64))
(/.f64 (*.f64 (sqrt.f64 (-.f64 (*.f64 (*.f64 (-.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) es) #s(literal -1 binary64))) x) k0)
(*.f64 (/.f64 x k0) #s(approx (sqrt (- (* (* (- (cos (+ phi phi)) 1) 1/2) es) -1)) (-.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) es) (-.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1 binary64))) #s(literal -1 binary64))))
(/.f64 #s(approx (* x (sqrt (- 1 (* es (* (sin phi) (sin phi)))))) x) k0)