void tridag(a, b, c, r, u, n) double *a, *b, *c, *r, *u;
int n;
{
  int j;
  double bet, gam[3844 + 1];
  if (b[1] == 0.0)
    problem("Error 1 in TRIDAG.");
  u[1] = r[1] / (bet = b[1]);
  for (j = 2; j <= n; j++) {
    gam[j] = c[j - 1] / bet;
    bet = b[j] - a[j] * gam[j];
    if (bet == 0.0)
      problem("Error 2 in TRIDAG.");
    u[j] = (r[j] - a[j] * u[j - 1]) / bet;
  }
  for (j = (n - 1); j >= 1; j--)
    u[j] -= gam[j + 1] * u[j + 1];
}


x86 assembly (clang, -O0):
<tridag>:
push   %rbp
mov    %rsp,%rbp
sub    $0x7860,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    115d <tridag+0x4d>
jp     115d <tridag+0x4d>
lea    0xeaa(%rip),%rdi
mov    $0x0,%al
call   1030 <problem@plt>
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
movsd  %xmm1,-0x38(%rbp)
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movl   $0x2,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1294 <tridag+0x184>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x38(%rbp),%xmm0
movslq -0x30(%rbp),%rax
movsd  %xmm0,-0x7860(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x30(%rbp),%rax
movsd  -0x7860(%rbp,%rax,8),%xmm2
ff ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1229 <tridag+0x119>
jp     1229 <tridag+0x119>
lea    0xdf1(%rip),%rdi
mov    $0x0,%al
call   1030 <problem@plt>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1188 <tridag+0x78>
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x1,-0x30(%rbp)
jl     1309 <tridag+0x1f9>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
cltq
movsd  -0x7860(%rbp,%rax,8),%xmm0
ff ff
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    129d <tridag+0x18d>
add    $0x7860,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tridag>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x7848,%rsp
mov    %r9d,%r14d
mov    %r8,%rbx
mov    %rcx,%r15
mov    %rdx,0x10(%rsp)
mov    %rsi,%r13
mov    %rdi,%rbp
movsd  0x8(%rsi),%xmm3
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm3
jne    115a <tridag+0x4a>
jp     115a <tridag+0x4a>
lea    0xecb(%rip),%rdi
xor    %eax,%eax
call   1030 <problem@plt>
movsd  0x8(%r13),%xmm3
movsd  0x8(%r15),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,0x8(%rbx)
cmp    $0x2,%r14d
jl     1287 <tridag+0x177>
mov    %r14,0x8(%rsp)
mov    %r14d,%eax
mov    %rax,(%rsp)
lea    -0x8(,%rax,8),%r14
ff
xor    %r12d,%r12d
jmp    11bb <tridag+0xab>
nopl   0x0(%rax)
xorpd  0xe68(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%r15,%r12,1),%xmm1
divsd  %xmm3,%xmm1
movsd  %xmm1,0x10(%rbx,%r12,1)
add    $0x8,%r12
movapd %xmm1,%xmm0
cmp    %r12,%r14
je     1221 <tridag+0x111>
mov    0x10(%rsp),%rax
movsd  0x8(%rax,%r12,1),%xmm2
divsd  %xmm3,%xmm2
movsd  %xmm2,0x30(%rsp,%r12,1)
movsd  0x10(%r13,%r12,1),%xmm3
movsd  0x10(%rbp,%r12,1),%xmm1
mulsd  %xmm1,%xmm2
subsd  %xmm2,%xmm3
ucomisd 0xe20(%rip),%xmm3
jne    1190 <tridag+0x80>
jp     1190 <tridag+0x80>
lea    0xe30(%rip),%rdi
xor    %eax,%eax
movsd  %xmm3,0x18(%rsp)
call   1030 <problem@plt>
movsd  0x18(%rsp),%xmm3
movsd  0x8(%rbx,%r12,1),%xmm0
movsd  0x10(%rbp,%r12,1),%xmm1
jmp    1190 <tridag+0x80>
mov    0x8(%rsp),%rsi
cmp    $0x2,%esi
jl     1287 <tridag+0x177>
xor    %ecx,%ecx
mov    (%rsp),%r11
mov    %r11,%rax
sub    $0x2,%rax
mov    $0x0,%r9d
cmovae %rax,%r9
add    $0x1,%r9
cmp    $0x4,%r9
jae    1299 <tridag+0x189>
mov    %r11,%rcx
add    $0x1,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    -0x2(%rcx),%eax
movsd  0x18(%rsp,%rcx,8),%xmm0
movsd  (%rbx,%rax,8),%xmm1
mulsd  -0x8(%rbx,%rcx,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
add    $0xffffffffffffffff,%rcx
cmp    $0x2,%rcx
ja     1260 <tridag+0x150>
add    $0x7848,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r11d,%eax
sub    $0x2,%eax
cmovb  %ecx,%eax
mov    %r11,%rdx
sub    $0x2,%rdx
cmovae %rdx,%rcx
add    $0xffffffff,%esi
cmp    %eax,%esi
jb     124c <tridag+0x13c>
mov    %rcx,%rax
shr    $0x20,%rax
jne    124c <tridag+0x13c>
lea    (%rbx,%rsi,8),%rax
shl    $0x3,%rcx
mov    %rax,%rdx
sub    %rcx,%rdx
cmp    %rax,%rdx
ja     124c <tridag+0x13c>
xor    %eax,%eax
test   %al,%al
jne    124c <tridag+0x13c>
lea    (%rbx,%r11,8),%rax
mov    %rax,%rdx
sub    %rcx,%rdx
cmp    %rax,%rdx
ja     124c <tridag+0x13c>
xor    %eax,%eax
test   %al,%al
jne    124c <tridag+0x13c>
cmp    $0x2,%r11
mov    $0x2,%eax
cmovb  %r11,%rax
lea    (%rax,%rsi,1),%rcx
sub    %r11,%rcx
lea    (%rbx,%rcx,8),%rcx
lea    (%rbx,%r11,8),%rdx
add    $0x8,%rdx
cmp    %rdx,%rcx
jae    1335 <tridag+0x225>
lea    (%rbx,%rsi,8),%rcx
add    $0x8,%rcx
lea    (%rbx,%rax,8),%rax
cmp    %rcx,%rax
jb     124c <tridag+0x13c>
mov    %r9,%r8
and    $0xfffffffffffffffc,%r8
mov    %r11,%rcx
sub    %r8,%rcx
mov    %r8,%r10
neg    %r10
mov    $0xffffffff,%edi
add    %r11,%rdi
lea    (%rsp,%r11,8),%rbp
add    $0x18,%rbp
lea    (%rbx,%r11,8),%rdx
add    $0xfffffffffffffff8,%rdx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%rdi,%rax,1),%esi
movupd -0x10(%rbp,%rax,8),%xmm0
movupd 0x0(%rbp,%rax,8),%xmm1
movupd -0x10(%rdx,%rax,8),%xmm2
mulpd  %xmm0,%xmm2
movupd (%rdx,%rax,8),%xmm0
mulpd  %xmm1,%xmm0
movupd -0x18(%rbx,%rsi,8),%xmm1
subpd  %xmm2,%xmm1
movupd -0x8(%rbx,%rsi,8),%xmm2
subpd  %xmm0,%xmm2
movupd %xmm2,-0x8(%rbx,%rsi,8)
movupd %xmm1,-0x18(%rbx,%rsi,8)
add    $0xfffffffffffffffc,%rax
cmp    %rax,%r10
jne    1370 <tridag+0x260>
cmp    %r8,%r9
jne    124f <tridag+0x13f>
jmp    1287 <tridag+0x177>

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

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