void mrqcof(x, y, sig, ndata, a, ma, lista, mfit, alpha, beta, chisq,
            funcs) float x[],
    y[], sig[], a[], **alpha, beta[], *chisq;
int ndata, ma, lista[], mfit;
void (*funcs)();
{
  int k, j, i;
  float ymod, wt, sig2i, dy, *dyda, *vector();
  void free_vector();
  dyda = vector(1, ma);
  for (j = 1; j <= mfit; j++) {
    for (k = 1; k <= j; k++)
      alpha[j][k] = 0.0;
    beta[j] = 0.0;
  }
  *chisq = 0.0;
  for (i = 1; i <= ndata; i++) {
    (*funcs)(x[i], a, &ymod, dyda, ma);
    sig2i = 1.0 / (sig[i] * sig[i]);
    dy = y[i] - ymod;
    for (j = 1; j <= mfit; j++) {
      wt = dyda[lista[j]] * sig2i;
      for (k = 1; k <= j; k++)
        alpha[j][k] += wt * dyda[lista[k]];
      beta[j] += dy * wt;
    }
    (*chisq) += dy * dy * sig2i;
  }
  for (j = 2; j <= mfit; j++)
    for (k = 1; k <= j - 1; k++)
      alpha[k][j] = alpha[j][k];
  free_vector(dyda, 1, ma);
}


x86 assembly (clang, -O0):
<mrqcof>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x2c(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1040 <vector@plt>
mov    %rax,-0x50(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     11d3 <mrqcof+0xb3>
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     11b5 <mrqcof+0x95>
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1183 <mrqcof+0x63>
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1170 <mrqcof+0x50>
mov    0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1356 <mrqcof+0x236>
mov    0x38(%rbp),%r8
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x28(%rbp),%rdi
mov    -0x50(%rbp),%rdx
mov    -0x2c(%rbp),%ecx
lea    -0x3c(%rbp),%rsi
mov    $0x1,%al
call   *%r8
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdc0(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     1325 <mrqcof+0x205>
mov    -0x50(%rbp),%rax
mov    0x10(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     12f3 <mrqcof+0x1d3>
movss  -0x40(%rbp),%xmm0
mov    -0x50(%rbp),%rax
mov    0x10(%rbp),%rcx
movslq -0x30(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    129d <mrqcof+0x17d>
movss  -0x48(%rbp),%xmm0
movss  -0x40(%rbp),%xmm2
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    126b <mrqcof+0x14b>
movss  -0x48(%rbp),%xmm0
mulss  -0x48(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
mov    0x30(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11e5 <mrqcof+0xc5>
movl   $0x2,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     13cc <mrqcof+0x2ac>
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     13b9 <mrqcof+0x299>
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1370 <mrqcof+0x250>
jmp    13be <mrqcof+0x29e>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    135d <mrqcof+0x23d>
mov    -0x50(%rbp),%rdi
mov    -0x2c(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1030 <free_vector@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mrqcof>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r8,0x40(%rsp)
mov    %rcx,0x8(%rsp)
mov    %rdx,0x38(%rsp)
mov    %rsi,0x30(%rsp)
mov    %rdi,0x28(%rsp)
mov    0x98(%rsp),%r14
mov    0x90(%rsp),%r12
mov    0x88(%rsp),%ebp
mov    $0x1,%edi
mov    %r9d,0x4(%rsp)
mov    %r9d,%esi
xor    %eax,%eax
call   1050 <vector@plt>
mov    %rax,%rbx
mov    %rbp,0x10(%rsp)
test   %ebp,%ebp
jle    1259 <mrqcof+0x129>
mov    0x88(%rsp),%eax
mov    %eax,%r13d
cmp    $0x1,%eax
jne    11bb <mrqcof+0x8b>
mov    $0x4,%r15d
mov    $0x1,%ebp
test   $0x1,%r13b
jne    123f <mrqcof+0x10f>
jmp    1259 <mrqcof+0x129>
mov    %r13,0x20(%rsp)
and    $0xfffffffe,%r13d
neg    %r13
mov    $0x1,%ebp
mov    $0x8,%r15d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    -0x4(%r15),%rdx
mov    (%r12,%rbp,8),%rdi
add    $0x4,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
movl   $0x0,(%r14,%rbp,4)
mov    0x8(%r12,%rbp,8),%rdi
add    $0x4,%rdi
xor    %esi,%esi
mov    %r15,%rdx
call   1040 <memset@plt>
movl   $0x0,0x4(%r14,%rbp,4)
add    $0x8,%r15
lea    0x2(,%rbp,1),%rax
add    %r13,%rax
add    $0x2,%rbp
cmp    $0x1,%rax
jne    11e0 <mrqcof+0xb0>
add    $0xfffffffffffffffc,%r15
mov    0x20(%rsp),%r13
test   $0x1,%r13b
je     1259 <mrqcof+0x129>
mov    (%r12,%rbp,8),%rdi
add    $0x4,%rdi
xor    %esi,%esi
mov    %r15,%rdx
call   1040 <memset@plt>
movl   $0x0,(%r14,%rbp,4)
mov    0xa0(%rsp),%rax
movl   $0x0,(%rax)
cmpl   $0x0,0x8(%rsp)
jle    13f8 <mrqcof+0x2c8>
mov    0x80(%rsp),%r15
mov    0x88(%rsp),%eax
lea    0x1(%rax),%r13d
mov    0x8(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x8(%rsp)
mov    $0x1,%r8d
jmp    12c7 <mrqcof+0x197>
nopw   0x0(%rax,%rax,1)
mulss  %xmm0,%xmm0
mulss  %xmm1,%xmm0
mov    0xa0(%rsp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
add    $0x1,%r8
cmp    0x8(%rsp),%r8
je     13f8 <mrqcof+0x2c8>
mov    0x28(%rsp),%rax
movss  (%rax,%r8,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x40(%rsp),%rdi
lea    0x1c(%rsp),%rsi
mov    %rbx,%rdx
mov    0x4(%rsp),%ecx
mov    $0x1,%al
mov    %r8,%rbp
call   *0xa8(%rsp)
mov    %rbp,%r8
mov    0x38(%rsp),%rax
movss  (%rax,%rbp,4),%xmm0
mulss  %xmm0,%xmm0
movss  0xcf4(%rip),%xmm1
divss  %xmm0,%xmm1
mov    0x30(%rsp),%rax
movss  (%rax,%rbp,4),%xmm0
subss  0x1c(%rsp),%xmm0
mov    0x88(%rsp),%eax
test   %eax,%eax
jle    12a0 <mrqcof+0x170>
mov    $0x1,%eax
xor    %ecx,%ecx
jmp    135d <mrqcof+0x22d>
nopl   0x0(%rax,%rax,1)
mulss  %xmm0,%xmm2
addss  (%r14,%rax,4),%xmm2
movss  %xmm2,(%r14,%rax,4)
add    $0x1,%rax
cmp    %r13,%rax
je     12a0 <mrqcof+0x170>
mov    %rcx,%rdi
add    $0x1,%rcx
movslq (%r15,%rax,4),%rdx
movss  (%rbx,%rdx,4),%xmm2
mulss  %xmm1,%xmm2
mov    (%r12,%rax,8),%rdx
mov    $0x1,%esi
test   %rdi,%rdi
je     13d3 <mrqcof+0x2a3>
mov    %rax,%rdi
and    $0xfffffffffffffffe,%rdi
neg    %rdi
mov    $0x1,%esi
xchg   %ax,%ax
movslq (%r15,%rsi,4),%rbp
movss  (%rbx,%rbp,4),%xmm3
mulss  %xmm2,%xmm3
addss  (%rdx,%rsi,4),%xmm3
movss  %xmm3,(%rdx,%rsi,4)
movslq 0x4(%r15,%rsi,4),%rbp
movss  (%rbx,%rbp,4),%xmm3
mulss  %xmm2,%xmm3
addss  0x4(%rdx,%rsi,4),%xmm3
movss  %xmm3,0x4(%rdx,%rsi,4)
lea    (%rdi,%rsi,1),%rbp
add    $0x2,%rbp
add    $0x2,%rsi
cmp    $0x1,%rbp
jne    1390 <mrqcof+0x260>
test   $0x1,%cl
je     1340 <mrqcof+0x210>
movslq (%r15,%rsi,4),%rdi
movss  (%rbx,%rdi,4),%xmm3
mulss  %xmm2,%xmm3
addss  (%rdx,%rsi,4),%xmm3
movss  %xmm3,(%rdx,%rsi,4)
jmp    1340 <mrqcof+0x210>
mov    0x88(%rsp),%eax
cmp    $0x2,%eax
jge    1426 <mrqcof+0x2f6>
mov    %rbx,%rdi
mov    $0x1,%esi
mov    0x4(%rsp),%edx
xor    %eax,%eax
call   1030 <free_vector@plt>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x10(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x10(%rsp)
mov    $0x2,%eax
mov    $0x1,%r9d
xor    %r8d,%r8d
jmp    145f <mrqcof+0x32f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rax
add    $0x1,%r9
cmp    0x10(%rsp),%rax
je     1404 <mrqcof+0x2d4>
mov    %r8,%rcx
add    $0x1,%r8
mov    (%r12,%rax,8),%rsi
mov    $0x1,%edi
cmp    $0x3,%rcx
jb     14cd <mrqcof+0x39d>
mov    %r9,%rcx
and    $0xfffffffffffffffc,%rcx
xor    %edi,%edi
xchg   %ax,%ax
movss  0x4(%rsi,%rdi,4),%xmm0
mov    0x8(%r12,%rdi,8),%rdx
movss  %xmm0,(%rdx,%rax,4)
movss  0x8(%rsi,%rdi,4),%xmm0
mov    0x10(%r12,%rdi,8),%rdx
movss  %xmm0,(%rdx,%rax,4)
movss  0xc(%rsi,%rdi,4),%xmm0
mov    0x18(%r12,%rdi,8),%rdx
movss  %xmm0,(%rdx,%rax,4)
movss  0x10(%rsi,%rdi,4),%xmm0
mov    0x20(%r12,%rdi,8),%rdx
movss  %xmm0,(%rdx,%rax,4)
add    $0x4,%rdi
cmp    %rdi,%rcx
jne    1480 <mrqcof+0x350>
add    $0x1,%rdi
test   $0x3,%r8b
je     1450 <mrqcof+0x320>
mov    %r9d,%ecx
and    $0x3,%ecx
lea    (%r12,%rdi,8),%rdx
lea    (%rsi,%rdi,4),%rsi
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
movss  (%rsi,%rdi,4),%xmm0
mov    (%rdx,%rdi,8),%rbp
movss  %xmm0,0x0(%rbp,%rax,4)
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    14f0 <mrqcof+0x3c0>
jmp    1450 <mrqcof+0x320>

LLVM IR (-O0):
define void @mrqcof(float* %0, float* %1, float* %2, i32 %3, float* %4, i32 %5, i32* %6, i32 %7, float** %8, float* %9, float* %10, void (...)* %11) {
%13 = alloca float*
%14 = alloca float*
%15 = alloca float*
%16 = alloca i32
%17 = alloca float*
%18 = alloca i32
%19 = alloca i32*
%20 = alloca i32
%21 = alloca float**
%22 = alloca float*
%23 = alloca float*
%24 = alloca void (...)*
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca float
%29 = alloca float
%30 = alloca float
%31 = alloca float
%32 = alloca float*
store float* %0, float** %13
store float* %1, float** %14
store float* %2, float** %15
store i32 %3, i32* %16
store float* %4, float** %17
store i32 %5, i32* %18
store i32* %6, i32** %19
store i32 %7, i32* %20
store float** %8, float*** %21
store float* %9, float** %22
store float* %10, float** %23
store void (...)* %11, void (...)** %24
%33 = load i32, i32* %18
%34 = call float* (i32, i32, ...) bitcast (float* (...)* @vector to float* (i32, i32, ...)*)(i32 1, i32 %33)
store float* %34, float** %32
store i32 1, i32* %26
br label %35
35: ; preds = %61, %12
%36 = load i32, i32* %26
%37 = load i32, i32* %20
%38 = icmp sle i32 %36, %37
br i1 %38, label %39, label %64
39: ; preds = %35
store i32 1, i32* %25
br label %40
40: ; preds = %53, %39
%41 = load i32, i32* %25
%42 = load i32, i32* %26
%43 = icmp sle i32 %41, %42
br i1 %43, label %44, label %56
44: ; preds = %40
%45 = load float**, float*** %21
%46 = load i32, i32* %26
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds float*, float** %45, i64 %47
%49 = load float*, float** %48
%50 = load i32, i32* %25
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %49, i64 %51
store float 0.000000e+00, float* %52
br label %53
53: ; preds = %44
%54 = load i32, i32* %25
%55 = add nsw i32 %54, 1
store i32 %55, i32* %25
br label %40
56: ; preds = %40
%57 = load float*, float** %22
%58 = load i32, i32* %26
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
store float 0.000000e+00, float* %60
br label %61
61: ; preds = %56
%62 = load i32, i32* %26
%63 = add nsw i32 %62, 1
store i32 %63, i32* %26
br label %35
64: ; preds = %35
%65 = load float*, float** %23
store float 0.000000e+00, float* %65
store i32 1, i32* %27
br label %66
66: ; preds = %167, %64
%67 = load i32, i32* %27
%68 = load i32, i32* %16
%69 = icmp sle i32 %67, %68
br i1 %69, label %70, label %170
70: ; preds = %66
%71 = load void (...)*, void (...)** %24
%72 = load float*, float** %13
%73 = load i32, i32* %27
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds float, float* %72, i64 %74
%76 = load float, float* %75
%77 = fpext float %76 to double
%78 = load float*, float** %17
%79 = load float*, float** %32
%80 = load i32, i32* %18
%81 = bitcast void (...)* %71 to void (double, float*, float*, float*, i32, ...)*
%82 = load float*, float** %15
%83 = load i32, i32* %27
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
%86 = load float, float* %85
%87 = load float*, float** %15
%88 = load i32, i32* %27
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %87, i64 %89
%91 = load float, float* %90
%92 = fmul float %86, %91
%93 = fpext float %92 to double
%94 = fdiv double 1.000000e+00, %93
%95 = fptrunc double %94 to float
store float %95, float* %30
%96 = load float*, float** %14
%97 = load i32, i32* %27
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds float, float* %96, i64 %98
%100 = load float, float* %99
%101 = load float, float* %28
%102 = fsub float %100, %101
store float %102, float* %31
store i32 1, i32* %26
br label %103
103: ; preds = %156, %70
%104 = load i32, i32* %26
%105 = load i32, i32* %20
%106 = icmp sle i32 %104, %105
br i1 %106, label %107, label %159
107: ; preds = %103
%108 = load float*, float** %32
%109 = load i32*, i32** %19
%110 = load i32, i32* %26
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds i32, i32* %109, i64 %111
%113 = load i32, i32* %112
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds float, float* %108, i64 %114
%116 = load float, float* %115
%117 = load float, float* %30
%118 = fmul float %116, %117
store float %118, float* %29
store i32 1, i32* %25
br label %119
119: ; preds = %144, %107
%120 = load i32, i32* %25
%121 = load i32, i32* %26
%122 = icmp sle i32 %120, %121
br i1 %122, label %123, label %147
123: ; preds = %119
%124 = load float, float* %29
%125 = load float*, float** %32
%126 = load i32*, i32** %19
%127 = load i32, i32* %25
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds i32, i32* %126, i64 %128
%130 = load i32, i32* %129
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds float, float* %125, i64 %131
%133 = load float, float* %132
%134 = load float**, float*** %21
%135 = load i32, i32* %26
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds float*, float** %134, i64 %136
%138 = load float*, float** %137
%139 = load i32, i32* %25
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds float, float* %138, i64 %140
%142 = load float, float* %141
%143 = call float @llvm.fmuladd.f32(float %124, float %133, float %142)
store float %143, float* %141
br label %144
144: ; preds = %123
%145 = load i32, i32* %25
%146 = add nsw i32 %145, 1
store i32 %146, i32* %25
br label %119
147: ; preds = %119
%148 = load float, float* %31
%149 = load float, float* %29
%150 = load float*, float** %22
%151 = load i32, i32* %26
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds float, float* %150, i64 %152
%154 = load float, float* %153
%155 = call float @llvm.fmuladd.f32(float %148, float %149, float %154)
store float %155, float* %153
br label %156
156: ; preds = %147
%157 = load i32, i32* %26
%158 = add nsw i32 %157, 1
store i32 %158, i32* %26
br label %103
159: ; preds = %103
%160 = load float, float* %31
%161 = load float, float* %31
%162 = fmul float %160, %161
%163 = load float, float* %30
%164 = load float*, float** %23
%165 = load float, float* %164
%166 = call float @llvm.fmuladd.f32(float %162, float %163, float %165)
store float %166, float* %164
br label %167
167: ; preds = %159
%168 = load i32, i32* %27
%169 = add nsw i32 %168, 1
store i32 %169, i32* %27
br label %66
170: ; preds = %66
store i32 2, i32* %26
br label %171
171: ; preds = %203, %170
%172 = load i32, i32* %26
%173 = load i32, i32* %20
%174 = icmp sle i32 %172, %173
br i1 %174, label %175, label %206
175: ; preds = %171
store i32 1, i32* %25
br label %176
176: ; preds = %199, %175
%177 = load i32, i32* %25
%178 = load i32, i32* %26
%179 = sub nsw i32 %178, 1
%180 = icmp sle i32 %177, %179
br i1 %180, label %181, label %202
181: ; preds = %176
%182 = load float**, float*** %21
%183 = load i32, i32* %26
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds float*, float** %182, i64 %184
%186 = load float*, float** %185
%187 = load i32, i32* %25
%188 = sext i32 %187 to i64
%189 = getelementptr inbounds float, float* %186, i64 %188
%190 = load float, float* %189
%191 = load float**, float*** %21
%192 = load i32, i32* %25
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds float*, float** %191, i64 %193
%195 = load float*, float** %194
%196 = load i32, i32* %26
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds float, float* %195, i64 %197
store float %190, float* %198
br label %199
199: ; preds = %181
%200 = load i32, i32* %25
%201 = add nsw i32 %200, 1
store i32 %201, i32* %25
br label %176
202: ; preds = %176
br label %203
203: ; preds = %202
%204 = load i32, i32* %26
%205 = add nsw i32 %204, 1
store i32 %205, i32* %26
br label %171
206: ; preds = %171
%207 = load float*, float** %32
%208 = load i32, i32* %18
ret void
}

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