(*.f64 esp (-.f64 angle phi0))
#s(approx (* esp (- angle phi0)) (*.f64 angle esp))
(if.f64 (<=.f64 (cos.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal 2 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 (+.f64 (*.f64 (sin.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64))) (sin.f64 #s(approx (- (* (PI) 1/2) (+ phi0 (/ y esp))) (*.f64 #s(literal -1 binary64) phi0)))) (*.f64 (cos.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64))) (sin.f64 #s(approx (+ phi0 (/ y esp)) phi0))))))
(if.f64 (<=.f64 (cos.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal 2 binary64)) (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 #s(literal 1 binary64) (/.f64 esp y))))) (atan2.f64 g (+.f64 (*.f64 (sin.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64))) (sin.f64 #s(approx (- (* (PI) 1/2) (+ phi0 (/ y esp))) (*.f64 #s(literal -1 binary64) phi0)))) (*.f64 (cos.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64))) (sin.f64 #s(approx (+ phi0 (/ y esp)) 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 #s(literal 1 binary64) (/.f64 esp y))))) (atan2.f64 g (sin.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) #s(approx (+ phi0 (/ y esp)) 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 (/.f64 #s(literal 1 binary64) esp) y)))) (atan2.f64 g (sin.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) #s(approx (+ phi0 (/ y esp)) phi0)))))
(if.f64 (<=.f64 (cos.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal -3602879701896397/18014398509481984 binary64)) (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))) (atan2.f64 g (cos.f64 #s(approx (+ phi0 (/ y esp)) phi0))))
(atan2.f64 g (cos.f64 #s(approx (+ phi0 (/ y esp)) phi0)))
(if.f64 (<=.f64 g #s(literal -7385903388887613/18014398509481984 binary64)) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64)))) (if.f64 (<=.f64 g #s(literal 8780918421921433/2787593149816327892691964784081045188247552 binary64)) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (+ (cos phi0) (* -1 (/ (* y (sin phi0)) esp))) (+.f64 #s(literal 1 binary64) (*.f64 phi0 (+.f64 (*.f64 #s(literal -1 binary64) (/.f64 y esp)) (*.f64 #s(literal -1/2 binary64) 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 -7385903388887613/18014398509481984 binary64)) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64)))) (if.f64 (<=.f64 g #s(literal 3356504041299755/799167628880894011233688890827050574271641124522232614619944181664095165137859998750798362384253944616915694367080095461234681773897801038410285056 binary64)) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (+ (cos phi0) (* -1 (/ (* y (sin phi0)) esp))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1 binary64) (/.f64 (*.f64 phi0 y) esp)))))) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 g #s(literal -7385903388887613/18014398509481984 binary64)) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64)))) (if.f64 (<=.f64 g #s(literal 3356504041299755/799167628880894011233688890827050574271641124522232614619944181664095165137859998750798362384253944616915694367080095461234681773897801038410285056 binary64)) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (+ (cos phi0) (* -1 (/ (* y (sin phi0)) esp))) (-.f64 #s(literal 1 binary64) (*.f64 (/.f64 y esp) phi0))))) (atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64))))))
(atan2.f64 g #s(approx (cos (+ phi0 (/ y esp))) #s(approx (cos (/ y esp)) #s(literal 1 binary64))))
(if.f64 (<=.f64 (+.f64 #s(literal 1 binary64) (*.f64 g g)) #s(literal 50000000000000004641673518601159954844517422622525385549225694063461714040984789960014820604544131271471563404911386848873613068925538235484773792943686604067961981752493137735453512646120016981018974140087018752579040234700800 binary64)) (asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))) (fabs.f64 (+.f64 (*.f64 (sin.f64 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)) (cos.f64 phi0)) (*.f64 (cos.f64 (*.f64 (/.f64 #s(literal 1 binary64) esp) y)) (sin.f64 phi0)))))) (asin.f64 #s(approx (sqrt (/ (pow (sin (+ (/ y esp) phi0)) 2) (+ 1 (* g g)))) (/.f64 (sqrt.f64 (pow.f64 (sin.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal 2 binary64))) g))))
(if.f64 (<=.f64 (+.f64 #s(literal 1 binary64) (*.f64 g g)) #s(literal 50000000000000000701959312789985260891230985285064568046915021472510652274325054054092066621782843422306142881889050953096494638431569844936383886042210844858380302841544704 binary64)) (asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))) (fabs.f64 (+.f64 (*.f64 (sin.f64 (/.f64 y esp)) (cos.f64 phi0)) (*.f64 (cos.f64 (/.f64 y esp)) (sin.f64 phi0)))))) (asin.f64 #s(approx (sqrt (/ (pow (sin (+ (/ y esp) phi0)) 2) (+ 1 (* g g)))) (/.f64 (sqrt.f64 (pow.f64 (sin.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal 2 binary64))) g))))
(if.f64 (<=.f64 g #s(literal 3800000000000000110810612650889027817164210252101892844867592768608318625454557010853888 binary64)) (asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))) (fabs.f64 (sin.f64 (+.f64 (*.f64 (/.f64 #s(literal 1 binary64) esp) y) phi0))))) (asin.f64 #s(approx (sqrt (/ (pow (sin (+ (/ y esp) phi0)) 2) (+ 1 (* g g)))) (/.f64 (sqrt.f64 (pow.f64 (sin.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal 2 binary64))) g))))
(if.f64 (<=.f64 g #s(literal 78999999999999996487205605809633643644636009416382715181688278112049243089992560308369520908837181858028073298178099137707861136933474594690282029056 binary64)) (asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))) (fabs.f64 (sin.f64 (+.f64 (*.f64 (/.f64 #s(literal 1 binary64) esp) y) phi0))))) (asin.f64 #s(approx (sqrt (/ (- 1/2 (* 1/2 (cos (* 2 (+ (/ y esp) phi0))))) (- (* g g) -1))) (/.f64 (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (+.f64 phi0 (/.f64 y esp))))))) g))))
(if.f64 (<=.f64 g #s(literal 46999999999999998917563339030654883119605688687423849686874405905413190203200343604726196596612913310699234154069380865516020118000250502183318148087808 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 (sqrt (/ (- 1/2 (* 1/2 (cos (* 2 (+ (/ y esp) phi0))))) (- (* g g) -1))) (/.f64 (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (+.f64 phi0 (/.f64 y esp))))))) g))))
(asin.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (/.f64 y esp) phi0))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))))
(if.f64 (<=.f64 (*.f64 g g) #s(literal 1730417938599789/23384026197294446691258957323460528314494920687616 binary64)) (asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 #s(approx (- -1 (* g g)) #s(literal -1 binary64))))) (fabs.f64 (sin.f64 (+.f64 (/.f64 y esp) phi0))))) (asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))) (fabs.f64 (sin.f64 #s(approx (+ (/ y esp) phi0) phi0))))))
(asin.f64 (*.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))) (fabs.f64 (sin.f64 #s(approx (+ (/ y esp) phi0) phi0)))))
(asin.f64 (*.f64 (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) #s(approx (cos (* 2 (+ (/ y esp) phi0))) #s(approx (+ (cos (* 2 phi0)) (* -2 (/ (* y (sin (* 2 phi0))) esp))) #s(approx (+ 1 (* -4 (/ (* phi0 y) esp))) (*.f64 y (+.f64 (*.f64 #s(literal -4 binary64) (/.f64 phi0 esp)) (/.f64 #s(literal 1 binary64) y))))))))) (/.f64 #s(literal 1 binary64) (sqrt.f64 (fabs.f64 (-.f64 (*.f64 g g) #s(literal -1 binary64)))))))
(asin.f64 (sqrt.f64 (/.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(approx (cos (* 2 (+ (/ y esp) phi0))) #s(approx (+ (cos (* 2 phi0)) (* -2 (/ (* y (sin (* 2 phi0))) esp))) (-.f64 (*.f64 (*.f64 phi0 (/.f64 y esp)) #s(literal -4 binary64)) #s(literal -1 binary64)))) #s(literal 1/2 binary64))) (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))))
(sinh.f64 (/.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 1 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp)) (*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (/.f64 (/.f64 (/.f64 (+.f64 x x) (*.f64 (*.f64 esp esp) esp)) esp) (*.f64 esp esp))) #s(literal 1/2 binary64)))
(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 1 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp)) (*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (/.f64 (/.f64 (+.f64 x x) (*.f64 (*.f64 esp esp) esp)) (*.f64 esp esp))) #s(literal 1/2 binary64)))
(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 1 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp)) (*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (/.f64 (/.f64 (+.f64 x x) (*.f64 esp esp)) (*.f64 esp esp))) #s(literal 1/2 binary64)))
(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 1 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp)) (*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (/.f64 (+.f64 x x) (*.f64 (*.f64 esp esp) esp))) #s(literal 1/2 binary64)))
(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 1 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp)) (*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (/.f64 (+.f64 x x) (*.f64 esp esp))) #s(literal 1/2 binary64)))
(if.f64 (<=.f64 (/.f64 x esp) #s(literal 200000 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp)) (if.f64 (<=.f64 (/.f64 x esp) #s(literal 499999999999999978930451175173142066076775939048257141926258886614516577002786239313118268535951812574041306445493431857101228510021003209840763187482937088894311771724997242528629131330872974013383816137806524948480039480659075272709232330533995834790894142764502740352844098034426819117056 binary64)) (*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (*.f64 esp (*.f64 (+.f64 x x) esp))) #s(literal 1/2 binary64)) #s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp))))
#s(approx (* (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) 1/2) (/.f64 x esp))
(*.f64 #s(approx (- (exp (/ x esp)) (/ 1 (exp (/ x esp)))) (+.f64 x x)) #s(literal 1/2 binary64))
(*.f64 #s(literal 2 binary64) (/.f64 (-.f64 (*.f64 (/.f64 (*.f64 #s(literal -1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))) (/.f64 (*.f64 #s(literal -1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64)))))) (*.f64 (/.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 3 binary64)) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))) (/.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 3 binary64)) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))))) (-.f64 (/.f64 (*.f64 #s(literal -1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))) (/.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 3 binary64)) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))))))
(*.f64 #s(literal 2 binary64) (+.f64 (/.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 3 binary64)) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64))))) (/.f64 (*.f64 #s(literal -1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (-.f64 (*.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (atan.f64 (exp.f64 (/.f64 y kR)))) (*.f64 #s(literal -1/16 binary64) (*.f64 (PI.f64) (PI.f64)))))))
(*.f64 #s(literal 2 binary64) (/.f64 (-.f64 (pow.f64 (atan.f64 (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 y) kR)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/16 binary64) (*.f64 (PI.f64) (PI.f64)))) (-.f64 (atan.f64 (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 y) kR)))) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))
(*.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 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (*.f64 y #s(approx (+ (* 1/2 (/ y (pow kR 2))) (/ 1 kR)) (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) (/.f64 y kR))) kR)))))) (/.f64 (PI.f64) #s(literal 4 binary64)))))
(if.f64 (<=.f64 (/.f64 y kR) #s(literal -4717453031026927/471745303102692669030194322623533611364615152552055053547025787530703020299018606510765840647207345942683546878845480612471323078740599079070179281429087005456174415872 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (*.f64 (-.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 (*.f64 (/.f64 y (*.f64 kR kR)) #s(literal 1/2 binary64)) y)) (/.f64 y kR)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal -1 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (/.f64 #s(literal 1 binary64) (/.f64 kR (+.f64 kR y))))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64))))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal -1 binary64)) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal -1/2 binary64))) (+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64))))
(+.f64 (*.f64 (atan.f64 #s(approx (exp (/ y kR)) (+.f64 #s(literal 1 binary64) (/.f64 y kR)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))
(*.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 (*.f64 (cos.f64 phipp) sinp0) (cos.f64 lampp))))
(asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 (*.f64 (cos.f64 phipp) (cos.f64 lampp)))))
(if.f64 (<=.f64 lampp #s(literal 4938706332868587/1496577676626844588240573268701473812127674924007424 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) (*.f64 sinp0 (cos.f64 phipp))))) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 #s(approx (* (cos phipp) (cos lampp)) (cos.f64 lampp))))))
(if.f64 (<=.f64 phipp #s(literal -1220 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 1357680363825023/1180591620717411303424 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 sinp0 #s(literal -5755292697852851/471745303102692669030194322623533611364615152552055053547025787530703020299018606510765840647207345942683546878845480612471323078740599079070179281429087005456174415872 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (* sinp0 (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0))) (*.f64 sinp0 (*.f64 (cos.f64 lampp) (cos.f64 phipp)))))) (if.f64 (<=.f64 sinp0 #s(literal 929962309538873/3957286423569672544968040336314355040562101585619389438625008897486355601704509874725062416979888719833562742816258229429621840548908019359480738481630290686665708354379186176 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) #s(approx (* sinp0 (cos phipp)) sinp0)))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (* sinp0 (+ (* (cos lampp) (cos phipp)) (/ (* cosp0 (sin phipp)) sinp0))) (*.f64 sinp0 (*.f64 (cos.f64 lampp) (cos.f64 phipp))))))))
(if.f64 (<=.f64 phipp #s(literal -1150 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) #s(approx (* sinp0 (cos phipp)) sinp0)))) (if.f64 (<=.f64 phipp #s(literal 1357680363825023/1180591620717411303424 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))) #s(approx (* sinp0 (cos phipp)) sinp0))))))
(if.f64 (<=.f64 cosp0 #s(literal -4263144513106179/1167984798111281975972139931059274579172666497855631342228273284582214442805421410945513679697247078343332431250840168271536308408672112127552681297848886832192510636636227827221215793215130566656 binary64)) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) #s(approx (* sinp0 (cos phipp)) sinp0)))) (if.f64 (<=.f64 cosp0 #s(literal 4883799633456925/813966605576154086191388158047028539690522256142078632860871314912602201882035164980643579207953404631705311351737456109117285358168867807661542486856790777111225290468019355281137486073219485440547487321275219968468108891293513341206528 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (*.f64 sinp0 (cos.f64 lampp))))) (asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) #s(approx (* sinp0 (* (cos phipp) (cos lampp))) #s(approx (* sinp0 (cos phipp)) sinp0))))))
(if.f64 (<=.f64 cosp0 #s(literal -4353782555323089/926336713898529563388567880069503262826159877325124512315660672063305037119488 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (+.f64 sinp0 (*.f64 cosp0 phipp))))) (if.f64 (<=.f64 cosp0 #s(literal 339999999999999991611392 binary64)) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (*.f64 sinp0 (cos.f64 lampp))))) (asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (+.f64 sinp0 (*.f64 cosp0 phipp)))))))
(asin.f64 #s(approx (+ (* cosp0 (sin phipp)) (* sinp0 (* (cos phipp) (cos lampp)))) #s(approx (+ (* cosp0 phipp) (* sinp0 (cos lampp))) (+.f64 sinp0 (*.f64 cosp0 phipp)))))
(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 ecc (+.f64 (*.f64 (*.f64 phip phip) #s(literal 1/120 binary64)) #s(literal -1/6 binary64))) (*.f64 phip phip)) ecc) phip))
#s(approx (* ecc (sin phip)) (-.f64 (*.f64 phip ecc) (*.f64 (*.f64 (*.f64 phip phip) phip) (*.f64 ecc #s(approx (- (* -1/120 (* phip phip)) -1/6) #s(literal 1/6 binary64))))))
#s(approx (* ecc (sin phip)) (*.f64 (+.f64 (*.f64 (*.f64 ecc #s(approx (+ (* (* phip phip) 1/120) -1/6) #s(literal -1/6 binary64))) (*.f64 phip phip)) ecc) phip))
#s(approx (* ecc (sin phip)) (*.f64 ecc phip))
(asin.f64 sinp0)
#s(approx (/ (/ (sqrt one_es) (- (* (- (* (cos (+ phi0 phi0)) 1/2) 1/2) (* ecc ecc)) -1)) (/ 1 k0)) (*.f64 k0 (sqrt.f64 one_es)))
(sqrt.f64 (-.f64 (*.f64 (*.f64 rone_es es) #s(approx (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) (-.f64 #s(literal 1 binary64) (*.f64 (+.f64 phi0 phi0) phi0)))) #s(literal -1 binary64)))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 es #s(approx (* (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) rone_es) rone_es))))
#s(approx (sqrt (+ 1 (* es (* (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0))) rone_es)))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) #s(approx (* es (* rone_es (pow (cos phi0) 4))) (*.f64 es rone_es)))))
(if.f64 (<=.f64 phip0 #s(literal 6592970859093183/42535295865117307932921825928971026432 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 (*.f64 ecc (+.f64 (*.f64 #s(literal 1/6 binary64) (*.f64 ecc (*.f64 hlf_e (+.f64 (*.f64 #s(literal -6 binary64) (*.f64 (sin.f64 phi0) (pow.f64 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0))) #s(literal 2 binary64)))) (+.f64 (*.f64 #s(literal 2 binary64) (pow.f64 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0))) #s(literal 3 binary64))) (*.f64 #s(literal 6 binary64) (*.f64 (pow.f64 (sin.f64 phi0) #s(literal 2 binary64)) (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))) (*.f64 #s(literal 1/2 binary64) (*.f64 hlf_e (+.f64 (*.f64 #s(literal -1 binary64) (pow.f64 (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0))) #s(literal 2 binary64))) (*.f64 #s(literal 2 binary64) (*.f64 (sin.f64 phi0) (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0)))))))))) (*.f64 hlf_e (-.f64 (sin.f64 phi0) (*.f64 #s(literal -1 binary64) (sin.f64 phi0))))))))) (-.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 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.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) phi0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.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))))))))))
(if.f64 (<=.f64 phip0 #s(literal 6592970859093183/42535295865117307932921825928971026432 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 #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 (+.f64 #s(approx (* ecc (+ (* 1/6 (* ecc (* hlf_e (+ (* -6 (* (sin phi0) (pow (- (sin phi0) (* -1 (sin phi0))) 2))) (+ (* 2 (pow (- (sin phi0) (* -1 (sin phi0))) 3)) (* 6 (* (pow (sin phi0) 2) (- (sin phi0) (* -1 (sin phi0)))))))))) (* 1/2 (* hlf_e (+ (* -1 (pow (- (sin phi0) (* -1 (sin phi0))) 2)) (* 2 (* (sin phi0) (- (sin phi0) (* -1 (sin phi0)))))))))) (*.f64 (*.f64 #s(literal 2/3 binary64) (*.f64 ecc ecc)) (*.f64 (*.f64 (*.f64 phi0 phi0) phi0) hlf_e))) (*.f64 (-.f64 (sin.f64 phi0) (neg.f64 (sin.f64 phi0))) hlf_e)) ecc)) c)) (-.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 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.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) phi0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.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))))))))))
#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 #s(approx (* hlf_e (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0)))))) (*.f64 (+.f64 #s(approx (* ecc (+ (* 1/6 (* ecc (* hlf_e (+ (* -6 (* (sin phi0) (pow (- (sin phi0) (* -1 (sin phi0))) 2))) (+ (* 2 (pow (- (sin phi0) (* -1 (sin phi0))) 3)) (* 6 (* (pow (sin phi0) 2) (- (sin phi0) (* -1 (sin phi0)))))))))) (* 1/2 (* hlf_e (+ (* -1 (pow (- (sin phi0) (* -1 (sin phi0))) 2)) (* 2 (* (sin phi0) (- (sin phi0) (* -1 (sin phi0)))))))))) (*.f64 (*.f64 #s(literal 2/3 binary64) (*.f64 ecc ecc)) (*.f64 (*.f64 (*.f64 phi0 phi0) phi0) hlf_e))) (*.f64 (-.f64 (sin.f64 phi0) (neg.f64 (sin.f64 phi0))) hlf_e)) ecc)) 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))))))))
#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 phi0 (+.f64 (*.f64 #s(literal 1/2 binary64) (*.f64 hlf_e (*.f64 phi0 (+.f64 (*.f64 #s(literal -1 binary64) (pow.f64 (-.f64 ecc (*.f64 #s(literal -1 binary64) ecc)) #s(literal 2 binary64))) (*.f64 #s(literal 2 binary64) (*.f64 ecc (-.f64 ecc (*.f64 #s(literal -1 binary64) ecc)))))))) (*.f64 hlf_e (-.f64 ecc (*.f64 #s(literal -1 binary64) ecc))))))))
#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)))))) #s(approx (* ecc (+ (* ecc (+ (* 1/6 (* ecc (* hlf_e (+ (* -6 (* (sin phi0) (pow (- (sin phi0) (* -1 (sin phi0))) 2))) (+ (* 2 (pow (- (sin phi0) (* -1 (sin phi0))) 3)) (* 6 (* (pow (sin phi0) 2) (- (sin phi0) (* -1 (sin phi0)))))))))) (* 1/2 (* hlf_e (+ (* -1 (pow (- (sin phi0) (* -1 (sin phi0))) 2)) (* 2 (* (sin phi0) (- (sin phi0) (* -1 (sin phi0)))))))))) (* hlf_e (- (sin phi0) (* -1 (sin phi0)))))) (*.f64 #s(literal 2 binary64) (*.f64 ecc (*.f64 hlf_e 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 hlf_e) #s(approx (log (/ (+ 1 (* ecc (sin phi0))) (- 1 (* ecc (sin phi0))))) (*.f64 (+.f64 ecc ecc) phi0))))
(*.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)) (pow (cos (* -1/4 (PI))) 2))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI))))))) (sqrt 2))) (* phipp (+ (* 1/6 (* phipp (+ (* -24 (/ (* (pow (cos (* -1/4 (PI))) 2) (* (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* -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)) (pow (cos (* -1/4 (PI))) 2)))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (+ (* 1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* 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)) (pow (cos (* -1/4 (PI))) 2))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 3)) (pow (sqrt 2) 3))))))) (* 1/2 (+ (* -4 (/ (* (pow (cos (* -1/4 (PI))) 2) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* -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)) (pow (cos (* -1/4 (PI))) 2))))) (sqrt 2)))))))) (/.f64 (*.f64 #s(literal 2 binary64) (*.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)))))) (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 (*.f64 #s(approx (+ (* 2 (/ (* (cos (* -1/4 (PI))) (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI))))))) (sqrt 2))) (* phipp (+ (* 1/6 (* phipp (+ (* -24 (/ (* (pow (cos (* -1/4 (PI))) 2) (* (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* -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)) (pow (cos (* -1/4 (PI))) 2)))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (+ (* 1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* 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)) (pow (cos (* -1/4 (PI))) 2))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 3)) (pow (sqrt 2) 3))))))) (* 1/2 (+ (* -4 (/ (* (pow (cos (* -1/4 (PI))) 2) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (pow (cos (* -1/4 (PI))) 2))) (* -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)) (pow (cos (* -1/4 (PI))) 2))))) (sqrt 2)))))))) (*.f64 (*.f64 (*.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)))))) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) #s(literal 2 binary64))) phipp) (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))))))) kR)
(if.f64 (<=.f64 phipp #s(literal -2894802230932905/7237005577332262213973186563042994240829374041602535252466099000494570602496 binary64)) (*.f64 kR (neg.f64 (log.f64 (/.f64 (fabs.f64 (cos.f64 (+.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp)))) (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))))))) (if.f64 (<=.f64 phipp #s(literal 4529838868451533/87112285931760246646623899502532662132736 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 phipp (/.f64 (/.f64 (-.f64 (/.f64 (*.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64))) (*.f64 (*.f64 phipp phipp) (*.f64 phipp phipp))) #s(literal 1/16 binary64)) (+.f64 (/.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)) (*.f64 phipp phipp)) #s(literal 1/4 binary64))) (-.f64 (*.f64 (/.f64 (PI.f64) phipp) #s(literal 1/4 binary64)) #s(literal 1/2 binary64))))))))))
(if.f64 (<=.f64 phipp #s(literal -2894802230932905/7237005577332262213973186563042994240829374041602535252466099000494570602496 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))))))) (if.f64 (<=.f64 phipp #s(literal 4529838868451533/87112285931760246646623899502532662132736 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 phipp (/.f64 (/.f64 (-.f64 (/.f64 (*.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64))) (*.f64 (*.f64 phipp phipp) (*.f64 phipp phipp))) #s(literal 1/16 binary64)) (+.f64 (/.f64 (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)) (*.f64 phipp phipp)) #s(literal 1/4 binary64))) (-.f64 (*.f64 (/.f64 (PI.f64) phipp) #s(literal 1/4 binary64)) #s(literal 1/2 binary64))))))))))
(if.f64 (<=.f64 phipp #s(literal -2894802230932905/7237005577332262213973186563042994240829374041602535252466099000494570602496 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))))))) (if.f64 (<=.f64 phipp #s(literal 6693195169205051/40564819207303340847894502572032 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 (/.f64 (-.f64 (*.f64 (*.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)))))))))
(if.f64 (<=.f64 phipp #s(literal -2894802230932905/7237005577332262213973186563042994240829374041602535252466099000494570602496 binary64)) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (if.f64 (<=.f64 phipp #s(literal 6693195169205051/40564819207303340847894502572032 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 (*.f64 (/.f64 (PI.f64) phipp) #s(literal -1/4 binary64)) phipp))))))))
(if.f64 (<=.f64 phipp #s(literal -2894802230932905/7237005577332262213973186563042994240829374041602535252466099000494570602496 binary64)) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (if.f64 (<=.f64 phipp #s(literal 6693195169205051/40564819207303340847894502572032 binary64)) (*.f64 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(*.f64 kR (log.f64 (/.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)))) (cos.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))
(*.f64 kR (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))))))
(*.f64 (log.f64 (fabs.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))) (*.f64 #s(literal 1/2 binary64) (PI.f64))))))) kR)
(if.f64 (<=.f64 phipp #s(literal -2894802230932905/7237005577332262213973186563042994240829374041602535252466099000494570602496 binary64)) (*.f64 kR (log.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -9/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 #s(approx kR #s(literal 0 binary64)) (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(*.f64 kR (log.f64 (tan.f64 #s(approx (+ (/ (PI) 4) (* 1/2 phipp)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(*.f64 kR lampp)
#s(approx (log (- (/ (* (cos (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))) (/ (* (sin (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))))) (+.f64 (log.f64 (*.f64 #s(literal 1/2 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (*.f64 phi (+.f64 (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (sqrt.f64 #s(literal 2 binary64)))) #s(approx (* phi (+ (* 1/6 (* phi (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (* (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (+ (* 1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4)))))) (* 1/2 (- (* 1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* 1/16 (/ (sqrt 2) (cos (* -1/4 (PI))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3)))))))) (sqrt 2))) (* 16 (/ (* (pow (cos (* -1/4 (PI))) 3) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 3)) (pow (sqrt 2) 3))))))) (* 1/2 (+ (* -4 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 2)) (pow (sqrt 2) 2))) (* 4 (/ (* (cos (* -1/4 (PI))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))))) (sqrt 2))))))) (*.f64 #s(literal 1/6 binary64) (*.f64 (pow.f64 phi #s(literal 2 binary64)) (+.f64 (*.f64 #s(literal -24 binary64) (/.f64 (*.f64 (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))) (*.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (-.f64 (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (pow.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64)))) (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64)))))))))) (pow.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)))) (+.f64 (*.f64 #s(literal 12 binary64) (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal -1/2 binary64) (+.f64 (*.f64 #s(literal -1/48 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (+.f64 (*.f64 #s(literal 1/16 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (pow.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64)) (sqrt.f64 #s(literal 2 binary64))) (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 4 binary64))))))) (*.f64 #s(literal 1/2 binary64) (-.f64 (*.f64 #s(literal 1/48 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (+.f64 (*.f64 #s(literal 1/16 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (*.f64 #s(literal 1/8 binary64) (/.f64 (*.f64 (pow.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64))))))))) (sqrt.f64 #s(literal 2 binary64)))) (*.f64 #s(literal 16 binary64) (/.f64 (*.f64 (pow.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 3 binary64)) (pow.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) #s(literal 3 binary64))) (pow.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 3 binary64)))))))))))))
#s(approx (log (- (/ (* (cos (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))) (/ (* (sin (* -1/2 phi)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* phi 1/2)))))) (+.f64 (log.f64 (*.f64 #s(literal 1/2 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (*.f64 phi #s(approx (+ (* 2 (/ (* (cos (* -1/4 (PI))) (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI))))))) (sqrt 2))) (* phi (+ (* 1/6 (* phi (+ (* -24 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (* (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI))))))))))) (pow (sqrt 2) 2))) (+ (* 12 (/ (* (cos (* -1/4 (PI))) (- (* -1/2 (+ (* -1/48 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (+ (* 1/16 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 3) (sqrt 2)) (pow (cos (* -1/4 (PI))) 4)))))) (* 1/2 (- (* 1/48 (/ (sqrt 2) (cos (* -1/4 (PI))))) (+ (* 1/16 (/ (sqrt 2) (cos (* -1/4 (PI))))) (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3)))))))) (sqrt 2))) (* 16 (/ (* (pow (cos (* -1/4 (PI))) 3) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 3)) (pow (sqrt 2) 3))))))) (* 1/2 (+ (* -4 (/ (* (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))) (pow (- (* -1/4 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))) (* -1/4 (/ (sqrt 2) (cos (* -1/4 (PI)))))) 2)) (pow (sqrt 2) 2))) (* 4 (/ (* (cos (* -1/4 (PI))) (- (* 1/8 (/ (* (pow (sin (* -1/4 (PI))) 2) (sqrt 2)) (pow (cos (* -1/4 (PI))) 3))) (* 1/8 (/ (* (sin (* -1/4 (PI))) (sqrt 2)) (+ 1/2 (* 1/2 (cos (* 2 (* -1/4 (PI)))))))))) (sqrt 2)))))))) (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (-.f64 (*.f64 #s(literal -1/4 binary64) (/.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sqrt.f64 #s(literal 2 binary64))) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal -1/2 binary64) (PI.f64))))))) (*.f64 #s(literal -1/4 binary64) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (sqrt.f64 #s(literal 2 binary64))))))))
(log.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 phi #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (+.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 binary64)))))))
(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))))
(if.f64 (<=.f64 phi #s(literal 8924260225606735/81129638414606681695789005144064 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 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))
(if.f64 (<=.f64 phi #s(literal 3010670175542045/79228162514264337593543950336 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 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi)))))
(log.f64 (tan.f64 (+.f64 (*.f64 phi #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (PI.f64)))))
(log.f64 (tan.f64 (+.f64 #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 (*.f64 (cos.f64 phip) sinp0) (cos.f64 lamp))))
(asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 (*.f64 (cos.f64 phip) (cos.f64 lamp)))))
(if.f64 (<=.f64 lamp #s(literal 5764607523034235/36028797018963968 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 (*.f64 #s(approx (cos lamp) (-.f64 (*.f64 (*.f64 (-.f64 (*.f64 (*.f64 lamp lamp) #s(literal 1/24 binary64)) #s(literal 1/2 binary64)) lamp) lamp) #s(literal -1 binary64))) (cos.f64 phip)) sinp0))) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 #s(approx (* (cos phip) (cos lamp)) (cos.f64 lamp))))))
(if.f64 (<=.f64 phip #s(literal -1450000 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 1357680363825023/1180591620717411303424 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 lamp #s(literal 3170534137668829/576460752303423488 binary64)) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) #s(approx (* sinp0 (* (cos phip) (cos lamp))) (*.f64 sinp0 (cos.f64 phip))))) (asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 #s(approx (* (cos phip) (cos lamp)) (cos.f64 lamp))))))
(if.f64 (<=.f64 phip #s(literal -14500000000000 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (/ 1 (/ 2 (+ (cos (- lamp phip)) (cos (+ lamp phip))))))) (*.f64 (*.f64 (neg.f64 sinp0) (cos.f64 phip)) (cos.f64 lamp)))) (if.f64 (<=.f64 phip #s(literal 419110025354681/295147905179352825856 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp))))) (asin.f64 #s(approx (- (* (sin phip) cosp0) (* (* sinp0 (cos lamp)) (cos phip))) (*.f64 (*.f64 (cos.f64 lamp) (neg.f64 sinp0)) (cos.f64 phip))))))
(if.f64 (<=.f64 phip #s(literal -14500000000000 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (/ 1 (/ 2 (+ (cos (- lamp phip)) (cos (+ lamp phip))))))) (*.f64 (*.f64 (neg.f64 sinp0) (cos.f64 phip)) (cos.f64 lamp)))) (if.f64 (<=.f64 phip #s(literal 419110025354681/295147905179352825856 binary64)) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp))))) (asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (/ 1 (/ 2 (+ (cos (- lamp phip)) (cos (+ lamp phip))))))) (*.f64 (*.f64 (neg.f64 sinp0) (cos.f64 phip)) (cos.f64 lamp))))))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (-.f64 (*.f64 cosp0 phip) (*.f64 sinp0 (cos.f64 lamp)))))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) (*.f64 (neg.f64 (cos.f64 lamp)) sinp0)))
(asin.f64 #s(approx (- (* cosp0 (sin phip)) (* sinp0 (* (cos phip) (cos lamp)))) #s(approx (* -1 (* sinp0 (cos lamp))) (neg.f64 sinp0))))
(if.f64 (<=.f64 phi #s(literal 7036417601633875/858099707516326214372737599885174152158679412517913176174307932398192897924707006515319955082681819372162038923935107254640248499964580476571753536389382144 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/4 binary64)) (*.f64 #s(literal -1/2 binary64) phi)))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 (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 (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi)))))))) K)))) (/.f64 (PI.f64) #s(literal 2 binary64))))
(-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 (*.f64 (-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/4 binary64)) (*.f64 #s(literal -1/2 binary64) phi)))) (*.f64 #s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi)) hlf_e)) c) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))
(if.f64 (<=.f64 K #s(literal 6880495847970215/286687326998758938951352611912760867599570623646035140467198604923365359511060601008752319138765710819328 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 (sqrt.f64 (exp.f64 (+.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 (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))
(if.f64 (<=.f64 K #s(literal 6880495847970215/286687326998758938951352611912760867599570623646035140467198604923365359511060601008752319138765710819328 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(approx 2 #s(literal 0 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))) (-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))
(-.f64 (*.f64 (atan.f64 (exp.f64 #s(approx (+ (* c (- (log (tan (+ (/ (PI) 4) (* 1/2 phi)))) (* hlf_e (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi)))))))) K) K))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 #s(approx (* ecc (+ (* 1/6 (* ecc (+ (* -6 (* (sin phi) (pow (- (sin phi) (* -1 (sin phi))) 2))) (+ (* 2 (pow (- (sin phi) (* -1 (sin phi))) 3)) (* 6 (* (pow (sin phi) 2) (- (sin phi) (* -1 (sin phi))))))))) (* 1/2 (+ (* -1 (pow (- (sin phi) (* -1 (sin phi))) 2)) (* 2 (* (sin phi) (- (sin phi) (* -1 (sin phi))))))))) (*.f64 (*.f64 ecc (*.f64 #s(literal 1/6 binary64) ecc)) #s(approx (- (- (* (* -6 (pow (* 2 (sin phi)) 2)) (sin phi)) (* -2 (pow (* 2 (sin phi)) 3))) (* (* -6 (* 2 (sin phi))) (- 1/2 (* 1/2 (cos (* 2 phi)))))) (-.f64 (*.f64 #s(literal -24 binary64) (pow.f64 (sin.f64 phi) #s(literal 3 binary64))) (+.f64 (*.f64 #s(literal -16 binary64) (pow.f64 (sin.f64 phi) #s(literal 3 binary64))) (*.f64 #s(literal -12 binary64) (*.f64 (sin.f64 phi) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) phi))))))))))) (*.f64 #s(literal 2 binary64) (sin.f64 phi))) ecc))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 ecc (-.f64 (sin.f64 phi) (*.f64 #s(literal -1 binary64) (sin.f64 phi)))))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) #s(approx (* (+ (* ecc (+ (* 1/6 (* ecc (+ (* -6 (* (sin phi) (pow (- (sin phi) (* -1 (sin phi))) 2))) (+ (* 2 (pow (- (sin phi) (* -1 (sin phi))) 3)) (* 6 (* (pow (sin phi) 2) (- (sin phi) (* -1 (sin phi))))))))) (* 1/2 (+ (* -1 (pow (- (sin phi) (* -1 (sin phi))) 2)) (* 2 (* (sin phi) (- (sin phi) (* -1 (sin phi))))))))) (* 2 (sin phi))) ecc) (*.f64 #s(literal 2 binary64) (*.f64 ecc (sin.f64 phi)))))
#s(approx (log (/ (+ 1 (* ecc (sin phi))) (- 1 (* ecc (sin phi))))) (*.f64 (+.f64 ecc ecc) phi))
(*.f64 #s(approx (sqrt (/ (- Up 1) (- -1 Up))) (+.f64 #s(literal 1 binary64) (*.f64 Up (-.f64 (*.f64 #s(literal 1/2 binary64) Up) #s(literal 1 binary64))))) Esc)
#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 (neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (cos.f64 (*.f64 BrA u))))) #s(literal 100000000000000000191567508573466873621595512726519201115280351459937932420398875596123614510818032353280 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 (/.f64 cosgam Vp) Sp) singam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 (*.f64 BrA u) #s(literal 500000000000000011872716179325552678704328963914341093736732494335118714771010286284088814108041647064672984566920058037896706584945040786718720 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 (-.f64 (*.f64 (/.f64 cosgam Vp) Sp) singam))) #s(approx (cos (* BrA u)) #s(literal 1 binary64))))))
(neg.f64 (*.f64 rB (atan2.f64 #s(approx (- (* Sp cosgam) (* Vp singam)) (*.f64 Vp (-.f64 (/.f64 (*.f64 Sp cosgam) Vp) singam))) #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)))))
(*.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))
(exp.f64 (neg.f64 (*.f64 BrA v)))
(if.f64 (<=.f64 (neg.f64 (*.f64 BrA v)) #s(literal -50 binary64)) #s(approx (exp (neg (* BrA v))) (/.f64 #s(literal 1 binary64) #s(approx (/ 2 (- (* -2 (* v BrA)) -2)) (/.f64 #s(literal -1 binary64) (*.f64 #s(approx BrA #s(literal 0 binary64)) v))))) (if.f64 (<=.f64 (neg.f64 (*.f64 BrA v)) #s(literal 9999999999999999538762658202121142272 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 4999999999999999930155298801282288858501320919063181937624830367941782926336371924532423207114480333393189640196327307696676586425126051668137976185307698505365345832344687589284519925536573169820811633035563360005510084776652009298228906344280973600585744230586460911069533464925641061001338333875010535424 binary64)) #s(approx (exp (neg (* BrA v))) (/.f64 (-.f64 #s(literal 1 binary64) (*.f64 (*.f64 (sqrt.f64 (*.f64 (*.f64 BrA BrA) (*.f64 BrA BrA))) 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 (*.f64 BrA v) #s(literal -inf.0 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) (if.f64 (<=.f64 (*.f64 BrA v) #s(literal -4999999999999999769381329101060571136 binary64)) #s(approx (exp (neg (* BrA v))) (/.f64 (-.f64 #s(literal 1 binary64) (*.f64 (*.f64 (*.f64 BrA BrA) v) v)) (+.f64 (*.f64 v BrA) #s(literal 1 binary64)))) (if.f64 (<=.f64 (*.f64 BrA v) #s(literal 1/2 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) #s(approx (exp (neg (* BrA v))) (/.f64 #s(literal 1 binary64) #s(approx (/ 2 (- (* -2 (* v BrA)) -2)) (/.f64 #s(literal -1 binary64) (*.f64 #s(approx BrA #s(literal 0 binary64)) v))))))))
(if.f64 (<=.f64 (*.f64 BrA v) #s(literal -inf.0 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) (if.f64 (<=.f64 (*.f64 BrA v) #s(literal -999999999999999993220948674361627976461708441944064 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)))) (if.f64 (<=.f64 (*.f64 BrA v) #s(literal 1/2 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) #s(approx (exp (neg (* BrA v))) (/.f64 #s(literal 1 binary64) #s(approx (/ 2 (- (* -2 (* v BrA)) -2)) (/.f64 #s(literal -1 binary64) (*.f64 #s(approx BrA #s(literal 0 binary64)) v))))))))
(if.f64 (<=.f64 (*.f64 BrA v) #s(literal 1/2 binary64)) #s(approx (exp (neg (* BrA v))) (-.f64 #s(literal 1 binary64) (*.f64 v BrA))) #s(approx (exp (neg (* BrA v))) (/.f64 #s(literal 1 binary64) #s(approx (/ 2 (- (* -2 (* v BrA)) -2)) (/.f64 #s(literal -1 binary64) (*.f64 #s(approx BrA #s(literal 0 binary64)) 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 4278419646001971/4503599627370496 binary64)) (*.f64 ArB (log.f64 (/.f64 (-.f64 (*.f64 (neg.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sqrt.f64 #s(literal 2 binary64))) #s(literal -1/2 binary64))) (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) #s(literal -1/2 binary64)))))) #s(approx (* ArB (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 (*.f64 #s(literal 2 binary64) (*.f64 (*.f64 ArB gamma0) (/.f64 (*.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #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 2 binary64) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 ArB (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))))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -3602879701896397/72057594037927936 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))) (*.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 #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 3022314549036573/151115727451828646838272 binary64)) #s(approx (* ArB (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 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 gamma0 ArB)) (*.f64 (*.f64 #s(literal -1/4 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (/.f64 (-.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64))))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) ArB))) (*.f64 (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))))) ArB)))
(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 #s(approx (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (*.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 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))) #s(approx (* ArB (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 (*.f64 #s(literal 2 binary64) (*.f64 (*.f64 ArB gamma0) (/.f64 (*.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #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 2 binary64) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 ArB (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))))))))))
(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 #s(approx (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (*.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 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))) #s(approx (* ArB (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 (*.f64 #s(literal 2 binary64) (*.f64 (*.f64 ArB gamma0) (/.f64 (*.f64 (-.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64))))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 #s(literal -1/4 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 ArB (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))))))))))
(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 #s(approx (/ (* (cos (* -1/2 gamma0)) (/ (sqrt 2) 2)) (cos (- (* -1/4 (PI)) (* gamma0 1/2)))) (*.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 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))) #s(approx (* ArB (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 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 gamma0 ArB)) (*.f64 (*.f64 #s(literal -1/4 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (/.f64 (-.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64))))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) ArB))))
(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 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0))))))) #s(approx (* ArB (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 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 gamma0 ArB)) (*.f64 (*.f64 #s(literal -1/4 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (/.f64 (-.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64))))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 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 3022314549036573/151115727451828646838272 binary64)) #s(approx (* ArB (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 (*.f64 (*.f64 #s(literal 2 binary64) (*.f64 gamma0 ArB)) (*.f64 (*.f64 #s(literal -1/4 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (/.f64 (-.f64 (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64))))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) (sqrt.f64 #s(literal 2 binary64))))) (*.f64 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64))))) ArB))) (*.f64 ArB (neg.f64 (log.f64 (/.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64))))) (fabs.f64 (sin.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 7737125245533627/38685626227668133590597632 binary64)) (*.f64 ArB (log.f64 (/.f64 (-.f64 #s(approx (* (neg (* (sqrt 2) 1/2)) (cos (* -1/2 gamma0))) (*.f64 #s(literal -1/2 binary64) (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))) (neg.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0))))))) (*.f64 ArB (neg.f64 (log.f64 (/.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64))))) (fabs.f64 (sin.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 2028240960365167/5070602400912917605986812821504 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* -1/2 gamma0))) (+.f64 #s(literal 1 binary64) (*.f64 gamma0 (+.f64 #s(literal 1/2 binary64) (*.f64 gamma0 (-.f64 (*.f64 #s(literal -1/48 binary64) gamma0) #s(literal 1/8 binary64)))))))))))) (*.f64 ArB (neg.f64 (log.f64 (/.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64))))) (fabs.f64 (sin.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 3022314549036573/151115727451828646838272 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* -1/2 gamma0))) (+.f64 #s(literal 1 binary64) (*.f64 gamma0 (+.f64 #s(literal 1/2 binary64) (*.f64 gamma0 (-.f64 (*.f64 #s(literal -1/48 binary64) gamma0) #s(literal 1/8 binary64)))))))))))) (*.f64 ArB (log.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))) (cos.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 3022314549036573/151115727451828646838272 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* -1/2 gamma0))) (+.f64 #s(literal 1 binary64) (*.f64 gamma0 (+.f64 #s(literal 1/2 binary64) (*.f64 gamma0 (-.f64 (*.f64 #s(literal -1/48 binary64) gamma0) #s(literal 1/8 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 3022314549036573/151115727451828646838272 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) #s(approx (* (* (sqrt 2) 1/2) (- (cos (* gamma0 1/2)) (sin (* -1/2 gamma0)))) (+.f64 (*.f64 #s(literal 1/4 binary64) (*.f64 gamma0 (sqrt.f64 #s(literal 2 binary64)))) (*.f64 #s(literal 1/2 binary64) (sqrt.f64 #s(literal 2 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 3022314549036573/151115727451828646838272 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(approx (- (* (neg (* (sqrt 2) 1/2)) (cos (* -1/2 gamma0))) (* (* (sin (* -1/2 gamma0)) (sqrt 2)) -1/2)) (+.f64 (*.f64 #s(literal -1/2 binary64) (sqrt.f64 #s(literal 2 binary64))) (*.f64 #s(literal -1/4 binary64) (*.f64 gamma0 (sqrt.f64 #s(literal 2 binary64)))))) (neg.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0))))))) (*.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 3022314549036573/151115727451828646838272 binary64)) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* gamma0 1/2)) (sin (* -1/2 gamma0))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 1/2 binary64) gamma0)))))))) (*.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 2028240960365167/5070602400912917605986812821504 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 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0)))))))
(*.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))))))
(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 (/.f64 (-.f64 (*.f64 (log.f64 (fabs.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (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 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))) (*.f64 (log.f64 (fabs.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))) (log.f64 (fabs.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))) (+.f64 (log.f64 (fabs.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (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))))))) ArB) #s(approx (* ArB (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 (*.f64 (*.f64 ArB (*.f64 (-.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(literal 1/4 binary64)) (sqrt.f64 #s(literal 2 binary64))))) (+.f64 gamma0 gamma0)) (*.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 1/2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) ArB))))
(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 (-.f64 (log.f64 (fabs.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 gamma0 #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))))))) #s(approx (* ArB (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 (*.f64 (*.f64 ArB (*.f64 (-.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(literal 1/4 binary64)) (sqrt.f64 #s(literal 2 binary64))))) (+.f64 gamma0 gamma0)) (*.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 1/2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) ArB))))
(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 (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (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) #s(approx (* ArB (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 (*.f64 (*.f64 ArB (*.f64 (-.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(literal 1/4 binary64)) (sqrt.f64 #s(literal 2 binary64))))) (+.f64 gamma0 gamma0)) (*.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 1/2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) ArB))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 5134103575202365/9007199254740992 binary64)) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 #s(approx (+ (* gamma0 -2) (PI)) (*.f64 gamma0 (-.f64 (/.f64 (PI.f64) gamma0) #s(literal 2 binary64)))) #s(literal -1/4 binary64))))) (log.f64 (fabs.f64 (cos.f64 (*.f64 #s(approx (+ (* gamma0 -2) (PI)) (*.f64 gamma0 (-.f64 (/.f64 (PI.f64) gamma0) #s(literal 2 binary64)))) #s(literal -1/4 binary64))))))) #s(approx (* ArB (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 (*.f64 (*.f64 ArB (*.f64 (-.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(literal 1/4 binary64)) (sqrt.f64 #s(literal 2 binary64))))) (+.f64 gamma0 gamma0)) (*.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 1/2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) ArB))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 5134103575202365/9007199254740992 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) #s(approx (* ArB (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 (*.f64 (*.f64 ArB (*.f64 (-.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sqrt.f64 #s(literal 2 binary64))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(literal 1/4 binary64)) (sqrt.f64 #s(literal 2 binary64))))) (+.f64 gamma0 gamma0)) (*.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (/.f64 #s(literal 1/2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) ArB))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 6124895493223875/36028797018963968 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 ArB (log.f64 (*.f64 (/.f64 #s(literal 1 binary64) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* -1/2 gamma0)) (sin (* gamma0 1/2))) (+.f64 #s(literal 1 binary64) (*.f64 #s(literal -1/2 binary64) gamma0))))))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 6124895493223875/36028797018963968 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 (log.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) #s(approx (- (cos (* -1/2 gamma0)) (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))))) ArB))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4503599627370495/4503599627370496 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 (log.f64 (/.f64 #s(approx (sin (- (* -1/2 gamma0) (* -1/4 (PI)))) (*.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))))) ArB))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4503599627370495/4503599627370496 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 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)))))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4503599627370495/4503599627370496 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 #s(approx ArB #s(literal 0 binary64)) (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))))
(*.f64 ArB (log.f64 (fabs.f64 (/.f64 #s(approx (sin (- (* -1/2 gamma0) (* -1/4 (PI)))) (*.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 #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))))))
#s(approx (- 1 (* es (* (sin phi0) (sin phi0)))) #s(literal 1 binary64))
(/.f64 (*.f64 (sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) one_es) es))) com) (*.f64 (cos.f64 phi0) (sqrt.f64 (-.f64 #s(literal 1 binary64) (*.f64 es (*.f64 (sin.f64 phi0) (sin.f64 phi0)))))))
(/.f64 #s(literal 1 binary64) (/.f64 (*.f64 (sqrt.f64 (-.f64 #s(literal 1 binary64) (*.f64 (*.f64 (-.f64 #s(literal 1 binary64) (cos.f64 (+.f64 phi0 phi0))) #s(literal 1/2 binary64)) es))) (cos.f64 phi0)) (*.f64 com (sqrt.f64 (-.f64 (/.f64 (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es) one_es) #s(literal -1 binary64))))))
#s(approx (/ (* (sqrt (/ (+ one_es (* (pow (cos phi0) 4) es)) one_es)) com) (* (cos phi0) (sqrt (- 1 (* es (* (sin phi0) (sin phi0))))))) #s(approx (+ (* com (sqrt (/ (+ es one_es) one_es))) (* (pow phi0 2) (- (* -1 (/ (* com es) (* one_es (sqrt (/ (+ es one_es) one_es))))) (* com (* (sqrt (/ (+ es one_es) one_es)) (- (* -1/2 es) 1/2)))))) (*.f64 (-.f64 (sqrt.f64 (/.f64 (+.f64 one_es es) one_es)) (*.f64 (+.f64 (*.f64 (sqrt.f64 (/.f64 (+.f64 one_es es) one_es)) (-.f64 (*.f64 #s(literal -1/2 binary64) es) #s(literal 1/2 binary64))) (/.f64 es (*.f64 (sqrt.f64 (/.f64 (+.f64 one_es es) one_es)) one_es))) (*.f64 phi0 phi0))) com)))
(/.f64 (*.f64 #s(literal 1 binary64) com) (*.f64 (sqrt.f64 (fabs.f64 (/.f64 one_es (+.f64 (*.f64 (*.f64 (*.f64 #s(approx (cos phi0) #s(literal 1 binary64)) #s(approx (cos phi0) #s(literal 1 binary64))) (*.f64 #s(approx (cos phi0) #s(literal 1 binary64)) #s(approx (cos phi0) #s(literal 1 binary64)))) es) one_es)))) (*.f64 #s(approx (cos phi0) #s(literal 1 binary64)) (sqrt.f64 (-.f64 #s(literal 1 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) #s(approx (* (cos (+ phi0 phi0)) 1/2) #s(literal 1/2 binary64))) es))))))
(/.f64 (*.f64 (sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 (/.f64 (*.f64 (*.f64 #s(approx (cos phi0) #s(literal 1 binary64)) #s(approx (cos phi0) #s(literal 1 binary64))) (*.f64 #s(approx (cos phi0) #s(literal 1 binary64)) #s(approx (cos phi0) #s(literal 1 binary64)))) one_es) es))) com) (*.f64 #s(approx (cos phi0) #s(literal 1 binary64)) (sqrt.f64 #s(approx (- 1 (* es (* (sin phi0) (sin phi0)))) #s(literal 1 binary64)))))
#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)))))
(if.f64 (<=.f64 one_es #s(literal 4837232111892845/460688772561223309599799143187044542348256984914116263229517370635452168260760357920669766257038423772151901248872539660616526444082616288154471954520592778765795328 binary64)) #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 (* es (* one_es (pow (cos phi0) 4))))) (* one_es (* (cos phi0) (sqrt (- 1 (* es (pow (sin phi0) 2))))))) #s(approx (/ (* com (sqrt (* es one_es))) one_es) (*.f64 com (sqrt.f64 (/.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 (cos phi0)) com)))
#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 (*.f64 (*.f64 (cos.f64 phi0) (cos.f64 phi0)) (*.f64 (cos.f64 phi0) (cos.f64 phi0)))) one_es)))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) one_es) es)))
#s(approx (sqrt (+ 1 (/ (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) one_es))) (-.f64 (sqrt.f64 (-.f64 (/.f64 es one_es) #s(literal -1 binary64))) (/.f64 (*.f64 (*.f64 phi0 phi0) es) (*.f64 (sqrt.f64 (-.f64 (/.f64 es one_es) #s(literal -1 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 #s(approx (* es (* (* (cos phi0) (cos phi0)) (* (cos phi0) (cos phi0)))) es) one_es)))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (*.f64 (-.f64 (*.f64 #s(literal -1/3 binary64) (*.f64 (*.f64 U ArB) U)) ArB) U))
#s(approx (* 1/2 (* ArB (log (/ (- 1 U) (+ 1 U))))) (*.f64 (neg.f64 U) ArB))
(*.f64 A lam)
(if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 49999999999999996962677625276823109300201436100586624765953857856616022815066169514216546287202538742184280590280810862893585968713180152651178994204334413874936507208410055205338551265812204529218599012742757995383198412754109163297745561348039749026730174593312862032038021904229799310374521740690718720 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (sin.f64 (-.f64 (*.f64 lam B) (*.f64 (PI.f64) #s(literal -1/2 binary64)))))) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 #s(approx (sin (* B lam)) (*.f64 lam (+.f64 B (*.f64 #s(literal -1/6 binary64) (*.f64 (pow.f64 B #s(literal 3 binary64)) (pow.f64 lam #s(literal 2 binary64))))))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))))
(if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 1999999999999999890644667736494890251419293140042495849331683229696 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) (if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 10000000000000000531660196626596490356033894575245100973356972987043891522292165594595004291349304909025721681812512093962950445138053653873169216309020403876699170397334223513449750683762833231235463783529148067211236930570359138156544 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (neg.f64 (+.f64 (PI.f64) (*.f64 lam B)))) singam)) (cos.f64 (*.f64 B lam)))) (*.f64 ArB (atan2.f64 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (/.f64 #s(approx (* singam (sin (* B lam))) (*.f64 B (*.f64 lam singam))) cosgam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 10000000000000000531660196626596490356033894575245100973356972987043891522292165594595004291349304909025721681812512093962950445138053653873169216309020403876699170397334223513449750683762833231235463783529148067211236930570359138156544 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 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (/.f64 #s(approx (* singam (sin (* B lam))) (*.f64 B (*.f64 lam singam))) cosgam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64)))))
(if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 1999999999999999890644667736494890251419293140042495849331683229696 binary64)) (*.f64 ArB (atan2.f64 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (/.f64 (*.f64 singam (sin.f64 (*.f64 B lam))) cosgam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64)))) (if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 249999999999999981127280380618811716294668050071597603343410047731917694218676725250216868646074079448025526803849324428504314495201949473268410963248002115317291743547418889035478194203043299371784807625853355592549187287252992 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (neg.f64 (+.f64 (PI.f64) (*.f64 lam B)))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))) (*.f64 ArB (atan2.f64 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (/.f64 #s(approx (* singam (sin (* B lam))) (*.f64 B (*.f64 lam singam))) cosgam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 1999999999999999890644667736494890251419293140042495849331683229696 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)))) (if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 249999999999999981127280380618811716294668050071597603343410047731917694218676725250216868646074079448025526803849324428504314495201949473268410963248002115317291743547418889035478194203043299371784807625853355592549187287252992 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (neg.f64 (+.f64 (PI.f64) (*.f64 lam B)))) singam)) #s(approx (cos (* B lam)) #s(literal 1 binary64)))) (*.f64 ArB (atan2.f64 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (/.f64 #s(approx (* singam (sin (* B lam))) (*.f64 B (*.f64 lam singam))) cosgam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64))))))
(if.f64 (<=.f64 (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))) #s(literal 200000000000000007778715510217677686261474498590404026668604764015382588578769793526159931215755402774652920623882426582706341222818875123308036734442537880708869172525233887089132911615311892438644481327104 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 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (/.f64 #s(approx (* singam (sin (* B lam))) (*.f64 B (*.f64 lam singam))) cosgam)))) #s(approx (cos (* B lam)) #s(literal 1 binary64)))))
(*.f64 ArB (atan2.f64 #s(approx (+ (* S cosgam) (* (sin (* B lam)) singam)) (*.f64 cosgam (+.f64 S (*.f64 #s(approx (/ (sin (* lam B)) cosgam) (/.f64 (*.f64 B lam) cosgam)) 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 (-.f64 u u0) cosrot) #s(literal -39999999999999997616291042021410332095933184403421192921799076575753209026647455352718401016207802277498189570060273431092509962742598192468558863886980588966057607168 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot (-.f64 u u0))) (if.f64 (<=.f64 (*.f64 (-.f64 u u0) cosrot) #s(literal 4999999999999999650599673463152198642336665750694884246307948430823614916415456951880981793447127233788614017024 binary64)) (-.f64 #s(approx (* (- u u0) cosrot) (*.f64 cosrot u)) (*.f64 v sinrot)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot (-.f64 u u0)))))
(if.f64 (<=.f64 (*.f64 (-.f64 u u0) cosrot) #s(literal -1018517988167243/509258994083621521567111422102344540262867098416484062659035112338595324940834176545849344 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot (-.f64 u u0))) (if.f64 (<=.f64 (*.f64 (-.f64 u u0) cosrot) #s(literal 40000000 binary64)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 (neg.f64 v) sinrot)) #s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 cosrot (-.f64 u u0)))))
#s(approx (- (* (- u u0) cosrot) (* v sinrot)) (*.f64 (neg.f64 v) sinrot))
(+.f64 (*.f64 v cosrot) (*.f64 (-.f64 u u0) sinrot))
(if.f64 (<=.f64 u #s(literal -140000000000000000466218943549635907608711084085526638444194611149466960436775348911123988960935926368320365421584610683565713532393619456 binary64)) (+.f64 (*.f64 v cosrot) #s(approx (* (- u u0) sinrot) (*.f64 sinrot u))) (if.f64 (<=.f64 u #s(literal -6084722881095501/81129638414606681695789005144064 binary64)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0))) (if.f64 (<=.f64 u #s(literal 872305872233851/83076749736557242056487941267521536 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (-.f64 (*.f64 cosrot v) (*.f64 sinrot u0))) (+.f64 (*.f64 v cosrot) #s(approx (* (- u u0) sinrot) (*.f64 sinrot u))))))
(if.f64 (<=.f64 (*.f64 v cosrot) #s(literal -348449143727041/348449143727040986586495598010130648530944 binary64)) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (-.f64 (*.f64 cosrot v) (*.f64 sinrot u0))) (if.f64 (<=.f64 (*.f64 v cosrot) #s(literal 100000000000000000026609864708367276537402401181200809098131977453489758916313088 binary64)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0))) #s(approx (+ (* v cosrot) (- (* sinrot u) (* u0 sinrot))) (-.f64 (*.f64 cosrot v) (*.f64 sinrot u0)))))
(if.f64 (<=.f64 (*.f64 v cosrot) #s(literal -5000000000000000235300672479527347945779800703933315382139354767449124765696 binary64)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 v #s(approx (+ cosrot (/ (* sinrot (- u u0)) v)) cosrot))) (if.f64 (<=.f64 (*.f64 v cosrot) #s(literal 100000000000000000026609864708367276537402401181200809098131977453489758916313088 binary64)) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 sinrot (-.f64 u u0))) #s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 v #s(approx (+ cosrot (/ (* sinrot (- u u0)) v)) cosrot)))))
#s(approx (+ (* v cosrot) (* (- u u0) sinrot)) (*.f64 v #s(approx (+ cosrot (/ (* sinrot (- u u0)) v)) cosrot)))
(sqrt.f64 (+.f64 (pow.f64 (sinh.f64 Ce) #s(literal 2 binary64)) #s(literal 1 binary64)))
(if.f64 (<=.f64 Ce #s(literal 480000000000000034529005388137557393408 binary64)) #s(approx (sqrt (+ (pow (sinh Ce) 2) 1)) (/.f64 (/.f64 (-.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 (*.f64 (*.f64 (*.f64 Ce Ce) #s(literal 1/4 binary64)) Ce) Ce) #s(literal -1 binary64)) (*.f64 (*.f64 #s(literal 1/16 binary64) (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce)) (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce)))) (-.f64 (*.f64 (*.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 (sqrt.f64 (*.f64 (*.f64 #s(literal 1/16 binary64) (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce)) (*.f64 (*.f64 (*.f64 Ce Ce) Ce) Ce)))))))
(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))
#s(approx (exp (* 2 Ce)) (+.f64 #s(literal 1 binary64) (*.f64 Ce (+.f64 #s(literal 2 binary64) (*.f64 Ce (-.f64 (*.f64 (*.f64 #s(literal 16/9 binary64) Ce) (/.f64 Ce (-.f64 (*.f64 #s(literal 4/3 binary64) Ce) #s(literal 2 binary64)))) #s(approx (/ -4 (- 2 (* 4/3 Ce))) (-.f64 (*.f64 Ce (-.f64 (*.f64 Ce (-.f64 (*.f64 #s(literal -16/27 binary64) Ce) #s(literal 8/9 binary64))) #s(literal 4/3 binary64))) #s(literal 2 binary64)))))))))
#s(approx (exp (* 2 Ce)) (+.f64 #s(literal 1 binary64) (*.f64 Ce (+.f64 #s(literal 2 binary64) (*.f64 Ce (-.f64 (*.f64 (*.f64 #s(literal 16/9 binary64) Ce) (/.f64 Ce (-.f64 (*.f64 #s(literal 4/3 binary64) Ce) #s(literal 2 binary64)))) #s(approx (/ -4 (- 2 (* 4/3 Ce))) (-.f64 (*.f64 Ce (-.f64 (*.f64 #s(literal -8/9 binary64) Ce) #s(literal 4/3 binary64))) #s(literal 2 binary64)))))))))
#s(approx (exp (* 2 Ce)) (+.f64 #s(literal 1 binary64) (+.f64 Ce (+.f64 Ce (/.f64 (*.f64 (*.f64 (-.f64 (*.f64 (*.f64 #s(literal 16/9 binary64) Ce) Ce) #s(literal 4 binary64)) Ce) Ce) (-.f64 (*.f64 #s(literal 4/3 binary64) Ce) #s(literal 2 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 (*.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 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64))) (*.f64 (cos.f64 Cn) (*.f64 (cos.f64 lam) (*.f64 (cos.f64 lam) (cos.f64 Cn)))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (pow.f64 (sin.f64 Cn) #s(literal 2 binary64)) (*.f64 (*.f64 (cos.f64 lam) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn))))) (cos.f64 lam)))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (pow.f64 (sin.f64 Cn) #s(literal 2 binary64)) (pow.f64 (*.f64 (cos.f64 lam) (cos.f64 Cn)) #s(literal 2 binary64)))))
(/.f64 #s(literal 2 binary64) (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 #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 (cos.f64 Cn) (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))))))) (sin.f64 lam))
(*.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 (/.f64 (cos.f64 Cn) #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))))))) (sin.f64 lam))
(if.f64 (<=.f64 lam #s(literal 8116567392432203/9223372036854775808 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) 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 (pow.f64 (cos.f64 lam) #s(literal 2 binary64))))))
(if.f64 (<=.f64 lam #s(literal 8116567392432203/9223372036854775808 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) lam)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (tan.f64 (asin.f64 (sin.f64 lam)))))
(if.f64 (<=.f64 lam #s(literal 8116567392432203/9223372036854775808 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) lam)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (sinh.f64 (atanh.f64 (sin.f64 lam)))))
(if.f64 (<=.f64 lam #s(literal 8116567392432203/9223372036854775808 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) lam)) #s(approx (* (/ (cos Cn) (sqrt (+ (* (- 1/2 (* -1/2 (cos (+ lam lam)))) (- 1/2 (* -1/2 (cos (+ Cn Cn))))) (- 1/2 (* 1/2 (cos (+ Cn Cn))))))) (sin lam)) (/.f64 (sin.f64 lam) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) lam))))))))
(if.f64 (<=.f64 lam #s(literal 2 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) 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 Cn)))))
(if.f64 (<=.f64 lam #s(literal 2 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) 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) (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal -1/2 binary64)))))))
(if.f64 (<=.f64 lam #s(literal 28000000000000000831033663644940973731051658669738426368 binary64)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) lam)) #s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (/.f64 #s(literal 1 binary64) (/.f64 (neg.f64 (fabs.f64 #s(approx (cos Cn) #s(literal 1 binary64)))) (neg.f64 (sin.f64 Cn))))))
#s(approx (/ (* (sin lam) (cos Cn)) (sqrt (+ (* (sin Cn) (sin Cn)) (* (* (cos Cn) (cos Cn)) (* (cos lam) (cos lam)))))) (*.f64 (cos.f64 Cn) 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)))
(if.f64 (<=.f64 (cos.f64 lam) #s(literal 281474976654361/281474976710656 binary64)) (atan2.f64 (sin.f64 Cn) #s(approx (* (cos Cn) (cos lam)) (cos.f64 lam))) (atan2.f64 (sin.f64 Cn) (/.f64 #s(literal 1 binary64) #s(approx (/ 1 (* (cos Cn) (cos lam))) (/.f64 #s(literal 1 binary64) (cos.f64 Cn))))))
(atan2.f64 (sin.f64 Cn) #s(approx (* (cos Cn) (cos lam)) (cos.f64 lam)))
(if.f64 (<=.f64 Cn #s(literal -2400000000000000067108864 binary64)) (atan2.f64 (sin.f64 Cn) (/.f64 #s(literal 1 binary64) #s(approx (/ 1 (* (cos Cn) (cos lam))) (/.f64 #s(literal 1 binary64) #s(approx (cos Cn) #s(literal 1 binary64)))))) (if.f64 (<=.f64 Cn #s(literal 16000000000000 binary64)) (atan2.f64 #s(approx (sin Cn) (*.f64 (-.f64 (*.f64 (*.f64 Cn Cn) #s(literal -1/6 binary64)) #s(literal -1 binary64)) Cn)) #s(approx (* (cos Cn) (cos lam)) (cos.f64 lam))) (atan2.f64 (sin.f64 Cn) (/.f64 #s(literal 1 binary64) #s(approx (/ 1 (* (cos Cn) (cos lam))) (/.f64 #s(literal 1 binary64) #s(approx (cos Cn) #s(literal 1 binary64))))))))
(atan2.f64 #s(approx (sin Cn) (*.f64 (-.f64 (*.f64 (*.f64 Cn Cn) #s(literal -1/6 binary64)) #s(literal -1 binary64)) Cn)) #s(approx (* (cos Cn) (cos lam)) (cos.f64 lam)))
(*.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 1488565707357403/372141426839350727961253789638658321589064376671906846864122981980487315514059736743009817965446945567110411062408283101969716033850703872 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 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 sinh_arg_i #s(literal -7241683770807669/43888992550349509466047490009497674160595141087458656560896015907649579054077383577321405596290902348906277802702976893959665470901957183225792829745965362396915989605680050116284721582335472197132100330098878361532751631431265351922342068003131649698083620772604076718737120590449149488485194416412164096 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi sinh_arg_i) cos_arg_r)) (if.f64 (<=.f64 sinh_arg_i #s(literal 8392312103782453/161390617380431786853494948250188242145606612051826469551916209783790476376052574664352834580008614464743948248296718336 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)) (*.f64 (-.f64 (*.f64 (*.f64 (/.f64 sin_arg_r hi) cosh_arg_i) hr) (*.f64 sinh_arg_i cos_arg_r)) hi))))
(if.f64 (<=.f64 hi #s(literal 1176225305839797/3794275180128377091639574036764685364535950857523710002444946112771297432041422848 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)) (if.f64 (<=.f64 hi #s(literal 58000000000000001685522420760352126685681153878680752122241818158390959365444175385176984151484393491410245566254478736654550257918172871923242562816406235623989395067041053724349210073407330046847625554034688 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)) (*.f64 hi (-.f64 (/.f64 (*.f64 cosh_arg_i (*.f64 hr sin_arg_r)) hi) (*.f64 cos_arg_r sinh_arg_i))))))
(if.f64 (<=.f64 sinh_arg_i #s(literal -7241683770807669/43888992550349509466047490009497674160595141087458656560896015907649579054077383577321405596290902348906277802702976893959665470901957183225792829745965362396915989605680050116284721582335472197132100330098878361532751631431265351922342068003131649698083620772604076718737120590449149488485194416412164096 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi sinh_arg_i) cos_arg_r)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)))
(if.f64 (<=.f64 hi #s(literal 1176225305839797/3794275180128377091639574036764685364535950857523710002444946112771297432041422848 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)))
(if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 20000000000000000057357570219907446497404120129229967567146859853820771307804544319366583914666449299233916626257196608020375872770963560895599534369611732108691868080208166641175396430819444098873307923634804982550384038403414239739984162143459594327374818907829826579083558912 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)) (*.f64 (*.f64 (neg.f64 hi) sinh_arg_i) cos_arg_r)))
(if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal -8920298079412249/178405961588244985132285746181186892047843328 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 sin_arg_r #s(approx (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr)) (*.f64 cosh_arg_i hr)))) (if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal 578960446186581/28948022309329048855892746252171976963317496166410141009864396001978282409984 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 (*.f64 (neg.f64 hi) sinh_arg_i) cos_arg_r)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 sin_arg_r #s(approx (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr)) (*.f64 cosh_arg_i hr))))))
(if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal -8920298079412249/178405961588244985132285746181186892047843328 binary64)) #s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 sin_arg_r #s(approx (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr)) (*.f64 cosh_arg_i hr)))) (if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) #s(literal 578960446186581/28948022309329048855892746252171976963317496166410141009864396001978282409984 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)) (*.f64 sin_arg_r #s(approx (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr)) (*.f64 cosh_arg_i hr))))))
#s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 sin_arg_r #s(approx (+ (* -1 (/ (* cos_arg_r (* hi sinh_arg_i)) sin_arg_r)) (* cosh_arg_i hr)) (*.f64 cosh_arg_i hr))))
#s(approx (- (* (* sin_arg_r cosh_arg_i) hr) (* (* cos_arg_r sinh_arg_i) hi)) (*.f64 cosh_arg_i (*.f64 hr sin_arg_r)))
(if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 5288447750321989/10576895500643977583230644928524336637254474927428499508554380724390492659780981533203027367035444557561459392400373732868096 binary64)) (+.f64 (*.f64 (*.f64 hi sin_arg_r) cosh_arg_i) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)) (if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 199999999999999984924696874707920970120897867915847050405221309697980696558932154585003938846536810050657940462325091296687310550613357744883467580356118956661470790120934939455989945801060127957611687906204227736000759240738169004268617951010459111545827258847272611683204755172652495528787968 binary64)) (+.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 cos_arg_r (*.f64 hr sinh_arg_i)))))
(if.f64 (<=.f64 cosh_arg_i #s(literal -2571913738095105/612360413832167804618481300175204950565278972827733245154176943827040045789568780701801476101110276210469073715640490142747228062965884323130340172286597179476547016660734615078156785793174374530940928 binary64)) (+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.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 (*.f64 cos_arg_r sinh_arg_i) #s(literal 199999999999999984924696874707920970120897867915847050405221309697980696558932154585003938846536810050657940462325091296687310550613357744883467580356118956661470790120934939455989945801060127957611687906204227736000759240738169004268617951010459111545827258847272611683204755172652495528787968 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 cos_arg_r (*.f64 hr sinh_arg_i))))
(if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 199999999999999984924696874707920970120897867915847050405221309697980696558932154585003938846536810050657940462325091296687310550613357744883467580356118956661470790120934939455989945801060127957611687906204227736000759240738169004268617951010459111545827258847272611683204755172652495528787968 binary64)) (+.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 cos_arg_r (*.f64 hr sinh_arg_i))))
(if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) #s(literal -49999999999999999737683287595902466157897305225341087810970847365954154269153922568421376 binary64)) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (*.f64 sin_arg_r #s(approx (+ (* cosh_arg_i hi) (/ (* cos_arg_r (* hr sinh_arg_i)) sin_arg_r)) (*.f64 cosh_arg_i hi)))) (if.f64 (<=.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) #s(literal 2000000000000000000000 binary64)) #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))) #s(approx (+ (* (* sin_arg_r cosh_arg_i) hi) (* (* cos_arg_r sinh_arg_i) hr)) (*.f64 sin_arg_r #s(approx (+ (* cosh_arg_i hi) (/ (* cos_arg_r (* hr sinh_arg_i)) sin_arg_r)) (*.f64 cosh_arg_i hi))))))
#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 (cos.f64 phi) 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))))
(/.f64 (*.f64 (cos.f64 phi) lam) #s(approx (sqrt (- 1 (* es (* (sin phi) (sin phi))))) (-.f64 (*.f64 (*.f64 #s(literal -1/2 binary64) es) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) #s(literal -1 binary64))))
(/.f64 (*.f64 (cos.f64 phi) lam) (sqrt.f64 (-.f64 (*.f64 (-.f64 #s(approx (* 1/2 (cos (+ phi phi))) #s(literal 1/2 binary64)) #s(literal 1/2 binary64)) es) #s(literal -1 binary64))))
(/.f64 #s(literal 1 binary64) #s(approx (/ (sqrt (- 1 (* es (* (sin phi) (sin phi))))) (* lam (cos phi))) (/.f64 #s(literal 1 binary64) (*.f64 lam (cos.f64 phi)))))
#s(approx (/ (* (cos phi) lam) (sqrt (- 1 (* es (* (sin phi) (sin phi)))))) (+.f64 lam (*.f64 #s(approx (* (* (- lam (* es lam)) -1/2) phi) (*.f64 #s(literal -1/2 binary64) (*.f64 lam #s(approx phi #s(literal 0 binary64))))) #s(approx phi #s(literal 0 binary64)))))
#s(approx (/ (* (cos phi) lam) (sqrt (- 1 (* es (* (sin phi) (sin phi)))))) (+.f64 lam (*.f64 #s(approx (* (* (- lam (* es lam)) -1/2) phi) (*.f64 #s(literal -1/2 binary64) (*.f64 lam phi))) 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 x (sqrt.f64 (-.f64 #s(literal 1 binary64) (*.f64 es (*.f64 (sin.f64 phi) (sin.f64 phi)))))) k0)
(/.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)