(*.f64 esp (-.f64 angle phi0))
(if.f64 (<=.f64 (cos.f64 (+.f64 phi0 (/.f64 y esp))) #s(literal -5404319552844595/18014398509481984 binary64)) (atan2.f64 g (+.f64 (*.f64 (sin.f64 phi0) (sin.f64 (+.f64 (/.f64 y esp) (PI.f64)))) (*.f64 (cos.f64 (/.f64 y esp)) (cos.f64 phi0)))) (atan2.f64 g (-.f64 (*.f64 (sin.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (/.f64 y esp))) (cos.f64 phi0)) (*.f64 (sin.f64 (/.f64 y esp)) (sin.f64 phi0)))))
(if.f64 (<=.f64 (/.f64 y esp) #s(literal -1000000000000000019884624838656 binary64)) (atan2.f64 g (sin.f64 (-.f64 (/.f64 (neg.f64 y) esp) (-.f64 phi0 (*.f64 #s(literal 1/2 binary64) (PI.f64)))))) (atan2.f64 g (-.f64 (*.f64 (cos.f64 (/.f64 y esp)) (cos.f64 phi0)) (*.f64 (sin.f64 (/.f64 y esp)) (sin.f64 phi0)))))
(if.f64 (<=.f64 (/.f64 y esp) #s(literal -1000000000000000019884624838656 binary64)) (atan2.f64 g (sin.f64 (-.f64 (/.f64 (neg.f64 y) esp) (-.f64 phi0 (*.f64 #s(literal 1/2 binary64) (PI.f64)))))) (atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp)))))
(atan2.f64 g (cos.f64 (+.f64 phi0 (*.f64 (/.f64 #s(literal 1 binary64) esp) y))))
(atan2.f64 g (cos.f64 (+.f64 phi0 (/.f64 y esp))))
(asin.f64 (/.f64 (fabs.f64 (+.f64 (*.f64 (sin.f64 phi0) (cos.f64 (/.f64 y esp))) (*.f64 (cos.f64 phi0) (sin.f64 (/.f64 y esp))))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))))
(asin.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (/.f64 y esp) phi0))) (sqrt.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 g g))))))
(sinh.f64 (/.f64 x esp))
(+.f64 (*.f64 (atan.f64 (/.f64 #s(literal 1 binary64) (sqrt.f64 (exp.f64 (*.f64 (/.f64 (neg.f64 y) kR) #s(literal 2 binary64)))))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))
(+.f64 (*.f64 (atan.f64 (/.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 y) kR)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal -1/2 binary64)))
(+.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 cosp0 (sin.f64 phipp)) (*.f64 (*.f64 (cos.f64 phipp) sinp0) (cos.f64 lampp))))
(asin.f64 (+.f64 (*.f64 cosp0 (sin.f64 phipp)) (*.f64 sinp0 (*.f64 (cos.f64 phipp) (cos.f64 lampp)))))
(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 (-.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))) (-.f64 (+.f64 (/.f64 (+.f64 (pow.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) phip)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip)))))) #s(literal 3 binary64)) (pow.f64 (log.f64 (/.f64 #s(literal 1 binary64) (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip)))))) #s(literal 3 binary64))) (+.f64 (*.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) phip)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip)))))) (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) phip)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip))))))) (-.f64 (*.f64 (log.f64 (/.f64 #s(literal 1 binary64) (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip)))))) (log.f64 (/.f64 #s(literal 1 binary64) (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip))))))) (*.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) phip)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip)))))) (log.f64 (/.f64 #s(literal 1 binary64) (fabs.f64 (cos.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))) rone_es) (cos.f64 phip))
(*.f64 (*.f64 (*.f64 (-.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))) (-.f64 (+.f64 (log.f64 (*.f64 (/.f64 #s(literal 1 binary64) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phip)) (sin.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))) rone_es) (cos.f64 phip))
(*.f64 (*.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 (+.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) phip)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phip)))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip))))) (-.f64 con (*.f64 (log.f64 (/.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 phip) ecc)) (-.f64 (*.f64 (sin.f64 phip) ecc) #s(literal 1 binary64)))) hlf_e))) rone_es)) (cos.f64 phip))
(*.f64 (-.f64 (+.f64 con (log.f64 (tan.f64 (/.f64 (-.f64 (*.f64 #s(literal 1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (*.f64 (*.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64))) (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64)))) (+.f64 (*.f64 #s(literal 1/16 binary64) (*.f64 (PI.f64) (PI.f64))) (+.f64 (*.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64))) (*.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64))))))))) (*.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 #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))) (-.f64 (+.f64 (log.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phip) (PI.f64))) (PI.f64)) (PI.f64)))) 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))) rone_es) (cos.f64 phip))
(*.f64 (*.f64 (*.f64 (-.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))) (-.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))) rone_es) (sin.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) phip)))
(*.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 (*.f64 rone_es (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)))))
(*.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 (*.f64 (-.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))) (-.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))) rone_es) (cos.f64 phip))
(/.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 K (log.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (+.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phip #s(literal 1/2 binary64))) (*.f64 (PI.f64) #s(literal 1/2 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 4506751413931401/40564819207303340847894502572032 binary64)) (*.f64 c (/.f64 (-.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 (log.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (+.f64 (-.f64 (*.f64 (-.f64 phip (PI.f64)) #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))) 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 2251799813685247/20282409603651670423947251286016 binary64)) (/.f64 (-.f64 K (log.f64 (/.f64 #s(literal 1 binary64) (tan.f64 (+.f64 (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (*.f64 phip #s(literal 1/2 binary64))) (*.f64 (PI.f64) #s(literal 1/2 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 2251799865151545/20282409603651670423947251286016 binary64)) (*.f64 c (/.f64 (-.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 phip #s(literal 1/2 binary64)) (*.f64 (PI.f64) #s(literal -1/4 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 2251799865151545/20282409603651670423947251286016 binary64)) (*.f64 c (/.f64 (-.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 phip #s(literal 1/2 binary64)) (*.f64 (PI.f64) #s(literal -1/4 binary64)))))) 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 es (*.f64 (pow.f64 (*.f64 (-.f64 (cos.f64 (+.f64 phi0 phi0)) #s(literal -1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64)) rone_es))))
(sqrt.f64 (-.f64 (*.f64 (*.f64 rone_es (pow.f64 (cos.f64 phi0) #s(literal 4 binary64))) es) #s(literal -1 binary64)))
(-.f64 (log.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) phip0)))) (*.f64 c (-.f64 (neg.f64 (log.f64 (/.f64 (fabs.f64 (sin.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))) (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) phi0) (*.f64 #s(literal 1/4 binary64) (PI.f64)))))))) (*.f64 hlf_e (log.f64 (/.f64 (+.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0))) (-.f64 #s(literal 1 binary64) (*.f64 ecc (sin.f64 phi0)))))))))
(-.f64 (log.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phip0)))) (*.f64 (-.f64 (log.f64 (fabs.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi0))) (sin.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.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 (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 kR (/.f64 (+.f64 (*.f64 (log.f64 (fabs.f64 (*.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (log.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 (log.f64 (fabs.f64 (*.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (neg.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))))) (+.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))))) (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))))))
(*.f64 kR (/.f64 (-.f64 (*.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))))) (log.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))))) (+.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))))) (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phipp #s(literal 1/2 binary64)))))))))
(*.f64 kR (*.f64 (log.f64 (fabs.f64 (*.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (/.f64 (log.f64 (fabs.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) (log.f64 (fabs.f64 (*.f64 (sin.f64 (*.f64 (-.f64 (*.f64 #s(literal -2 binary64) phipp) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (*.f64 #s(literal 1/4 binary64) (PI.f64))))))))))
(*.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)
(*.f64 kR (log.f64 (tan.f64 (+.f64 (-.f64 (*.f64 phipp #s(literal 1/2 binary64)) (-.f64 (-.f64 (-.f64 (-.f64 (*.f64 #s(literal -9/4 binary64) (PI.f64)) (PI.f64)) (PI.f64)) (PI.f64)) (PI.f64))) (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 -9/4 binary64) (PI.f64))))))
(*.f64 kR lampp)
(-.f64 (+.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (log.f64 (fabs.f64 (-.f64 (sin.f64 (*.f64 #s(literal -1/2 binary64) phi)) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))
(+.f64 (log.f64 (fabs.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi))))) (-.f64 (log.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(-.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) phi)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))
(-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 phi #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 binary64)))))))
(log.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 phi #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 phi #s(literal 1/2 binary64))))))
(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 cosp0 (sin.f64 phip)) (*.f64 sinp0 (*.f64 (cos.f64 phip) (cos.f64 lamp)))))
(if.f64 (<=.f64 phi #s(literal 3852638599766875/40131652080904949243476790488282231640246122763238325954424140190648896440865179612073261537762363061729301215028215161995082338334532195000669973530974432754174985283877903733762083113741475809259744657408 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 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 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 4 binary64)) (PI.f64)))) (/.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 (*.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)) #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 (*.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 (PI.f64) #s(literal -1/2 binary64)))))
(if.f64 (<=.f64 phi #s(literal 4975428362386363/76545051729020975577310162521900618820659871603466655644272117978380005723696097587725184512638784526308634214455061267843403507870735540391292521535824647434568377082591826884769598224146796816367616 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 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 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) phi))))) c)))) #s(literal 4 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 4975428362386363/76545051729020975577310162521900618820659871603466655644272117978380005723696097587725184512638784526308634214455061267843403507870735540391292521535824647434568377082591826884769598224146796816367616 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 3827252586451049/38272525864510487788655081260950309410329935801733327822136058989190002861848048793862592256319392263154317107227530633921701753935367770195646260767912323717284188541295913442384799112073398408183808 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 (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 (PI.f64) (PI.f64)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (*.f64 #s(literal -1/2 binary64) phi) (PI.f64)))))) 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 cp (cos.f64 phipp)) (sin.f64 lamp)))
(-.f64 (log.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 phi) ecc)))) (*.f64 #s(literal 1/2 binary64) (log.f64 (pow.f64 (-.f64 #s(literal 1 binary64) (*.f64 (sin.f64 phi) ecc)) #s(literal 2 binary64)))))
(-.f64 (log.f64 (fabs.f64 (-.f64 #s(literal -1 binary64) (*.f64 (sin.f64 phi) ecc)))) (log.f64 (fabs.f64 (-.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)))))
(*.f64 (sqrt.f64 (fabs.f64 (/.f64 (-.f64 #s(literal 1 binary64) Up) (-.f64 #s(literal -1 binary64) Up)))) Esc)
(*.f64 (neg.f64 rB) (atan2.f64 (-.f64 (*.f64 cosgam Sp) (*.f64 singam Vp)) (*.f64 (+.f64 (*.f64 (-.f64 (cos.f64 (*.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 u BrA)) #s(literal 2 binary64))) #s(literal 1 binary64)) (/.f64 #s(literal -1/2 binary64) (cos.f64 (*.f64 u BrA)))) (cos.f64 (*.f64 u BrA))) #s(literal 1/2 binary64))))
(*.f64 (neg.f64 rB) (atan2.f64 (-.f64 (*.f64 cosgam Sp) (*.f64 singam Vp)) (*.f64 (-.f64 (cos.f64 (*.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 u BrA)) #s(literal 2 binary64))) #s(literal 1 binary64)) (/.f64 #s(literal -1/2 binary64) (cos.f64 (*.f64 u BrA))))))
(neg.f64 (*.f64 rB (atan2.f64 (-.f64 (*.f64 Sp cosgam) (*.f64 Vp singam)) (sin.f64 (-.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 u BrA))))))
(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 Tp #s(literal 46000000000000001107925088917691237645105463841205082554960431831992090594522407266654667193450496 binary64)) (/.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) Tp) (+.f64 (*.f64 (/.f64 cosgam Tp) Vp) (*.f64 (/.f64 singam Tp) Sp)))
(if.f64 (<=.f64 Tp #s(literal 3602879701896397/72057594037927936 binary64)) (/.f64 (+.f64 (*.f64 Vp cosgam) (*.f64 Sp singam)) Tp) (+.f64 (*.f64 (/.f64 Vp Tp) cosgam) (*.f64 (/.f64 singam Tp) Sp)))
(/.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 (-.f64 Qp #s(literal 1 binary64)) (/.f64 (-.f64 #s(literal 1 binary64) Qp) 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 (log.f64 (/.f64 #s(literal 1 binary64) (/.f64 (-.f64 (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (cos.f64 (*.f64 (PI.f64) #s(literal -3/4 binary64)))) (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 (PI.f64) #s(literal -3/4 binary64))))) (*.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))))))))
(*.f64 ArB (neg.f64 (log.f64 (/.f64 (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.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)) (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)))))))
(*.f64 ArB (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (/.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (cos.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))))
(if.f64 (<=.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 9007199074597007/9007199254740992 binary64)) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))) (log.f64 (fabs.f64 (cos.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))) (*.f64 ArB (log.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (+.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))))
(if.f64 (<=.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 9007199074597007/9007199254740992 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (*.f64 ArB (log.f64 (/.f64 (sin.f64 (*.f64 (-.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))) (sin.f64 (+.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 gamma0 #s(literal 1/2 binary64)))))))))
(if.f64 (<=.f64 (tan.f64 (+.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 4503599627370495/4503599627370496 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.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 4503599627370495/4503599627370496 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))))) (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 gamma0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))))
(*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (-.f64 (*.f64 #s(literal 5/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) gamma0))))))
(*.f64 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 (log.f64 (fabs.f64 (-.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (/.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))) (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))))) (log.f64 (fabs.f64 (*.f64 (/.f64 (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))) (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (*.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 ArB (log.f64 (-.f64 (/.f64 #s(literal 1 binary64) (*.f64 (/.f64 #s(literal -1 binary64) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))) (/.f64 (*.f64 (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))) (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64))) (cos.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))))
(*.f64 ArB (log.f64 (/.f64 (-.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0))) (*.f64 (*.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 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))))))
(*.f64 ArB (log.f64 (*.f64 (/.f64 #s(literal 1 binary64) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64))))))))
(*.f64 ArB (log.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal -1/2 binary64)) (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (/.f64 #s(literal -1 binary64) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64))))))))
(*.f64 (log.f64 (/.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (-.f64 (cos.f64 (*.f64 #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 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (/.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))))) ArB)
(*.f64 ArB (log.f64 (*.f64 (/.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) gamma0)) (sin.f64 (*.f64 gamma0 #s(literal 1/2 binary64)))) (cos.f64 (*.f64 (+.f64 (*.f64 #s(literal -2 binary64) gamma0) (PI.f64)) #s(literal 1/4 binary64)))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)))))
(if.f64 (<=.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))) #s(literal 9007199254740987/9007199254740992 binary64)) (*.f64 ArB (log.f64 (fabs.f64 (tan.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64)))))) (*.f64 ArB (log.f64 (*.f64 (/.f64 #s(literal -1 binary64) (neg.f64 (cos.f64 (/.f64 (-.f64 (*.f64 #s(literal 1/64 binary64) (*.f64 (*.f64 (PI.f64) (PI.f64)) (PI.f64))) (*.f64 (*.f64 (-.f64 (PI.f64) (*.f64 gamma0 #s(literal 1/2 binary64))) (-.f64 (PI.f64) (*.f64 gamma0 #s(literal 1/2 binary64)))) (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64)))) (+.f64 (*.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64)) (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (-.f64 (PI.f64) (*.f64 gamma0 #s(literal 1/2 binary64))))) (*.f64 #s(literal 1/16 binary64) (*.f64 (PI.f64) (PI.f64)))))))) (sin.f64 (-.f64 (-.f64 (-.f64 (*.f64 gamma0 #s(literal 1/2 binary64)) (PI.f64)) (PI.f64)) (*.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 -20240225/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) ArB)) (*.f64 ArB (-.f64 (log.f64 (fabs.f64 (sin.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))))) (log.f64 (fabs.f64 (cos.f64 (*.f64 (+.f64 (*.f64 gamma0 #s(literal -2 binary64)) (PI.f64)) #s(literal -1/4 binary64))))))))
(if.f64 (<=.f64 (*.f64 ArB (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -20240225/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) ArB)) (*.f64 (log.f64 (fabs.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.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) gamma0)))))) ArB))
(if.f64 (<=.f64 (*.f64 ArB (log.f64 (tan.f64 (-.f64 (/.f64 (PI.f64) #s(literal 4 binary64)) (*.f64 #s(literal 1/2 binary64) gamma0))))) #s(literal -20240225/202402253307310618352495346718917307049556649764142118356901358027430339567995346891960383701437124495187077864316811911389808737385793476867013399940738509921517424276566361364466907742093216341239767678472745068562007483424692698618103355649159556340810056512358769552333414615230502532186327508646006263307707741093494784 binary64)) (log.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal -1/2 binary64) gamma0) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) 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 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 (/.f64 com (cos.f64 phi0)) (sqrt.f64 (/.f64 (-.f64 (/.f64 #s(literal 1 binary64) (/.f64 one_es (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es))) #s(literal -1 binary64)) (-.f64 #s(literal 1 binary64) (*.f64 (*.f64 (-.f64 #s(literal 1 binary64) (cos.f64 (+.f64 phi0 phi0))) #s(literal 1/2 binary64)) es)))))
(*.f64 (/.f64 com (cos.f64 phi0)) (sqrt.f64 (/.f64 (fabs.f64 (+.f64 (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es) one_es)) (*.f64 (-.f64 #s(literal 1 binary64) (*.f64 (*.f64 #s(literal 1/2 binary64) es) (-.f64 #s(literal 1 binary64) (cos.f64 (+.f64 phi0 phi0))))) (fabs.f64 one_es)))))
(sqrt.f64 (/.f64 (+.f64 one_es (*.f64 (pow.f64 (cos.f64 phi0) #s(literal 4 binary64)) es)) one_es))
(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 (log.f64 (/.f64 (-.f64 #s(literal 1 binary64) U) (+.f64 #s(literal 1 binary64) U)))))
(*.f64 A lam)
(if.f64 (<=.f64 B #s(literal -150000000000000005706058209634146259469505575867375116725271140279175986656890458835150470724443609261810727281002205504459950368998149529861458487918866276066675366113915585151323551800165366268582300726006810521341414866944 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 (sin.f64 (-.f64 (PI.f64) (*.f64 lam B))) singam) (*.f64 cosgam S)) (cos.f64 (*.f64 B lam)))) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (sin.f64 (-.f64 (*.f64 lam B) (*.f64 (PI.f64) #s(literal -1/2 binary64)))))))
(if.f64 (<=.f64 (*.f64 B lam) #s(literal -999999999999999945575230987042816 binary64)) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 (sin.f64 (-.f64 (PI.f64) (*.f64 lam B))) singam) (*.f64 cosgam S)) (cos.f64 (*.f64 B lam)))) (*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (sin.f64 (*.f64 B lam)) singam)) (cos.f64 (*.f64 B lam)))))
(*.f64 ArB (atan2.f64 (+.f64 (*.f64 S cosgam) (*.f64 (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 (-.f64 W #s(literal 1 binary64)) (/.f64 (-.f64 #s(literal 1 binary64) W) W)))
(*.f64 #s(literal 1/2 binary64) (-.f64 W (/.f64 #s(literal 1 binary64) W)))
(/.f64 (-.f64 (*.f64 S singam) (*.f64 (sin.f64 (*.f64 B lam)) cosgam)) T)
(asin.f64 (/.f64 (sin.f64 phi0) alpha))
(*.f64 (*.f64 n0 (/.f64 k0 (sin.f64 #s(literal 27401669256311/20000000000000 binary64)))) (sin.f64 (+.f64 #s(literal 27401669256311/20000000000000 binary64) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))
(*.f64 (*.f64 #s(literal -3 binary64) (*.f64 n0 k0)) (*.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) #s(literal -27401669256311/20000000000000 binary64))) #s(literal 1/3 binary64)))
(*.f64 (*.f64 #s(literal -3 binary64) n0) (*.f64 k0 (*.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) #s(literal -27401669256311/20000000000000 binary64))) #s(literal 1/3 binary64))))
(*.f64 (*.f64 (tan.f64 (-.f64 #s(literal -27401669256311/20000000000000 binary64) (*.f64 #s(literal 1/2 binary64) (PI.f64)))) k0) n0)
(*.f64 (*.f64 n0 k0) (tan.f64 (-.f64 #s(literal -27401669256311/20000000000000 binary64) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))
(/.f64 (*.f64 k0 n0) (tan.f64 #s(literal 27401669256311/20000000000000 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 phi0 #s(literal 8338720222318261/94758184344525691842589080106353915726128296943157752144717531617800961467674370503593652882607817257720198406807316479868870852301929589321550737002025216015896910157522577243058183937475491017166931103132108688408987234729984 binary64)) (*.f64 (*.f64 g (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64)))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (neg.f64 alpha))) (if.f64 (<=.f64 phi0 #s(literal 5000000000 binary64)) (/.f64 (*.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) g) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) alpha)) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (neg.f64 alpha))))))
(if.f64 (<=.f64 phi0 #s(literal 8338720222318261/94758184344525691842589080106353915726128296943157752144717531617800961467674370503593652882607817257720198406807316479868870852301929589321550737002025216015896910157522577243058183937475491017166931103132108688408987234729984 binary64)) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))))) (neg.f64 alpha)))) (if.f64 (<=.f64 phi0 #s(literal 5000000000 binary64)) (/.f64 (*.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) g) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) alpha)) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (sin.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0)) (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 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 9007199254740091/9007199254740992 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 (sin.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (/.f64 #s(literal 1 binary64) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) u0))) (*.f64 (sin.f64 (*.f64 u0 #s(literal 1/2 binary64))) (sin.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64))))))) g) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha)) (/.f64 (*.f64 (tan.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 4 binary64) (+.f64 (+.f64 u0 u0) (PI.f64))))) 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 phi0 #s(literal 8338720222318261/94758184344525691842589080106353915726128296943157752144717531617800961467674370503593652882607817257720198406807316479868870852301929589321550737002025216015896910157522577243058183937475491017166931103132108688408987234729984 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 phi0 #s(literal 5000000000 binary64)) (/.f64 (*.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) g) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) alpha)) (*.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))))))
(if.f64 (<=.f64 phi0 #s(literal 8338720222318261/94758184344525691842589080106353915726128296943157752144717531617800961467674370503593652882607817257720198406807316479868870852301929589321550737002025216015896910157522577243058183937475491017166931103132108688408987234729984 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)))) (if.f64 (<=.f64 phi0 #s(literal 5000000000 binary64)) (/.f64 (*.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) g) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) alpha)) (*.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))))))
(if.f64 (<=.f64 alpha #s(literal 2476979795053773/2251799813685248 binary64)) (/.f64 (*.f64 (tan.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 4 binary64) (+.f64 (+.f64 u0 u0) (PI.f64))))) 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)) (/.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 3602879701896397/4503599627370496 binary64)) (/.f64 (*.f64 (sin.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g) (*.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (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))))) (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 (tan.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 4 binary64) (+.f64 (+.f64 u0 u0) (PI.f64))))) 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 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) alpha) #s(literal 3602879701896397/4503599627370496 binary64)) (/.f64 (*.f64 (sin.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g) (*.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (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))))) (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 #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 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 4458563631096791/4503599627370496 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))) g) (/.f64 (*.f64 (/.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) g) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) alpha)))
(if.f64 (<=.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 4458563631096791/4503599627370496 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))) g) (/.f64 (*.f64 (tan.f64 (-.f64 (*.f64 u0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) g) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) alpha)))
(if.f64 (<=.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 4458563631096791/4503599627370496 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))) g) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (tan.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 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 4458563631096791/4503599627370496 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))) g) (/.f64 (*.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) g) (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha)))
(if.f64 (<=.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi0) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 4458563631096791/4503599627370496 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 alpha)) (tan.f64 (*.f64 (+.f64 (+.f64 u0 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 (+.f64 u0 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 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))) g)
(*.f64 (*.f64 g (pow.f64 (tan.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))))
(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 #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 (+.f64 Xr Yr) #s(literal 304747143/200000000000000000000 binary64)) (-.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 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr)) (*.f64 #s(literal 212097503/125000000000000000000000000 binary64) Yr))) (*.f64 (*.f64 (*.f64 #s(literal -3689471323/250000000000000000000000000000000 binary64) (*.f64 Xr Yr)) (+.f64 Xr Yr)) (-.f64 Xr Yr))) (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 4616990181/250000000000000000000 binary64) Xr) Yr) (-.f64 (*.f64 Xr #s(literal -4668270147/10000000000000000 binary64)) #s(literal -9827991/390625000 binary64))) (*.f64 #s(literal -149230739/1250000000000000 binary64) Yr))) (*.f64 #s(literal 4165541759/500000000000000000000000000000000 binary64) (+.f64 (*.f64 (-.f64 (*.f64 Yr Yr) (*.f64 #s(literal 6 binary64) (*.f64 Xr Xr))) (*.f64 Yr Yr)) (*.f64 (*.f64 (*.f64 Xr Xr) Xr) Xr))))
(-.f64 #s(literal 2946529277/100000000000 binary64) (-.f64 (-.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 (-.f64 (*.f64 (*.f64 (+.f64 Yr Xr) #s(literal 4616990181/500000000000000000000 binary64)) (-.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 #s(literal -881662847/200000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr)))) (*.f64 (*.f64 (*.f64 #s(literal 4165541759/125000000000000000000000000000000 binary64) (*.f64 Yr Xr)) (+.f64 Yr Xr)) (-.f64 Xr Yr)))) (*.f64 (+.f64 (*.f64 (-.f64 (*.f64 Yr Yr) (*.f64 #s(literal 6 binary64) (*.f64 Xr Xr))) (*.f64 Yr Yr)) (*.f64 (*.f64 (*.f64 Xr Xr) Xr) Xr)) #s(literal -3689471323/1000000000000000000000000000000000 binary64))))
(asin.f64 (-.f64 (*.f64 (sin.f64 s) (cos.f64 ad)) (*.f64 (*.f64 (cos.f64 s) (sin.f64 ad)) (cos.f64 d))))
(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 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (+.f64 (/.f64 (*.f64 (/.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 2 binary64)) (cos.f64 (*.f64 u #s(literal 1/2 binary64)))) (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 u u) (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 1/2 binary64)) (sin.f64 (*.f64 u #s(literal 1/2 binary64)))) (cos.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/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 1/4 binary64))))
(if.f64 (<=.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u) (/.f64 (PI.f64) #s(literal 4 binary64)))) #s(literal 7421932185906577/18014398509481984 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (/.f64 (cos.f64 (+.f64 (*.f64 (PI.f64) #s(literal 1/2 binary64)) (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64)))) (neg.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/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 (pow.f64 (PI.f64) #s(literal 2/3 binary64)) (*.f64 (cbrt.f64 (PI.f64)) #s(literal 1/4 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 (+.f64 u u) (PI.f64)) #s(literal 1/4 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 (+.f64 u u) (PI.f64)) #s(literal 1/4 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))
(if.f64 (<=.f64 k #s(literal 54000 binary64)) (*.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 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) (/.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 (*.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 (/.f64 (sin.f64 (+.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64))) (neg.f64 (cos.f64 (*.f64 (+.f64 (+.f64 u u) (PI.f64)) #s(literal 1/4 binary64))))) (/.f64 #s(literal 1 binary64) alpha)))) (*.f64 #s(literal 1/4 binary64) (PI.f64))) #s(literal 2 binary64)))
(if.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))) #s(literal 1 binary64)) (*.f64 (-.f64 (atan.f64 (*.f64 (pow.f64 (/.f64 #s(literal 1 binary64) (/.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))) (*.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 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (tan.f64 (+.f64 (+.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (-.f64 (*.f64 (PI.f64) #s(literal -1/4 binary64)) (PI.f64))) (PI.f64)) (PI.f64))) (/.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 #s(literal 1 binary64) (/.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))) (*.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 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 (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 -5043456793138493/2521728396569246669585858566409191283525103313309788586748690777871726193375821479130513040312634601011624191379636224 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 2043740476963553/10218702384817765435680628290748613458265350453429542612493041881278524886369096016818984783322294789577433327842265575649138882500575199542984559607218336872038429045509558663769793133795138494375185186532064890845853749530218856391110938974453986086436459043203870933208875495579361330830770176 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 -5043456793138493/2521728396569246669585858566409191283525103313309788586748690777871726193375821479130513040312634601011624191379636224 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 2043740476963553/10218702384817765435680628290748613458265350453429542612493041881278524886369096016818984783322294789577433327842265575649138882500575199542984559607218336872038429045509558663769793133795138494375185186532064890845853749530218856391110938974453986086436459043203870933208875495579361330830770176 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 -7872201966280717/15744403932561434696684473303452629045213679255131528440951130063136634306810047014785327192773139116009306758441243430342744218096217082060889571263281690386187633395165356521866664817226721079737670210248565328244806179188238434160900023542852296724603729870848 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))))
(asin.f64 (+.f64 (*.f64 (cos.f64 ad) (sin.f64 u)) (*.f64 (*.f64 (cos.f64 deltav) (cos.f64 u)) (sin.f64 ad))))
(/.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 rho0 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) #s(literal 27401669256311/20000000000000 binary64)) (/.f64 (PI.f64) #s(literal 4 binary64)))) n)) (pow.f64 (*.f64 (*.f64 (cos.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (+.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) s)) (sin.f64 (*.f64 s #s(literal 1/2 binary64)))))) (/.f64 #s(literal 1 binary64) (+.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (sin.f64 (+.f64 (*.f64 (fabs.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 s #s(literal 1/2 binary64)))) #s(literal 2 binary64)) (*.f64 (PI.f64) #s(literal 1/2 binary64)))))))) 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))))) (cos.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))))) (neg.f64 n))) rho0)
(if.f64 (<=.f64 (/.f64 (*.f64 rho0 (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) #s(literal 27401669256311/20000000000000 binary64)) (/.f64 (PI.f64) #s(literal 4 binary64)))) n)) (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) s) (/.f64 (PI.f64) #s(literal 4 binary64)))) n)) #s(literal 500000000000000004246810716844851488074434962299380307947499551351398452953088 binary64)) (*.f64 (*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 s s) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 n))) rho0) (*.f64 (pow.f64 (tan.f64 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 4 binary64) (+.f64 (+.f64 s s) (PI.f64))))) (neg.f64 n)) (*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) 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 (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 (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 (/.f64 #s(literal 1 binary64) (/.f64 #s(literal 2 binary64) (-.f64 #s(literal 1 binary64) (cos.f64 (+.f64 Cn Cn))))) (*.f64 (*.f64 (cos.f64 Cn) (cos.f64 Cn)) (*.f64 (cos.f64 lam) (cos.f64 lam))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (*.f64 (*.f64 (cos.f64 Cn) (cos.f64 Cn)) (*.f64 (cos.f64 lam) (cos.f64 lam))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (pow.f64 (*.f64 (cos.f64 lam) (cos.f64 Cn)) #s(literal 2 binary64)))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (+.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 lam lam)))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn))))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (-.f64 (+.f64 (pow.f64 (*.f64 (-.f64 (cos.f64 (+.f64 lam lam)) #s(literal -1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64)) #s(literal 1/2 binary64)) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 lam lam))))))
(/.f64 #s(literal 2 binary64) (sqrt.f64 (-.f64 #s(literal 1/2 binary64) (-.f64 (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 lam lam))) (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 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn)))) (pow.f64 (*.f64 (-.f64 (cos.f64 (+.f64 lam lam)) #s(literal -1 binary64)) #s(literal 1/2 binary64)) #s(literal 2 binary64)))))
(*.f64 (cos.f64 Cn) (/.f64 (sin.f64 lam) (sqrt.f64 (+.f64 (pow.f64 (*.f64 (cos.f64 lam) (cos.f64 Cn)) #s(literal 2 binary64)) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn))))))))
(*.f64 (cos.f64 Cn) (/.f64 (sin.f64 lam) (sqrt.f64 (+.f64 (*.f64 (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 lam lam)))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal -1/2 binary64) (cos.f64 (+.f64 Cn Cn))))) (-.f64 #s(literal 1/2 binary64) (*.f64 #s(literal 1/2 binary64) (cos.f64 (+.f64 Cn Cn))))))))
(*.f64 (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)))))))
(*.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)))))))
(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 (*.f64 cos_arg_r sinh_arg_i) #s(literal -8997827589086393/899782758908639276562107701537196371773939423660383326620151114522367516134297574063808137220778171429984182126704179024641653211098859937801703036173032771026944 binary64)) (-.f64 (*.f64 (*.f64 hr cosh_arg_i) sin_arg_r) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) (if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 5593535921577141/13983839803942852150595109342614667231797724205161451430391332862456574221363722734055314582922181823940116947860839707067576592166539734717014172871857704892828551908192974659309961553456532605588398441519406859119106715628289835749686907766833152 binary64)) (-.f64 (*.f64 (*.f64 hr sin_arg_r) cosh_arg_i) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)) (if.f64 (<=.f64 (*.f64 cos_arg_r sinh_arg_i) #s(literal 100000000000000000744898050207431989144199493858315387235964254131263985246781616026371987637390705840846560260278464628372543383280977318309056924111623883709653889736043921408 binary64)) (-.f64 (*.f64 (*.f64 sin_arg_r cosh_arg_i) hr) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hi)) (-.f64 (*.f64 (*.f64 hr sin_arg_r) cosh_arg_i) (*.f64 (*.f64 hi cos_arg_r) sinh_arg_i)))))
(if.f64 (<=.f64 cosh_arg_i #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 cosh_arg_i #s(literal -5902958103587057/295147905179352825856 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 cos_arg_r) sinh_arg_i)))
(if.f64 (<=.f64 cos_arg_r #s(literal -5415987671873539/19342813113834066795298816 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 hi cos_arg_r) sinh_arg_i)))
(if.f64 (<=.f64 sin_arg_r #s(literal -999999999999999966484112715463900049825186092620125502979674597309179755437379230686511104 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 2 binary64)) (+.f64 (*.f64 (*.f64 hi cosh_arg_i) sin_arg_r) (*.f64 (*.f64 hr sinh_arg_i) cos_arg_r)) (+.f64 (*.f64 (*.f64 hi sin_arg_r) cosh_arg_i) (*.f64 (*.f64 cos_arg_r sinh_arg_i) hr)))
(if.f64 (<=.f64 sinh_arg_i #s(literal 6156563468186637/615656346818663737691860001564743965704370926101022604186692084441339402679643915803347910232576806887603562348544 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 lam (/.f64 (cos.f64 phi) (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))))