#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 get_length(double *, int);
scalar_t__ isnormal(double);
double scalar_product(double *, double *, int);
int /*<<< orphan*/ swap_maxlen_to_head(double **, int, int);
void gram_schmidt_orthogonalization(double **vector, int m, int n) {
  if (m <= 1)
    return;
  for (int i = 0; i < m; i++) {
    swap_maxlen_to_head(vector + i, m - i, n);
    for (int j = i + 1; j < m; j++) {
      double alpha = scalar_product(vector[i], vector[j], n);
      double length = get_length(vector[i], n);
      if (isnormal(length)) {
        alpha /= length * length;
      }
      for (int k = 0; k < n; k++) {
        vector[j][k] -= alpha * vector[i][k];
      }
    }
  }
}


x86 assembly (clang, -O0):
<gram_schmidt_orthogonalization>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
cmpl   $0x1,-0xc(%rbp)
jg     1161 <gram_schmidt_orthogonalization+0x21>
jmp    12aa <gram_schmidt_orthogonalization+0x16a>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12aa <gram_schmidt_orthogonalization+0x16a>
mov    -0x8(%rbp),%rdi
movslq -0x14(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0xc(%rbp),%esi
sub    -0x14(%rbp),%esi
mov    -0x10(%rbp),%edx
call   1060 <swap_maxlen_to_head@plt>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1297 <gram_schmidt_orthogonalization+0x157>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
mov    -0x10(%rbp),%edx
call   1050 <scalar_product@plt>
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x10(%rbp),%esi
call   1030 <get_length@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1040 <isnormal@plt>
cmp    $0x0,%rax
je     1210 <gram_schmidt_orthogonalization+0xd0>
movsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1284 <gram_schmidt_orthogonalization+0x144>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1217 <gram_schmidt_orthogonalization+0xd7>
jmp    1289 <gram_schmidt_orthogonalization+0x149>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    119a <gram_schmidt_orthogonalization+0x5a>
jmp    129c <gram_schmidt_orthogonalization+0x15c>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1168 <gram_schmidt_orthogonalization+0x28>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gram_schmidt_orthogonalization>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %esi,0xc(%rsp)
cmp    $0x2,%esi
jl     13e8 <gram_schmidt_orthogonalization+0x2a8>
mov    %edx,%ebx
mov    %rdi,%r15
mov    0xc(%rsp),%r12d
mov    %edx,%ebp
mov    %ebp,%eax
and    $0xfffffffc,%eax
mov    %rax,0x38(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,0x30(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,0x28(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x18(%rsp)
mov    %rbp,%rax
neg    %rax
mov    %rax,0x40(%rsp)
mov    $0x1,%eax
mov    %rax,0x10(%rsp)
xor    %r14d,%r14d
jmp    11c4 <gram_schmidt_orthogonalization+0x84>
addq   $0x1,0x10(%rsp)
mov    0x20(%rsp),%r14
cmp    %r12,%r14
je     13e8 <gram_schmidt_orthogonalization+0x2a8>
lea    (%r15,%r14,8),%r13
mov    0xc(%rsp),%esi
sub    %r14d,%esi
mov    %r13,%rdi
mov    %ebx,%edx
call   1060 <swap_maxlen_to_head@plt>
add    $0x1,%r14
mov    %r14,0x20(%rsp)
cmp    %r12,%r14
jae    11b0 <gram_schmidt_orthogonalization+0x70>
mov    0x10(%rsp),%r14
jmp    11f9 <gram_schmidt_orthogonalization+0xb9>
xchg   %ax,%ax
add    $0x1,%r14
cmp    %r12,%r14
je     11b0 <gram_schmidt_orthogonalization+0x70>
mov    0x0(%r13),%rdi
mov    (%r15,%r14,8),%rsi
mov    %ebx,%edx
call   1050 <scalar_product@plt>
movapd %xmm0,0x50(%rsp)
mov    0x0(%r13),%rdi
mov    %ebx,%esi
call   1030 <get_length@plt>
movsd  %xmm0,0x48(%rsp)
call   1040 <isnormal@plt>
test   %rax,%rax
je     1250 <gram_schmidt_orthogonalization+0x110>
movsd  0x48(%rsp),%xmm0
mulsd  %xmm0,%xmm0
movapd 0x50(%rsp),%xmm6
divsd  %xmm0,%xmm6
test   %ebx,%ebx
jg     125a <gram_schmidt_orthogonalization+0x11a>
jmp    11f0 <gram_schmidt_orthogonalization+0xb0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd 0x50(%rsp),%xmm6
test   %ebx,%ebx
jle    11f0 <gram_schmidt_orthogonalization+0xb0>
mov    0x0(%r13),%rax
mov    (%r15,%r14,8),%rcx
xorpd  0xd96(%rip),%xmm6
cmp    $0x4,%ebx
jb     1281 <gram_schmidt_orthogonalization+0x141>
lea    (%rax,%rbp,8),%rdx
cmp    %rdx,%rcx
jae    12f7 <gram_schmidt_orthogonalization+0x1b7>
lea    (%rcx,%rbp,8),%rdx
cmp    %rdx,%rax
jae    12f7 <gram_schmidt_orthogonalization+0x1b7>
xor    %esi,%esi
mov    %rsi,%rdx
test   $0x1,%bpl
je     12a6 <gram_schmidt_orthogonalization+0x166>
movsd  (%rax,%rsi,8),%xmm0
mulsd  %xmm6,%xmm0
addsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rcx,%rsi,8)
mov    %rsi,%rdx
or     $0x1,%rdx
not    %rsi
cmp    0x40(%rsp),%rsi
je     11f0 <gram_schmidt_orthogonalization+0xb0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rax,%rdx,8),%xmm0
mulsd  %xmm6,%xmm0
addsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
movsd  0x8(%rax,%rdx,8),%xmm0
mulsd  %xmm6,%xmm0
addsd  0x8(%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rbp
jne    12c0 <gram_schmidt_orthogonalization+0x180>
jmp    11f0 <gram_schmidt_orthogonalization+0xb0>
movapd %xmm6,%xmm0
unpcklpd %xmm6,%xmm0
cmpq   $0x0,0x30(%rsp)
je     13dd <gram_schmidt_orthogonalization+0x29d>
mov    0x18(%rsp),%rsi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rax,%rdx,8),%xmm1
movupd 0x10(%rax,%rdx,8),%xmm2
mulpd  %xmm0,%xmm1
movupd (%rcx,%rdx,8),%xmm3
addpd  %xmm1,%xmm3
movupd 0x10(%rcx,%rdx,8),%xmm1
movupd 0x20(%rcx,%rdx,8),%xmm4
movupd 0x30(%rcx,%rdx,8),%xmm5
mulpd  %xmm0,%xmm2
addpd  %xmm1,%xmm2
movupd %xmm3,(%rcx,%rdx,8)
movupd %xmm2,0x10(%rcx,%rdx,8)
movupd 0x20(%rax,%rdx,8),%xmm1
movupd 0x30(%rax,%rdx,8),%xmm2
mulpd  %xmm0,%xmm1
addpd  %xmm4,%xmm1
mulpd  %xmm0,%xmm2
addpd  %xmm5,%xmm2
movupd %xmm1,0x20(%rcx,%rdx,8)
movupd %xmm2,0x30(%rcx,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1320 <gram_schmidt_orthogonalization+0x1e0>
testb  $0x1,0x28(%rsp)
je     13c7 <gram_schmidt_orthogonalization+0x287>
movupd (%rax,%rdx,8),%xmm1
movupd 0x10(%rax,%rdx,8),%xmm2
movupd (%rcx,%rdx,8),%xmm3
movupd 0x10(%rcx,%rdx,8),%xmm4
mulpd  %xmm0,%xmm1
addpd  %xmm3,%xmm1
mulpd  %xmm0,%xmm2
addpd  %xmm4,%xmm2
movupd %xmm1,(%rcx,%rdx,8)
movupd %xmm2,0x10(%rcx,%rdx,8)
mov    0x38(%rsp),%rdx
mov    %rdx,%rsi
cmp    %rbp,%rdx
je     11f0 <gram_schmidt_orthogonalization+0xb0>
jmp    1283 <gram_schmidt_orthogonalization+0x143>
xor    %edx,%edx
testb  $0x1,0x28(%rsp)
jne    1396 <gram_schmidt_orthogonalization+0x256>
jmp    13c7 <gram_schmidt_orthogonalization+0x287>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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