#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 exp(int);
void reaction_rate(double r[], double PR[], double TR[], double Clock) {
  double yr, alpha;
  yr = Clock / 60 / 24 / 365;
  if (TR[0] < 180)
    alpha = exp(-yr / 4.481420117724551);
  else
    alpha = exp(-(TR[0] - 180) / 50) * exp(-yr / 0.621);
  TR[0] = TR[0] + 273.15;
  r[0] = 0.1036 * exp(-3674 / TR[0]) * PR[0] * PR[2] * PR[6] *
         (1 + 1.7 * PR[5]) /
         ((1 + 0.583 * PR[0] * (1 + 1.7 * PR[5])) * (1 + 6.8 * PR[6]));
  r[1] = 1.9365e5 * exp(-10116 / TR[0]) * PR[0] * (1 + 0.68 * PR[5]) /
         (1 + 0.76 * PR[0] * (1 + 0.68 * PR[5]));
  r[0] = alpha * r[0];
  r[1] = alpha * r[1];
}


x86 assembly (clang, -O0):
<reaction_rate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  0xee2(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xece(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xeba(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xe9d(%rip),%xmm0
ucomisd (%rax),%xmm0
jbe    11a8 <reaction_rate+0x98>
movsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xe9f(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
movsd  %xmm0,-0x30(%rbp)
jmp    122f <reaction_rate+0x11f>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xe48(%rip),%xmm1
subsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xe4d(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xe0f(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0x38(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xe4d(%rip),%xmm0
addsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movsd  0xe2d(%rip),%xmm0
divsd  (%rax),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
movaps %xmm0,%xmm1
movsd  0xe0d(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x30(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm1
movsd  0xdce(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xd92(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  0xdb6(%rip),%xmm1
mulsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm2
movsd  0xd99(%rip),%xmm3
mulsd  %xmm3,%xmm2
movsd  0xd5d(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
movsd  0xd4d(%rip),%xmm2
addsd  %xmm2,%xmm1
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm2
movsd  0xd60(%rip),%xmm3
mulsd  %xmm3,%xmm2
movsd  0xd2c(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movsd  0xd2c(%rip),%xmm0
divsd  (%rax),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
movaps %xmm0,%xmm1
movsd  0xd0c(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm1
movsd  0xcdf(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xccb(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  0xcc7(%rip),%xmm1
mulsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm2
movsd  0xcaa(%rip),%xmm3
mulsd  %xmm3,%xmm2
movsd  0xc96(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
movsd  0xc86(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x8(%rax)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<reaction_rate>:
push   %r15
push   %r14
push   %rbx
sub    $0x20,%rsp
mov    %rdx,%r15
mov    %rsi,%rbx
mov    %rdi,%r14
divsd  0xed6(%rip),%xmm0
divsd  0xed6(%rip),%xmm0
divsd  0xed6(%rip),%xmm0
movsd  (%rdx),%xmm2
movsd  0xed2(%rip),%xmm1
ucomisd %xmm2,%xmm1
jbe    115f <reaction_rate+0x4f>
divsd  0xee4(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
jmp    119f <reaction_rate+0x8f>
addsd  0xeb9(%rip),%xmm2
divsd  0xeb9(%rip),%xmm2
cvttsd2si %xmm2,%edi
movsd  %xmm0,0x18(%rsp)
call   1030 <exp@plt>
movsd  %xmm0,(%rsp)
movsd  0x18(%rsp),%xmm0
divsd  0xe9f(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
mulsd  (%rsp),%xmm0
movapd %xmm0,(%rsp)
movsd  (%r15),%xmm0
addsd  0xe8f(%rip),%xmm0
movsd  %xmm0,(%r15)
movsd  0xe8a(%rip),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
call   1030 <exp@plt>
mulsd  0xe7d(%rip),%xmm0
movsd  (%rbx),%xmm1
movsd  0x28(%rbx),%xmm2
mulsd  %xmm1,%xmm0
mulsd  0x10(%rbx),%xmm0
movsd  0x30(%rbx),%xmm3
mulsd  %xmm3,%xmm0
mulsd  0xe62(%rip),%xmm2
movsd  0xe62(%rip),%xmm4
addsd  %xmm4,%xmm2
mulsd  %xmm2,%xmm0
mulsd  0xe5a(%rip),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm4,%xmm1
mulsd  0xe52(%rip),%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm1,%xmm3
divsd  %xmm3,%xmm0
movsd  %xmm0,(%r14)
movsd  0xe41(%rip),%xmm0
divsd  (%r15),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <exp@plt>
mulsd  0xe33(%rip),%xmm0
movsd  (%rbx),%xmm1
movsd  0x28(%rbx),%xmm2
mulsd  %xmm1,%xmm0
mulsd  0xe26(%rip),%xmm2
movsd  0xdf6(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
mulsd  0xe16(%rip),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
movapd (%rsp),%xmm2
unpcklpd %xmm2,%xmm2
movsd  (%r14),%xmm1
unpcklpd %xmm0,%xmm1
mulpd  %xmm2,%xmm1
movupd %xmm1,(%r14)
add    $0x20,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @reaction_rate(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 double
store double* %0, double** %5
store double* %1, double** %6
store double* %2, double** %7
store double %3, double* %8
%11 = load double, double* %8
%12 = fdiv double %11, 6.000000e+01
%13 = fdiv double %12, 2.400000e+01
%14 = fdiv double %13, 3.650000e+02
store double %14, double* %9
%15 = load double*, double** %7
%16 = getelementptr inbounds double, double* %15, i64 0
%17 = load double, double* %16
%18 = fcmp olt double %17, 1.800000e+02
br i1 %18, label %19, label %25
19: ; preds = %4
%20 = load double, double* %9
%21 = fneg double %20
%22 = fdiv double %21, 0x4011ECF965350DBD
%23 = fptosi double %22 to i32
%24 = call double @exp(i32 %23)
store double %24, double* %10
br label %40
25: ; preds = %4
%26 = load double*, double** %7
%27 = getelementptr inbounds double, double* %26, i64 0
%28 = load double, double* %27
%29 = fsub double %28, 1.800000e+02
%30 = fneg double %29
%31 = fdiv double %30, 5.000000e+01
%32 = fptosi double %31 to i32
%33 = call double @exp(i32 %32)
%34 = load double, double* %9
%35 = fneg double %34
%36 = fdiv double %35, 6.210000e-01
%37 = fptosi double %36 to i32
%38 = call double @exp(i32 %37)
%39 = fmul double %33, %38
store double %39, double* %10
br label %40
40: ; preds = %25, %19
%41 = load double*, double** %7
%42 = getelementptr inbounds double, double* %41, i64 0
%43 = load double, double* %42
%44 = fadd double %43, 2.731500e+02
%45 = load double*, double** %7
%46 = getelementptr inbounds double, double* %45, i64 0
store double %44, double* %46
%47 = load double*, double** %7
%48 = getelementptr inbounds double, double* %47, i64 0
%49 = load double, double* %48
%50 = fdiv double -3.674000e+03, %49
%51 = fptosi double %50 to i32
%52 = call double @exp(i32 %51)
%53 = fmul double 1.036000e-01, %52
%54 = load double*, double** %6
%55 = getelementptr inbounds double, double* %54, i64 0
%56 = load double, double* %55
%57 = fmul double %53, %56
%58 = load double*, double** %6
%59 = getelementptr inbounds double, double* %58, i64 2
%60 = load double, double* %59
%61 = fmul double %57, %60
%62 = load double*, double** %6
%63 = getelementptr inbounds double, double* %62, i64 6
%64 = load double, double* %63
%65 = fmul double %61, %64
%66 = load double*, double** %6
%67 = getelementptr inbounds double, double* %66, i64 5
%68 = load double, double* %67
%69 = call double @llvm.fmuladd.f64(double 1.700000e+00, double %68, double 1.000000e+00)
%70 = fmul double %65, %69
%71 = load double*, double** %6
%72 = getelementptr inbounds double, double* %71, i64 0
%73 = load double, double* %72
%74 = fmul double 5.830000e-01, %73
%75 = load double*, double** %6
%76 = getelementptr inbounds double, double* %75, i64 5
%77 = load double, double* %76
%78 = call double @llvm.fmuladd.f64(double 1.700000e+00, double %77, double 1.000000e+00)
%79 = call double @llvm.fmuladd.f64(double %74, double %78, double 1.000000e+00)
%80 = load double*, double** %6
%81 = getelementptr inbounds double, double* %80, i64 6
%82 = load double, double* %81
%83 = call double @llvm.fmuladd.f64(double 6.800000e+00, double %82, double 1.000000e+00)
%84 = fmul double %79, %83
%85 = fdiv double %70, %84
%86 = load double*, double** %5
%87 = getelementptr inbounds double, double* %86, i64 0
store double %85, double* %87
%88 = load double*, double** %7
%89 = getelementptr inbounds double, double* %88, i64 0
%90 = load double, double* %89
%91 = fdiv double -1.011600e+04, %90
%92 = fptosi double %91 to i32
%93 = call double @exp(i32 %92)
%94 = fmul double 1.936500e+05, %93
%95 = load double*, double** %6
%96 = getelementptr inbounds double, double* %95, i64 0
%97 = load double, double* %96
%98 = fmul double %94, %97
%99 = load double*, double** %6
%100 = getelementptr inbounds double, double* %99, i64 5
%101 = load double, double* %100
%102 = call double @llvm.fmuladd.f64(double 6.800000e-01, double %101, double 1.000000e+00)
%103 = fmul double %98, %102
%104 = load double*, double** %6
%105 = getelementptr inbounds double, double* %104, i64 0
%106 = load double, double* %105
%107 = fmul double 7.600000e-01, %106
%108 = load double*, double** %6
%109 = getelementptr inbounds double, double* %108, i64 5
%110 = load double, double* %109
%111 = call double @llvm.fmuladd.f64(double 6.800000e-01, double %110, double 1.000000e+00)
%112 = call double @llvm.fmuladd.f64(double %107, double %111, double 1.000000e+00)
%113 = fdiv double %103, %112
%114 = load double*, double** %5
%115 = getelementptr inbounds double, double* %114, i64 1
store double %113, double* %115
%116 = load double, double* %10
%117 = load double*, double** %5
%118 = getelementptr inbounds double, double* %117, i64 0
%119 = load double, double* %118
%120 = fmul double %116, %119
%121 = load double*, double** %5
%122 = getelementptr inbounds double, double* %121, i64 0
store double %120, double* %122
%123 = load double, double* %10
%124 = load double*, double** %5
%125 = getelementptr inbounds double, double* %124, i64 1
%126 = load double, double* %125
%127 = fmul double %123, %126
%128 = load double*, double** %5
%129 = getelementptr inbounds double, double* %128, i64 1
store double %127, double* %129
ret void
}

LLVM IR (-O2):
define void @reaction_rate(double* nocapture %0, double* nocapture %1, double* nocapture %2, double %3) {
%5 = fdiv double %3, 6.000000e+01
%6 = fdiv double %5, 2.400000e+01
%7 = fdiv double %6, 3.650000e+02
%8 = load double, double* %2
%9 = fcmp olt double %8, 1.800000e+02
br i1 %9, label %10, label %14
10: ; preds = %4
%11 = fdiv double %7, 0xC011ECF965350DBD
%12 = fptosi double %11 to i32
%13 = tail call double @exp(i32 %12)
br label %23
14: ; preds = %4
%15 = fadd double %8, -1.800000e+02
%16 = fdiv double %15, -5.000000e+01
%17 = fptosi double %16 to i32
%18 = tail call double @exp(i32 %17)
%19 = fdiv double %7, -6.210000e-01
%20 = fptosi double %19 to i32
%21 = tail call double @exp(i32 %20)
%22 = fmul double %18, %21
br label %23
23: ; preds = %14, %10
%24 = phi double [ %13, %10 ], [ %22, %14 ]
%25 = load double, double* %2
%26 = fadd double %25, 2.731500e+02
store double %26, double* %2
%27 = fdiv double -3.674000e+03, %26
%28 = fptosi double %27 to i32
%29 = tail call double @exp(i32 %28)
%30 = fmul double %29, 1.036000e-01
%31 = load double, double* %1
%32 = fmul double %30, %31
%33 = getelementptr inbounds double, double* %1, i64 2
%34 = load double, double* %33
%35 = fmul double %32, %34
%36 = getelementptr inbounds double, double* %1, i64 6
%37 = load double, double* %36
%38 = fmul double %35, %37
%39 = getelementptr inbounds double, double* %1, i64 5
%40 = load double, double* %39
%41 = tail call double @llvm.fmuladd.f64(double %40, double 1.700000e+00, double 1.000000e+00)
%42 = fmul double %38, %41
%43 = fmul double %31, 5.830000e-01
%44 = tail call double @llvm.fmuladd.f64(double %43, double %41, double 1.000000e+00)
%45 = tail call double @llvm.fmuladd.f64(double %37, double 6.800000e+00, double 1.000000e+00)
%46 = fmul double %45, %44
%47 = fdiv double %42, %46
store double %47, double* %0
%48 = load double, double* %2
%49 = fdiv double -1.011600e+04, %48
%50 = fptosi double %49 to i32
%51 = tail call double @exp(i32 %50)
%52 = fmul double %51, 1.936500e+05
%53 = load double, double* %1
%54 = fmul double %52, %53
%55 = load double, double* %39
%56 = tail call double @llvm.fmuladd.f64(double %55, double 6.800000e-01, double 1.000000e+00)
%57 = fmul double %54, %56
%58 = fmul double %53, 7.600000e-01
%59 = tail call double @llvm.fmuladd.f64(double %58, double %56, double 1.000000e+00)
%60 = fdiv double %57, %59
%61 = load double, double* %0
%62 = insertelement <2 x double> poison, double %24, i64 0
%63 = shufflevector <2 x double> %62, <2 x double> poison, <2 x i32> zeroinitializer
%64 = insertelement <2 x double> poison, double %61, i64 0
%65 = insertelement <2 x double> %64, double %60, i64 1
%66 = fmul <2 x double> %63, %65
%67 = bitcast double* %0 to <2 x double>*
store <2 x double> %66, <2 x double>* %67
ret void
}