#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 atan2(double, double);
int /*<<< orphan*/ matCalcCross(double *, double *, double *);
double matCalcDot(double *, double *);
double matCalcLen(double *);
double matCalcTorsion(double *v1, double *v2, double *v3, double *v4) {
  double d1[3], d2[3], d3[3], d4[3], c1[3], c2[3], c3[3], l1, l2, l3, r;
  d1[0] = v1[0] - v2[0];
  d1[1] = v1[1] - v2[1];
  d1[2] = v1[2] - v2[2];
  d2[0] = v3[0] - v2[0];
  d2[1] = v3[1] - v2[1];
  d2[2] = v3[2] - v2[2];
  d3[0] = v2[0] - v3[0];
  d3[1] = v2[1] - v3[1];
  d3[2] = v2[2] - v3[2];
  d4[0] = v4[0] - v3[0];
  d4[1] = v4[1] - v3[1];
  d4[2] = v4[2] - v3[2];
  matCalcCross(d1, d2, c1);
  matCalcCross(d3, d4, c2);
  matCalcCross(c1, d3, c3);
  l1 = matCalcLen(c1);
  l2 = matCalcLen(c2);
  l3 = matCalcLen(c3);
  if (l1 == 0.0 || l2 == 0.0 || l3 == 0.0) {
    return 0.0;
  }
  r = atan2(matCalcDot(c3, c2) / (l3 * l2), matCalcDot(c1, c2) / (l1 * l2));
  return r * 180.0 / M_PI;
}


x86 assembly (clang, -O0):
<matCalcTorsion>:
push   %rbp
mov    %rsp,%rbp
sub    $0x130,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x28(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x28(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
lea    -0x40(%rbp),%rdi
lea    -0x60(%rbp),%rsi
lea    -0xc0(%rbp),%rdx
call   1030 <matCalcCross@plt>
lea    -0x80(%rbp),%rdi
lea    -0xa0(%rbp),%rsi
lea    -0xe0(%rbp),%rdx
call   1030 <matCalcCross@plt>
lea    -0xc0(%rbp),%rdi
lea    -0x80(%rbp),%rsi
lea    -0x100(%rbp),%rdx
call   1030 <matCalcCross@plt>
lea    -0xc0(%rbp),%rdi
call   1050 <matCalcLen@plt>
movsd  %xmm0,-0x108(%rbp)
ff
lea    -0xe0(%rbp),%rdi
call   1050 <matCalcLen@plt>
movsd  %xmm0,-0x110(%rbp)
ff
lea    -0x100(%rbp),%rdi
call   1050 <matCalcLen@plt>
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x108(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    130e <matCalcTorsion+0x1ce>
jp     130e <matCalcTorsion+0x1ce>
jmp    1349 <matCalcTorsion+0x209>
movsd  -0x110(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    132e <matCalcTorsion+0x1ee>
jp     132e <matCalcTorsion+0x1ee>
jmp    1349 <matCalcTorsion+0x209>
movsd  -0x118(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1356 <matCalcTorsion+0x216>
jp     1356 <matCalcTorsion+0x216>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    13e4 <matCalcTorsion+0x2a4>
lea    -0x100(%rbp),%rdi
lea    -0xe0(%rbp),%rsi
call   1060 <matCalcDot@plt>
movsd  -0x118(%rbp),%xmm1
ff
mulsd  -0x110(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
lea    -0xc0(%rbp),%rdi
lea    -0xe0(%rbp),%rsi
call   1060 <matCalcDot@plt>
movaps %xmm0,%xmm1
movsd  -0x128(%rbp),%xmm0
ff
movsd  -0x108(%rbp),%xmm2
ff
mulsd  -0x110(%rbp),%xmm2
ff
divsd  %xmm2,%xmm1
call   1040 <atan2@plt>
movsd  %xmm0,-0x120(%rbp)
ff
movsd  0xc34(%rip),%xmm0
mulsd  -0x120(%rbp),%xmm0
ff
mov    0x2c0d(%rip),%rax
divsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x130,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matCalcTorsion>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0xf8,%rsp
movupd (%rdi),%xmm0
movupd (%rsi),%xmm1
subpd  %xmm1,%xmm0
movapd %xmm0,0x80(%rsp)
movsd  0x10(%rdi),%xmm0
movsd  0x10(%rsi),%xmm2
subsd  %xmm2,%xmm0
movsd  %xmm0,0x90(%rsp)
movupd (%rdx),%xmm0
movapd %xmm0,%xmm3
subpd  %xmm1,%xmm3
movapd %xmm3,0x60(%rsp)
movsd  0x10(%rdx),%xmm3
movapd %xmm3,%xmm4
subsd  %xmm2,%xmm4
movsd  %xmm4,0x70(%rsp)
subpd  %xmm0,%xmm1
movapd %xmm1,0x40(%rsp)
subsd  %xmm3,%xmm2
movsd  %xmm2,0x50(%rsp)
movupd (%rcx),%xmm1
subpd  %xmm0,%xmm1
movapd %xmm1,0x20(%rsp)
movsd  0x10(%rcx),%xmm0
subsd  %xmm3,%xmm0
movsd  %xmm0,0x30(%rsp)
lea    0x80(%rsp),%rdi
lea    0x60(%rsp),%rsi
lea    0xe0(%rsp),%rbx
mov    %rbx,%rdx
call   1030 <matCalcCross@plt>
lea    0x40(%rsp),%r14
lea    0x20(%rsp),%rsi
lea    0xc0(%rsp),%r15
mov    %r14,%rdi
mov    %r15,%rdx
call   1030 <matCalcCross@plt>
lea    0xa0(%rsp),%r12
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r12,%rdx
call   1030 <matCalcCross@plt>
mov    %rbx,%rdi
call   1050 <matCalcLen@plt>
movsd  %xmm0,0x18(%rsp)
mov    %r15,%rdi
call   1050 <matCalcLen@plt>
movsd  %xmm0,0x10(%rsp)
mov    %r12,%rdi
call   1050 <matCalcLen@plt>
movapd %xmm0,%xmm3
xorpd  %xmm0,%xmm0
movsd  0x18(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jne    125e <matCalcTorsion+0x11e>
jnp    12ee <matCalcTorsion+0x1ae>
xorpd  %xmm1,%xmm1
movsd  0x10(%rsp),%xmm2
ucomisd %xmm1,%xmm2
jne    1270 <matCalcTorsion+0x130>
jnp    12ee <matCalcTorsion+0x1ae>
ucomisd %xmm1,%xmm3
jne    1278 <matCalcTorsion+0x138>
jnp    12ee <matCalcTorsion+0x1ae>
lea    0xa0(%rsp),%rdi
lea    0xc0(%rsp),%rbx
mov    %rbx,%rsi
movsd  %xmm3,0x8(%rsp)
call   1060 <matCalcDot@plt>
movsd  0x8(%rsp),%xmm1
mulsd  0x10(%rsp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0xe0(%rsp),%rdi
mov    %rbx,%rsi
call   1060 <matCalcDot@plt>
movapd %xmm0,%xmm1
movsd  0x18(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
divsd  %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
call   1040 <atan2@plt>
mulsd  0xd1d(%rip),%xmm0
mov    0x2cfe(%rip),%rax
divsd  (%rax),%xmm0
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define double @matCalcTorsion(double* nocapture %0, double* nocapture %1, double* nocapture %2, double* nocapture %3) {
%5 = alloca [3 x double]
%6 = alloca [3 x double]
%7 = alloca [3 x double]
%8 = alloca [3 x double]
%9 = alloca [3 x double]
%10 = alloca [3 x double]
%11 = alloca [3 x double]
%12 = bitcast [3 x double]* %5 to i8*
%13 = bitcast [3 x double]* %6 to i8*
%14 = bitcast [3 x double]* %7 to i8*
%15 = bitcast [3 x double]* %8 to i8*
%16 = bitcast [3 x double]* %9 to i8*
%17 = bitcast [3 x double]* %10 to i8*
%18 = bitcast [3 x double]* %11 to i8*
%19 = getelementptr inbounds [3 x double], [3 x double]* %5, i64 0, i64 0
%20 = bitcast double* %0 to <2 x double>*
%21 = load <2 x double>, <2 x double>* %20
%22 = bitcast double* %1 to <2 x double>*
%23 = load <2 x double>, <2 x double>* %22
%24 = fsub <2 x double> %21, %23
%25 = bitcast [3 x double]* %5 to <2 x double>*
store <2 x double> %24, <2 x double>* %25
%26 = getelementptr inbounds double, double* %0, i64 2
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %1, i64 2
%29 = load double, double* %28
%30 = fsub double %27, %29
%31 = getelementptr inbounds [3 x double], [3 x double]* %5, i64 0, i64 2
store double %30, double* %31
%32 = getelementptr inbounds [3 x double], [3 x double]* %6, i64 0, i64 0
%33 = bitcast double* %2 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = fsub <2 x double> %34, %23
%36 = bitcast [3 x double]* %6 to <2 x double>*
store <2 x double> %35, <2 x double>* %36
%37 = getelementptr inbounds double, double* %2, i64 2
%38 = load double, double* %37
%39 = fsub double %38, %29
%40 = getelementptr inbounds [3 x double], [3 x double]* %6, i64 0, i64 2
store double %39, double* %40
%41 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 0
%42 = fsub <2 x double> %23, %34
%43 = bitcast [3 x double]* %7 to <2 x double>*
store <2 x double> %42, <2 x double>* %43
%44 = fsub double %29, %38
%45 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 2
store double %44, double* %45
%46 = getelementptr inbounds [3 x double], [3 x double]* %8, i64 0, i64 0
%47 = bitcast double* %3 to <2 x double>*
%48 = load <2 x double>, <2 x double>* %47
%49 = fsub <2 x double> %48, %34
%50 = bitcast [3 x double]* %8 to <2 x double>*
store <2 x double> %49, <2 x double>* %50
%51 = getelementptr inbounds double, double* %3, i64 2
%52 = load double, double* %51
%53 = fsub double %52, %38
%54 = getelementptr inbounds [3 x double], [3 x double]* %8, i64 0, i64 2
store double %53, double* %54
%55 = getelementptr inbounds [3 x double], [3 x double]* %9, i64 0, i64 0
%56 = call i32 @matCalcCross(double* %19, double* %32, double* %55)
%57 = getelementptr inbounds [3 x double], [3 x double]* %10, i64 0, i64 0
%58 = call i32 @matCalcCross(double* %41, double* %46, double* %57)
%59 = getelementptr inbounds [3 x double], [3 x double]* %11, i64 0, i64 0
%60 = call i32 @matCalcCross(double* %55, double* %41, double* %59)
%61 = call double @matCalcLen(double* %55)
%62 = call double @matCalcLen(double* %57)
%63 = call double @matCalcLen(double* %59)
%64 = fcmp oeq double %61, 0.000000e+00
%65 = fcmp oeq double %62, 0.000000e+00
%66 = select i1 %64, i1 true, i1 %65
%67 = fcmp oeq double %63, 0.000000e+00
%68 = select i1 %66, i1 true, i1 %67
br i1 %68, label %80, label %69
69: ; preds = %4
%70 = call double @matCalcDot(double* %59, double* %57)
%71 = fmul double %62, %63
%72 = fdiv double %70, %71
%73 = call double @matCalcDot(double* %55, double* %57)
%74 = fmul double %61, %62
%75 = fdiv double %73, %74
%76 = call double @atan2(double %72, double %75)
%77 = fmul double %76, 1.800000e+02
%78 = load double, double* @M_PI
%79 = fdiv double %77, %78
br label %80
80: ; preds = %4, %69
%81 = phi double [ %79, %69 ], [ 0.000000e+00, %4 ]
ret double %81
}