#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 */
int /*<<< orphan*/ dpofa(double *, int *, int *, int *);
int formt(int *m, double *wt, double *sy, double *ss, int *col, double *theta,
          int *info) {
  int wt_dim1, wt_offset, sy_dim1, sy_offset, ss_dim1, ss_offset, i__1, i__2,
      i__3;
  int i__, j, k, k1;
  double ddum;
  ss_dim1 = *m;
  ss_offset = 1 + ss_dim1;
  ss -= ss_offset;
  sy_dim1 = *m;
  sy_offset = 1 + sy_dim1;
  sy -= sy_offset;
  wt_dim1 = *m;
  wt_offset = 1 + wt_dim1;
  wt -= wt_offset;
  i__1 = *col;
  for (j = 1; j <= i__1; ++j) {
    wt[j * wt_dim1 + 1] = *theta * ss[j * ss_dim1 + 1];
  }
  i__1 = *col;
  for (i__ = 2; i__ <= i__1; ++i__) {
    i__2 = *col;
    for (j = i__; j <= i__2; ++j) {
      k1 = ((i__) <= (j) ? (i__) : (j)) - 1;
      ddum = 0.;
      i__3 = k1;
      for (k = 1; k <= i__3; ++k) {
        ddum +=
            sy[i__ + k * sy_dim1] * sy[j + k * sy_dim1] / sy[k + k * sy_dim1];
      }
      wt[i__ + j * wt_dim1] = ddum + *theta * ss[i__ + j * ss_dim1];
    }
  }
  dpofa(&wt[wt_offset], m, col, info);
  if (*info != 0) {
    *info = -3;
  }
  return 0;
}


x86 assembly (clang, -O0):
<formt>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%rax
movslq %ecx,%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x1,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jg     1219 <formt+0x109>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x5c(%rbp),%ecx
imul   -0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    11cb <formt+0xbb>
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x2,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jg     135d <formt+0x24d>
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x50(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x50(%rbp),%eax
jg     134a <formt+0x23a>
mov    -0x58(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jg     1267 <formt+0x157>
mov    -0x58(%rbp),%eax
mov    %eax,-0x74(%rbp)
jmp    126d <formt+0x15d>
mov    -0x5c(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x74(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x64(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x54(%rbp)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x54(%rbp),%eax
jg     12f7 <formt+0x1e7>
mov    -0x18(%rbp),%rax
mov    -0x58(%rbp),%ecx
mov    -0x60(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x5c(%rbp),%ecx
mov    -0x60(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x60(%rbp),%ecx
mov    -0x60(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    128b <formt+0x17b>
movsd  -0x70(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x58(%rbp),%ecx
mov    -0x5c(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%ecx
mov    -0x5c(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1244 <formt+0x134>
jmp    134f <formt+0x23f>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1229 <formt+0x119>
mov    -0x10(%rbp),%rdi
movslq -0x38(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x8(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    0x10(%rbp),%rcx
call   1030 <dpofa@plt>
mov    0x10(%rbp),%rax
cmpl   $0x0,(%rax)
je     1394 <formt+0x284>
mov    0x10(%rbp),%rax
movl   $0xfffffffd,(%rax)
xor    %eax,%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<formt>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,(%rsp)
mov    %rdi,0x8(%rsp)
mov    %r8,0x10(%rsp)
mov    (%r8),%r15d
test   %r15d,%r15d
jle    1210 <formt+0x100>
mov    0x8(%rsp),%rax
movslq (%rax),%r14
mov    %r14,%r8
not    %r8
mov    $0x1,%r10d
cmp    $0x1,%r15d
je     11e1 <formt+0xd1>
lea    (%rcx,%r14,8),%rsi
add    $0x8,%rsi
mov    %r14,%r12
shl    $0x4,%r12
mov    (%rsp),%rax
lea    (%rax,%r14,8),%rbp
add    $0x8,%rbp
lea    (%rax,%r12,1),%rdi
add    $0x8,%rdi
lea    (%rcx,%r12,1),%rax
add    $0x8,%rax
mov    %r15d,%r13d
and    $0xfffffffe,%r13d
neg    %r13
mov    $0x1,%r10d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%r9),%xmm0
mulsd  (%rsi,%r8,8),%xmm0
movsd  %xmm0,0x0(%rbp,%r8,8)
movsd  (%r9),%xmm0
mulsd  (%rax,%r8,8),%xmm0
movsd  %xmm0,(%rdi,%r8,8)
add    %r12,%rsi
add    %r12,%rbp
add    %r12,%rdi
add    %r12,%rax
lea    (%r10,%r13,1),%rbx
add    $0x2,%rbx
add    $0x2,%r10
cmp    $0x1,%rbx
jne    11a0 <formt+0x90>
lea    (%rcx,%r8,8),%r12
mov    (%rsp),%rax
lea    (%rax,%r8,8),%r13
test   $0x1,%r15b
je     120a <formt+0xfa>
movsd  (%r9),%xmm0
imul   %r14,%r10
mulsd  0x8(%r12,%r10,8),%xmm0
movsd  %xmm0,0x8(%r13,%r10,8)
cmp    $0x2,%r15d
jge    1247 <formt+0x137>
mov    (%rsp),%rdi
mov    0x8(%rsp),%rsi
mov    0x10(%rsp),%rdx
mov    0x50(%rsp),%rbx
mov    %rbx,%rcx
call   1030 <dpofa@plt>
cmpl   $0x0,(%rbx)
je     1236 <formt+0x126>
movl   $0xfffffffd,(%rbx)
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    (%rdx,%r8,8),%rax
add    $0x1,%r15
add    %r14,%r8
lea    (%rdx,%r8,8),%rcx
add    $0x10,%rcx
mov    $0x2,%r8d
jmp    127d <formt+0x16d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r8
add    $0x8,%rcx
cmp    %r15,%r8
je     1210 <formt+0x100>
mov    %rcx,%rdx
mov    %r8,%r10
jmp    12bd <formt+0x1ad>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r9),%xmm1
mov    %r10,%rsi
imul   %r14,%rsi
add    %r8,%rsi
mulsd  (%r12,%rsi,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rsi,8)
add    $0x1,%r10
add    $0x8,%rdx
cmp    %r15,%r10
je     1270 <formt+0x160>
cmp    %r10,%r8
mov    %r8,%rsi
cmova  %r10,%rsi
xorpd  %xmm0,%xmm0
cmp    $0x2,%esi
jl     1290 <formt+0x180>
mov    %esi,%ebp
xorpd  %xmm0,%xmm0
mov    $0x1,%esi
xor    %edi,%edi
mov    %r14d,%ebx
movslq %ebx,%rbx
movsd  (%rcx,%rdi,8),%xmm1
mulsd  (%rdx,%rdi,8),%xmm1
lea    (%rsi,%rbx,1),%r11
divsd  (%rax,%r11,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rsi
add    %r14d,%ebx
add    %r14,%rdi
cmp    %rbp,%rsi
jb     12e0 <formt+0x1d0>
jmp    1290 <formt+0x180>

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

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