#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 *dvector(int, int);
int /*<<< orphan*/ free_dvector(double *, int, int);
void solveWithSvd(double u[][100], double w[], double v[][100], int m, int n,
                  double b[], double x[]) {
  int jj, j, i;
  double s, *tmp;
  tmp = dvector(1, n);
  for (j = 1; j <= n; j++) {
    s = 0.0;
    if (w[j]) {
      for (i = 1; i <= m; i++)
        s += u[i][j] * b[i];
      s /= w[j];
    }
    tmp[j] = s;
  }
  for (j = 1; j <= n; j++) {
    s = 0.0;
    for (jj = 1; jj <= n; jj++)
      s += v[j][jj] * tmp[jj];
    x[j] = s;
  }
  free_dvector(tmp, 1, n);
}


x86 assembly (clang, -O0):
<solveWithSvd>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x20(%rbp),%esi
mov    $0x1,%edi
call   1030 <dvector@plt>
mov    %rax,-0x48(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     122a <solveWithSvd+0x10a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1194 <solveWithSvd+0x74>
jp     1194 <solveWithSvd+0x74>
jmp    120a <solveWithSvd+0xea>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     11ef <solveWithSvd+0xcf>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x320,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    119b <solveWithSvd+0x7b>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x40(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    115b <solveWithSvd+0x3b>
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     12c0 <solveWithSvd+0x1a0>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     12a0 <solveWithSvd+0x180>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0x320,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    124c <solveWithSvd+0x12c>
movsd  -0x40(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1231 <solveWithSvd+0x111>
mov    -0x48(%rbp),%rdi
mov    -0x20(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_dvector@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<solveWithSvd>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,%rbx
mov    %r8d,%ebp
mov    %ecx,%r12d
mov    %rdx,%r15
mov    %rsi,%r13
mov    %rdi,%r14
mov    0x50(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    $0x1,%edi
mov    %r8d,%esi
call   1030 <dvector@plt>
mov    %rbp,%rdx
test   %edx,%edx
jle    139c <solveWithSvd+0x27c>
lea    0x1(%r12),%r8d
mov    %rdx,0x8(%rsp)
lea    0x1(%rdx),%r9d
lea    -0x1(%r8),%r10
add    $0xfffffffffffffffe,%r8
mov    %r10d,%edi
and    $0x3,%edi
and    $0xfffffffffffffffc,%r10
lea    0xc88(%r14),%r11
add    $0x8,%r14
mov    $0x1,%ecx
jmp    11be <solveWithSvd+0x9e>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rcx,8)
add    $0x1,%rcx
add    $0x8,%r11
add    $0x8,%r14
cmp    %r9,%rcx
je     12a3 <solveWithSvd+0x183>
movsd  0x0(%r13,%rcx,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11d1 <solveWithSvd+0xb1>
jnp    11a4 <solveWithSvd+0x84>
xorpd  %xmm1,%xmm1
test   %r12d,%r12d
jle    11a0 <solveWithSvd+0x80>
cmp    $0x3,%r8
jae    11ec <solveWithSvd+0xcc>
mov    $0x1,%ebp
test   %rdi,%rdi
jne    1261 <solveWithSvd+0x141>
jmp    11a0 <solveWithSvd+0x80>
mov    %r11,%rdx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x960(%rdx),%xmm2
ff
mulsd  0x8(%rbx,%rbp,8),%xmm2
addsd  %xmm1,%xmm2
movsd  -0x640(%rdx),%xmm1
ff
mulsd  0x10(%rbx,%rbp,8),%xmm1
addsd  %xmm2,%xmm1
movsd  -0x320(%rdx),%xmm2
ff
mulsd  0x18(%rbx,%rbp,8),%xmm2
addsd  %xmm1,%xmm2
movsd  (%rdx),%xmm1
mulsd  0x20(%rbx,%rbp,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x4,%rbp
add    $0xc80,%rdx
cmp    %rbp,%r10
jne    1200 <solveWithSvd+0xe0>
add    $0x1,%rbp
test   %rdi,%rdi
je     11a0 <solveWithSvd+0x80>
lea    (%rbx,%rbp,8),%rdx
imul   $0x320,%rbp,%rbp
add    %r14,%rbp
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp),%xmm2
mulsd  (%rdx,%rsi,8),%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rsi
add    $0x320,%rbp
cmp    %rsi,%rdi
jne    1280 <solveWithSvd+0x160>
jmp    11a0 <solveWithSvd+0x80>
mov    0x8(%rsp),%rdx
test   %edx,%edx
jle    139c <solveWithSvd+0x27c>
lea    -0x1(%r9),%r11
lea    -0x2(%r9),%r8
mov    %r11d,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%r11
lea    0x340(%r15),%rdi
add    $0x320,%r15
mov    $0x1,%r10d
jmp    1306 <solveWithSvd+0x1e6>
nopl   0x0(%rax,%rax,1)
mov    0x10(%rsp),%rcx
movsd  %xmm0,(%rcx,%r10,8)
add    $0x1,%r10
add    $0x320,%rdi
add    $0x320,%r15
cmp    %r9,%r10
je     139c <solveWithSvd+0x27c>
xorpd  %xmm0,%xmm0
mov    $0x1,%ebx
cmp    $0x3,%r8
jb     136c <solveWithSvd+0x24c>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movsd  -0x18(%rdi,%rbx,8),%xmm1
movsd  -0x10(%rdi,%rbx,8),%xmm2
mulsd  0x8(%rax,%rbx,8),%xmm1
mulsd  0x10(%rax,%rbx,8),%xmm2
addsd  %xmm0,%xmm1
movsd  -0x8(%rdi,%rbx,8),%xmm3
mulsd  0x18(%rax,%rbx,8),%xmm3
addsd  %xmm1,%xmm2
movsd  (%rdi,%rbx,8),%xmm0
mulsd  0x20(%rax,%rbx,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rbx
cmp    %rbx,%r11
jne    1320 <solveWithSvd+0x200>
add    $0x1,%rbx
test   %rsi,%rsi
je     12e0 <solveWithSvd+0x1c0>
lea    (%rax,%rbx,8),%rcx
lea    (%r15,%rbx,8),%rbx
xor    %ebp,%ebp
nop
movsd  (%rbx,%rbp,8),%xmm1
mulsd  (%rcx,%rbp,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbp
cmp    %rbp,%rsi
jne    1380 <solveWithSvd+0x260>
jmp    12e0 <solveWithSvd+0x1c0>
mov    %rax,%rdi
mov    $0x1,%esi
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_dvector@plt>

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

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