(*.f64 esp (-.f64 angle phi0))
(if.f64 (<=.f64 (cos.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal -7566047373982433/18014398509481984 binary64)) (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 #s(literal 1 binary64) (/.f64 esp y))))) (atan2.f64 g (+.f64 (*.f64 (sin.f64 phi0) (cos.f64 (-.f64 (/.f64 y esp) (*.f64 (PI.f64) #s(literal -1/2 binary64))))) (*.f64 (sin.f64 (+.f64 (/.f64 (neg.f64 y) esp) (*.f64 #s(literal 1/2 binary64) (PI.f64)))) (cos.f64 phi0)))))
(if.f64 (<=.f64 (cos.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal -3422735716801577/9007199254740992 binary64)) (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 #s(literal 1 binary64) (/.f64 esp y))))) (atan2.f64 g (+.f64 (*.f64 (sin.f64 phi0) (cos.f64 (-.f64 (/.f64 y esp) (*.f64 (PI.f64) #s(literal -1/2 binary64))))) (*.f64 (cos.f64 (/.f64 y esp)) (cos.f64 phi0)))))
(atan2.f64 g (+.f64 (*.f64 (sin.f64 phi0) (cos.f64 (*.f64 (+.f64 (PI.f64) (*.f64 (/.f64 #s(literal 2 binary64) esp) y)) #s(literal 1/2 binary64)))) (*.f64 (cos.f64 (/.f64 y esp)) (cos.f64 phi0))))
(atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 #s(literal 1 binary64) (/.f64 esp y)))))
(atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp))))
(asin.f64 (/.f64 (sqrt.f64 (pow.f64 (+.f64 (*.f64 (sin.f64 (-.f64 (/.f64 (neg.f64 y) esp) (*.f64 #s(literal -1/2 binary64) (PI.f64)))) (sin.f64 phi0)) (*.f64 (sin.f64 (/.f64 y esp)) (cos.f64 phi0))) #s(literal 2 binary64))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))))
(asin.f64 (sqrt.f64 (/.f64 (pow.f64 (+.f64 (*.f64 (sin.f64 (/.f64 y esp)) (cos.f64 phi0)) (*.f64 (cos.f64 (/.f64 y esp)) (sin.f64 phi0))) #s(literal 2 binary64)) (+.f64 #s(literal 1 binary64) (*.f64 g g)))))
(asin.f64 (sqrt.f64 (/.f64 (pow.f64 (sin.f64 (+.f64 (/.f64 y esp) phi0)) #s(literal 2 binary64)) (+.f64 #s(literal 1 binary64) (*.f64 g g)))))
(asin.f64 (/.f64 (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (+.f64 (/.f64 y esp) phi0)))))) (sqrt.f64 (fabs.f64 (-.f64 (*.f64 g g) #s(literal -1 binary64))))))
(asin.f64 (sqrt.f64 (*.f64 (/.f64 #s(literal 1 binary64) (-.f64 (*.f64 g g) #s(literal -1 binary64))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (+.f64 (/.f64 y esp) phi0))))))))
(asin.f64 (sqrt.f64 (/.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (+.f64 (/.f64 y esp) phi0))))) (-.f64 (*.f64 g g) #s(literal -1 binary64)))))
(sinh.f64 (/.f64 x esp))
(/.f64 (-.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (sqrt.f64 (pow.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal 1/2 binary64) (PI.f64))) #s(literal 4 binary64)))) (-.f64 #s(literal 0 binary64) (*.f64 #s(literal -1/2 binary64) (PI.f64))))
(/.f64 (-.f64 (neg.f64 (pow.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64))) (neg.f64 (pow.f64 (-.f64 (atan.f64 (exp.f64 (/.f64 y kR))) (*.f64 #s(literal 1/2 binary64) (PI.f64))) #s(literal 2 binary64)))) (+.f64 #s(literal 0 binary64) (*.f64 #s(literal -1/2 binary64) (PI.f64))))
(/.f64 (*.f64 (-.f64 (*.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 #s(literal 0 binary64) (*.f64 #s(literal -1/2 binary64) (PI.f64)))) (-.f64 #s(literal 0 binary64) (*.f64 #s(literal -1/2 binary64) (PI.f64))))
(+.f64 (*.f64 (atan.f64 (exp.f64 (/.f64 y kR))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))
(-.f64 phip delp)
(asin.f64 (+.f64 (*.f64 (*.f64 sinp0 (cos.f64 lampp)) (cos.f64 phipp)) (*.f64 (sin.f64 phipp) cosp0)))
(asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 (*.f64 (cos.f64 phipp) (cos.f64 lampp)))))
(asin.f64 (/.f64 (*.f64 (cos.f64 phipp) (sin.f64 lampp)) (cos.f64 phip)))
(asin.f64 (*.f64 (/.f64 (sin.f64 lampp) (cos.f64 phip)) (cos.f64 phipp)))
(asin.f64 (*.f64 (/.f64 (cos.f64 phipp) (cos.f64 phip)) (sin.f64 lampp)))
(*.f64 ecc (sin.f64 phip))
(*.f64 (-.f64 (+.f64 con (log.f64 (*.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip))) (/.f64 #s(literal 1 binary64) (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phip))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phip))))))) (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phip phip)) #s(literal 1/2 binary64))) (*.f64 ecc ecc))) (*.f64 (cos.f64 phip) rone_es)))
(*.f64 (*.f64 (-.f64 (+.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip)))) con) (*.f64 (log.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phip) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phip) ecc)))) hlf_e)) (*.f64 (cos.f64 phip) (-.f64 #s(literal 1 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phip phip)))) (*.f64 ecc ecc))))) rone_es)
(/.f64 (-.f64 K (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phip #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phip #s(literal 1/2 binary64)))))))) c)
(/.f64 (-.f64 K (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 #s(literal 1/2 binary64) (-.f64 phip (PI.f64))))) #s(literal -1 binary64)))) c)
(if.f64 (<=.f64 (-.f64 K (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip))))) #s(literal 2251799813685247/20282409603651670423947251286016 binary64)) (/.f64 (+.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 #s(literal 1/2 binary64) (-.f64 phip (PI.f64)))))) K) c) (if.f64 (<=.f64 (-.f64 K (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip))))) #s(literal 1/9007199254740992 binary64)) (/.f64 (*.f64 c (-.f64 K (log.f64 (tan.f64 (-.f64 (*.f64 phip #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (*.f64 c c)) (/.f64 (-.f64 K (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64)))))) c)))
(/.f64 (-.f64 K (log.f64 (tan.f64 (-.f64 (*.f64 phip #s(literal 1/2 binary64)) (*.f64 (PI.f64) #s(literal -1/4 binary64)))))) c)
(/.f64 (sin.f64 phi0) c)
(asin.f64 sinp0)
(/.f64 (*.f64 (sqrt.f64 one_es) k0) (-.f64 (*.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi0 phi0))) #s(literal 1/2 binary64)) (*.f64 ecc ecc)) #s(literal -1 binary64)))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es) rone_es)))
(-.f64 (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (/.f64 #s(literal 1 binary64) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (neg.f64 (log.f64 (/.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (fabs.f64 (-.f64 (*.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/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))))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0)))))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (log.f64 (neg.f64 (/.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (/.f64 (neg.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0))))) (*.f64 (-.f64 (log.f64 (fabs.f64 (/.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) #s(literal 1 binary64))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))))) (*.f64 (log.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 phi0) ecc)) (-.f64 (*.f64 (sin.f64 phi0) ecc) #s(literal 1 binary64)))) hlf_e)) 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) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64)) (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)))))))))
(-.f64 (log.f64 (-.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (cos.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64))) (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)))))))))
(-.f64 (log.f64 (neg.f64 (/.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0)))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0)))))) (*.f64 c (-.f64 (log.f64 (/.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64))) (PI.f64))) (neg.f64 (neg.f64 (neg.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0)))))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0)))) (neg.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0)))))) (*.f64 c (-.f64 (log.f64 (/.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64))) (PI.f64))) (cos.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64))) (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)))))))))
(-.f64 (log.f64 (neg.f64 (/.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip0))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phip0)))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0)))))) (*.f64 c (-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64)) (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)))))))))
(-.f64 (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip0) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (/.f64 (PI.f64) #s(literal 2 binary64))))))) (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi0) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi0) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0))))) c))
(-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (log.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64))) (PI.f64)) (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)))))))))
(-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0)))) (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi0) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi0) ecc))) (neg.f64 hlf_e)) (tan.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (PI.f64)) (PI.f64))) (PI.f64))))) c))
(-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0)))) (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi0) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi0) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0))))) c))
(*.f64 (/.f64 (-.f64 (pow.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))))) #s(literal 2 binary64)) (pow.f64 (neg.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))))) #s(literal 2 binary64))) (log.f64 (/.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64)))) (/.f64 #s(literal 1 binary64) (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))))))) kR)
(*.f64 (log.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))) kR)
(if.f64 (<=.f64 phipp #s(literal -6211381958656337/365375409332725729550921208179070754913983135744 binary64)) (*.f64 kR (log.f64 (tan.f64 (/.f64 (-.f64 (*.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) phipp) phipp) (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))) (*.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)))) (*.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))))) (if.f64 (<=.f64 phipp #s(literal 8958978968711217/279968092772225526319680285071055534765205687154331191862498637620473983897520118172609686658950889472 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) kR)) (*.f64 kR (log.f64 (tan.f64 (+.f64 (*.f64 phipp #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (PI.f64))))))))
(if.f64 (<=.f64 kR #s(literal 4776913109852041/20769187434139310514121985316880384 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))) kR)) (*.f64 kR (log.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 phipp #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))) (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64))))))))
(if.f64 (<=.f64 kR #s(literal 4776913109852041/20769187434139310514121985316880384 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))) kR)) (*.f64 kR (log.f64 (tan.f64 (/.f64 (-.f64 (*.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) phipp) phipp) (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))) (*.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64)))) (+.f64 (-.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) phipp) phipp) (*.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 phipp #s(literal 1/2 binary64)) #s(literal 1/4 binary64)) (PI.f64)))) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/16 binary64))))))))
(if.f64 (<=.f64 kR #s(literal 4776913109852041/20769187434139310514121985316880384 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phipp))) kR)) (*.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 (tan.f64 (/.f64 (-.f64 (*.f64 (*.f64 #s(literal 1/4 binary64) phipp) phipp) (*.f64 #s(literal 1/16 binary64) (*.f64 (PI.f64) (PI.f64)))) (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))))
(*.f64 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 (tan.f64 (+.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 13/4 binary64) (PI.f64))))))
(*.f64 kR lampp)
(neg.f64 (log.f64 (/.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 2 binary64))) #s(literal -1 binary64)) #s(literal 1/2 binary64)) (fabs.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (fabs.f64 (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(neg.f64 (log.f64 (/.f64 (*.f64 (+.f64 (cos.f64 (*.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) #s(literal 2 binary64))) #s(literal 1 binary64)) #s(literal 1/2 binary64)) (fabs.f64 (*.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi)))) (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))))
(neg.f64 (log.f64 (/.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64)))) (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phi) (PI.f64)) #s(literal -1/4 binary64))))))
(log.f64 (tan.f64 (+.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -9/4 binary64) (PI.f64)) (PI.f64))) (PI.f64))))
(log.f64 (tan.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -9/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)))))
(log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))
(asin.f64 (-.f64 (*.f64 (sin.f64 phip) cosp0) (*.f64 (*.f64 sinp0 (cos.f64 lamp)) (cos.f64 phip))))
(asin.f64 (-.f64 (*.f64 cosp0 (sin.f64 phip)) (*.f64 sinp0 (*.f64 (cos.f64 phip) (cos.f64 lamp)))))
(asin.f64 (-.f64 (*.f64 (sin.f64 phip) cosp0) (*.f64 (*.f64 (*.f64 (cos.f64 lamp) sinp0) (+.f64 phip phip)) #s(literal 1/2 binary64))))
(if.f64 (<=.f64 phi #s(literal -8892599035402499/1218164251424999885044172798484398538859528357199375940858488307151618586345803262808201883235251282403163114528926083522932396233150386755822248412039081677441409712494559128733848706936256706044099949184902297359210699740674359368218295451933620701603467350388034693385228573748989263872 binary64)) (/.f64 (-.f64 (*.f64 (pow.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) #s(literal 4 binary64)) (*.f64 (*.f64 (PI.f64) (PI.f64)) #s(literal 1/4 binary64))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (+.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64))) (*.f64 #s(literal 3 binary64) (PI.f64))) (PI.f64)) (PI.f64))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(if.f64 (<=.f64 phi #s(literal -8892599035402499/1218164251424999885044172798484398538859528357199375940858488307151618586345803262808201883235251282403163114528926083522932396233150386755822248412039081677441409712494559128733848706936256706044099949184902297359210699740674359368218295451933620701603467350388034693385228573748989263872 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (+.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64))) (*.f64 #s(literal 3 binary64) (PI.f64))) (PI.f64)) (PI.f64))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(if.f64 (<=.f64 phi #s(literal -8892599035402499/1218164251424999885044172798484398538859528357199375940858488307151618586345803262808201883235251282403163114528926083522932396233150386755822248412039081677441409712494559128733848706936256706044099949184902297359210699740674359368218295451933620701603467350388034693385228573748989263872 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64))) (*.f64 #s(literal 3 binary64) (PI.f64))) (PI.f64))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(if.f64 (<=.f64 phi #s(literal -8892599035402499/1218164251424999885044172798484398538859528357199375940858488307151618586345803262808201883235251282403163114528926083522932396233150386755822248412039081677441409712494559128733848706936256706044099949184902297359210699740674359368218295451933620701603467350388034693385228573748989263872 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64))) (*.f64 #s(literal 3 binary64) (PI.f64)))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(if.f64 (<=.f64 phi #s(literal -8892599035402499/1218164251424999885044172798484398538859528357199375940858488307151618586345803262808201883235251282403163114528926083522932396233150386755822248412039081677441409712494559128733848706936256706044099949184902297359210699740674359368218295451933620701603467350388034693385228573748989263872 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64))) (+.f64 (PI.f64) (PI.f64)))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(if.f64 (<=.f64 phi #s(literal -6635262925566485/3016028602530220424421062271578867838525126125874252890898307532480811172246863646685486789069132871854112712507200418582406770249989531589552959798883475511096677545853948017425342019529539590198394174553047792145043723038206885734470238620775335376572758801465612196212014773714996912064944617579072807182486470656 binary64)) (-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))) (-.f64 (*.f64 (atan.f64 (*.f64 (pow.f64 (*.f64 (tan.f64 (+.f64 (+.f64 (PI.f64) (PI.f64)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64))))) (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 phi) ecc)) (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64))) (neg.f64 hlf_e))) c) (exp.f64 K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64))))
(-.f64 (*.f64 (atan.f64 (exp.f64 (+.f64 K (*.f64 (log.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))
(-.f64 (*.f64 (atan.f64 (*.f64 (pow.f64 (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 phi) ecc)) (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64))) (neg.f64 hlf_e)) (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi)))) c) (exp.f64 K))) #s(literal 2 binary64)) (*.f64 #s(literal 1/2 binary64) (PI.f64)))
(asin.f64 (*.f64 (/.f64 (sin.f64 lamp) (cos.f64 phipp)) cp))
(asin.f64 (*.f64 (/.f64 cp (cos.f64 phipp)) (sin.f64 lamp)))
(-.f64 (log.f64 (fabs.f64 (-.f64 #s(literal 1 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi phi)))) (*.f64 ecc ecc))))) (-.f64 (log.f64 (fabs.f64 (-.f64 #s(literal 1 binary64) (pow.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64))) #s(literal 2 binary64))))) (log.f64 (fabs.f64 (-.f64 (-.f64 (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64)) #s(literal -1 binary64)) (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))))))))
(+.f64 (log.f64 (/.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) #s(literal 1 binary64))) (fabs.f64 (-.f64 #s(literal 1 binary64) (pow.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64))) #s(literal 2 binary64)))))) (log.f64 (fabs.f64 (-.f64 (-.f64 (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64)) #s(literal -1 binary64)) (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64))))))))
(log.f64 (/.f64 (fabs.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) #s(literal 1 binary64))) (/.f64 (fabs.f64 (-.f64 #s(literal 1 binary64) (pow.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64))) #s(literal 2 binary64)))) (fabs.f64 (-.f64 (-.f64 (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64)) #s(literal -1 binary64)) (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))))))))
(neg.f64 (log.f64 (/.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (/.f64 (-.f64 (*.f64 #s(literal 4 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64))) (*.f64 ecc ecc))) #s(literal 1 binary64)) (-.f64 (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64)) #s(literal -1 binary64)))) (fabs.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) #s(literal 1 binary64))))))
(+.f64 (log.f64 (fabs.f64 (-.f64 #s(literal 1 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi phi)))) (*.f64 ecc ecc))))) (log.f64 (pow.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64)) #s(literal -2 binary64))))
(neg.f64 (log.f64 (/.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) (-.f64 (*.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 2 binary64)) #s(literal 1 binary64))) (fabs.f64 (-.f64 (*.f64 (*.f64 ecc ecc) (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1/2 binary64)))) #s(literal 1 binary64))))))
(log.f64 (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi phi)))) (*.f64 ecc ecc))) (pow.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64)) #s(literal -2 binary64))))
(neg.f64 (log.f64 (fabs.f64 (/.f64 (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc)) (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal -1 binary64))))))
(log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi)))))
(if.f64 (<=.f64 phi #s(literal 6100 binary64)) (-.f64 (*.f64 (*.f64 phi phi) (+.f64 phi phi)) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64))))) (-.f64 (/.f64 (+.f64 phi phi) (-.f64 #s(literal 0 binary64) (*.f64 phi phi))) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64))))))
(if.f64 (<=.f64 phi #s(literal 6100 binary64)) (-.f64 (*.f64 (*.f64 phi phi) (+.f64 phi phi)) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64))))) (-.f64 (*.f64 (+.f64 phi phi) (/.f64 #s(literal 1 binary64) (-.f64 #s(literal 0 binary64) (*.f64 phi phi)))) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64))))))
(-.f64 (*.f64 (*.f64 phi phi) (+.f64 phi phi)) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64)))))
(-.f64 (*.f64 (+.f64 phi phi) phi) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64)))))
(-.f64 (+.f64 phi phi) (log.f64 (fabs.f64 (-.f64 (*.f64 (sin.f64 phi) ecc) #s(literal 1 binary64)))))
(/.f64 Esc (sqrt.f64 (/.f64 (+.f64 #s(literal 1 binary64) Up) (-.f64 #s(literal 1 binary64) Up))))
(*.f64 (sqrt.f64 (/.f64 (-.f64 Up #s(literal 1 binary64)) (-.f64 #s(literal -1 binary64) Up))) Esc)
(neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (cos.f64 (*.f64 BrA u)))))
(exp.f64 (neg.f64 (*.f64 BrA v)))
(if.f64 (<=.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) #s(literal -1000000000000000044140518902895287779286391397382581274563006173283444396083023609274483667691850832398819698877547611031397112968428705874685599733334034192471780653571870045215197739635249206690814463183771858052833032509915549602573975010166573043840478561173504 binary64)) (-.f64 (*.f64 (/.f64 Vp Tp) cosgam) (*.f64 (neg.f64 singam) (/.f64 Sp Tp))) (if.f64 (<=.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) #s(literal 9999999999999999860310597602564577717002641838126363875249660735883565852672743849064846414228960666786379280392654615393353172850252103336275952370615397010730691664689375178569039851073146339641623266071126720011020169553304018596457812688561947201171488461172921822139066929851282122002676667750021070848 binary64)) (/.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) Tp) (+.f64 (*.f64 Sp (/.f64 singam Tp)) (*.f64 (/.f64 cosgam Tp) Vp))))
(if.f64 (<=.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) #s(literal -1000000000000000038893577551088388431307372492952020133343023820076912942893848967630799656078777013873264603119412132913531706114094375616540183672212689403544345862626169435445664558076559462193222406635520 binary64)) (+.f64 (*.f64 Sp (/.f64 singam Tp)) (*.f64 (/.f64 cosgam Tp) Vp)) (if.f64 (<=.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) #s(literal 9999999999999999860310597602564577717002641838126363875249660735883565852672743849064846414228960666786379280392654615393353172850252103336275952370615397010730691664689375178569039851073146339641623266071126720011020169553304018596457812688561947201171488461172921822139066929851282122002676667750021070848 binary64)) (/.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) Tp) (+.f64 (*.f64 Sp (/.f64 singam Tp)) (*.f64 (/.f64 cosgam Tp) Vp))))
(/.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) Tp)
(*.f64 #s(literal 1/2 binary64) (+.f64 Qp (/.f64 #s(literal 1 binary64) Qp)))
(*.f64 #s(literal 1/2 binary64) (-.f64 Qp (/.f64 #s(literal 1 binary64) Qp)))
(*.f64 ArB (log.f64 (-.f64 (/.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.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 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)))))))))
(*.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 #s(literal -1/2 binary64) gamma0)))))) (log.f64 (fabs.f64 (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (sin.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64))) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))))))))
(if.f64 (<=.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0)))) #s(literal -944473296573929/9444732965739290427392 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 #s(literal -1/2 binary64) gamma0)))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))))) (+.f64 (*.f64 (log.f64 (*.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) ArB) (*.f64 (neg.f64 (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) ArB)))
(*.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 #s(literal -1/2 binary64) gamma0)))))) (log.f64 (fabs.f64 (sin.f64 (+.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))))))
(+.f64 (*.f64 (log.f64 (*.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) ArB) (*.f64 (neg.f64 (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) ArB))
(*.f64 ArB (log.f64 (/.f64 (-.f64 (*.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sqrt.f64 #s(literal 2 binary64))) #s(literal -1/2 binary64)) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))) (neg.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(*.f64 (log.f64 (*.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0))))) ArB)
(*.f64 ArB (log.f64 (/.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 #s(literal -1/2 binary64) gamma0)))) (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))
(if.f64 (<=.f64 (*.f64 ArB (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -10120113/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (log.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (PI.f64))) (PI.f64))) ArB)) (*.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 (*.f64 ArB (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -10120113/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (log.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (PI.f64))) (PI.f64))) ArB)) (*.f64 ArB (log.f64 (*.f64 (/.f64 #s(literal -1 binary64) (cos.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0)))) (sin.f64 (*.f64 (-.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))))
(if.f64 (<=.f64 (*.f64 ArB (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -10120113/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (log.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (+.f64 (*.f64 #s(literal 13/4 binary64) (PI.f64)) (PI.f64))) (PI.f64))) ArB)) (*.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/20282409603651670423947251286016 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) ArB)) (*.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/20282409603651670423947251286016 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0))) ArB)) (*.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 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0))))))
(*.f64 ArB (log.f64 (tan.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))
(*.f64 ArB (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 9/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0)))))
(*.f64 ArB (-.f64 (+.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (log.f64 (fabs.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))))))) (log.f64 (sqrt.f64 (*.f64 (+.f64 (cos.f64 (*.f64 (-.f64 (PI.f64) (+.f64 gamma0 gamma0)) #s(literal 1/2 binary64))) #s(literal 1 binary64)) #s(literal 1/2 binary64))))))
(*.f64 ArB (-.f64 (+.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (log.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)))))) (log.f64 (fabs.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 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 (sqrt.f64 (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (*.f64 #s(literal 2 binary64) (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64))))))))))
(*.f64 ArB (neg.f64 (log.f64 (/.f64 (fabs.f64 (cos.f64 (*.f64 (-.f64 (PI.f64) (+.f64 gamma0 gamma0)) #s(literal -1/4 binary64)))) (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 (/.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))))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))) ArB)
(*.f64 (log.f64 (*.f64 (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal -1/4 binary64)))) (-.f64 (cos.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 gamma0 #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 (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64))))))) (*.f64 ArB (log.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (sin.f64 (+.f64 (*.f64 (-.f64 (PI.f64) (+.f64 gamma0 gamma0)) #s(literal -1/4 binary64)) (*.f64 #s(literal 1/2 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 1125899906831365/1125899906842624 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64))))))) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) ArB)))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 1125899906831365/1125899906842624 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64))))))) (log.f64 (pow.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))) ArB)))
(*.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 (tan.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) ArB)
(if.f64 (<.f64 phi0 #s(literal 0 binary64)) (neg.f64 (fabs.f64 (*.f64 ArB (atan.f64 (/.f64 (sqrt.f64 (-.f64 (*.f64 D D) #s(literal 1 binary64))) (cos.f64 alpha_c)))))) (fabs.f64 (*.f64 ArB (atan.f64 (/.f64 (sqrt.f64 (-.f64 (*.f64 D D) #s(literal 1 binary64))) (cos.f64 alpha_c))))))
(-.f64 (*.f64 (*.f64 (-.f64 (cos.f64 (+.f64 phi0 phi0)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) es) #s(literal -1 binary64))
(/.f64 com (*.f64 (sqrt.f64 (*.f64 (fabs.f64 (/.f64 one_es (+.f64 (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es) one_es))) (-.f64 #s(literal 1 binary64) (*.f64 (*.f64 #s(literal 1/2 binary64) es) (-.f64 #s(literal 1 binary64) (cos.f64 (+.f64 phi0 phi0))))))) (cos.f64 phi0)))
(sqrt.f64 (-.f64 (/.f64 (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es) one_es) #s(literal -1 binary64)))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) one_es) es)))
(*.f64 #s(literal 1/2 binary64) (*.f64 ArB (-.f64 (*.f64 #s(literal 1/2 binary64) (log.f64 (fabs.f64 (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 U U)) (-.f64 #s(literal 1 binary64) (*.f64 U U)))))) (log.f64 (+.f64 (*.f64 U (+.f64 #s(literal 2 binary64) U)) #s(literal 1 binary64))))))
(*.f64 #s(literal 1/2 binary64) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (-.f64 (*.f64 U U) #s(literal 1 binary64)))) (log.f64 (+.f64 (*.f64 U (+.f64 #s(literal 2 binary64) U)) #s(literal 1 binary64))))))
(*.f64 #s(literal 1/2 binary64) (*.f64 ArB (neg.f64 (log.f64 (/.f64 (-.f64 (-.f64 #s(literal 1 binary64) U) (*.f64 (-.f64 U #s(literal 1 binary64)) U)) (*.f64 (-.f64 U #s(literal 1 binary64)) (-.f64 U #s(literal 1 binary64))))))))
(*.f64 #s(literal 1/2 binary64) (*.f64 ArB (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (-.f64 U #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) U))))))
(*.f64 #s(literal 1/2 binary64) (*.f64 ArB (neg.f64 (log.f64 (/.f64 (-.f64 #s(literal -1 binary64) U) (-.f64 U #s(literal 1 binary64)))))))
(*.f64 #s(literal 1/2 binary64) (*.f64 ArB (log.f64 (/.f64 (-.f64 #s(literal 1 binary64) U) (+.f64 #s(literal 1 binary64) U)))))
(*.f64 A lam)
(*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (/.f64 (*.f64 (pow.f64 (cos.f64 (*.f64 lam B)) #s(literal 3 binary64)) #s(literal 2 binary64)) (+.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) (*.f64 lam B)) (*.f64 lam B))) (sin.f64 (-.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) (*.f64 lam B)) (*.f64 lam B)))))))
(*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (sin.f64 (-.f64 (*.f64 (neg.f64 B) lam) (*.f64 (PI.f64) #s(literal -1/2 binary64))))))
(*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam))))
(-.f64 (*.f64 (-.f64 u u0) cosrot) (*.f64 v sinrot))
(+.f64 (*.f64 v cosrot) (*.f64 (-.f64 u u0) sinrot))
(*.f64 #s(literal 1/2 binary64) (+.f64 W (/.f64 #s(literal 1 binary64) W)))
(*.f64 #s(literal 1/2 binary64) (-.f64 W (/.f64 #s(literal 1 binary64) W)))
(if.f64 (<=.f64 (-.f64 (*.f64 S singam) (*.f64 (sin.f64 (*.f64 B lam)) cosgam)) #s(literal 5000000000000000233009035874103487842025429049746884307104902290093413906615431497863838561070978561605169882979927449326586308330034490456803110487463217202937150636836581109497436029475276191632298678857801213921774796800 binary64)) (/.f64 (-.f64 (*.f64 S singam) (*.f64 (sin.f64 (*.f64 B lam)) cosgam)) T) (-.f64 (*.f64 S (/.f64 singam T)) (*.f64 (/.f64 cosgam T) (sin.f64 (*.f64 lam B)))))
(/.f64 (-.f64 (*.f64 S singam) (*.f64 (sin.f64 (*.f64 B lam)) cosgam)) T)
(asin.f64 (/.f64 (sin.f64 phi0) alpha))
(/.f64 (*.f64 (*.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) #s(literal -27401669256311/20000000000000 binary64))) n0) k0) (sin.f64 #s(literal 27401669256311/20000000000000 binary64)))
(/.f64 (*.f64 (*.f64 n0 k0) (sin.f64 (+.f64 #s(literal 27401669256311/20000000000000 binary64) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (sin.f64 #s(literal 27401669256311/20000000000000 binary64)))
(/.f64 (*.f64 (*.f64 #s(literal 3 binary64) n0) (*.f64 k0 #s(literal 1/3 binary64))) (tan.f64 (+.f64 (+.f64 (-.f64 (PI.f64) (-.f64 #s(literal -27401669256311/20000000000000 binary64) (PI.f64))) (PI.f64)) (PI.f64))))
(*.f64 (/.f64 k0 (tan.f64 (+.f64 (+.f64 (-.f64 (PI.f64) (-.f64 #s(literal -27401669256311/20000000000000 binary64) (PI.f64))) (PI.f64)) (PI.f64)))) n0)
(*.f64 (*.f64 n0 k0) (tan.f64 (-.f64 #s(literal -27401669256311/20000000000000 binary64) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))
(*.f64 (/.f64 n0 (tan.f64 #s(literal 27401669256311/20000000000000 binary64))) k0)
(*.f64 (/.f64 k0 (tan.f64 #s(literal 27401669256311/20000000000000 binary64))) n0)
(/.f64 (sqrt.f64 (/.f64 (-.f64 (*.f64 es es) #s(literal 1 binary64)) (-.f64 #s(literal -1 binary64) es))) (-.f64 (*.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi0 phi0))) #s(literal 1/2 binary64)) es) #s(literal -1 binary64)))
(/.f64 (sqrt.f64 (-.f64 #s(literal 1 binary64) es)) (-.f64 (*.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 phi0 phi0))) #s(literal 1/2 binary64)) es) #s(literal -1 binary64)))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2251799802426249/2251799813685248 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 (*.f64 (+.f64 (/.f64 (/.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u0)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) (/.f64 (*.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2251799802426249/2251799813685248 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal +inf.0 binary64)) (/.f64 (*.f64 (+.f64 (/.f64 (/.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u0)) (/.f64 #s(literal 2 binary64) (sqrt.f64 #s(literal 2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) (/.f64 (*.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (tan.f64 (*.f64 #s(literal -1/4 binary64) (-.f64 (-.f64 (neg.f64 (PI.f64)) u0) u0))) alpha))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2251799802426249/2251799813685248 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal +inf.0 binary64)) (/.f64 (*.f64 (+.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) (/.f64 (*.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (tan.f64 (*.f64 #s(literal -1/4 binary64) (-.f64 (-.f64 (neg.f64 (PI.f64)) u0) u0))) alpha))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2251799802426249/2251799813685248 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 (*.f64 (+.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (+.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))))))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) (/.f64 (*.f64 (+.f64 (/.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 u0 #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 u0 #s(literal 1/2 binary64)))))) g) (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha))))
(if.f64 (<=.f64 alpha #s(literal 1170935903116329/4503599627370496 binary64)) (/.f64 (/.f64 #s(literal 1 binary64) (/.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (sin.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g))) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phi0 #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phi0 #s(literal 1/2 binary64)))))) alpha)) (if.f64 (<=.f64 alpha #s(literal 104999999999999998921267926997898560112690709440242490226450292851480670657516514628969468579740626803766033680988594369320572941999352591811666873294015856401631865208832 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (/.f64 (*.f64 (+.f64 (/.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 u0 #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 u0 #s(literal 1/2 binary64)))))) g) (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha))))
(if.f64 (<=.f64 alpha #s(literal 170 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)))) alpha) (*.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g))) (if.f64 (<=.f64 alpha #s(literal 104999999999999998921267926997898560112690709440242490226450292851480670657516514628969468579740626803766033680988594369320572941999352591811666873294015856401631865208832 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g)))
(if.f64 (<=.f64 alpha #s(literal 170 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)))) alpha) (*.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g))) (if.f64 (<=.f64 alpha #s(literal 104999999999999998921267926997898560112690709440242490226450292851480670657516514628969468579740626803766033680988594369320572941999352591811666873294015856401631865208832 binary64)) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) alpha)) (/.f64 (*.f64 (+.f64 (/.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 2 binary64) (*.f64 (sqrt.f64 #s(literal 2 binary64)) (cos.f64 (*.f64 u0 #s(literal 1/2 binary64)))))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 u0 #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 u0 #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 u0 #s(literal 1/2 binary64)))))) g) (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha))))
(if.f64 (<=.f64 alpha #s(literal 65999999999999997607243397563642606242302454441349835381681422336 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)))) alpha) (*.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g))) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g))
(if.f64 (<=.f64 alpha #s(literal 65999999999999997607243397563642606242302454441349835381681422336 binary64)) (*.f64 (*.f64 g (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))))) (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)))) (neg.f64 alpha))) (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g))
(if.f64 (<=.f64 (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) #s(literal -inf.0 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) #s(literal 200000000000000015842876501691535308251363838339942186816779868466887151795034205545089069114411529059504325666588836124813676426230104197677563914641752713707086241642983763505789334141041164451549418938435594261010114368138763297090749546488747114934452621501484084432923307385290752 binary64)) (/.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g) (/.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (cos.f64 (+.f64 (neg.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))))) (/.f64 (/.f64 (*.f64 (sin.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g) (sin.f64 (+.f64 (neg.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 1/2 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g) (/.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (cos.f64 (+.f64 (neg.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))))) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (/.f64 (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64))) (neg.f64 (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))))))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 1/2 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 (*.f64 (/.f64 (sin.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (/.f64 (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64))) (neg.f64 (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))))))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 1/2 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g) (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) #s(literal -1 binary64)))) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (/.f64 (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64))) (neg.f64 (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))))))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 1/2 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g) (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) #s(literal -1 binary64)))) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (/.f64 (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64))) (neg.f64 (cos.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha))))))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 1/2 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 #s(literal 1 binary64) (*.f64 (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g) (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) #s(literal -1 binary64)))) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g)))
(if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 1/2 binary64)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (if.f64 (<=.f64 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 2 binary64)) (/.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g)) (*.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g)))
(if.f64 (<=.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 1 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 (PI.f64) u0) u0) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (/.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) g)))
(if.f64 (<=.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 1 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 (PI.f64) u0) u0) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g))
(*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 u0 (+.f64 u0 (PI.f64))) #s(literal 1/4 binary64)))) g)
(pow.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 phi0)))) (*.f64 (*.f64 #s(literal 1/2 binary64) alpha) evar))
(sqrt.f64 (/.f64 (+.f64 (*.f64 (-.f64 es #s(literal 1 binary64)) (-.f64 es #s(literal 1 binary64))) (*.f64 (-.f64 #s(literal 1 binary64) es) (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es))) (*.f64 (-.f64 es #s(literal 1 binary64)) (-.f64 es #s(literal 1 binary64)))))
(sqrt.f64 (/.f64 (-.f64 (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es) (-.f64 es #s(literal 1 binary64))) (-.f64 #s(literal 1 binary64) es)))
(sqrt.f64 (+.f64 #s(literal 1 binary64) (/.f64 (*.f64 es (pow.f64 (cos.f64 phi0) #s(literal 4 binary64))) (-.f64 #s(literal 1 binary64) es))))
(-.f64 (-.f64 (-.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 4616990181/250000000000000000000 binary64) Xr) Yr) (-.f64 (+.f64 (*.f64 Xr #s(literal -4668270147/10000000000000000 binary64)) #s(literal 9827991/390625000 binary64)) (*.f64 #s(literal -149230739/1250000000000000 binary64) Yr))) (*.f64 (*.f64 #s(literal -304747143/200000000000000000000 binary64) (+.f64 Xr Yr)) (-.f64 Xr Yr))) (*.f64 (*.f64 #s(literal -881662847/200000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr))))) (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr)))) (-.f64 (*.f64 (*.f64 (*.f64 #s(literal -3689471323/250000000000000000000000000000000 binary64) (*.f64 Xr Yr)) (+.f64 Xr Yr)) (-.f64 Xr Yr)) (*.f64 (+.f64 (*.f64 (*.f64 (-.f64 (*.f64 Yr Yr) (*.f64 #s(literal 6 binary64) (*.f64 Xr Xr))) Yr) Yr) (*.f64 (*.f64 Xr Xr) (*.f64 Xr Xr))) #s(literal -4165541759/500000000000000000000000000000000 binary64))))
(-.f64 (-.f64 (-.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 4616990181/250000000000000000000 binary64) Xr) Yr) (-.f64 (+.f64 (*.f64 Xr #s(literal -4668270147/10000000000000000 binary64)) #s(literal 9827991/390625000 binary64)) (*.f64 #s(literal -149230739/1250000000000000 binary64) Yr))) (*.f64 (*.f64 #s(literal -304747143/200000000000000000000 binary64) (+.f64 Xr Yr)) (-.f64 Xr Yr))) (*.f64 (*.f64 #s(literal -881662847/200000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr))))) (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr)))) (-.f64 (*.f64 (*.f64 (*.f64 #s(literal -3689471323/250000000000000000000000000000000 binary64) (*.f64 Xr Yr)) (+.f64 Xr Yr)) (-.f64 Xr Yr)) (*.f64 (+.f64 (*.f64 (*.f64 Yr Yr) (+.f64 (*.f64 Yr Yr) (*.f64 #s(literal -6 binary64) (*.f64 Xr Xr)))) (*.f64 (*.f64 Xr Xr) (*.f64 Xr Xr))) #s(literal -4165541759/500000000000000000000000000000000 binary64))))
(-.f64 (-.f64 (-.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 4616990181/250000000000000000000 binary64) Xr) Yr) (-.f64 (+.f64 (*.f64 Xr #s(literal -4668270147/10000000000000000 binary64)) #s(literal 9827991/390625000 binary64)) (*.f64 #s(literal -149230739/1250000000000000 binary64) Yr))) (*.f64 (*.f64 #s(literal -304747143/200000000000000000000 binary64) (+.f64 Xr Yr)) (-.f64 Xr Yr))) (*.f64 (*.f64 #s(literal -881662847/200000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr))))) (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr)))) (-.f64 (*.f64 (*.f64 (*.f64 #s(literal -3689471323/250000000000000000000000000000000 binary64) (*.f64 Xr Yr)) (+.f64 Xr Yr)) (-.f64 Xr Yr)) (*.f64 (+.f64 (*.f64 Yr (*.f64 Yr (-.f64 (*.f64 Yr Yr) (*.f64 #s(literal 6 binary64) (*.f64 Xr Xr))))) (*.f64 (*.f64 (*.f64 Xr Xr) Xr) Xr)) #s(literal -4165541759/500000000000000000000000000000000 binary64))))
(-.f64 (-.f64 (-.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 4616990181/250000000000000000000 binary64) Xr) Yr) (-.f64 (+.f64 (*.f64 Xr #s(literal -4668270147/10000000000000000 binary64)) #s(literal 9827991/390625000 binary64)) (*.f64 #s(literal -149230739/1250000000000000 binary64) Yr))) (*.f64 #s(literal -304747143/200000000000000000000 binary64) (*.f64 (+.f64 Xr Yr) (-.f64 Xr Yr)))) (*.f64 (*.f64 #s(literal -881662847/200000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr))))) (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr)))) (-.f64 (*.f64 (*.f64 (*.f64 #s(literal -3689471323/250000000000000000000000000000000 binary64) (*.f64 Xr Yr)) (+.f64 Xr Yr)) (-.f64 Xr Yr)) (*.f64 (+.f64 (*.f64 (*.f64 Yr Yr) (+.f64 (*.f64 Yr Yr) (*.f64 #s(literal -6 binary64) (*.f64 Xr Xr)))) (*.f64 (*.f64 (*.f64 Xr Xr) Xr) Xr)) #s(literal -4165541759/500000000000000000000000000000000 binary64))))
(-.f64 (-.f64 (-.f64 (-.f64 #s(literal 2946529277/100000000000 binary64) (-.f64 (*.f64 #s(literal -149230739/1250000000000000 binary64) Xr) (*.f64 Yr (+.f64 #s(literal 4668270147/10000000000000000 binary64) (*.f64 #s(literal -304747143/100000000000000000000 binary64) Xr))))) (*.f64 (*.f64 #s(literal -4616990181/500000000000000000000 binary64) (+.f64 Yr Xr)) (-.f64 Xr Yr))) (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr))))) (-.f64 (*.f64 (*.f64 #s(literal 881662847/200000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr))) (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 Yr Yr) (+.f64 (*.f64 Yr Yr) (*.f64 #s(literal -6 binary64) (*.f64 Xr Xr)))) (pow.f64 Xr #s(literal 4 binary64))) #s(literal -3689471323/1000000000000000000000000000000000 binary64)) (*.f64 (*.f64 (*.f64 #s(literal -4165541759/125000000000000000000000000000000 binary64) Xr) Yr) (*.f64 (-.f64 Xr Yr) (+.f64 Yr Xr))))))
(-.f64 (-.f64 (-.f64 #s(literal 2946529277/100000000000 binary64) (+.f64 (-.f64 (*.f64 #s(literal -149230739/1250000000000000 binary64) Xr) (*.f64 (+.f64 (*.f64 #s(literal -304747143/100000000000000000000 binary64) Xr) #s(literal 4668270147/10000000000000000 binary64)) Yr)) (*.f64 (*.f64 #s(literal -4616990181/500000000000000000000 binary64) (+.f64 Yr Xr)) (-.f64 Xr Yr)))) (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr))))) (-.f64 (*.f64 (*.f64 #s(literal 881662847/200000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr))) (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 Yr Yr) (+.f64 (*.f64 Yr Yr) (*.f64 #s(literal -6 binary64) (*.f64 Xr Xr)))) (*.f64 (*.f64 (*.f64 Xr Xr) Xr) Xr)) #s(literal -3689471323/1000000000000000000000000000000000 binary64)) (*.f64 (*.f64 (*.f64 #s(literal -4165541759/125000000000000000000000000000000 binary64) Xr) Yr) (*.f64 (-.f64 Xr Yr) (+.f64 Yr Xr))))))
(asin.f64 (-.f64 (*.f64 (sin.f64 s) (cos.f64 ad)) (*.f64 (sin.f64 ad) (*.f64 (cos.f64 d) (cos.f64 s)))))
(*.f64 (-.f64 (atan.f64 (*.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) (pow.f64 (/.f64 rho0 rho) (/.f64 #s(literal 1 binary64) n)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (*.f64 (pow.f64 (/.f64 (*.f64 (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u)) (sin.f64 (*.f64 u #s(literal 1/2 binary64)))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 u #s(literal 1/2 binary64)))))) (/.f64 #s(literal 1 binary64) alpha)) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u)) (sin.f64 (*.f64 u #s(literal 1/2 binary64))))) (sin.f64 (+.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.f64 (PI.f64) #s(literal 1/2 binary64))))) (/.f64 #s(literal 1 binary64) alpha)) (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u)) (sin.f64 (*.f64 u #s(literal 1/2 binary64))))) (cos.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (/.f64 #s(literal 1 binary64) alpha)) (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u) (/.f64 (PI.f64) #s(literal 4 binary64)))) (/.f64 #s(literal 1 binary64) alpha))) (pow.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 fi))) (-.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 fi)))) (*.f64 #s(literal 1/2 binary64) evar)))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 0 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (/.f64 #s(literal 1 binary64) (pow.f64 k (/.f64 #s(literal 1 binary64) alpha))) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u) (/.f64 (PI.f64) #s(literal 4 binary64)))) (/.f64 #s(literal 1 binary64) alpha))) (pow.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 fi))) (-.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 fi)))) (*.f64 #s(literal 1/2 binary64) evar)))) (/.f64 (PI.f64) #s(literal 4 binary64)))) (*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (*.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)) (PI.f64)) (PI.f64))) (/.f64 #s(literal 1 binary64) alpha)) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)))
(if.f64 (<=.f64 (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u) (/.f64 (PI.f64) #s(literal 4 binary64)))) (/.f64 #s(literal 1 binary64) alpha))) (pow.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 fi))) (-.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 fi)))) (*.f64 #s(literal 1/2 binary64) evar)))) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 0 binary64)) (*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64))) (/.f64 #s(literal 1 binary64) alpha)) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (*.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)) (PI.f64)) (PI.f64))) (/.f64 #s(literal 1 binary64) alpha)) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)))
(*.f64 (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 (/.f64 (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 fi) evar))) (*.f64 evar #s(literal 1/2 binary64))) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))) (pow.f64 (tan.f64 (+.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64))) (/.f64 #s(literal 1 binary64) alpha)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 fi) evar)) (-.f64 (*.f64 (sin.f64 fi) evar) #s(literal 1 binary64))) (*.f64 evar #s(literal 1/2 binary64))) (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64))) (/.f64 #s(literal 1 binary64) alpha)) (pow.f64 k (/.f64 #s(literal -1 binary64) alpha))))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64))
(asin.f64 (/.f64 (*.f64 (cos.f64 s) (sin.f64 d)) (cos.f64 u)))
(asin.f64 (*.f64 (/.f64 (sin.f64 d) (cos.f64 u)) (cos.f64 s)))
(asin.f64 (*.f64 (/.f64 (cos.f64 s) (cos.f64 u)) (sin.f64 d)))
(if.f64 (<=.f64 phi #s(literal -7365128867766319/409173825987017733751648712103449894027080255755383098685411421012016724550584319360408761540738019643860835515945008876152157068235674131666065948672 binary64)) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (if.f64 (<=.f64 phi #s(literal 5363123171977039/107262463439540776796592199985646769019834926564739147021788491549774112240588375814414994385335227421520254865491888406830031062495572559571469192048672768 binary64)) (*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi)) (sin.f64 (*.f64 phi #s(literal 1/2 binary64))))) (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) alpha) (/.f64 k gfi))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)) (PI.f64)) (PI.f64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(if.f64 (<=.f64 phi #s(literal -7365128867766319/409173825987017733751648712103449894027080255755383098685411421012016724550584319360408761540738019643860835515945008876152157068235674131666065948672 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))) (if.f64 (<=.f64 phi #s(literal 5363123171977039/107262463439540776796592199985646769019834926564739147021788491549774112240588375814414994385335227421520254865491888406830031062495572559571469192048672768 binary64)) (*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi)) (sin.f64 (*.f64 phi #s(literal 1/2 binary64))))) (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) alpha) (/.f64 k gfi))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)) (PI.f64)) (PI.f64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))
(if.f64 (<=.f64 phi #s(literal -7365128867766319/409173825987017733751648712103449894027080255755383098685411421012016724550584319360408761540738019643860835515945008876152157068235674131666065948672 binary64)) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 k (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 phi #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phi #s(literal 1/2 binary64))))) (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (sin.f64 (*.f64 phi #s(literal 1/2 binary64)))) (cos.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phi #s(literal 1/2 binary64)))))) alpha)) gfi)) (/.f64 (PI.f64) #s(literal 4 binary64)))))
(if.f64 (<=.f64 phi #s(literal 1320736327839163/66036816391958155794042473114564575814855915754515190265995088628334237819590607434302156985957148486834299374157141084860293444937835692156317372307867447113548335499875635521091439709553305131628029113526398550518503536593945217354799612345879559561192322612217249792 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)) (PI.f64)) (PI.f64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))
(*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha) k) gfi)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))
(*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (/.f64 k gfi) (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi phi) (PI.f64)) #s(literal 1/4 binary64))) alpha))) (*.f64 #s(literal 1/4 binary64) (PI.f64))))
(asin.f64 (+.f64 (*.f64 (cos.f64 ad) (sin.f64 u)) (*.f64 (*.f64 (cos.f64 deltav) (sin.f64 ad)) (cos.f64 u))))
(/.f64 (*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) rho0) (pow.f64 (+.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 s #s(literal 1/2 binary64)))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) s))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 #s(literal -1/2 binary64) s))))) (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 s #s(literal 1/2 binary64)))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) s))) (*.f64 (sin.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (sin.f64 (*.f64 #s(literal -1/2 binary64) s)))))) n))
(*.f64 (*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) s)) (sin.f64 (*.f64 s #s(literal 1/2 binary64))))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) s)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (sin.f64 (*.f64 s #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (neg.f64 n))) rho0)
(*.f64 (*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) (pow.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) s)) (sin.f64 (*.f64 s #s(literal 1/2 binary64))))) (cos.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (neg.f64 n))) rho0)
(*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 s s) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 n)) rho0))
(pow.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 phi))) (-.f64 #s(literal 1 binary64) (*.f64 evar (sin.f64 phi)))) (*.f64 (*.f64 #s(literal 1/2 binary64) alpha) evar))
(asin.f64 (/.f64 (*.f64 (cos.f64 u) (sin.f64 deltav)) (cos.f64 s)))
(asin.f64 (*.f64 (/.f64 (sin.f64 deltav) (cos.f64 s)) (cos.f64 u)))
(asin.f64 (*.f64 (/.f64 (cos.f64 u) (cos.f64 s)) (sin.f64 deltav)))
(cosh.f64 (asinh.f64 (sinh.f64 Ce)))
(sqrt.f64 (+.f64 (pow.f64 (sinh.f64 Ce) #s(literal 2 binary64)) #s(literal 1 binary64)))
(*.f64 (exp.f64 (*.f64 #s(literal -2 binary64) Ce)) #s(literal 1/2 binary64))
(exp.f64 (+.f64 Ce Ce))
(cosh.f64 (+.f64 Ce Ce))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (pow.f64 (sin.f64 Cn) #s(literal 2 binary64)) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (pow.f64 (cos.f64 lam) #s(literal 2 binary64))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (*.f64 (sin.f64 Cn) (sin.f64 Cn)) (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 #s(literal 1/2 binary64) (-.f64 (*.f64 (cos.f64 (+.f64 lam lam)) #s(literal 1/2 binary64)) (pow.f64 (*.f64 (+.f64 (cos.f64 (+.f64 lam lam)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) #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))) (pow.f64 (*.f64 (+.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64))))))
(/.f64 (*.f64 (sin.f64 lam) (cos.f64 Cn)) (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 Cn Cn)))) (*.f64 (cos.f64 lam) (cos.f64 lam))))))
(/.f64 (*.f64 (sin.f64 lam) (cos.f64 Cn)) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (pow.f64 (*.f64 (cos.f64 lam) (cos.f64 Cn)) #s(literal 2 binary64)))))
(*.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))))))))
(if.f64 (<=.f64 Cn #s(literal 8669969714643489/36893488147419103232 binary64)) (/.f64 (*.f64 (sin.f64 lam) (cos.f64 lam)) (sqrt.f64 (-.f64 (pow.f64 (cos.f64 lam) #s(literal 4 binary64)) (-.f64 (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64)) #s(literal 1/2 binary64))))) (*.f64 (cos.f64 Cn) (*.f64 (sin.f64 lam) (/.f64 #s(literal 1 binary64) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 (cos.f64 (+.f64 lam lam)) #s(literal 1/2 binary64))) (pow.f64 (*.f64 (-.f64 (cos.f64 (+.f64 lam lam)) #s(literal -1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64))))))))
(if.f64 (<=.f64 Cn #s(literal 8669969714643489/36893488147419103232 binary64)) (/.f64 (*.f64 (sin.f64 lam) (cos.f64 lam)) (sqrt.f64 (-.f64 (pow.f64 (cos.f64 lam) #s(literal 4 binary64)) (-.f64 (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64)) #s(literal 1/2 binary64))))) (*.f64 (cos.f64 Cn) (/.f64 (sin.f64 lam) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (-.f64 (*.f64 (cos.f64 (+.f64 lam lam)) #s(literal 1/2 binary64)) (pow.f64 (*.f64 (+.f64 (cos.f64 (+.f64 lam lam)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64))))))))
(if.f64 (<=.f64 Cn #s(literal 8116567392432203/73786976294838206464 binary64)) (/.f64 (*.f64 (sin.f64 lam) (cos.f64 lam)) (sqrt.f64 (-.f64 (pow.f64 (cos.f64 lam) #s(literal 4 binary64)) (-.f64 (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64)) #s(literal 1/2 binary64))))) (*.f64 (cos.f64 Cn) (/.f64 (sin.f64 lam) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (-.f64 (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64)) (pow.f64 (*.f64 (+.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64))))))))
(/.f64 (*.f64 (sin.f64 lam) (cos.f64 lam)) (sqrt.f64 (-.f64 (pow.f64 (cos.f64 lam) #s(literal 4 binary64)) (-.f64 (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64)) #s(literal 1/2 binary64)))))
(*.f64 (sin.f64 lam) (/.f64 (cos.f64 lam) (sqrt.f64 (-.f64 (pow.f64 (cos.f64 lam) #s(literal 4 binary64)) (-.f64 (*.f64 (cos.f64 (+.f64 Cn Cn)) #s(literal 1/2 binary64)) #s(literal 1/2 binary64))))))
(atan2.f64 (sin.f64 Cn) (*.f64 (cos.f64 Cn) (cos.f64 lam)))
(atan2.f64 (sin.f64 Cn) (-.f64 (*.f64 (cos.f64 (+.f64 lam lam)) #s(literal 1/2 binary64)) #s(literal -1/2 binary64)))
(*.f64 (+.f64 cosh_arg_i cosh_arg_i) cos_arg_r)
(+.f64 cos_arg_r cos_arg_r)
(*.f64 (*.f64 #s(literal -2 binary64) sin_arg_r) sinh_arg_i)
(if.f64 (<=.f64 hr #s(literal 3602879701896397/18014398509481984 binary64)) (-.f64 (*.f64 (*.f64 hr sin_arg_r) cosh_arg_i) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi sinh_arg_i) cos_arg_r)))
(if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 4999999999999999916795901115958608572801861375087352681835038072302342087505062772657389384729693708756186941721815525337672536740822868667327551851630428160 binary64)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)))
(if.f64 (<=.f64 sin_arg_r #s(literal -56000000000000001560727714672506550224046618630900641843237521303187488768 binary64)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) (-.f64 (*.f64 (*.f64 hr sin_arg_r) cosh_arg_i) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)))
(-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi))
(if.f64 (<=.f64 cosh_arg_i #s(literal -5483051187099521/14059105607947488696282932836518693308967803494693489478439861164411992439598399594747002144074658928593502845729752797260025831423419686528151609940203637047296 binary64)) (+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 hr sinh_arg_i) cos_arg_r)) (+.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hi) (*.f64 (*.f64 hr cos_arg_r) sinh_arg_i)))
(if.f64 (<=.f64 sinh_arg_i #s(literal 6338253001141147/1267650600228229401496703205376 binary64)) (+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 hr sinh_arg_i) cos_arg_r)) (+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)))
(+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr))
(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 (*.f64 (-.f64 (cos.f64 (+.f64 phi phi)) #s(literal 1 binary64)) #s(literal 1/2 binary64)) es) #s(literal -1 binary64))
(*.f64 (/.f64 x k0) (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))))