void compute_abc_nobc(double *dx, double *dfactor, double *delj, double *MInt,
                      double *V, double dt, int N, double *a, double *b,
                      double *c) {
  int ii;
  double atemp, ctemp;
  a[0] = 0;
  c[N - 1] = 0;
  for (ii = 0; ii < N; ii++)
    b[ii] = 1. / dt;
  for (ii = 0; ii < N - 1; ii++) {
    atemp = MInt[ii] * delj[ii] + V[ii] / (2 * dx[ii]);
    a[ii + 1] = -dfactor[ii + 1] * atemp;
    b[ii] += dfactor[ii] * atemp;
    ctemp = -MInt[ii] * (1 - delj[ii]) + V[ii + 1] / (2 * dx[ii]);
    b[ii + 1] += dfactor[ii + 1] * ctemp;
    c[ii] = -dfactor[ii] * ctemp;
  }
}


x86 assembly (clang, -O0):
<compute_abc_nobc>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
mov    0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1188 <compute_abc_nobc+0x88>
movsd  0xea0(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1154 <compute_abc_nobc+0x54>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1344 <compute_abc_nobc+0x244>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  0xe29(%rip),%xmm3
mulsd  (%rax,%rcx,8),%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x40(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  0xd7b(%rip),%xmm2
subsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  0xd4c(%rip),%xmm3
mulsd  (%rax,%rcx,8),%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
mov    0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x48(%rbp),%xmm0
mov    0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118f <compute_abc_nobc+0x8f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_abc_nobc>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    0x40(%rsp),%r10
mov    0x30(%rsp),%r11
movq   $0x0,(%r11)
movslq %r9d,%rax
movq   $0x0,-0x8(%r10,%rax,8)
test   %r9d,%r9d
jle    1335 <compute_abc_nobc+0x235>
mov    0x38(%rsp),%r12
movsd  0xec4(%rip),%xmm1
movapd %xmm1,%xmm2
divsd  %xmm0,%xmm2
mov    %r9d,%r14d
cmp    $0x4,%r9d
jae    1154 <compute_abc_nobc+0x54>
xor    %ebx,%ebx
jmp    1260 <compute_abc_nobc+0x160>
mov    %r14d,%ebx
and    $0xfffffffc,%ebx
movapd %xmm2,%xmm0
unpcklpd %xmm2,%xmm0
lea    -0x4(%rbx),%rax
mov    %rax,%r13
shr    $0x2,%r13
add    $0x1,%r13
mov    %r13d,%r15d
and    $0x7,%r15d
cmp    $0x1c,%rax
jae    1185 <compute_abc_nobc+0x85>
xor    %eax,%eax
jmp    1225 <compute_abc_nobc+0x125>
and    $0xfffffffffffffff8,%r13
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movupd %xmm0,(%r12,%rax,8)
movupd %xmm0,0x10(%r12,%rax,8)
movupd %xmm0,0x20(%r12,%rax,8)
movupd %xmm0,0x30(%r12,%rax,8)
movupd %xmm0,0x40(%r12,%rax,8)
movupd %xmm0,0x50(%r12,%rax,8)
movupd %xmm0,0x60(%r12,%rax,8)
movupd %xmm0,0x70(%r12,%rax,8)
movupd %xmm0,0x80(%r12,%rax,8)
movupd %xmm0,0x90(%r12,%rax,8)
movupd %xmm0,0xa0(%r12,%rax,8)
movupd %xmm0,0xb0(%r12,%rax,8)
movupd %xmm0,0xc0(%r12,%rax,8)
movupd %xmm0,0xd0(%r12,%rax,8)
movupd %xmm0,0xe0(%r12,%rax,8)
movupd %xmm0,0xf0(%r12,%rax,8)
add    $0x20,%rax
add    $0xfffffffffffffff8,%r13
jne    1190 <compute_abc_nobc+0x90>
test   %r15,%r15
je     1257 <compute_abc_nobc+0x157>
lea    (%r12,%rax,8),%r13
add    $0x10,%r13
shl    $0x5,%r15
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%r13,%rax,1)
movupd %xmm0,0x0(%r13,%rax,1)
add    $0x20,%rax
cmp    %rax,%r15
jne    1240 <compute_abc_nobc+0x140>
cmp    %r14,%rbx
je     126f <compute_abc_nobc+0x16f>
nopl   0x0(%rax)
movsd  %xmm2,(%r12,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1260 <compute_abc_nobc+0x160>
cmp    $0x2,%r9d
jl     1335 <compute_abc_nobc+0x235>
add    $0xffffffff,%r9d
mov    %r9d,%r9d
shl    $0x3,%r9
xor    %ebx,%ebx
movapd 0xd82(%rip),%xmm0
xchg   %ax,%ax
movsd  (%rcx,%rbx,1),%xmm2
movsd  (%r8,%rbx,1),%xmm3
movsd  (%rdi,%rbx,1),%xmm4
addsd  %xmm4,%xmm4
divsd  %xmm4,%xmm3
mulsd  (%rdx,%rbx,1),%xmm2
addsd  %xmm3,%xmm2
movsd  0x8(%rsi,%rbx,1),%xmm3
xorpd  %xmm0,%xmm3
mulsd  %xmm2,%xmm3
movsd  %xmm3,0x8(%r11,%rbx,1)
mulsd  (%rsi,%rbx,1),%xmm2
addsd  (%r12,%rbx,1),%xmm2
movsd  %xmm2,(%r12,%rbx,1)
movapd %xmm1,%xmm2
subsd  (%rdx,%rbx,1),%xmm2
movsd  0x8(%r8,%rbx,1),%xmm3
movsd  (%rdi,%rbx,1),%xmm4
addsd  %xmm4,%xmm4
divsd  %xmm4,%xmm3
mulsd  (%rcx,%rbx,1),%xmm2
subsd  %xmm2,%xmm3
movsd  0x8(%rsi,%rbx,1),%xmm2
mulsd  %xmm3,%xmm2
addsd  0x8(%r12,%rbx,1),%xmm2
movsd  %xmm2,0x8(%r12,%rbx,1)
movsd  (%rsi,%rbx,1),%xmm2
xorpd  %xmm0,%xmm2
mulsd  %xmm3,%xmm2
movsd  %xmm2,(%r10,%rbx,1)
add    $0x8,%rbx
cmp    %rbx,%r9
jne    1290 <compute_abc_nobc+0x190>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @compute_abc_nobc(double* nocapture %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, double* nocapture %4, double %5, i32 %6, double* nocapture writeonly %7, double* nocapture %8, double* nocapture writeonly %9) {
store double 0.000000e+00, double* %7
%11 = add nsw i32 %6, -1
%12 = sext i32 %11 to i64
%13 = getelementptr inbounds double, double* %9, i64 %12
store double 0.000000e+00, double* %13
%14 = icmp sgt i32 %6, 0
br i1 %14, label %15, label %146
15: ; preds = %10
%16 = fdiv double 1.000000e+00, %5
%17 = zext i32 %6 to i64
%18 = icmp ult i32 %6, 4
br i1 %18, label %92, label %19
19: ; preds = %15
%20 = and i64 %17, 4294967292
%21 = insertelement <2 x double> poison, double %16, i64 0
%22 = shufflevector <2 x double> %21, <2 x double> poison, <2 x i32> zeroinitializer
%23 = insertelement <2 x double> poison, double %16, i64 0
%24 = shufflevector <2 x double> %23, <2 x double> poison, <2 x i32> zeroinitializer
%25 = add nsw i64 %20, -4
%26 = lshr exact i64 %25, 2
%27 = add nuw nsw i64 %26, 1
%28 = and i64 %27, 7
%29 = icmp ult i64 %25, 28
br i1 %29, label %77, label %30
30: ; preds = %19
%31 = and i64 %27, 9223372036854775800
br label %32
32: ; preds = %32, %30
%33 = phi i64 [ 0, %30 ], [ %74, %32 ]
%34 = phi i64 [ 0, %30 ], [ %75, %32 ]
%35 = getelementptr inbounds double, double* %8, i64 %33
%36 = bitcast double* %35 to <2 x double>*
store <2 x double> %22, <2 x double>* %36
%37 = getelementptr inbounds double, double* %35, i64 2
%38 = bitcast double* %37 to <2 x double>*
store <2 x double> %24, <2 x double>* %38
%39 = or i64 %33, 4
%40 = getelementptr inbounds double, double* %8, i64 %39
%41 = bitcast double* %40 to <2 x double>*
store <2 x double> %22, <2 x double>* %41
%42 = getelementptr inbounds double, double* %40, i64 2
%43 = bitcast double* %42 to <2 x double>*
store <2 x double> %24, <2 x double>* %43
%44 = or i64 %33, 8
%45 = getelementptr inbounds double, double* %8, i64 %44
%46 = bitcast double* %45 to <2 x double>*
store <2 x double> %22, <2 x double>* %46
%47 = getelementptr inbounds double, double* %45, i64 2
%48 = bitcast double* %47 to <2 x double>*
store <2 x double> %24, <2 x double>* %48
%49 = or i64 %33, 12
%50 = getelementptr inbounds double, double* %8, i64 %49
%51 = bitcast double* %50 to <2 x double>*
store <2 x double> %22, <2 x double>* %51
%52 = getelementptr inbounds double, double* %50, i64 2
%53 = bitcast double* %52 to <2 x double>*
store <2 x double> %24, <2 x double>* %53
%54 = or i64 %33, 16
%55 = getelementptr inbounds double, double* %8, i64 %54
%56 = bitcast double* %55 to <2 x double>*
store <2 x double> %22, <2 x double>* %56
%57 = getelementptr inbounds double, double* %55, i64 2
%58 = bitcast double* %57 to <2 x double>*
store <2 x double> %24, <2 x double>* %58
%59 = or i64 %33, 20
%60 = getelementptr inbounds double, double* %8, i64 %59
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %22, <2 x double>* %61
%62 = getelementptr inbounds double, double* %60, i64 2
%63 = bitcast double* %62 to <2 x double>*
store <2 x double> %24, <2 x double>* %63
%64 = or i64 %33, 24
%65 = getelementptr inbounds double, double* %8, i64 %64
%66 = bitcast double* %65 to <2 x double>*
store <2 x double> %22, <2 x double>* %66
%67 = getelementptr inbounds double, double* %65, i64 2
%68 = bitcast double* %67 to <2 x double>*
store <2 x double> %24, <2 x double>* %68
%69 = or i64 %33, 28
%70 = getelementptr inbounds double, double* %8, i64 %69
%71 = bitcast double* %70 to <2 x double>*
store <2 x double> %22, <2 x double>* %71
%72 = getelementptr inbounds double, double* %70, i64 2
%73 = bitcast double* %72 to <2 x double>*
store <2 x double> %24, <2 x double>* %73
%74 = add nuw i64 %33, 32
%75 = add i64 %34, 8
%76 = icmp eq i64 %75, %31
br i1 %76, label %77, label %32
77: ; preds = %32, %19
%78 = phi i64 [ 0, %19 ], [ %74, %32 ]
%79 = icmp eq i64 %28, 0
br i1 %79, label %90, label %80
80: ; preds = %77, %80
%81 = phi i64 [ %87, %80 ], [ %78, %77 ]
%82 = phi i64 [ %88, %80 ], [ 0, %77 ]
%83 = getelementptr inbounds double, double* %8, i64 %81
%84 = bitcast double* %83 to <2 x double>*
store <2 x double> %22, <2 x double>* %84
%85 = getelementptr inbounds double, double* %83, i64 2
%86 = bitcast double* %85 to <2 x double>*
store <2 x double> %24, <2 x double>* %86
%87 = add nuw i64 %81, 4
%88 = add i64 %82, 1
%89 = icmp eq i64 %88, %28
br i1 %89, label %90, label %80
90: ; preds = %80, %77
%91 = icmp eq i64 %20, %17
br i1 %91, label %94, label %92
92: ; preds = %15, %90
%93 = phi i64 [ 0, %15 ], [ %20, %90 ]
br label %98
94: ; preds = %98, %90
%95 = icmp sgt i32 %6, 1
br i1 %95, label %96, label %146
96: ; preds = %94
%97 = zext i32 %11 to i64
br label %103
98: ; preds = %92, %98
%99 = phi i64 [ %101, %98 ], [ %93, %92 ]
%100 = getelementptr inbounds double, double* %8, i64 %99
store double %16, double* %100
%101 = add nuw nsw i64 %99, 1
%102 = icmp eq i64 %101, %17
br i1 %102, label %94, label %98
103: ; preds = %96, %103
%104 = phi i64 [ 0, %96 ], [ %116, %103 ]
%105 = getelementptr inbounds double, double* %3, i64 %104
%106 = load double, double* %105
%107 = getelementptr inbounds double, double* %2, i64 %104
%108 = load double, double* %107
%109 = getelementptr inbounds double, double* %4, i64 %104
%110 = load double, double* %109
%111 = getelementptr inbounds double, double* %0, i64 %104
%112 = load double, double* %111
%113 = fmul double %112, 2.000000e+00
%114 = fdiv double %110, %113
%115 = tail call double @llvm.fmuladd.f64(double %106, double %108, double %114)
%116 = add nuw nsw i64 %104, 1
%117 = getelementptr inbounds double, double* %1, i64 %116
%118 = load double, double* %117
%119 = fneg double %118
%120 = fmul double %115, %119
%121 = getelementptr inbounds double, double* %7, i64 %116
store double %120, double* %121
%122 = getelementptr inbounds double, double* %1, i64 %104
%123 = load double, double* %122
%124 = getelementptr inbounds double, double* %8, i64 %104
%125 = load double, double* %124
%126 = tail call double @llvm.fmuladd.f64(double %123, double %115, double %125)
store double %126, double* %124
%127 = load double, double* %105
%128 = fneg double %127
%129 = load double, double* %107
%130 = fsub double 1.000000e+00, %129
%131 = getelementptr inbounds double, double* %4, i64 %116
%132 = load double, double* %131
%133 = load double, double* %111
%134 = fmul double %133, 2.000000e+00
%135 = fdiv double %132, %134
%136 = tail call double @llvm.fmuladd.f64(double %128, double %130, double %135)
%137 = load double, double* %117
%138 = getelementptr inbounds double, double* %8, i64 %116
%139 = load double, double* %138
%140 = tail call double @llvm.fmuladd.f64(double %137, double %136, double %139)
store double %140, double* %138
%141 = load double, double* %122
%142 = fneg double %141
%143 = fmul double %136, %142
%144 = getelementptr inbounds double, double* %9, i64 %104
store double %143, double* %144
%145 = icmp eq i64 %116, %97
br i1 %145, label %146, label %103
146: ; preds = %103, %10, %94
ret void
}