#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 fabs(double);
double round(double);
void nearest_integer_shift(double *mu, double *b, int k, int l, int n,
                           double lc) {
  int i, j;
  double r = round(mu[k * n + l]);
  if (fabs(mu[k * n + l]) > lc) {
    for (i = 0; i < n; i++)
      b[i * n + k] -= r * b[i * n + l];
    for (j = 0; j < l - 1; j++)
      mu[k * n + j] -= r * mu[l * n + j];
    mu[k * n + l] -= r;
  }
}


x86 assembly (clang, -O0):
<nearest_integer_shift>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <round@plt>
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe94(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    1290 <nearest_integer_shift+0x180>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11f1 <nearest_integer_shift+0xe1>
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%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    1182 <nearest_integer_shift+0x72>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    126c <nearest_integer_shift+0x15c>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%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    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11f8 <nearest_integer_shift+0xe8>
movsd  -0x38(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nearest_integer_shift>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
movsd  %xmm0,(%rsp)
mov    %r8d,%r13d
mov    %ecx,%r14d
mov    %edx,%ebx
mov    %rsi,0x8(%rsp)
mov    %rdi,%rbp
mov    %r8d,%r12d
imul   %edx,%r12d
lea    (%r12,%r14,1),%eax
movslq %eax,%r15
movsd  (%rdi,%r15,8),%xmm0
movapd %xmm0,0x10(%rsp)
call   1030 <round@plt>
movapd 0x10(%rsp),%xmm1
andpd  0xea0(%rip),%xmm1
ucomisd (%rsp),%xmm1
jbe    1442 <nearest_integer_shift+0x332>
mov    %r12,%r10
mov    %r15,%r11
test   %r13d,%r13d
mov    0x8(%rsp),%r12
mov    %r14,%r15
jle    122c <nearest_integer_shift+0x11c>
mov    %rbp,(%rsp)
movapd 0xe82(%rip),%xmm1
xorpd  %xmm0,%xmm1
mov    %r13d,%r8d
movslq %r15d,%r9
movslq %ebx,%rcx
mov    %r13d,%r14d
cmp    $0x1,%r13d
jne    11a8 <nearest_integer_shift+0x98>
xor    %esi,%esi
jmp    1203 <nearest_integer_shift+0xf3>
mov    %r8d,%edi
and    $0xfffffffe,%edi
lea    (%rcx,%r8,1),%rbx
mov    %r8,%rdx
shl    $0x4,%rdx
lea    (%r9,%r8,1),%rbp
xor    %esi,%esi
mov    %r12,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rax,%r9,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rax,%rcx,8),%xmm2
movsd  %xmm2,(%rax,%rcx,8)
movsd  (%rax,%rbp,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rax,%rbx,8),%xmm2
movsd  %xmm2,(%rax,%rbx,8)
add    $0x2,%rsi
add    %rdx,%rax
cmp    %rsi,%rdi
jne    11d0 <nearest_integer_shift+0xc0>
test   $0x1,%r8b
mov    (%rsp),%rbp
mov    %r14d,%r13d
je     122c <nearest_integer_shift+0x11c>
imul   %r8,%rsi
add    %rsi,%r9
mulsd  (%r12,%r9,8),%xmm1
add    %rcx,%rsi
addsd  (%r12,%rsi,8),%xmm1
movsd  %xmm1,(%r12,%rsi,8)
cmp    $0x2,%r15d
jl     1430 <nearest_integer_shift+0x320>
imul   %r15d,%r13d
add    $0xffffffff,%r15d
movapd 0xdca(%rip),%xmm1
xorpd  %xmm0,%xmm1
movslq %r13d,%r8
movslq %r10d,%rdx
cmp    $0x4,%r15d
jb     129c <nearest_integer_shift+0x18c>
lea    0x0(,%rdx,8),%rax
add    %rbp,%rax
lea    (%r8,%r15,1),%rcx
lea    0x0(,%rcx,8),%rcx
add    %rbp,%rcx
cmp    %rcx,%rax
jae    1327 <nearest_integer_shift+0x217>
lea    (%rdx,%r15,1),%rax
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
lea    0x0(,%r8,8),%rcx
add    %rbp,%rcx
cmp    %rax,%rcx
jae    1327 <nearest_integer_shift+0x217>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
test   $0x1,%r15b
je     12cc <nearest_integer_shift+0x1bc>
lea    (%rax,%r8,1),%rsi
lea    (%rax,%rdx,1),%rdi
movsd  0x0(%rbp,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x0(%rbp,%rdi,8),%xmm2
movsd  %xmm2,0x0(%rbp,%rdi,8)
or     $0x1,%rax
add    %r15,%rcx
je     1430 <nearest_integer_shift+0x320>
lea    0x8(,%rdx,8),%rcx
add    %rbp,%rcx
lea    0x8(,%r8,8),%rdx
add    %rbp,%rdx
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rdx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  -0x8(%rcx,%rax,8),%xmm2
movsd  %xmm2,-0x8(%rcx,%rax,8)
movsd  (%rdx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rcx,%rax,8),%xmm2
movsd  %xmm2,(%rcx,%rax,8)
add    $0x2,%rax
cmp    %rax,%r15
jne    12f0 <nearest_integer_shift+0x1e0>
jmp    1430 <nearest_integer_shift+0x320>
mov    %rbp,%r9
mov    %r15d,%eax
and    $0xfffffffc,%eax
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
lea    -0x4(%rax),%rcx
mov    %rcx,%rbp
shr    $0x2,%rbp
add    $0x1,%rbp
test   %rcx,%rcx
je     13e1 <nearest_integer_shift+0x2d1>
mov    %rbp,%rdi
and    $0xfffffffffffffffe,%rdi
mov    %r9,%rcx
lea    (%r9,%r8,8),%rbx
add    $0x30,%rbx
lea    (%r9,%rdx,8),%rcx
add    $0x30,%rcx
xor    %esi,%esi
nopl   0x0(%rax)
movupd -0x30(%rbx,%rsi,8),%xmm3
movupd -0x20(%rbx,%rsi,8),%xmm4
mulpd  %xmm2,%xmm3
movupd -0x30(%rcx,%rsi,8),%xmm5
addpd  %xmm3,%xmm5
movupd -0x20(%rcx,%rsi,8),%xmm3
movupd -0x10(%rcx,%rsi,8),%xmm6
movupd (%rcx,%rsi,8),%xmm7
mulpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm5,-0x30(%rcx,%rsi,8)
movupd %xmm4,-0x20(%rcx,%rsi,8)
movupd -0x10(%rbx,%rsi,8),%xmm3
movupd (%rbx,%rsi,8),%xmm4
mulpd  %xmm2,%xmm3
addpd  %xmm6,%xmm3
mulpd  %xmm2,%xmm4
addpd  %xmm7,%xmm4
movupd %xmm3,-0x10(%rcx,%rsi,8)
movupd %xmm4,(%rcx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1370 <nearest_integer_shift+0x260>
jmp    13e3 <nearest_integer_shift+0x2d3>
xor    %esi,%esi
test   $0x1,%bpl
mov    %r9,%rbp
je     1427 <nearest_integer_shift+0x317>
lea    (%rsi,%r8,1),%rcx
movupd 0x0(%rbp,%rcx,8),%xmm3
movupd 0x10(%rbp,%rcx,8),%xmm4
add    %rdx,%rsi
movupd 0x0(%rbp,%rsi,8),%xmm5
movupd 0x10(%rbp,%rsi,8),%xmm6
mulpd  %xmm2,%xmm3
addpd  %xmm5,%xmm3
mulpd  %xmm2,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm3,0x0(%rbp,%rsi,8)
movupd %xmm4,0x10(%rbp,%rsi,8)
cmp    %r15,%rax
jne    129e <nearest_integer_shift+0x18e>
movsd  0x0(%rbp,%r11,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp,%r11,8)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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