#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 exp(double);
double fabs(double);
int /*<<< orphan*/ gamma_wilson_3(double *, double, double *, double **,
                                  double **, double, double);
double pow(double, int);
int /*<<< orphan*/ printf(char *);
void bubpt_column(double TT[], double YY[], double PP, double XX[], double AA[],
                  double BB[], double CC[], double VVIJ[][7], double AAIJ[][7],
                  double R, double T_ref) {
  int i, j, count;
  double sumY, f, fslope, TTT, abs_sumY1;
  double aa[3], bb[3], cc[3], ps[3];
  double xx[3] = {0., 0., 0.};
  double yy[3] = {0., 0., 0.};
  double gg[3] = {0., 0., 0.};
  double vvij[3][3], aaij[3][3];
  double bbb[3] = {0., 0., 0.};
  for (i = 4; i < 7; i++) {
    aa[i - 4] = AA[i];
    bb[i - 4] = BB[i];
    cc[i - 4] = CC[i];
    xx[i - 4] = XX[i];
    yy[i - 4] = YY[i];
    for (j = 4; j < 7; j++) {
      vvij[i - 4][j - 4] = VVIJ[i][j];
      aaij[i - 4][j - 4] = AAIJ[i][j];
    }
  }
  TTT = TT[0];
  count = 0;
  while (1) {
    for (i = 0; i < 3; i++)
      ps[i] = exp(aa[i] + bb[i] / (TTT + cc[i]));
    gamma_wilson_3(gg, TTT, xx, vvij, aaij, R, T_ref);
    sumY = 0.;
    fslope = 0.;
    for (i = 0; i < 3; i++) {
      yy[i] = ps[i] * xx[i] * gg[i] / PP;
      sumY += yy[i];
    }
    abs_sumY1 = fabs(sumY - 1);
    if (abs_sumY1 > 1e-10) {
      f = sumY * PP - PP;
      for (j = 0; j < 3; j++) {
        bbb[j] = bb[j] * xx[j] * ps[j] * gg[j] / (pow((TTT + cc[j]), 2.));
        fslope += -bbb[j];
      }
      TTT = TTT - f / fslope;
    } else
      break;
    if (count++ > 100) {
      printf("bubpt_column.c: Breaking further calculation due to infinite "
             "loop \n");
      break;
    }
  }
  for (i = 4; i < 7; i++)
    YY[i] = yy[i - 4];
  TT[0] = TTT;
}


x86 assembly (clang, -O0):
<bubpt_column>:
push   %rbp
mov    %rsp,%rbp
sub    $0x230,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movsd  %xmm1,-0x40(%rbp)
movsd  %xmm2,-0x48(%rbp)
lea    -0x120(%rbp),%rdi
xor    %esi,%esi
mov    $0x18,%edx
call   1040 <memset@plt>
lea    -0x140(%rbp),%rdi
xor    %esi,%esi
mov    $0x18,%edx
call   1040 <memset@plt>
lea    -0x160(%rbp),%rdi
xor    %esi,%esi
mov    $0x18,%edx
call   1040 <memset@plt>
lea    -0x220(%rbp),%rdi
xor    %esi,%esi
mov    $0x18,%edx
call   1040 <memset@plt>
movl   $0x4,-0x4c(%rbp)
cmpl   $0x7,-0x4c(%rbp)
jge    1329 <bubpt_column+0x1d9>
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
cltq
movsd  %xmm0,-0xa0(%rbp,%rax,8)
ff ff
mov    -0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
cltq
movsd  %xmm0,-0xc0(%rbp,%rax,8)
ff ff
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
cltq
movsd  %xmm0,-0xe0(%rbp,%rax,8)
ff ff
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
cltq
movsd  %xmm0,-0x120(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
cltq
movsd  %xmm0,-0x140(%rbp,%rax,8)
ff ff
movl   $0x4,-0x50(%rbp)
cmpl   $0x7,-0x50(%rbp)
jge    1316 <bubpt_column+0x1c6>
mov    0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
imul   $0x38,%rcx,%rcx
add    %rcx,%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
movslq %eax,%rcx
lea    -0x1b0(%rbp),%rax
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
mov    -0x50(%rbp),%ecx
sub    $0x4,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
imul   $0x38,%rcx,%rcx
add    %rcx,%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
movslq %eax,%rcx
lea    -0x200(%rbp),%rax
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
mov    -0x50(%rbp),%ecx
sub    $0x4,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1284 <bubpt_column+0x134>
jmp    131b <bubpt_column+0x1cb>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11dd <bubpt_column+0x8d>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x54(%rbp)
movl   $0x0,-0x4c(%rbp)
cmpl   $0x3,-0x4c(%rbp)
jge    13a2 <bubpt_column+0x252>
movslq -0x4c(%rbp),%rax
movsd  -0xa0(%rbp,%rax,8),%xmm0
ff ff
movslq -0x4c(%rbp),%rax
movsd  -0xc0(%rbp,%rax,8),%xmm1
ff ff
movsd  -0x78(%rbp),%xmm2
movslq -0x4c(%rbp),%rax
addsd  -0xe0(%rbp,%rax,8),%xmm2
ff ff
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
call   1050 <exp@plt>
movslq -0x4c(%rbp),%rax
movsd  %xmm0,-0x100(%rbp,%rax,8)
ff ff
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1344 <bubpt_column+0x1f4>
lea    -0x160(%rbp),%rdi
movsd  -0x78(%rbp),%xmm0
lea    -0x120(%rbp),%rsi
lea    -0x1b0(%rbp),%rdx
lea    -0x200(%rbp),%rcx
movsd  -0x40(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm2
call   1070 <gamma_wilson_3@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x4c(%rbp)
cmpl   $0x3,-0x4c(%rbp)
jge    1451 <bubpt_column+0x301>
movslq -0x4c(%rbp),%rax
movsd  -0x100(%rbp,%rax,8),%xmm0
ff ff
movslq -0x4c(%rbp),%rax
mulsd  -0x120(%rbp,%rax,8),%xmm0
ff ff
movslq -0x4c(%rbp),%rax
mulsd  -0x160(%rbp,%rax,8),%xmm0
ff ff
divsd  -0x18(%rbp),%xmm0
movslq -0x4c(%rbp),%rax
movsd  %xmm0,-0x140(%rbp,%rax,8)
ff ff
movslq -0x4c(%rbp),%rax
movsd  -0x140(%rbp,%rax,8),%xmm0
ff ff
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    13e9 <bubpt_column+0x299>
movsd  -0x60(%rbp),%xmm0
movsd  0xbaa(%rip),%xmm1
subsd  %xmm1,%xmm0
movaps 0xba7(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  0xb81(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    159a <bubpt_column+0x44a>
movsd  -0x60(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x18(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x50(%rbp)
cmpl   $0x3,-0x50(%rbp)
jge    157d <bubpt_column+0x42d>
movslq -0x50(%rbp),%rax
movsd  -0xc0(%rbp,%rax,8),%xmm0
ff ff
movslq -0x50(%rbp),%rax
mulsd  -0x120(%rbp,%rax,8),%xmm0
ff ff
movslq -0x50(%rbp),%rax
mulsd  -0x100(%rbp,%rax,8),%xmm0
ff ff
movslq -0x50(%rbp),%rax
mulsd  -0x160(%rbp,%rax,8),%xmm0
ff ff
movsd  %xmm0,-0x228(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movslq -0x50(%rbp),%rax
addsd  -0xe0(%rbp,%rax,8),%xmm0
ff ff
mov    $0x2,%edi
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x228(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movslq -0x50(%rbp),%rax
movsd  %xmm0,-0x220(%rbp,%rax,8)
ff ff
movslq -0x50(%rbp),%rax
movsd  -0x220(%rbp,%rax,8),%xmm0
ff ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    14c3 <bubpt_column+0x373>
movsd  -0x78(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
divsd  -0x70(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    159f <bubpt_column+0x44f>
jmp    15c9 <bubpt_column+0x479>
mov    -0x54(%rbp),%eax
mov    %eax,%ecx
add    $0x1,%ecx
mov    %ecx,-0x54(%rbp)
cmp    $0x64,%eax
jle    15c4 <bubpt_column+0x474>
lea    0xa66(%rip),%rdi
call   1030 <printf@plt>
jmp    15c9 <bubpt_column+0x479>
jmp    133d <bubpt_column+0x1ed>
movl   $0x4,-0x4c(%rbp)
cmpl   $0x7,-0x4c(%rbp)
jge    1606 <bubpt_column+0x4b6>
mov    -0x4c(%rbp),%eax
sub    $0x4,%eax
cltq
movsd  -0x140(%rbp,%rax,8),%xmm0
ff ff
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    15d0 <bubpt_column+0x480>
movsd  -0x78(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x230,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bubpt_column>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x248,%rsp
movsd  %xmm2,0xe0(%rsp)
movsd  %xmm1,0xd8(%rsp)
mov    %rsi,0x98(%rsp)
mov    %rdi,%r15
mov    0x288(%rsp),%rax
mov    0x280(%rsp),%rsi
xorps  %xmm1,%xmm1
movaps %xmm1,0x40(%rsp)
movq   $0x0,0x50(%rsp)
movsd  0x20(%rcx),%xmm1
movsd  %xmm1,0xd0(%rsp)
movsd  0x28(%rcx),%xmm1
movsd  %xmm1,0xc8(%rsp)
movsd  0x30(%rcx),%xmm1
movsd  %xmm1,0xc0(%rsp)
movsd  0x20(%r8),%xmm2
movupd 0x28(%r8),%xmm3
movsd  0x20(%r9),%xmm1
movsd  %xmm1,0xb8(%rsp)
movsd  0x28(%r9),%xmm1
movsd  %xmm1,0xb0(%rsp)
movsd  0x30(%r9),%xmm1
movsd  %xmm1,0xa8(%rsp)
mov    0x30(%rdx),%rcx
mov    %rcx,0x70(%rsp)
movups 0x20(%rdx),%xmm1
movaps %xmm1,0x60(%rsp)
mov    0x110(%rsi),%rcx
mov    %rcx,0x210(%rsp)
movups 0x100(%rsi),%xmm1
movaps %xmm1,0x200(%rsp)
mov    0x110(%rax),%rcx
mov    %rcx,0x1c0(%rsp)
movups 0x100(%rax),%xmm1
movaps %xmm1,0x1b0(%rsp)
movups 0x138(%rsi),%xmm1
movups %xmm1,0x218(%rsp)
mov    0x148(%rsi),%rcx
mov    %rcx,0x228(%rsp)
movups 0x138(%rax),%xmm1
movups %xmm1,0x1c8(%rsp)
mov    0x148(%rax),%rcx
mov    %rcx,0x1d8(%rsp)
movups 0x170(%rsi),%xmm1
movaps %xmm1,0x230(%rsp)
mov    0x180(%rsi),%rcx
mov    %rcx,0x240(%rsp)
movups 0x170(%rax),%xmm1
movaps %xmm1,0x1e0(%rsp)
mov    0x180(%rax),%rax
mov    %rax,0x1f0(%rsp)
movapd 0xd40(%rip),%xmm1
movapd %xmm0,%xmm4
xorpd  %xmm1,%xmm4
movapd %xmm4,0x140(%rsp)
movapd %xmm2,0x190(%rsp)
xorpd  %xmm1,%xmm2
movapd %xmm2,0x130(%rsp)
xorpd  %xmm3,%xmm1
movapd %xmm1,0x150(%rsp)
movapd %xmm0,0xf0(%rsp)
unpcklpd %xmm0,%xmm0
movapd %xmm0,0x170(%rsp)
movapd %xmm3,0x180(%rsp)
unpckhpd %xmm3,%xmm3
movapd %xmm3,0x160(%rsp)
mov    $0x66,%r14d
lea    0x40(%rsp),%r12
lea    0x60(%rsp),%r13
lea    0x200(%rsp),%rbx
lea    0x1b0(%rsp),%rbp
movsd  (%rdi),%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm0,0x8(%rsp)
movapd %xmm0,%xmm1
addsd  0xb8(%rsp),%xmm1
movapd 0x190(%rsp),%xmm0
movsd  %xmm1,0x30(%rsp)
divsd  %xmm1,%xmm0
addsd  0xd0(%rsp),%xmm0
call   1040 <exp@plt>
movapd %xmm0,0x10(%rsp)
movsd  0x8(%rsp),%xmm1
addsd  0xb0(%rsp),%xmm1
movapd 0x180(%rsp),%xmm0
movsd  %xmm1,0x80(%rsp)
divsd  %xmm1,%xmm0
addsd  0xc8(%rsp),%xmm0
call   1040 <exp@plt>
movapd %xmm0,0x110(%rsp)
movsd  0x8(%rsp),%xmm1
addsd  0xa8(%rsp),%xmm1
movapd 0x160(%rsp),%xmm0
movsd  %xmm1,0xe8(%rsp)
divsd  %xmm1,%xmm0
addsd  0xc0(%rsp),%xmm0
call   1040 <exp@plt>
movapd %xmm0,0x100(%rsp)
mov    %r12,%rdi
movsd  0x8(%rsp),%xmm0
mov    %r13,%rsi
mov    %rbx,%rdx
mov    %rbp,%rcx
movsd  0xd8(%rsp),%xmm1
movsd  0xe0(%rsp),%xmm2
call   1060 <gamma_wilson_3@plt>
movapd 0x60(%rsp),%xmm3
movapd 0x10(%rsp),%xmm1
unpcklpd 0x110(%rsp),%xmm1
mulpd  %xmm3,%xmm1
movapd 0x40(%rsp),%xmm5
mulpd  %xmm5,%xmm1
divpd  0x170(%rsp),%xmm1
movapd %xmm1,%xmm0
addsd  0xbcd(%rip),%xmm0
movapd %xmm1,%xmm4
unpckhpd %xmm1,%xmm4
addsd  %xmm0,%xmm4
movsd  0x70(%rsp),%xmm2
mulsd  0x100(%rsp),%xmm2
mulsd  0x50(%rsp),%xmm2
divsd  0xf0(%rsp),%xmm2
addsd  %xmm2,%xmm4
movapd %xmm4,%xmm0
addsd  0xb83(%rip),%xmm0
andpd  0xb6b(%rip),%xmm0
ucomisd 0xb7b(%rip),%xmm0
jbe    161f <bubpt_column+0x4df>
movsd  %xmm2,0xa0(%rsp)
movapd %xmm1,0x120(%rsp)
movsd  0x30(%rsp),%xmm0
mov    $0x2,%edi
movapd %xmm3,0x30(%rsp)
movapd %xmm4,0x1a0(%rsp)
movapd %xmm5,0x20(%rsp)
call   1050 <pow@plt>
movapd 0x30(%rsp),%xmm1
mulsd  0x130(%rsp),%xmm1
mulsd  0x10(%rsp),%xmm1
mulsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0xb1f(%rip),%xmm1
movapd %xmm1,0x30(%rsp)
movups 0x68(%rsp),%xmm0
movaps %xmm0,0x20(%rsp)
movups 0x48(%rsp),%xmm0
movaps %xmm0,0x10(%rsp)
movsd  0x80(%rsp),%xmm0
mov    $0x2,%edi
call   1050 <pow@plt>
movapd 0x20(%rsp),%xmm1
movhpd 0x70(%rsp),%xmm1
movaps %xmm0,0x80(%rsp)
mulpd  0x150(%rsp),%xmm1
movapd %xmm1,0x20(%rsp)
movaps 0x10(%rsp),%xmm0
movhps 0x50(%rsp),%xmm0
movaps %xmm0,0x10(%rsp)
movsd  0xe8(%rsp),%xmm0
mov    $0x2,%edi
call   1050 <pow@plt>
movapd 0x110(%rsp),%xmm1
unpcklpd 0x100(%rsp),%xmm1
mulpd  0x20(%rsp),%xmm1
mulpd  0x10(%rsp),%xmm1
movapd 0x80(%rsp),%xmm2
unpcklpd %xmm0,%xmm2
divpd  %xmm2,%xmm1
movapd 0x30(%rsp),%xmm0
addsd  %xmm1,%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm0,%xmm1
movapd 0x1a0(%rsp),%xmm0
mulsd  0xf0(%rsp),%xmm0
addsd  0x140(%rsp),%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
subsd  %xmm1,%xmm0
add    $0xffffffff,%r14d
jne    1350 <bubpt_column+0x210>
movsd  %xmm0,0x8(%rsp)
lea    0xa30(%rip),%rdi
call   1030 <puts@plt>
movapd 0x120(%rsp),%xmm1
movsd  0xa0(%rsp),%xmm2
mov    0x98(%rsp),%rax
movupd %xmm1,0x20(%rax)
movsd  %xmm2,0x30(%rax)
movsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%r15)
add    $0x248,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [68 x i8] c"bubpt_column.c: Breaking further calculation due to infinite loop \0A\00"
define void @bubpt_column(double* %0, double* %1, double %2, double* %3, double* %4, double* %5, double* %6, [7 x double]* %7, [7 x double]* %8, double %9, double %10) {
%12 = alloca double*
%13 = alloca double*
%14 = alloca double
%15 = alloca double*
%16 = alloca double*
%17 = alloca double*
%18 = alloca double*
%19 = alloca [7 x double]*
%20 = alloca [7 x double]*
%21 = alloca double
%22 = alloca double
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca [3 x double]
%32 = alloca [3 x double]
%33 = alloca [3 x double]
%34 = alloca [3 x double]
%35 = alloca [3 x double]
%36 = alloca [3 x double]
%37 = alloca [3 x double]
%38 = alloca [3 x [3 x double]]
%39 = alloca [3 x [3 x double]]
%40 = alloca [3 x double]
store double* %0, double** %12
store double* %1, double** %13
store double %2, double* %14
store double* %3, double** %15
store double* %4, double** %16
store double* %5, double** %17
store double* %6, double** %18
store [7 x double]* %7, [7 x double]** %19
store [7 x double]* %8, [7 x double]** %20
store double %9, double* %21
store double %10, double* %22
%41 = bitcast [3 x double]* %35 to i8*
%42 = bitcast [3 x double]* %36 to i8*
%43 = bitcast [3 x double]* %37 to i8*
%44 = bitcast [3 x double]* %40 to i8*
store i32 4, i32* %23
br label %45
45: ; preds = %134, %11
%46 = load i32, i32* %23
%47 = icmp slt i32 %46, 7
br i1 %47, label %48, label %137
48: ; preds = %45
%49 = load double*, double** %16
%50 = load i32, i32* %23
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load double, double* %52
%54 = load i32, i32* %23
%55 = sub nsw i32 %54, 4
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds [3 x double], [3 x double]* %31, i64 0, i64 %56
store double %53, double* %57
%58 = load double*, double** %17
%59 = load i32, i32* %23
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load double, double* %61
%63 = load i32, i32* %23
%64 = sub nsw i32 %63, 4
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds [3 x double], [3 x double]* %32, i64 0, i64 %65
store double %62, double* %66
%67 = load double*, double** %18
%68 = load i32, i32* %23
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
%71 = load double, double* %70
%72 = load i32, i32* %23
%73 = sub nsw i32 %72, 4
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 0, i64 %74
store double %71, double* %75
%76 = load double*, double** %15
%77 = load i32, i32* %23
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %76, i64 %78
%80 = load double, double* %79
%81 = load i32, i32* %23
%82 = sub nsw i32 %81, 4
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds [3 x double], [3 x double]* %35, i64 0, i64 %83
store double %80, double* %84
%85 = load double*, double** %13
%86 = load i32, i32* %23
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
%89 = load double, double* %88
%90 = load i32, i32* %23
%91 = sub nsw i32 %90, 4
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds [3 x double], [3 x double]* %36, i64 0, i64 %92
store double %89, double* %93
store i32 4, i32* %24
br label %94
94: ; preds = %130, %48
%95 = load i32, i32* %24
%96 = icmp slt i32 %95, 7
br i1 %96, label %97, label %133
97: ; preds = %94
%98 = load [7 x double]*, [7 x double]** %19
%99 = load i32, i32* %23
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds [7 x double], [7 x double]* %98, i64 %100
%102 = load i32, i32* %24
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds [7 x double], [7 x double]* %101, i64 0, i64 %103
%105 = load double, double* %104
%106 = load i32, i32* %23
%107 = sub nsw i32 %106, 4
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %38, i64 0, i64 %108
%110 = load i32, i32* %24
%111 = sub nsw i32 %110, 4
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds [3 x double], [3 x double]* %109, i64 0, i64 %112
store double %105, double* %113
%114 = load [7 x double]*, [7 x double]** %20
%115 = load i32, i32* %23
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds [7 x double], [7 x double]* %114, i64 %116
%118 = load i32, i32* %24
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds [7 x double], [7 x double]* %117, i64 0, i64 %119
%121 = load double, double* %120
%122 = load i32, i32* %23
%123 = sub nsw i32 %122, 4
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 %124
%126 = load i32, i32* %24
%127 = sub nsw i32 %126, 4
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds [3 x double], [3 x double]* %125, i64 0, i64 %128
store double %121, double* %129
br label %130
130: ; preds = %97
%131 = load i32, i32* %24
%132 = add nsw i32 %131, 1
store i32 %132, i32* %24
br label %94
133: ; preds = %94
br label %134
134: ; preds = %133
%135 = load i32, i32* %23
%136 = add nsw i32 %135, 1
store i32 %136, i32* %23
br label %45
137: ; preds = %45
%138 = load double*, double** %12
%139 = getelementptr inbounds double, double* %138, i64 0
%140 = load double, double* %139
store double %140, double* %29
store i32 0, i32* %25
br label %141
141: ; preds = %137, %281
store i32 0, i32* %23
br label %142
142: ; preds = %166, %141
%143 = load i32, i32* %23
%144 = icmp slt i32 %143, 3
br i1 %144, label %145, label %169
145: ; preds = %142
%146 = load i32, i32* %23
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds [3 x double], [3 x double]* %31, i64 0, i64 %147
%149 = load double, double* %148
%150 = load i32, i32* %23
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds [3 x double], [3 x double]* %32, i64 0, i64 %151
%153 = load double, double* %152
%154 = load double, double* %29
%155 = load i32, i32* %23
%156 = sext i32 %155 to i64
%157 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 0, i64 %156
%158 = load double, double* %157
%159 = fadd double %154, %158
%160 = fdiv double %153, %159
%161 = fadd double %149, %160
%162 = call double @exp(double %161)
%163 = load i32, i32* %23
%164 = sext i32 %163 to i64
%165 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 %164
store double %162, double* %165
br label %166
166: ; preds = %145
%167 = load i32, i32* %23
%168 = add nsw i32 %167, 1
store i32 %168, i32* %23
br label %142
169: ; preds = %142
%170 = getelementptr inbounds [3 x double], [3 x double]* %37, i64 0, i64 0
%171 = load double, double* %29
%172 = getelementptr inbounds [3 x double], [3 x double]* %35, i64 0, i64 0
%173 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %38, i64 0, i64 0
%174 = bitcast [3 x double]* %173 to double**
%175 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %39, i64 0, i64 0
%176 = bitcast [3 x double]* %175 to double**
%177 = load double, double* %21
%178 = load double, double* %22
%179 = call i32 @gamma_wilson_3(double* %170, double %171, double* %172, double** %174, double** %176, double %177, double %178)
store double 0.000000e+00, double* %26
store double 0.000000e+00, double* %28
store i32 0, i32* %23
br label %180
180: ; preds = %209, %169
%181 = load i32, i32* %23
%182 = icmp slt i32 %181, 3
br i1 %182, label %183, label %212
183: ; preds = %180
%184 = load i32, i32* %23
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 %185
%187 = load double, double* %186
%188 = load i32, i32* %23
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds [3 x double], [3 x double]* %35, i64 0, i64 %189
%191 = load double, double* %190
%192 = fmul double %187, %191
%193 = load i32, i32* %23
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds [3 x double], [3 x double]* %37, i64 0, i64 %194
%196 = load double, double* %195
%197 = fmul double %192, %196
%198 = load double, double* %14
%199 = fdiv double %197, %198
%200 = load i32, i32* %23
%201 = sext i32 %200 to i64
%202 = getelementptr inbounds [3 x double], [3 x double]* %36, i64 0, i64 %201
store double %199, double* %202
%203 = load i32, i32* %23
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds [3 x double], [3 x double]* %36, i64 0, i64 %204
%206 = load double, double* %205
%207 = load double, double* %26
%208 = fadd double %207, %206
store double %208, double* %26
br label %209
209: ; preds = %183
%210 = load i32, i32* %23
%211 = add nsw i32 %210, 1
store i32 %211, i32* %23
br label %180
212: ; preds = %180
%213 = load double, double* %26
%214 = fsub double %213, 1.000000e+00
%215 = call double @llvm.fabs.f64(double %214)
store double %215, double* %30
%216 = load double, double* %30
%217 = fcmp ogt double %216, 1.000000e-10
br i1 %217, label %218, label %274
218: ; preds = %212
%219 = load double, double* %26
%220 = load double, double* %14
%221 = load double, double* %14
%222 = fneg double %221
%223 = call double @llvm.fmuladd.f64(double %219, double %220, double %222)
store double %223, double* %27
store i32 0, i32* %24
br label %224
224: ; preds = %265, %218
%225 = load i32, i32* %24
%226 = icmp slt i32 %225, 3
br i1 %226, label %227, label %268
227: ; preds = %224
%228 = load i32, i32* %24
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds [3 x double], [3 x double]* %32, i64 0, i64 %229
%231 = load double, double* %230
%232 = load i32, i32* %24
%233 = sext i32 %232 to i64
%234 = getelementptr inbounds [3 x double], [3 x double]* %35, i64 0, i64 %233
%235 = load double, double* %234
%236 = fmul double %231, %235
%237 = load i32, i32* %24
%238 = sext i32 %237 to i64
%239 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 %238
%240 = load double, double* %239
%241 = fmul double %236, %240
%242 = load i32, i32* %24
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds [3 x double], [3 x double]* %37, i64 0, i64 %243
%245 = load double, double* %244
%246 = fmul double %241, %245
%247 = load double, double* %29
%248 = load i32, i32* %24
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 0, i64 %249
%251 = load double, double* %250
%252 = fadd double %247, %251
%253 = call double @pow(double %252, i32 2)
%254 = fdiv double %246, %253
%255 = load i32, i32* %24
%256 = sext i32 %255 to i64
%257 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 %256
store double %254, double* %257
%258 = load i32, i32* %24
%259 = sext i32 %258 to i64
%260 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 %259
%261 = load double, double* %260
%262 = fneg double %261
%263 = load double, double* %28
%264 = fadd double %263, %262
store double %264, double* %28
br label %265
265: ; preds = %227
%266 = load i32, i32* %24
%267 = add nsw i32 %266, 1
store i32 %267, i32* %24
br label %224
268: ; preds = %224
%269 = load double, double* %29
%270 = load double, double* %27
%271 = load double, double* %28
%272 = fdiv double %270, %271
%273 = fsub double %269, %272
store double %273, double* %29
br label %275
274: ; preds = %212
br label %282
275: ; preds = %268
%276 = load i32, i32* %25
%277 = add nsw i32 %276, 1
store i32 %277, i32* %25
%278 = icmp sgt i32 %276, 100
br i1 %278, label %279, label %281
279: ; preds = %275
%280 = call i32 @printf(i8* getelementptr inbounds ([68 x i8], [68 x i8]* @.str, i64 0, i64 0))
br label %282
281: ; preds = %275
br label %141
282: ; preds = %279, %274
store i32 4, i32* %23
br label %283
283: ; preds = %296, %282
%284 = load i32, i32* %23
%285 = icmp slt i32 %284, 7
br i1 %285, label %286, label %299
286: ; preds = %283
%287 = load i32, i32* %23
%288 = sub nsw i32 %287, 4
%289 = sext i32 %288 to i64
%290 = getelementptr inbounds [3 x double], [3 x double]* %36, i64 0, i64 %289
%291 = load double, double* %290
%292 = load double*, double** %13
%293 = load i32, i32* %23
%294 = sext i32 %293 to i64
%295 = getelementptr inbounds double, double* %292, i64 %294
store double %291, double* %295
br label %296
296: ; preds = %286
%297 = load i32, i32* %23
%298 = add nsw i32 %297, 1
store i32 %298, i32* %23
br label %283
299: ; preds = %283
%300 = load double, double* %29
%301 = load double*, double** %12
%302 = getelementptr inbounds double, double* %301, i64 0
store double %300, double* %302
ret void
}

LLVM IR (-O2):
@str = [67 x i8] c"bubpt_column.c: Breaking further calculation due to infinite loop \00"
define void @bubpt_column(double* nocapture %0, double* nocapture %1, double %2, double* nocapture %3, double* nocapture %4, double* nocapture %5, double* nocapture %6, [7 x double]* nocapture %7, [7 x double]* nocapture %8, double %9, double %10) {
%12 = alloca [3 x double]
%13 = bitcast [3 x double]* %12 to i8*
%14 = alloca [3 x double]
%15 = alloca [3 x [3 x double]]
%16 = alloca [3 x [3 x double]]
%17 = bitcast [3 x double]* %14 to i8*
%18 = bitcast [3 x [3 x double]]* %15 to i8*
%19 = bitcast [3 x [3 x double]]* %16 to i8*
%20 = getelementptr inbounds double, double* %4, i64 4
%21 = load double, double* %20
%22 = getelementptr inbounds double, double* %4, i64 5
%23 = load double, double* %22
%24 = getelementptr inbounds double, double* %4, i64 6
%25 = load double, double* %24
%26 = getelementptr inbounds double, double* %5, i64 4
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %5, i64 5
%29 = bitcast double* %28 to <2 x double>*
%30 = load <2 x double>, <2 x double>* %29
%31 = getelementptr inbounds double, double* %6, i64 4
%32 = load double, double* %31
%33 = getelementptr inbounds double, double* %6, i64 5
%34 = load double, double* %33
%35 = getelementptr inbounds double, double* %6, i64 6
%36 = load double, double* %35
%37 = getelementptr double, double* %3, i64 4
%38 = bitcast double* %37 to i8*
%39 = bitcast [3 x [3 x double]]* %16 to i8*
%40 = getelementptr [7 x double], [7 x double]* %8, i64 4, i64 4
%41 = bitcast double* %40 to i8*
%42 = bitcast [3 x [3 x double]]* %15 to i8*
%43 = getelementptr [7 x double], [7 x double]* %7, i64 4, i64 4
%44 = bitcast double* %43 to i8*
%45 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %16, i64 0, i64 1, i64 0
%46 = bitcast double* %45 to i8*
%47 = getelementptr [7 x double], [7 x double]* %8, i64 5, i64 4
%48 = bitcast double* %47 to i8*
%49 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %15, i64 0, i64 1, i64 0
%50 = bitcast double* %49 to i8*
%51 = getelementptr [7 x double], [7 x double]* %7, i64 5, i64 4
%52 = bitcast double* %51 to i8*
%53 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %16, i64 0, i64 2, i64 0
%54 = bitcast double* %53 to i8*
%55 = getelementptr [7 x double], [7 x double]* %8, i64 6, i64 4
%56 = bitcast double* %55 to i8*
%57 = getelementptr inbounds [3 x [3 x double]], [3 x [3 x double]]* %15, i64 0, i64 2, i64 0
%58 = bitcast double* %57 to i8*
%59 = getelementptr [7 x double], [7 x double]* %7, i64 6, i64 4
%60 = bitcast double* %59 to i8*
%61 = load double, double* %0
%62 = getelementptr inbounds [3 x double], [3 x double]* %14, i64 0, i64 0
%63 = getelementptr inbounds [3 x double], [3 x double]* %12, i64 0, i64 0
%64 = bitcast [3 x [3 x double]]* %15 to double**
%65 = bitcast [3 x [3 x double]]* %16 to double**
%66 = fneg double %2
%67 = getelementptr inbounds [3 x double], [3 x double]* %12, i64 0, i64 1
%68 = getelementptr inbounds [3 x double], [3 x double]* %14, i64 0, i64 1
%69 = getelementptr inbounds [3 x double], [3 x double]* %12, i64 0, i64 2
%70 = getelementptr inbounds [3 x double], [3 x double]* %14, i64 0, i64 2
%71 = fneg double %27
%72 = fneg <2 x double> %30
%73 = insertelement <2 x double> poison, double %2, i64 0
%74 = shufflevector <2 x double> %73, <2 x double> poison, <2 x i32> zeroinitializer
%75 = extractelement <2 x double> %30, i64 0
%76 = extractelement <2 x double> %30, i64 1
%77 = bitcast [3 x double]* %12 to <2 x double>*
%78 = bitcast [3 x double]* %14 to <2 x double>*
%79 = bitcast double* %67 to <2 x double>*
%80 = bitcast double* %68 to <2 x double>*
br label %81
81: ; preds = %117, %11
%82 = phi i32 [ 0, %11 ], [ %149, %117 ]
%83 = phi double [ %61, %11 ], [ %148, %117 ]
%84 = fadd double %83, %32
%85 = fdiv double %27, %84
%86 = fadd double %21, %85
%87 = call double @exp(double %86)
%88 = fadd double %83, %34
%89 = fdiv double %75, %88
%90 = fadd double %23, %89
%91 = call double @exp(double %90)
%92 = fadd double %83, %36
%93 = fdiv double %76, %92
%94 = fadd double %25, %93
%95 = call double @exp(double %94)
%96 = call i32 @gamma_wilson_3(double* %62, double %83, double* %63, double** %64, double** %65, double %9, double %10)
%97 = load <2 x double>, <2 x double>* %77
%98 = insertelement <2 x double> poison, double %87, i64 0
%99 = insertelement <2 x double> %98, double %91, i64 1
%100 = fmul <2 x double> %99, %97
%101 = load <2 x double>, <2 x double>* %78
%102 = fmul <2 x double> %100, %101
%103 = fdiv <2 x double> %102, %74
%104 = extractelement <2 x double> %103, i64 0
%105 = fadd double %104, 0.000000e+00
%106 = extractelement <2 x double> %103, i64 1
%107 = fadd double %105, %106
%108 = load double, double* %69
%109 = fmul double %95, %108
%110 = load double, double* %70
%111 = fmul double %109, %110
%112 = fdiv double %111, %2
%113 = fadd double %107, %112
%114 = fadd double %113, -1.000000e+00
%115 = call double @llvm.fabs.f64(double %114)
%116 = fcmp ogt double %115, 1.000000e-10
br i1 %116, label %117, label %153
117: ; preds = %81
%118 = call double @pow(double %84, i32 2)
%119 = extractelement <2 x double> %97, i64 0
%120 = fmul double %119, %71
%121 = fmul double %120, %87
%122 = extractelement <2 x double> %101, i64 0
%123 = fmul double %121, %122
%124 = fdiv double %123, %118
%125 = fadd double %124, 0.000000e+00
%126 = load <2 x double>, <2 x double>* %79
%127 = load <2 x double>, <2 x double>* %80
%128 = call double @pow(double %88, i32 2)
%129 = load double, double* %69
%130 = load double, double* %70
%131 = call double @pow(double %92, i32 2)
%132 = insertelement <2 x double> %126, double %129, i64 1
%133 = fmul <2 x double> %132, %72
%134 = insertelement <2 x double> poison, double %91, i64 0
%135 = insertelement <2 x double> %134, double %95, i64 1
%136 = fmul <2 x double> %133, %135
%137 = insertelement <2 x double> %127, double %130, i64 1
%138 = fmul <2 x double> %136, %137
%139 = insertelement <2 x double> poison, double %128, i64 0
%140 = insertelement <2 x double> %139, double %131, i64 1
%141 = fdiv <2 x double> %138, %140
%142 = extractelement <2 x double> %141, i64 0
%143 = fadd double %125, %142
%144 = extractelement <2 x double> %141, i64 1
%145 = fadd double %143, %144
%146 = call double @llvm.fmuladd.f64(double %113, double %2, double %66)
%147 = fdiv double %146, %145
%148 = fsub double %83, %147
%149 = add nuw nsw i32 %82, 1
%150 = icmp eq i32 %149, 102
br i1 %150, label %151, label %81
151: ; preds = %117
%152 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([67 x i8], [67 x i8]* @str, i64 0, i64 0))
br label %153
153: ; preds = %81, %151
%154 = phi double [ %148, %151 ], [ %83, %81 ]
%155 = getelementptr inbounds double, double* %1, i64 6
%156 = getelementptr inbounds double, double* %1, i64 4
%157 = bitcast double* %156 to <2 x double>*
store <2 x double> %103, <2 x double>* %157
store double %112, double* %155
store double %154, double* %0
ret void
}