#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double M_PI;
double asin(double);
double atan2(double, double);
int /*<<< orphan*/ geocentric2geodetic_vec(double, double, double, double,
                                           double, double *, double *,
                                           double *);
int /*<<< orphan*/ geodetic2geocentric(double, double, double *, double *);
int /*<<< orphan*/ pomme_external_geocentric(int *, double, double, double,
                                             double, double, double, double,
                                             double, double, double *, double *,
                                             double *);
int /*<<< orphan*/ pomme_internal_geocentric(int, double, double, double,
                                             double, double *, double *,
                                             double *, double *, double *,
                                             double *);
double sqrt(double);
void pomme_geodetic(int control[17], double fday, double geod_alt,
                    double geod_lat, double lon, double est, double ist,
                    double imf_by, double f107, double Em, double *geod_x,
                    double *geod_y, double *geod_z, double *h, double *f,
                    double *incl, double *decl) {
  double theta, phi, lat_rad, delta, x, y, z, r, rrel;
  double mx, my, mz, cx, cy, cz, ex, ey, ez, fyear;
  geodetic2geocentric(geod_lat * M_PI / 180.0, geod_alt, &lat_rad, &r);
  theta = M_PI / 2.0 - lat_rad;
  phi = lon * M_PI / 180.0;
  fyear = 2000.0 + fday / 365.25;
  pomme_internal_geocentric(control[0], fyear, r, theta, phi, &mx, &my, &mz,
                            &cx, &cy, &cz);
  rrel = r / 6371.2;
  pomme_external_geocentric(control, fday, rrel, theta, phi, est, ist, imf_by,
                            f107, Em, &ex, &ey, &ez);
  x = mx + cx + ex;
  y = my + cy + ey;
  z = mz + cz + ez;
  *f = sqrt(x * x + y * y + z * z);
  delta = (90.0 - geod_lat) * M_PI / 180.0;
  geocentric2geodetic_vec(theta, delta, x, y, z, geod_x, geod_y, geod_z);
  *h = sqrt(*geod_x * *geod_x + *geod_y * *geod_y);
  *incl = 0.0;
  if (*f > 0.0)
    *incl = asin(*geod_z / *f);
  *incl *= 180.0 / M_PI;
  *decl = atan2(*geod_y, *geod_x);
  *decl *= 180.0 / M_PI;
  return;
}


x86 assembly (clang, -O0):
<pomme_geodetic>:
push   %rbp
mov    %rsp,%rbp
sub    $0x110,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
movsd  0x10(%rbp),%xmm8
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
movsd  %xmm4,-0x30(%rbp)
movsd  %xmm5,-0x38(%rbp)
movsd  %xmm6,-0x40(%rbp)
movsd  %xmm7,-0x48(%rbp)
mov    %rsi,-0x50(%rbp)
mov    %rdx,-0x58(%rbp)
mov    %rcx,-0x60(%rbp)
mov    %r8,-0x68(%rbp)
mov    %r9,-0x70(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    0x2e13(%rip),%rax
mulsd  (%rax),%xmm0
movsd  0xe1f(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x18(%rbp),%xmm1
lea    -0x88(%rbp),%rdi
lea    -0xb0(%rbp),%rsi
call   1090 <geodetic2geocentric@plt>
mov    0x2de4(%rip),%rax
movsd  (%rax),%xmm0
movsd  0xe18(%rip),%xmm1
divsd  %xmm1,%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x78(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    0x2dbb(%rip),%rax
mulsd  (%rax),%xmm0
movsd  0xdc7(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x10(%rbp),%xmm1
movsd  0xdd1(%rip),%xmm0
divsd  %xmm0,%xmm1
movsd  0xdbd(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    (%rax),%edi
movsd  -0x108(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm1
ff
movsd  -0x78(%rbp),%xmm2
movsd  -0x80(%rbp),%xmm3
lea    -0xc0(%rbp),%rsi
lea    -0xc8(%rbp),%rdx
lea    -0xd0(%rbp),%rcx
lea    -0xd8(%rbp),%r8
lea    -0xe0(%rbp),%r9
lea    -0xe8(%rbp),%rax
mov    %rax,(%rsp)
call   1050 <pomme_internal_geocentric@plt>
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0xd46(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x8(%rbp),%rdi
movsd  -0x10(%rbp),%xmm0
movsd  -0xb8(%rbp),%xmm1
ff
movsd  -0x78(%rbp),%xmm2
movsd  -0x80(%rbp),%xmm3
movsd  -0x30(%rbp),%xmm4
movsd  -0x38(%rbp),%xmm5
movsd  -0x40(%rbp),%xmm6
movsd  -0x48(%rbp),%xmm7
movsd  0x10(%rbp),%xmm8
lea    -0xf0(%rbp),%rsi
lea    -0xf8(%rbp),%rdx
lea    -0x100(%rbp),%rcx
movsd  %xmm8,(%rsp)
call   1040 <pomme_external_geocentric@plt>
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0xd8(%rbp),%xmm0
ff
addsd  -0xf0(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
addsd  -0xe0(%rbp),%xmm0
ff
addsd  -0xf8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
addsd  -0xe8(%rbp),%xmm0
ff
addsd  -0x100(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x98(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm0
ff
mulsd  -0xa0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1080 <sqrt@plt>
mov    -0x70(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  0xc28(%rip),%xmm0
subsd  -0x20(%rbp),%xmm0
mov    0x2bfc(%rip),%rax
mulsd  (%rax),%xmm0
movsd  0xc08(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm3
ff
movsd  -0xa8(%rbp),%xmm4
ff
mov    -0x50(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x60(%rbp),%rdx
call   1070 <geocentric2geodetic_vec@plt>
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x58(%rbp),%rax
mulsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1080 <sqrt@plt>
mov    -0x68(%rbp),%rax
movsd  %xmm0,(%rax)
mov    0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x70(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    14ac <pomme_geodetic+0x33c>
mov    -0x60(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x70(%rbp),%rax
divsd  (%rax),%xmm0
call   1060 <asin@plt>
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  0xb4c(%rip),%xmm0
mov    0x2b2d(%rip),%rax
divsd  (%rax),%xmm0
mov    0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm1
call   1030 <atan2@plt>
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  0xb10(%rip),%xmm0
mov    0x2af1(%rip),%rax
divsd  (%rax),%xmm0
mov    0x20(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
add    $0x110,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pomme_geodetic>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r9,%r12
mov    %r8,0x50(%rsp)
mov    %rcx,%r15
mov    %rdx,%r14
mov    %rsi,%rbx
movsd  %xmm7,0x18(%rsp)
movsd  %xmm6,0x40(%rsp)
movsd  %xmm5,0x38(%rsp)
movsd  %xmm4,0x30(%rsp)
movsd  %xmm3,0x8(%rsp)
movsd  %xmm0,0x10(%rsp)
mov    %rdi,%rbp
mov    0x2e28(%rip),%r13
movsd  0x0(%r13),%xmm0
movsd  %xmm2,0x48(%rsp)
mulsd  %xmm2,%xmm0
divsd  0xe28(%rip),%xmm0
lea    0xa0(%rsp),%rdi
lea    0x28(%rsp),%rsi
call   1090 <geodetic2geocentric@plt>
movsd  0x0(%r13),%xmm0
movsd  0xe10(%rip),%xmm1
mulsd  %xmm0,%xmm1
subsd  0xa0(%rsp),%xmm1
movapd %xmm1,%xmm2
movsd  %xmm1,0x20(%rsp)
movsd  0x8(%rsp),%xmm1
mulsd  %xmm0,%xmm1
divsd  0xddf(%rip),%xmm1
movapd %xmm1,%xmm3
movsd  %xmm1,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
divsd  0xdd7(%rip),%xmm0
addsd  0xdd7(%rip),%xmm0
mov    0x0(%rbp),%edi
movsd  0x28(%rsp),%xmm1
lea    0x70(%rsp),%rax
mov    %rax,(%rsp)
lea    0x98(%rsp),%rsi
lea    0x90(%rsp),%rdx
lea    0x88(%rsp),%rcx
lea    0x80(%rsp),%r8
lea    0x78(%rsp),%r9
call   1050 <pomme_internal_geocentric@plt>
movsd  0x28(%rsp),%xmm1
divsd  0xd95(%rip),%xmm1
movsd  0xe0(%rsp),%xmm0
movsd  %xmm0,(%rsp)
lea    0x68(%rsp),%rsi
lea    0x60(%rsp),%rdx
lea    0x58(%rsp),%rcx
mov    %rbp,%rdi
movsd  0x10(%rsp),%xmm0
movsd  0x20(%rsp),%xmm2
movsd  0x8(%rsp),%xmm3
movsd  0x30(%rsp),%xmm4
movsd  0x38(%rsp),%xmm5
movsd  0x40(%rsp),%xmm6
movsd  0x18(%rsp),%xmm7
call   1040 <pomme_external_geocentric@plt>
movsd  0x98(%rsp),%xmm2
addsd  0x80(%rsp),%xmm2
addsd  0x68(%rsp),%xmm2
movsd  0x90(%rsp),%xmm3
addsd  0x78(%rsp),%xmm3
addsd  0x60(%rsp),%xmm3
movsd  0x88(%rsp),%xmm4
addsd  0x70(%rsp),%xmm4
addsd  0x58(%rsp),%xmm4
movapd %xmm3,%xmm0
mulsd  %xmm3,%xmm0
movapd %xmm2,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movapd %xmm4,%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     134c <pomme_geodetic+0x1dc>
sqrtsd %xmm0,%xmm0
jmp    1375 <pomme_geodetic+0x205>
movsd  %xmm2,0x8(%rsp)
movsd  %xmm3,0x10(%rsp)
movsd  %xmm4,0x18(%rsp)
call   1080 <sqrt@plt>
movsd  0x18(%rsp),%xmm4
movsd  0x10(%rsp),%xmm3
movsd  0x8(%rsp),%xmm2
movsd  %xmm0,(%r12)
movsd  0xca5(%rip),%xmm1
subsd  0x48(%rsp),%xmm1
mulsd  0x0(%r13),%xmm1
divsd  0xc69(%rip),%xmm1
movsd  0x20(%rsp),%xmm0
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15,%rdx
call   1070 <geocentric2geodetic_vec@plt>
movsd  (%rbx),%xmm0
movsd  (%r14),%xmm1
mulsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xc68(%rip),%xmm0
jb     13d0 <pomme_geodetic+0x260>
sqrtsd %xmm0,%xmm0
jmp    13d5 <pomme_geodetic+0x265>
call   1080 <sqrt@plt>
mov    0xf0(%rsp),%rbp
mov    0x50(%rsp),%rax
movsd  %xmm0,(%rax)
mov    0xe8(%rsp),%rax
movq   $0x0,(%rax)
movsd  (%r12),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jbe    141f <pomme_geodetic+0x2af>
movsd  (%r15),%xmm0
divsd  %xmm1,%xmm0
call   1060 <asin@plt>
mov    0xe8(%rsp),%rax
movsd  %xmm0,(%rax)
movsd  0xbd9(%rip),%xmm1
divsd  0x0(%r13),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax)
movsd  (%r14),%xmm0
movsd  (%rbx),%xmm1
call   1030 <atan2@plt>
movsd  %xmm0,0x0(%rbp)
movsd  0xbb0(%rip),%xmm1
divsd  0x0(%r13),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp)
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define void @pomme_geodetic(i32* %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8, double %9, double* %10, double* %11, double* %12, double* %13, double* %14, double* %15, double* %16) {
%18 = alloca i32*
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double*
%29 = alloca double*
%30 = alloca double*
%31 = alloca double*
%32 = alloca double*
%33 = alloca double*
%34 = alloca double*
%35 = alloca double
%36 = alloca double
%37 = alloca double
%38 = alloca double
%39 = alloca double
%40 = alloca double
%41 = alloca double
%42 = alloca double
%43 = alloca double
%44 = alloca double
%45 = alloca double
%46 = alloca double
%47 = alloca double
%48 = alloca double
%49 = alloca double
%50 = alloca double
%51 = alloca double
%52 = alloca double
%53 = alloca double
store i32* %0, i32** %18
store double %1, double* %19
store double %2, double* %20
store double %3, double* %21
store double %4, double* %22
store double %5, double* %23
store double %6, double* %24
store double %7, double* %25
store double %8, double* %26
store double %9, double* %27
store double* %10, double** %28
store double* %11, double** %29
store double* %12, double** %30
store double* %13, double** %31
store double* %14, double** %32
store double* %15, double** %33
store double* %16, double** %34
%54 = load double, double* %21
%55 = load double, double* @M_PI
%56 = fmul double %54, %55
%57 = fdiv double %56, 1.800000e+02
%58 = load double, double* %20
%59 = call i32 @geodetic2geocentric(double %57, double %58, double* %37, double* %42)
%60 = load double, double* @M_PI
%61 = fdiv double %60, 2.000000e+00
%62 = load double, double* %37
%63 = fsub double %61, %62
store double %63, double* %35
%64 = load double, double* %22
%65 = load double, double* @M_PI
%66 = fmul double %64, %65
%67 = fdiv double %66, 1.800000e+02
store double %67, double* %36
%68 = load double, double* %19
%69 = fdiv double %68, 3.652500e+02
%70 = fadd double 2.000000e+03, %69
store double %70, double* %53
%71 = load i32*, i32** %18
%72 = getelementptr inbounds i32, i32* %71, i64 0
%73 = load i32, i32* %72
%74 = load double, double* %53
%75 = load double, double* %42
%76 = load double, double* %35
%77 = load double, double* %36
%78 = call i32 @pomme_internal_geocentric(i32 %73, double %74, double %75, double %76, double %77, double* %44, double* %45, double* %46, double* %47, double* %48, double* %49)
%79 = load double, double* %42
%80 = fdiv double %79, 6.371200e+03
store double %80, double* %43
%81 = load i32*, i32** %18
%82 = load double, double* %19
%83 = load double, double* %43
%84 = load double, double* %35
%85 = load double, double* %36
%86 = load double, double* %23
%87 = load double, double* %24
%88 = load double, double* %25
%89 = load double, double* %26
%90 = load double, double* %27
%91 = call i32 @pomme_external_geocentric(i32* %81, double %82, double %83, double %84, double %85, double %86, double %87, double %88, double %89, double %90, double* %50, double* %51, double* %52)
%92 = load double, double* %44
%93 = load double, double* %47
%94 = fadd double %92, %93
%95 = load double, double* %50
%96 = fadd double %94, %95
store double %96, double* %39
%97 = load double, double* %45
%98 = load double, double* %48
%99 = fadd double %97, %98
%100 = load double, double* %51
%101 = fadd double %99, %100
store double %101, double* %40
%102 = load double, double* %46
%103 = load double, double* %49
%104 = fadd double %102, %103
%105 = load double, double* %52
%106 = fadd double %104, %105
store double %106, double* %41
%107 = load double, double* %39
%108 = load double, double* %39
%109 = load double, double* %40
%110 = load double, double* %40
%111 = fmul double %109, %110
%112 = call double @llvm.fmuladd.f64(double %107, double %108, double %111)
%113 = load double, double* %41
%114 = load double, double* %41
%115 = call double @llvm.fmuladd.f64(double %113, double %114, double %112)
%116 = call double @sqrt(double %115)
%117 = load double*, double** %32
store double %116, double* %117
%118 = load double, double* %21
%119 = fsub double 9.000000e+01, %118
%120 = load double, double* @M_PI
%121 = fmul double %119, %120
%122 = fdiv double %121, 1.800000e+02
store double %122, double* %38
%123 = load double, double* %35
%124 = load double, double* %38
%125 = load double, double* %39
%126 = load double, double* %40
%127 = load double, double* %41
%128 = load double*, double** %28
%129 = load double*, double** %29
%130 = load double*, double** %30
%131 = call i32 @geocentric2geodetic_vec(double %123, double %124, double %125, double %126, double %127, double* %128, double* %129, double* %130)
%132 = load double*, double** %28
%133 = load double, double* %132
%134 = load double*, double** %28
%135 = load double, double* %134
%136 = load double*, double** %29
%137 = load double, double* %136
%138 = load double*, double** %29
%139 = load double, double* %138
%140 = fmul double %137, %139
%141 = call double @llvm.fmuladd.f64(double %133, double %135, double %140)
%142 = call double @sqrt(double %141)
%143 = load double*, double** %31
store double %142, double* %143
%144 = load double*, double** %33
store double 0.000000e+00, double* %144
%145 = load double*, double** %32
%146 = load double, double* %145
%147 = fcmp ogt double %146, 0.000000e+00
br i1 %147, label %148, label %156
148: ; preds = %17
%149 = load double*, double** %30
%150 = load double, double* %149
%151 = load double*, double** %32
%152 = load double, double* %151
%153 = fdiv double %150, %152
%154 = call double @asin(double %153)
%155 = load double*, double** %33
store double %154, double* %155
br label %156
156: ; preds = %148, %17
%157 = load double, double* @M_PI
%158 = fdiv double 1.800000e+02, %157
%159 = load double*, double** %33
%160 = load double, double* %159
%161 = fmul double %160, %158
store double %161, double* %159
%162 = load double*, double** %29
%163 = load double, double* %162
%164 = load double*, double** %28
%165 = load double, double* %164
%166 = call double @atan2(double %163, double %165)
%167 = load double*, double** %34
store double %166, double* %167
%168 = load double, double* @M_PI
%169 = fdiv double 1.800000e+02, %168
%170 = load double*, double** %34
%171 = load double, double* %170
%172 = fmul double %171, %169
store double %172, double* %170
ret void
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define void @pomme_geodetic(i32* %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8, double %9, double* %10, double* %11, double* %12, double* nocapture writeonly %13, double* nocapture %14, double* nocapture %15, double* nocapture writeonly %16) {
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = bitcast double* %18 to i8*
%30 = bitcast double* %19 to i8*
%31 = bitcast double* %20 to i8*
%32 = bitcast double* %21 to i8*
%33 = bitcast double* %22 to i8*
%34 = bitcast double* %23 to i8*
%35 = bitcast double* %24 to i8*
%36 = bitcast double* %25 to i8*
%37 = bitcast double* %26 to i8*
%38 = bitcast double* %27 to i8*
%39 = bitcast double* %28 to i8*
%40 = load double, double* @M_PI
%41 = fmul double %40, %3
%42 = fdiv double %41, 1.800000e+02
%43 = call i32 @geodetic2geocentric(double %42, double %2, double* %18, double* %19)
%44 = load double, double* @M_PI
%45 = fmul double %44, 5.000000e-01
%46 = load double, double* %18
%47 = fsub double %45, %46
%48 = fmul double %44, %4
%49 = fdiv double %48, 1.800000e+02
%50 = fdiv double %1, 3.652500e+02
%51 = fadd double %50, 2.000000e+03
%52 = load i32, i32* %0
%53 = load double, double* %19
%54 = call i32 @pomme_internal_geocentric(i32 %52, double %51, double %53, double %47, double %49, double* %20, double* %21, double* %22, double* %23, double* %24, double* %25)
%55 = load double, double* %19
%56 = fdiv double %55, 6.371200e+03
%57 = call i32 @pomme_external_geocentric(i32* %0, double %1, double %56, double %47, double %49, double %5, double %6, double %7, double %8, double %9, double* %26, double* %27, double* %28)
%58 = load double, double* %20
%59 = load double, double* %23
%60 = fadd double %58, %59
%61 = load double, double* %26
%62 = fadd double %60, %61
%63 = load double, double* %21
%64 = load double, double* %24
%65 = fadd double %63, %64
%66 = load double, double* %27
%67 = fadd double %65, %66
%68 = load double, double* %22
%69 = load double, double* %25
%70 = fadd double %68, %69
%71 = load double, double* %28
%72 = fadd double %70, %71
%73 = fmul double %67, %67
%74 = call double @llvm.fmuladd.f64(double %62, double %62, double %73)
%75 = call double @llvm.fmuladd.f64(double %72, double %72, double %74)
%76 = call double @sqrt(double %75)
store double %76, double* %14
%77 = fsub double 9.000000e+01, %3
%78 = load double, double* @M_PI
%79 = fmul double %77, %78
%80 = fdiv double %79, 1.800000e+02
%81 = call i32 @geocentric2geodetic_vec(double %47, double %80, double %62, double %67, double %72, double* %10, double* %11, double* %12)
%82 = load double, double* %10
%83 = load double, double* %11
%84 = fmul double %83, %83
%85 = call double @llvm.fmuladd.f64(double %82, double %82, double %84)
%86 = call double @sqrt(double %85)
store double %86, double* %13
store double 0.000000e+00, double* %15
%87 = load double, double* %14
%88 = fcmp ogt double %87, 0.000000e+00
br i1 %88, label %89, label %93
89: ; preds = %17
%90 = load double, double* %12
%91 = fdiv double %90, %87
%92 = call double @asin(double %91)
store double %92, double* %15
br label %93
93: ; preds = %89, %17
%94 = phi double [ %92, %89 ], [ 0.000000e+00, %17 ]
%95 = load double, double* @M_PI
%96 = fdiv double 1.800000e+02, %95
%97 = fmul double %94, %96
store double %97, double* %15
%98 = load double, double* %11
%99 = load double, double* %10
%100 = call double @atan2(double %98, double %99)
store double %100, double* %16
%101 = load double, double* @M_PI
%102 = fdiv double 1.800000e+02, %101
%103 = fmul double %100, %102
store double %103, double* %16
ret void
}