void disVelEvaluate(const int nPoints, const int nSpace, const double self_a,
                    const double *x, const double *u, double *m, double *dm,
                    double *f, double *df, double *a, double *da, double *phi,
                    double *dphi) {
  int k, I, J;
  const int nSpace2 = nSpace * nSpace;
  for (k = 0; k < nPoints; k++) {
    m[k] = u[k];
    phi[k] = 0.0;
    dm[k] = dphi[k] = 1.0;
    f[k * nSpace] = u[k];
    f[k * nSpace + 1] = 0.0;
    df[k * nSpace] = 1.0;
    df[k * nSpace + 1] = 0.0;
    if (x[k * 3 + 1] > 0.5) {
      f[k * nSpace] *= 0.25;
      df[k * nSpace] *= 0.25;
    }
    for (I = 0; I < nSpace; I++) {
      a[k * nSpace2 + I * nSpace + I] = self_a;
      for (J = 0; J < nSpace; J++) {
        da[k * nSpace2 + I * nSpace + J] = 0.0;
      }
    }
  }
}


x86 assembly (clang, -O0):
<disVelEvaluate>:
push   %rbp
mov    %rsp,%rbp
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    131a <disVelEvaluate+0x21a>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x38(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xe7a(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xe65(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xe11(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xdd0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    127a <disVelEvaluate+0x17a>
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xdc0(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xda0(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1307 <disVelEvaluate+0x207>
movsd  -0x10(%rbp),%xmm0
mov    0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x40(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x8(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12f4 <disVelEvaluate+0x1f4>
mov    0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x40(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x8(%rbp),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12b8 <disVelEvaluate+0x1b8>
jmp    12f9 <disVelEvaluate+0x1f9>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1281 <disVelEvaluate+0x181>
jmp    130c <disVelEvaluate+0x20c>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1148 <disVelEvaluate+0x48>
pop    %rbp
ret

x86 assembly (clang, -O2):
<disVelEvaluate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9,0x60(%rsp)
mov    %rdx,0x58(%rsp)
test   %edi,%edi
jle    1413 <disVelEvaluate+0x303>
mov    %rcx,%r9
mov    %esi,%r13d
mov    0xc0(%rsp),%rdx
mov    0xb8(%rsp),%r12
mov    0xb0(%rsp),%rbp
mov    %esi,%r14d
imul   %r13d,%r14d
mov    %esi,%eax
lea    0x0(,%rax,8),%r11
movslq %esi,%rcx
mov    %edi,%edi
mov    %rdi,0x48(%rsp)
mov    %rcx,0x50(%rsp)
add    %ecx,%ecx
mov    %ecx,0x1c(%rsp)
mov    %rax,0x40(%rsp)
and    $0xfffffffe,%eax
neg    %rax
mov    %rax,0x68(%rsp)
xor    %edi,%edi
mov    %esi,%ebx
xor    %r10d,%r10d
movsd  %xmm0,0x10(%rsp)
mov    %r14d,0x4(%rsp)
mov    %r11,0x8(%rsp)
mov    %r8,0x30(%rsp)
mov    %r9,0x28(%rsp)
mov    %r13,0x70(%rsp)
jmp    11c5 <disVelEvaluate+0xb5>
nopl   0x0(%rax)
add    $0x1,%r10
add    %r14d,%edi
add    %r14d,%ebx
cmp    0x48(%rsp),%r10
je     1413 <disVelEvaluate+0x303>
movsd  (%r9,%r10,8),%xmm1
movsd  %xmm1,(%r8,%r10,8)
mov    0xd0(%rsp),%rax
movq   $0x0,(%rax,%r10,8)
mov    0xd8(%rsp),%rax
movabs $0x3ff0000000000000,%rcx
mov    %rcx,(%rax,%r10,8)
mov    0x60(%rsp),%rax
mov    %rcx,(%rax,%r10,8)
movsd  (%r9,%r10,8),%xmm1
mov    %r10,%rax
imul   0x50(%rsp),%rax
movsd  %xmm1,0x0(%rbp,%rax,8)
movq   $0x0,0x8(%rbp,%rax,8)
movsd  0xdda(%rip),%xmm1
movups %xmm1,(%r12,%rax,8)
lea    (%r10,%r10,2),%rcx
mov    0x58(%rsp),%rsi
movsd  0x8(%rsi,%rcx,8),%xmm1
ucomisd 0xdc6(%rip),%xmm1
jbe    126c <disVelEvaluate+0x15c>
movsd  0x0(%rbp,%rax,8),%xmm1
movsd  0xdbe(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  %xmm1,0x0(%rbp,%rax,8)
movsd  (%r12,%rax,8),%xmm1
mulsd  %xmm2,%xmm1
movsd  %xmm1,(%r12,%rax,8)
test   %r13d,%r13d
jle    11b0 <disVelEvaluate+0xa0>
cmp    $0x1,%r13d
mov    %rdi,0x20(%rsp)
jne    12a0 <disVelEvaluate+0x190>
xor    %r15d,%r15d
testb  $0x1,0x40(%rsp)
je     11b0 <disVelEvaluate+0xa0>
jmp    138f <disVelEvaluate+0x27f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r10,0x38(%rsp)
mov    $0x1,%r12d
xor    %r15d,%r15d
mov    %ebx,0x18(%rsp)
mov    %edi,%r14d
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movslq %ebx,%rbx
lea    (%r12,%rbx,1),%r13
movslq %r14d,%r14
lea    (%r14,%rbp,1),%rax
movsd  %xmm0,(%rdx,%rax,8)
mov    0xc8(%rsp),%rax
lea    (%rax,%r14,8),%rdi
xor    %esi,%esi
mov    %r11,%rdx
call   1030 <memset@plt>
movsd  0x10(%rsp),%xmm0
mov    0xc0(%rsp),%rax
movsd  %xmm0,(%rax,%r13,8)
mov    0x70(%rsp),%r13
lea    0x0(%r13,%r14,1),%eax
cltq
mov    0xc8(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
mov    0x8(%rsp),%rdx
call   1030 <memset@plt>
mov    0x1c(%rsp),%eax
mov    0x8(%rsp),%r11
mov    0xc0(%rsp),%rdx
movsd  0x10(%rsp),%xmm0
add    $0x2,%rbp
add    %eax,%r14d
add    $0x2,%r12
add    %eax,%ebx
add    $0xfffffffffffffffe,%r15
cmp    %r15,0x68(%rsp)
jne    12c0 <disVelEvaluate+0x1b0>
neg    %r15
mov    0x30(%rsp),%r8
mov    0x28(%rsp),%r9
mov    0xb8(%rsp),%r12
mov    0xb0(%rsp),%rbp
mov    0x4(%rsp),%r14d
mov    0x20(%rsp),%rdi
mov    0x18(%rsp),%ebx
mov    0x38(%rsp),%r10
testb  $0x1,0x40(%rsp)
je     11b0 <disVelEvaluate+0xa0>
mov    %r14d,%eax
imul   %r10d,%eax
mov    %r10d,%ecx
imul   %r13d,%ecx
add    %r15d,%ecx
imul   %r13d,%ecx
movslq %ecx,%rcx
add    %r15,%rcx
movsd  %xmm0,(%rdx,%rcx,8)
imul   %r13d,%r15d
add    %eax,%r15d
movslq %r15d,%rax
mov    0xc8(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
mov    %r11,%rdx
mov    %r10,%rbp
mov    %r9,%r15
mov    %r8,%r12
call   1030 <memset@plt>
mov    %r12,%r8
mov    %r15,%r9
mov    %rbp,%r10
mov    0x20(%rsp),%rdi
mov    0x8(%rsp),%r11
mov    0x4(%rsp),%r14d
mov    0xb0(%rsp),%rbp
mov    0xb8(%rsp),%r12
mov    0xc0(%rsp),%rdx
movsd  0x10(%rsp),%xmm0
jmp    11b0 <disVelEvaluate+0xa0>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @disVelEvaluate(i32 %0, i32 %1, double %2, double* nocapture %3, double* nocapture %4, double* nocapture writeonly %5, double* nocapture writeonly %6, double* nocapture %7, double* nocapture %8, double* nocapture writeonly %9, double* nocapture writeonly %10, double* nocapture writeonly %11, double* nocapture writeonly %12) {
%14 = icmp sgt i32 %0, 0
br i1 %14, label %15, label %102
15: ; preds = %13
%16 = icmp sgt i32 %1, 0
%17 = mul i32 %1, %1
%18 = zext i32 %1 to i64
%19 = shl nuw nsw i64 %18, 3
%20 = sext i32 %1 to i64
%21 = zext i32 %0 to i64
%22 = and i64 %18, 1
%23 = icmp eq i32 %1, 1
%24 = and i64 %18, 4294967294
%25 = icmp eq i64 %22, 0
br label %26
26: ; preds = %15, %99
%27 = phi i64 [ 0, %15 ], [ %100, %99 ]
%28 = trunc i64 %27 to i32
%29 = mul i32 %17, %28
%30 = getelementptr inbounds double, double* %4, i64 %27
%31 = load double, double* %30
%32 = getelementptr inbounds double, double* %5, i64 %27
store double %31, double* %32
%33 = getelementptr inbounds double, double* %11, i64 %27
store double 0.000000e+00, double* %33
%34 = getelementptr inbounds double, double* %12, i64 %27
store double 1.000000e+00, double* %34
%35 = getelementptr inbounds double, double* %6, i64 %27
store double 1.000000e+00, double* %35
%36 = load double, double* %30
%37 = mul nsw i64 %27, %20
%38 = getelementptr inbounds double, double* %7, i64 %37
store double %36, double* %38
%39 = add nsw i64 %37, 1
%40 = getelementptr inbounds double, double* %7, i64 %39
store double 0.000000e+00, double* %40
%41 = getelementptr inbounds double, double* %8, i64 %37
%42 = bitcast double* %41 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 0.000000e+00>, <2 x double>* %42
%43 = mul nuw nsw i64 %27, 3
%44 = add nuw nsw i64 %43, 1
%45 = getelementptr inbounds double, double* %3, i64 %44
%46 = load double, double* %45
%47 = fcmp ogt double %46, 5.000000e-01
br i1 %47, label %48, label %53
48: ; preds = %26
%49 = load double, double* %38
%50 = fmul double %49, 2.500000e-01
store double %50, double* %38
%51 = load double, double* %41
%52 = fmul double %51, 2.500000e-01
store double %52, double* %41
br label %53
53: ; preds = %48, %26
br i1 %16, label %54, label %99
54: ; preds = %53
%55 = mul i32 %28, %1
br i1 %23, label %85, label %56
56: ; preds = %54, %56
%57 = phi i64 [ %82, %56 ], [ 0, %54 ]
%58 = phi i64 [ %83, %56 ], [ 0, %54 ]
%59 = trunc i64 %57 to i32
%60 = add i32 %55, %59
%61 = mul i32 %60, %1
%62 = sext i32 %61 to i64
%63 = add nsw i64 %57, %62
%64 = getelementptr inbounds double, double* %9, i64 %63
store double %2, double* %64
%65 = mul i32 %59, %1
%66 = add i32 %29, %65
%67 = sext i32 %66 to i64
%68 = getelementptr double, double* %10, i64 %67
%69 = bitcast double* %68 to i8*
%70 = or i64 %57, 1
%71 = trunc i64 %70 to i32
%72 = add i32 %55, %71
%73 = mul i32 %72, %1
%74 = sext i32 %73 to i64
%75 = add nsw i64 %70, %74
%76 = getelementptr inbounds double, double* %9, i64 %75
store double %2, double* %76
%77 = mul i32 %71, %1
%78 = add i32 %29, %77
%79 = sext i32 %78 to i64
%80 = getelementptr double, double* %10, i64 %79
%81 = bitcast double* %80 to i8*
%82 = add nuw nsw i64 %57, 2
%83 = add i64 %58, 2
%84 = icmp eq i64 %83, %24
br i1 %84, label %85, label %56
85: ; preds = %56, %54
%86 = phi i64 [ 0, %54 ], [ %82, %56 ]
br i1 %25, label %99, label %87
87: ; preds = %85
%88 = trunc i64 %86 to i32
%89 = add i32 %55, %88
%90 = mul i32 %89, %1
%91 = sext i32 %90 to i64
%92 = add nsw i64 %86, %91
%93 = getelementptr inbounds double, double* %9, i64 %92
store double %2, double* %93
%94 = mul i32 %88, %1
%95 = add i32 %29, %94
%96 = sext i32 %95 to i64
%97 = getelementptr double, double* %10, i64 %96
%98 = bitcast double* %97 to i8*
br label %99
99: ; preds = %87, %85, %53
%100 = add nuw nsw i64 %27, 1
%101 = icmp eq i64 %100, %21
br i1 %101, label %102, label %26
102: ; preds = %99, %13
ret void
}