void lu_sustitucion(double **a, double *o, int n, double *b, double *x) {
  int i, j;
  double sum;
  for (i = 1; i < n; i++) {
    sum = b[i];
    for (j = 0; j < i - 1; j++) {
      sum = sum - a[i][j] * b[j];
    }
    b[i] = sum;
  }
  x[n - 1] = b[n - 1] / a[n - 1][n - 1];
  for (i = n - 2; i >= 0; i--) {
    sum = 0;
    for (j = i + 1; j < n; j++) {
      sum = sum + a[i][j] * x[j];
    }
    x[i] = (b[i] - sum) / a[i][i];
  }
}


x86 assembly (clang, -O0):
<lu_sustitucion>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11cd <lu_sustitucion+0xcd>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11ad <lu_sustitucion+0xad>
movsd  -0x38(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1143 <lu_sustitucion+0x43>
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    111e <lu_sustitucion+0x1e>
mov    -0x20(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     12c4 <lu_sustitucion+0x1c4>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1282 <lu_sustitucion+0x182>
movsd  -0x38(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1234 <lu_sustitucion+0x134>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1219 <lu_sustitucion+0x119>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lu_sustitucion>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    -0x1(%rdx),%r9d
cmp    $0x1,%edx
jle    120d <lu_sustitucion+0x10d>
mov    %r9d,%r10d
mov    $0x1,%r14d
xor    %r11d,%r11d
jmp    1147 <lu_sustitucion+0x47>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm0,(%rcx,%r14,8)
add    $0x1,%r14
add    $0x1,%r11
cmp    %r10,%r11
je     1229 <lu_sustitucion+0x129>
movsd  (%rcx,%r14,8),%xmm0
cmp    $0x2,%r14
jb     1130 <lu_sustitucion+0x30>
lea    -0x1(%r11),%rax
mov    (%rdi,%r14,8),%rsi
cmp    $0x3,%rax
jae    1170 <lu_sustitucion+0x70>
xor    %eax,%eax
jmp    11c7 <lu_sustitucion+0xc7>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r11,%rbx
and    $0xfffffffffffffffc,%rbx
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%rsi,%rax,8),%xmm1
movsd  0x8(%rsi,%rax,8),%xmm2
mulsd  (%rcx,%rax,8),%xmm1
mulsd  0x8(%rcx,%rax,8),%xmm2
subsd  %xmm1,%xmm0
movsd  0x10(%rsi,%rax,8),%xmm1
mulsd  0x10(%rcx,%rax,8),%xmm1
subsd  %xmm2,%xmm0
movsd  0x18(%rsi,%rax,8),%xmm2
mulsd  0x18(%rcx,%rax,8),%xmm2
subsd  %xmm1,%xmm0
subsd  %xmm2,%xmm0
add    $0x4,%rax
cmp    %rax,%rbx
jne    1180 <lu_sustitucion+0x80>
test   $0x3,%r11b
je     1130 <lu_sustitucion+0x30>
mov    %r11d,%ebx
and    $0x3,%ebx
lea    (%rcx,%rax,8),%r15
lea    (%rsi,%rax,8),%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rsi,8),%xmm1
mulsd  (%r15,%rsi,8),%xmm1
subsd  %xmm1,%xmm0
add    $0x1,%rsi
cmp    %rsi,%rbx
jne    11f0 <lu_sustitucion+0xf0>
jmp    1130 <lu_sustitucion+0x30>
movslq %r9d,%rax
movsd  (%rcx,%rax,8),%xmm0
mov    (%rdi,%rax,8),%rcx
divsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%r8,%rax,8)
jmp    1370 <lu_sustitucion+0x270>
movslq %r9d,%rax
movsd  (%rcx,%rax,8),%xmm0
mov    (%rdi,%rax,8),%rsi
divsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%r8,%rax,8)
cmp    $0x2,%edx
jl     1370 <lu_sustitucion+0x270>
mov    %edx,%r15d
add    $0xfffffffe,%edx
mov    $0x1,%r9d
sub    %r15,%r9
lea    -0x1(%r15),%r10
lea    (%r8,%r15,8),%r12
add    $0xfffffffffffffff8,%r12
xor    %r11d,%r11d
mov    $0x1,%r13b
mov    %r15,%r14
jmp    12a8 <lu_sustitucion+0x1a8>
nop
movsd  (%rcx,%rdx,8),%xmm1
subsd  %xmm0,%xmm1
divsd  (%rax,%rdx,8),%xmm1
movsd  %xmm1,(%r8,%rdx,8)
lea    -0x1(%rdx),%rax
add    $0x1,%r11d
add    $0x1,%r9
add    $0xfffffffffffffff8,%r12
add    $0xffffffffffffffff,%r10
add    $0x1,%r13b
test   %rdx,%rdx
mov    %rax,%rdx
jle    1370 <lu_sustitucion+0x270>
movzbl %r13b,%r13d
and    $0x3,%r13d
add    $0xffffffffffffffff,%r14
lea    0x1(%rdx),%rsi
mov    (%rdi,%rdx,8),%rax
xorpd  %xmm0,%xmm0
cmp    %r15,%rsi
jge    1270 <lu_sustitucion+0x170>
lea    0x1(%r11),%esi
test   $0x3,%sil
je     1306 <lu_sustitucion+0x206>
lea    (%rax,%r10,8),%rbx
xorpd  %xmm0,%xmm0
xor    %esi,%esi
nopl   0x0(%rax)
movsd  (%rbx,%rsi,8),%xmm1
mulsd  (%r12,%rsi,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rsi
cmp    %esi,%r13d
jne    12e0 <lu_sustitucion+0x1e0>
sub    %r9,%rsi
cmp    $0x3,%r11d
jae    1320 <lu_sustitucion+0x220>
jmp    1270 <lu_sustitucion+0x170>
xorpd  %xmm0,%xmm0
mov    %r14,%rsi
cmp    $0x3,%r11d
jb     1270 <lu_sustitucion+0x170>
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rsi,8),%xmm1
movsd  0x8(%rax,%rsi,8),%xmm2
mulsd  (%r8,%rsi,8),%xmm1
mulsd  0x8(%r8,%rsi,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x10(%rax,%rsi,8),%xmm3
mulsd  0x10(%r8,%rsi,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x18(%rax,%rsi,8),%xmm0
mulsd  0x18(%r8,%rsi,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rsi
cmp    %esi,%r15d
jne    1320 <lu_sustitucion+0x220>
jmp    1270 <lu_sustitucion+0x170>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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