(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 #s(literal 1/2 binary64) (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 #s(literal -1/3 binary64) (*.f64 #s(literal -3 binary64) (*.f64 (*.f64 k0 (tan.f64 (-.f64 #s(literal -27401669256311/20000000000000 binary64) (*.f64 #s(literal 1/2 binary64) (PI.f64))))) n0)))
(*.f64 (*.f64 (neg.f64 n0) (tan.f64 (-.f64 (*.f64 #s(literal 1/2 binary64) (PI.f64)) #s(literal -27401669256311/20000000000000 binary64)))) k0)
(*.f64 (*.f64 k0 (tan.f64 (-.f64 #s(literal -27401669256311/20000000000000 binary64) (*.f64 #s(literal 1/2 binary64) (PI.f64))))) n0)
(*.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 (*.f64 es es) es) #s(literal 1 binary64)) (-.f64 #s(literal -1 binary64) (*.f64 (-.f64 es #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 2558465639041233/492012622892544834271389790732894657662927476722860263779722814473019822087813969212041474774160597375290836201288857198210756815506783814402799101977552824568363543598917391308333275538335033741802194070267666507650193099632451067528125735714134272643866558464 binary64)) (*.f64 (*.f64 g (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (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 phi0 #s(literal 1/2 binary64))) (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 20499999999999999986131350259736876692371668992 binary64)) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (tan.f64 (+.f64 (+.f64 (PI.f64) (PI.f64)) (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)))) (neg.f64 alpha)))) (/.f64 (/.f64 (*.f64 (sin.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))) g) (cos.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) (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)))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) #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)))) (-.f64 (*.f64 (cos.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64))) (cos.f64 (*.f64 #s(literal -1/2 binary64) phi0))) (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) (sin.f64 (*.f64 phi0 #s(literal 1/2 binary64)))) #s(literal 1/2 binary64))))) alpha))))
(if.f64 (<=.f64 phi0 #s(literal 1196574698874669/31488807865122869393368946606905258090427358510263056881902260126273268613620094029570654385546278232018613516882486860685488436192434164121779142526563380772375266790330713043733329634453442159475340420497130656489612358376476868321800047085704593449207459741696 binary64)) (*.f64 (*.f64 g (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (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 phi0 #s(literal 1/2 binary64))) (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 2605242419472011/3064991081731777716716694054300618367237478244367204352 binary64)) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (tan.f64 (+.f64 (+.f64 (PI.f64) (PI.f64)) (+.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)))) (neg.f64 alpha)))) (*.f64 (*.f64 g (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0)))) (pow.f64 (/.f64 (*.f64 (-.f64 (cos.f64 (*.f64 phi0 #s(literal 1/2 binary64))) (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 phi0 #s(literal 1/2 binary64))) (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 alpha #s(literal 3872340515953109/4722366482869645213696 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 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))))) (if.f64 (<=.f64 alpha #s(literal 639999999999999961860656672342565313534930950454739086748785225212051344426359285643494416259324836439971033120964374897480159403881571079496941822191689423456921714688 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 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 binary64))) g)))
(if.f64 (<=.f64 alpha #s(literal 489999999999999985024337226896815131773586071661677202250276896638568910765107229839901806701940341473280 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 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))))) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 binary64))) g))
(if.f64 (<=.f64 alpha #s(literal 489999999999999985024337226896815131773586071661677202250276896638568910765107229839901806701940341473280 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)))) (neg.f64 alpha)) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))))) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 binary64))) g))
(if.f64 (<=.f64 phi0 #s(literal -4499999999999999969494269034627325894077687988748288 binary64)) (*.f64 (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))) (cos.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal -1/4 binary64)))) (neg.f64 alpha)) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))))) (if.f64 (<=.f64 phi0 #s(literal 1196574698874669/31488807865122869393368946606905258090427358510263056881902260126273268613620094029570654385546278232018613516882486860685488436192434164121779142526563380772375266790330713043733329634453442159475340420497130656489612358376476868321800047085704593449207459741696 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)))) (neg.f64 alpha)) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))))) (if.f64 (<=.f64 phi0 #s(literal 50000000000000002683081102196736 binary64)) (/.f64 (*.f64 (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64))) 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 (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))) (cos.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal -1/4 binary64)))) (neg.f64 alpha)) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (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 0 binary64)) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 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 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))))) (/.f64 (*.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) u0) (/.f64 (PI.f64) #s(literal 4 binary64)))) g) (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))))
(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 0 binary64)) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 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 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 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 0 binary64)) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 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 +inf.0 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))))) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (tan.f64 (+.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64)) (PI.f64)) (+.f64 phi0 phi0))) alpha) (*.f64 g (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) (*.f64 #s(literal -1/2 binary64) u0))))))))
(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 0 binary64)) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 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 +inf.0 binary64)) (/.f64 #s(literal 1 binary64) (/.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64))) alpha) (*.f64 g (tan.f64 (*.f64 (+.f64 (+.f64 u0 u0) (PI.f64)) #s(literal 1/4 binary64)))))) (*.f64 (*.f64 (pow.f64 (tan.f64 (-.f64 (+.f64 phi0 phi0) (*.f64 (*.f64 #s(literal 3 binary64) (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)))
(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 0 binary64)) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 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 +inf.0 binary64)) (*.f64 g (*.f64 (tan.f64 (*.f64 (+.f64 (*.f64 u0 #s(literal 2 binary64)) (PI.f64)) #s(literal 1/4 binary64))) (pow.f64 (tan.f64 (-.f64 (*.f64 phi0 #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (neg.f64 alpha)))) (*.f64 (*.f64 (pow.f64 (tan.f64 (-.f64 (+.f64 phi0 phi0) (*.f64 (*.f64 #s(literal 3 binary64) (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)))
(if.f64 (<=.f64 alpha #s(literal 15000000000000000000000 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 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 binary64))) g))
(if.f64 (<=.f64 alpha #s(literal 15000000000000000000000 binary64)) (*.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 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 binary64))) g))
(if.f64 (<=.f64 alpha #s(literal 15000000000000000000000 binary64)) (*.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 phi0 phi0) (PI.f64)) #s(literal 1/4 binary64)))) g) (*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 binary64))) g))
(*.f64 (pow.f64 (tan.f64 (*.f64 (*.f64 #s(literal 3 binary64) (PI.f64)) #s(literal 1/4 binary64))) (+.f64 (neg.f64 alpha) #s(literal 1 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 #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 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr))) (+.f64 (*.f64 (+.f64 (*.f64 Yr (*.f64 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 #s(literal 3689471323/250000000000000000000000000000000 binary64) Xr) Yr) (*.f64 (-.f64 Xr Yr) (+.f64 Xr Yr))))))
(-.f64 (-.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 (*.f64 #s(literal -149230739/1250000000000000 binary64) Yr) (+.f64 (*.f64 (*.f64 #s(literal -881662847/200000000000000000000000000 binary64) Xr) (-.f64 (*.f64 Xr Xr) (*.f64 #s(literal 3 binary64) (*.f64 Yr Yr)))) (*.f64 (*.f64 #s(literal -304747143/200000000000000000000 binary64) (+.f64 Xr Yr)) (-.f64 Xr Yr))))) (-.f64 (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr))) (+.f64 (*.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 #s(literal 3689471323/250000000000000000000000000000000 binary64) Xr) Yr) (*.f64 (-.f64 Xr Yr) (+.f64 Xr Yr))))))
(-.f64 (+.f64 (-.f64 (*.f64 Xr #s(literal -4668270147/10000000000000000 binary64)) #s(literal -9827991/390625000 binary64)) (*.f64 (*.f64 #s(literal 4616990181/250000000000000000000 binary64) Xr) Yr)) (+.f64 (+.f64 (+.f64 (*.f64 #s(literal -149230739/1250000000000000 binary64) 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 -304747143/200000000000000000000 binary64) (+.f64 Xr Yr)) (-.f64 Xr Yr))) (-.f64 (*.f64 (*.f64 #s(literal -212097503/125000000000000000000000000 binary64) Yr) (-.f64 (*.f64 #s(literal 3 binary64) (*.f64 Xr Xr)) (*.f64 Yr Yr))) (+.f64 (*.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 -4165541759/500000000000000000000000000000000 binary64)) (*.f64 (*.f64 (*.f64 #s(literal 3689471323/250000000000000000000000000000000 binary64) Xr) Yr) (*.f64 (+.f64 Xr Yr) (-.f64 Xr Yr)))))))
(-.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 (-.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 Yr (*.f64 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))))))
(-.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 Yr (*.f64 Yr (-.f64 (*.f64 Yr Yr) (*.f64 #s(literal 6 binary64) (*.f64 Xr Xr))))) (*.f64 (*.f64 Xr Xr) (*.f64 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 (*.f64 (pow.f64 (/.f64 (+.f64 (*.f64 (*.f64 (sqrt.f64 #s(literal 2 binary64)) #s(literal 1/2 binary64)) (sin.f64 (*.f64 u #s(literal 1/2 binary64)))) (*.f64 (cos.f64 (*.f64 #s(literal -1/2 binary64) u)) (*.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) u))) (*.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))) (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 #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 u #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (*.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 #s(literal -1/4 binary64) (PI.f64)) (*.f64 u #s(literal 1/2 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 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (/.f64 (sin.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (PI.f64)))) (neg.f64 (sin.f64 (+.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (*.f64 u #s(literal 1/2 binary64))) (*.f64 (PI.f64) #s(literal 1/2 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 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 k #s(literal 45000000000000000000 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 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (*.f64 (pow.f64 k (/.f64 #s(literal -1 binary64) alpha)) (pow.f64 (/.f64 (sin.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (PI.f64)))) (neg.f64 (cos.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (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 1/4 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 -4951760157141521/4951760157141521099596496896 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (*.f64 (/.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 u #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (/.f64 #s(literal 1 binary64) alpha)) (pow.f64 k (/.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)))) (*.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 (*.f64 (atan.f64 (*.f64 (*.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))) (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))))) #s(literal 4 binary64)) (PI.f64)) #s(literal 1/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 -1094449103198247/1886981212410770676120777290494134445458460610208220214188103150122812081196074426043063362588829383770734187515381922449885292314962396316280717125716348021824697663488 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 k (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) alpha)) gfi)) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) (if.f64 (<=.f64 phi #s(literal 4713081767787343/673297395398191808926846705008656469190443494761366254015779955592797715750026345230549316643531321757053820013590811867266487328023526894480360290643440748794268629904440908568532961009193307810725536222180016128 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 k (pow.f64 (tan.f64 (-.f64 (+.f64 (PI.f64) (*.f64 phi #s(literal 1/2 binary64))) (-.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (PI.f64)) (PI.f64)))) alpha)) gfi)) (*.f64 (PI.f64) #s(literal 1/4 binary64))))))
(if.f64 (<=.f64 phi #s(literal 3974446316289815/794889263257962974796277498092801308291525640763748664903194643469338087775424965801409745320266996710649718116931109481559848982586784968419475084821084743272680947722675151641735826243378403750534655587182832000457137589153821622272 binary64)) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 k (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) alpha)) gfi)) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) (*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 k (pow.f64 (tan.f64 (-.f64 (+.f64 (PI.f64) (*.f64 phi #s(literal 1/2 binary64))) (-.f64 (-.f64 (*.f64 #s(literal -1/4 binary64) (PI.f64)) (PI.f64)) (PI.f64)))) alpha)) gfi)) (*.f64 (PI.f64) #s(literal 1/4 binary64)))))
(*.f64 #s(literal 2 binary64) (-.f64 (atan.f64 (/.f64 (*.f64 k (pow.f64 (tan.f64 (+.f64 (*.f64 #s(literal 1/2 binary64) phi) (*.f64 (PI.f64) #s(literal 1/4 binary64)))) alpha)) gfi)) (*.f64 (PI.f64) #s(literal 1/4 binary64))))
(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 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) (*.f64 rho0 (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))))
(/.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 (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))))) n))
(*.f64 (pow.f64 (tan.f64 (-.f64 (*.f64 #s(literal 1/4 binary64) (PI.f64)) #s(literal -27401669256311/40000000000000 binary64))) n) (*.f64 rho0 (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))))
(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 +inf.0 binary64)) (*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 s s) (PI.f64)) #s(literal 1/4 binary64))) (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) rho0) (pow.f64 (/.f64 (sin.f64 (+.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64))) (PI.f64))) (neg.f64 (cos.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))))) 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 (sin.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 #s(literal -1/4 binary64) (PI.f64)))) (sin.f64 (-.f64 (*.f64 s #s(literal 1/2 binary64)) (*.f64 (PI.f64) #s(literal -3/4 binary64))))) n))
(*.f64 (pow.f64 (tan.f64 (*.f64 (+.f64 (+.f64 s s) (PI.f64)) #s(literal 1/4 binary64))) (neg.f64 n)) (*.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 #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)))