#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 paw_Def_Integr(double, double *, double, int);
int paw_N_LogGrid();
double *paw_alloc_LogGrid();
int /*<<< orphan*/ paw_dealloc_LogGrid(double *);
double paw_log_amesh_LogGrid();
double *paw_r_LogGrid();
double paw_get_kinetic_energy(int num_states, int *l, double *fill,
                              double **psi, double **psi_prime) {
  int i;
  int k;
  int Ngrid;
  double ekin;
  double ekin_total;
  double log_amesh;
  double *f;
  double *r;
  Ngrid = paw_N_LogGrid();
  log_amesh = paw_log_amesh_LogGrid();
  r = paw_r_LogGrid();
  f = paw_alloc_LogGrid();
  ekin_total = 0.0;
  for (i = 0; i <= num_states - 1; i++) {
    for (k = 0; k <= Ngrid - 1; k++) {
      f[k] = 0.5 * psi_prime[i][k] / (r[k] * log_amesh) * psi_prime[i][k] /
                 (r[k] * log_amesh) +
             0.5 * l[i] * (l[i] + 1) / (r[k] * r[k]) * psi[i][k] * psi[i][k];
    }
    ekin = paw_Def_Integr(0.0, f, 0.0, Ngrid - 1);
    ekin_total = ekin_total + ekin * fill[i];
  }
  paw_dealloc_LogGrid(f);
  return ekin_total;
}


x86 assembly (clang, -O0):
<paw_get_kinetic_energy>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    $0x0,%al
call   1030 <paw_N_LogGrid@plt>
mov    %eax,-0x34(%rbp)
mov    $0x0,%al
call   1040 <paw_log_amesh_LogGrid@plt>
movsd  %xmm0,-0x50(%rbp)
mov    $0x0,%al
call   1050 <paw_r_LogGrid@plt>
mov    %rax,-0x60(%rbp)
mov    $0x0,%al
call   1070 <paw_alloc_LogGrid@plt>
mov    %rax,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1323 <paw_get_kinetic_energy+0x1c3>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     12d7 <paw_get_kinetic_energy+0x177>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  0xe09(%rip),%xmm1
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x60(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x50(%rbp),%xmm0
divsd  %xmm0,%xmm1
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x60(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x50(%rbp),%xmm0
divsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm2
movsd  0xdae(%rip),%xmm0
mulsd  %xmm2,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
mulsd  %xmm2,%xmm0
mov    -0x60(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x60(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm2
divsd  %xmm2,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x58(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11ce <paw_get_kinetic_energy+0x6e>
mov    -0x58(%rbp),%rdi
mov    -0x34(%rbp),%esi
sub    $0x1,%esi
xorps  %xmm1,%xmm1
movaps %xmm1,%xmm0
call   1060 <paw_Def_Integr@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b6 <paw_get_kinetic_energy+0x56>
mov    -0x58(%rbp),%rdi
call   1080 <paw_dealloc_LogGrid@plt>
movsd  -0x48(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<paw_get_kinetic_energy>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r8,0x28(%rsp)
mov    %rcx,0x38(%rsp)
mov    %rdx,0x48(%rsp)
mov    %rsi,0x30(%rsp)
mov    %edi,%ebx
xor    %eax,%eax
call   1030 <paw_N_LogGrid@plt>
mov    %eax,%r15d
xor    %eax,%eax
call   1040 <paw_log_amesh_LogGrid@plt>
movapd %xmm0,0x50(%rsp)
xor    %eax,%eax
call   1050 <paw_r_LogGrid@plt>
mov    %rax,%r14
xor    %eax,%eax
call   1070 <paw_alloc_LogGrid@plt>
mov    %rax,%rbp
test   %ebx,%ebx
jle    13e8 <paw_get_kinetic_energy+0x288>
lea    -0x1(%r15),%eax
mov    %eax,0x1c(%rsp)
mov    %ebx,%eax
mov    %rax,0x40(%rsp)
mov    %r15d,%r13d
lea    0x0(,%r13,8),%rcx
add    %rbp,%rcx
lea    (%r14,%r13,8),%rax
cmp    %rax,%rbp
setb   %al
mov    %rcx,0x20(%rsp)
cmp    %rcx,%r14
setb   %cl
and    %al,%cl
mov    %cl,0xf(%rsp)
mov    %r13d,%r12d
and    $0xfffffffe,%r12d
movapd 0x50(%rsp),%xmm6
movapd %xmm6,%xmm8
unpcklpd %xmm6,%xmm8
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %r15,%rbx
xor    %r15d,%r15d
movsd  0xde2(%rip),%xmm7
movapd 0xde9(%rip),%xmm9
movapd %xmm8,0x60(%rsp)
jmp    128c <paw_get_kinetic_energy+0x12c>
xorpd  %xmm0,%xmm0
xorpd  %xmm1,%xmm1
mov    %rbp,%rdi
mov    0x1c(%rsp),%esi
call   1060 <paw_Def_Integr@plt>
movapd 0xdc3(%rip),%xmm9
movsd  0xdab(%rip),%xmm7
movapd 0x60(%rsp),%xmm8
movapd 0x50(%rsp),%xmm6
mov    0x48(%rsp),%rax
mulsd  (%rax,%r15,8),%xmm0
movsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
add    $0x1,%r15
cmp    0x40(%rsp),%r15
je     13f2 <paw_get_kinetic_energy+0x292>
test   %ebx,%ebx
jle    1230 <paw_get_kinetic_energy+0xd0>
mov    0x30(%rsp),%rax
mov    (%rax,%r15,4),%ecx
lea    0x1(%rcx),%edx
mov    0x28(%rsp),%rax
mov    (%rax,%r15,8),%rax
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
mulsd  %xmm7,%xmm1
mulsd  %xmm1,%xmm0
mov    0x38(%rsp),%rcx
mov    (%rcx,%r15,8),%rcx
cmp    $0x1,%ebx
jne    12d0 <paw_get_kinetic_energy+0x170>
xor    %edx,%edx
jmp    1390 <paw_get_kinetic_energy+0x230>
lea    (%rax,%r13,8),%rdx
lea    (%rcx,%r13,8),%r8
cmp    %rdx,%rbp
setb   %dil
mov    0x20(%rsp),%rsi
cmp    %rsi,%rax
setb   %dl
and    %dil,%dl
or     0xf(%rsp),%dl
cmp    %r8,%rbp
setb   %dil
cmp    %rsi,%rcx
setb   %sil
test   %dl,%dl
jne    1382 <paw_get_kinetic_energy+0x222>
mov    $0x0,%edx
and    %sil,%dil
jne    1390 <paw_get_kinetic_energy+0x230>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movupd (%rax,%rdx,8),%xmm2
movapd %xmm2,%xmm3
mulpd  %xmm9,%xmm3
movupd (%r14,%rdx,8),%xmm4
movapd %xmm8,%xmm5
mulpd  %xmm4,%xmm5
divpd  %xmm5,%xmm3
mulpd  %xmm2,%xmm3
divpd  %xmm5,%xmm3
mulpd  %xmm4,%xmm4
movapd %xmm1,%xmm2
divpd  %xmm4,%xmm2
movupd (%rcx,%rdx,8),%xmm4
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm2
addpd  %xmm3,%xmm2
movupd %xmm2,0x0(%rbp,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%r12
jne    1320 <paw_get_kinetic_energy+0x1c0>
mov    %r12,%rdx
cmp    %r13,%r12
jne    1390 <paw_get_kinetic_energy+0x230>
jmp    1230 <paw_get_kinetic_energy+0xd0>
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rax,%rdx,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm7,%xmm2
movsd  (%r14,%rdx,8),%xmm3
movapd %xmm6,%xmm4
mulsd  %xmm3,%xmm4
divsd  %xmm4,%xmm2
mulsd  %xmm1,%xmm2
divsd  %xmm4,%xmm2
mulsd  %xmm3,%xmm3
movapd %xmm0,%xmm1
divsd  %xmm3,%xmm1
movsd  (%rcx,%rdx,8),%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x0(%rbp,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r13
jne    1390 <paw_get_kinetic_energy+0x230>
jmp    1230 <paw_get_kinetic_energy+0xd0>
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %rbp,%rdi
call   1080 <paw_dealloc_LogGrid@plt>
movsd  0x10(%rsp),%xmm0
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @paw_get_kinetic_energy(i32 %0, i32* %1, double* %2, double** %3, double** %4) {
%6 = alloca i32
%7 = alloca i32*
%8 = alloca double*
%9 = alloca double**
%10 = alloca double**
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double*
%18 = alloca double*
store i32 %0, i32* %6
store i32* %1, i32** %7
store double* %2, double** %8
store double** %3, double*** %9
store double** %4, double*** %10
%19 = call i32 (...) @paw_N_LogGrid()
store i32 %19, i32* %13
%20 = call double (...) @paw_log_amesh_LogGrid()
store double %20, double* %16
%21 = call double* (...) @paw_r_LogGrid()
store double* %21, double** %18
%22 = call double* (...) @paw_alloc_LogGrid()
store double* %22, double** %17
store double 0.000000e+00, double* %15
store i32 0, i32* %11
br label %23
23: ; preds = %138, %5
%24 = load i32, i32* %11
%25 = load i32, i32* %6
%26 = sub nsw i32 %25, 1
%27 = icmp sle i32 %24, %26
br i1 %27, label %28, label %141
28: ; preds = %23
store i32 0, i32* %12
br label %29
29: ; preds = %122, %28
%30 = load i32, i32* %12
%31 = load i32, i32* %13
%32 = sub nsw i32 %31, 1
%33 = icmp sle i32 %30, %32
br i1 %33, label %34, label %125
34: ; preds = %29
%35 = load double**, double*** %10
%36 = load i32, i32* %11
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double*, double** %35, i64 %37
%39 = load double*, double** %38
%40 = load i32, i32* %12
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %39, i64 %41
%43 = load double, double* %42
%44 = fmul double 5.000000e-01, %43
%45 = load double*, double** %18
%46 = load i32, i32* %12
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %45, i64 %47
%49 = load double, double* %48
%50 = load double, double* %16
%51 = fmul double %49, %50
%52 = fdiv double %44, %51
%53 = load double**, double*** %10
%54 = load i32, i32* %11
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double*, double** %53, i64 %55
%57 = load double*, double** %56
%58 = load i32, i32* %12
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %57, i64 %59
%61 = load double, double* %60
%62 = fmul double %52, %61
%63 = load double*, double** %18
%64 = load i32, i32* %12
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = load double, double* %16
%69 = fmul double %67, %68
%70 = fdiv double %62, %69
%71 = load i32*, i32** %7
%72 = load i32, i32* %11
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds i32, i32* %71, i64 %73
%75 = load i32, i32* %74
%76 = sitofp i32 %75 to double
%77 = fmul double 5.000000e-01, %76
%78 = load i32*, i32** %7
%79 = load i32, i32* %11
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds i32, i32* %78, i64 %80
%82 = load i32, i32* %81
%83 = add nsw i32 %82, 1
%84 = sitofp i32 %83 to double
%85 = fmul double %77, %84
%86 = load double*, double** %18
%87 = load i32, i32* %12
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %86, i64 %88
%90 = load double, double* %89
%91 = load double*, double** %18
%92 = load i32, i32* %12
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %91, i64 %93
%95 = load double, double* %94
%96 = fmul double %90, %95
%97 = fdiv double %85, %96
%98 = load double**, double*** %9
%99 = load i32, i32* %11
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double*, double** %98, i64 %100
%102 = load double*, double** %101
%103 = load i32, i32* %12
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %102, i64 %104
%106 = load double, double* %105
%107 = fmul double %97, %106
%108 = load double**, double*** %9
%109 = load i32, i32* %11
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double*, double** %108, i64 %110
%112 = load double*, double** %111
%113 = load i32, i32* %12
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %112, i64 %114
%116 = load double, double* %115
%117 = call double @llvm.fmuladd.f64(double %107, double %116, double %70)
%118 = load double*, double** %17
%119 = load i32, i32* %12
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %118, i64 %120
store double %117, double* %121
br label %122
122: ; preds = %34
%123 = load i32, i32* %12
%124 = add nsw i32 %123, 1
store i32 %124, i32* %12
br label %29
125: ; preds = %29
%126 = load double*, double** %17
%127 = load i32, i32* %13
%128 = sub nsw i32 %127, 1
%129 = call double @paw_Def_Integr(double 0.000000e+00, double* %126, double 0.000000e+00, i32 %128)
store double %129, double* %14
%130 = load double, double* %15
%131 = load double, double* %14
%132 = load double*, double** %8
%133 = load i32, i32* %11
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = call double @llvm.fmuladd.f64(double %131, double %136, double %130)
store double %137, double* %15
br label %138
138: ; preds = %125
%139 = load i32, i32* %11
%140 = add nsw i32 %139, 1
store i32 %140, i32* %11
br label %23
141: ; preds = %23
%142 = load double*, double** %17
%143 = call i32 @paw_dealloc_LogGrid(double* %142)
%144 = load double, double* %15
ret double %144
}

LLVM IR (-O2):
define double @paw_get_kinetic_energy(i32 %0, i32* nocapture %1, double* nocapture %2, double** nocapture %3, double** nocapture %4) {
%6 = tail call i32 (...) @paw_N_LogGrid()
%7 = tail call double (...) @paw_log_amesh_LogGrid()
%8 = tail call double* (...) @paw_r_LogGrid()
%9 = tail call double* (...) @paw_alloc_LogGrid()
%10 = icmp sgt i32 %0, 0
br i1 %10, label %11, label %109
11: ; preds = %5
%12 = icmp sgt i32 %6, 0
%13 = add nsw i32 %6, -1
%14 = zext i32 %0 to i64
%15 = zext i32 %6 to i64
%16 = getelementptr double, double* %9, i64 %15
%17 = getelementptr double, double* %8, i64 %15
%18 = icmp eq i32 %6, 1
%19 = icmp ult double* %9, %17
%20 = icmp ult double* %8, %16
%21 = and i1 %19, %20
%22 = and i64 %15, 4294967294
%23 = insertelement <2 x double> poison, double %7, i64 0
%24 = shufflevector <2 x double> %23, <2 x double> poison, <2 x i32> zeroinitializer
%25 = icmp eq i64 %22, %15
br label %26
26: ; preds = %11, %102
%27 = phi i64 [ 0, %11 ], [ %107, %102 ]
%28 = phi double [ 0.000000e+00, %11 ], [ %106, %102 ]
br i1 %12, label %29, label %102
29: ; preds = %26
%30 = getelementptr inbounds double*, double** %4, i64 %27
%31 = load double*, double** %30
%32 = getelementptr inbounds i32, i32* %1, i64 %27
%33 = load i32, i32* %32
%34 = sitofp i32 %33 to double
%35 = fmul double %34, 5.000000e-01
%36 = add nsw i32 %33, 1
%37 = sitofp i32 %36 to double
%38 = fmul double %35, %37
%39 = getelementptr inbounds double*, double** %3, i64 %27
%40 = load double*, double** %39
br i1 %18, label %80, label %41
41: ; preds = %29
%42 = getelementptr double, double* %31, i64 %15
%43 = getelementptr double, double* %40, i64 %15
%44 = icmp ult double* %9, %42
%45 = icmp ult double* %31, %16
%46 = and i1 %44, %45
%47 = or i1 %46, %21
%48 = icmp ult double* %9, %43
%49 = icmp ult double* %40, %16
%50 = and i1 %48, %49
%51 = or i1 %47, %50
br i1 %51, label %80, label %52
52: ; preds = %41
%53 = insertelement <2 x double> poison, double %38, i64 0
%54 = shufflevector <2 x double> %53, <2 x double> poison, <2 x i32> zeroinitializer
br label %55
55: ; preds = %55, %52
%56 = phi i64 [ 0, %52 ], [ %77, %55 ]
%57 = getelementptr inbounds double, double* %31, i64 %56
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = fmul <2 x double> %59, <double 5.000000e-01, double 5.000000e-01>
%61 = getelementptr inbounds double, double* %8, i64 %56
%62 = bitcast double* %61 to <2 x double>*
%63 = load <2 x double>, <2 x double>* %62
%64 = fmul <2 x double> %24, %63
%65 = fdiv <2 x double> %60, %64
%66 = fmul <2 x double> %59, %65
%67 = fdiv <2 x double> %66, %64
%68 = fmul <2 x double> %63, %63
%69 = fdiv <2 x double> %54, %68
%70 = getelementptr inbounds double, double* %40, i64 %56
%71 = bitcast double* %70 to <2 x double>*
%72 = load <2 x double>, <2 x double>* %71
%73 = fmul <2 x double> %72, %69
%74 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %73, <2 x double> %72, <2 x double> %67)
%75 = getelementptr inbounds double, double* %9, i64 %56
%76 = bitcast double* %75 to <2 x double>*
store <2 x double> %74, <2 x double>* %76
%77 = add nuw i64 %56, 2
%78 = icmp eq i64 %77, %22
br i1 %78, label %79, label %55
79: ; preds = %55
br i1 %25, label %102, label %80
80: ; preds = %41, %29, %79
%81 = phi i64 [ 0, %41 ], [ 0, %29 ], [ %22, %79 ]
br label %82
82: ; preds = %80, %82
%83 = phi i64 [ %100, %82 ], [ %81, %80 ]
%84 = getelementptr inbounds double, double* %31, i64 %83
%85 = load double, double* %84
%86 = fmul double %85, 5.000000e-01
%87 = getelementptr inbounds double, double* %8, i64 %83
%88 = load double, double* %87
%89 = fmul double %7, %88
%90 = fdiv double %86, %89
%91 = fmul double %85, %90
%92 = fdiv double %91, %89
%93 = fmul double %88, %88
%94 = fdiv double %38, %93
%95 = getelementptr inbounds double, double* %40, i64 %83
%96 = load double, double* %95
%97 = fmul double %96, %94
%98 = tail call double @llvm.fmuladd.f64(double %97, double %96, double %92)
%99 = getelementptr inbounds double, double* %9, i64 %83
store double %98, double* %99
%100 = add nuw nsw i64 %83, 1
%101 = icmp eq i64 %100, %15
br i1 %101, label %102, label %82
102: ; preds = %82, %79, %26
%103 = tail call double @paw_Def_Integr(double 0.000000e+00, double* %9, double 0.000000e+00, i32 %13)
%104 = getelementptr inbounds double, double* %2, i64 %27
%105 = load double, double* %104
%106 = tail call double @llvm.fmuladd.f64(double %103, double %105, double %28)
%107 = add nuw nsw i64 %27, 1
%108 = icmp eq i64 %107, %14
br i1 %108, label %109, label %26
109: ; preds = %102, %5
%110 = phi double [ 0.000000e+00, %5 ], [ %106, %102 ]
%111 = tail call i32 @paw_dealloc_LogGrid(double* %9)
ret double %110
}