#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 */
typedef scalar_t__ int32_t;

/* Variables and functions */
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(scalar_t__);
int /*<<< orphan*/ memcpy(double *, double *, scalar_t__);
int /*<<< orphan*/ uav(scalar_t__, scalar_t__, double *, double *, double *,
                       double *, double *, double *, int);
void pinv(double *A, int32_t m, int32_t n, double *Y, int32_t *size) {
  int32_t ka = ((((m) > (n)) ? (m) : (n)) + 1);
  double *u =
      (double *)malloc(((m * m) + (m * n) + (n * n) + ka * 3) * sizeof(double));
  double *a = u + m * m;
  double *v = a + m * n;
  double *s = v + n * n;
  double *e = s + ka;
  double *w = e + ka;
  size[0] = n;
  size[1] = m;
  int32_t i, j, k, l;
  memcpy(a, A, m * n * sizeof(double));
  uav(m, n, u, a, v, s, e, w, 100);
  j = n;
  if (m < n)
    j = m;
  j = j - 1;
  k = 0;
  while ((k <= j) && (a[k * n + k] != 0.0))
    k = k + 1;
  k = k - 1;
  double sum;
  for (i = 0; i <= n - 1; i++)
    for (j = 0; j <= m - 1; j++) {
      sum = 0.0;
      for (l = 0; l <= k; l++)
        sum += v[l * n + i] * u[j * m + l] / a[l * n + l];
      Y[i * m + j] = sum;
    }
  free(u);
}


x86 assembly (clang, -O0):
<pinv>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
jle    117d <pinv+0x3d>
mov    -0x10(%rbp),%rax
mov    %rax,-0x90(%rbp)
jmp    1188 <pinv+0x48>
mov    -0x18(%rbp),%rax
mov    %rax,-0x90(%rbp)
mov    -0x90(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rdi
imul   -0x10(%rbp),%rdi
mov    -0x10(%rbp),%rax
imul   -0x18(%rbp),%rax
add    %rax,%rdi
mov    -0x18(%rbp),%rax
imul   -0x18(%rbp),%rax
add    %rax,%rdi
imul   $0x3,-0x30(%rbp),%rax
add    %rax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
imul   -0x10(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
imul   -0x18(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
mov    -0x18(%rbp),%rcx
imul   -0x18(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
mov    -0x30(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x18(%rbp),%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x10(%rbp),%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,0x8(%rax)
mov    -0x40(%rbp),%rdi
mov    -0x8(%rbp),%rsi
mov    -0x10(%rbp),%rdx
imul   -0x18(%rbp),%rdx
shl    $0x3,%rdx
call   1040 <memcpy@plt>
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x38(%rbp),%rdx
mov    -0x40(%rbp),%rcx
mov    -0x48(%rbp),%r8
mov    -0x50(%rbp),%r9
mov    -0x58(%rbp),%r10
mov    -0x60(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
movl   $0x64,0x10(%rsp)
call   1060 <uav@plt>
mov    -0x18(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    12c0 <pinv+0x180>
mov    -0x10(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x70(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x70(%rbp)
movq   $0x0,-0x78(%rbp)
mov    -0x78(%rbp),%rcx
xor    %eax,%eax
cmp    -0x70(%rbp),%rcx
mov    %al,-0x91(%rbp)
jg     1315 <pinv+0x1d5>
mov    -0x40(%rbp),%rax
mov    -0x78(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
setne  %al
setp   %cl
or     %cl,%al
mov    %al,-0x91(%rbp)
mov    -0x91(%rbp),%al
test   $0x1,%al
jne    1328 <pinv+0x1e8>
jmp    1339 <pinv+0x1f9>
mov    -0x78(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x78(%rbp)
jmp    12d4 <pinv+0x194>
mov    -0x78(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x78(%rbp)
movq   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%rax
mov    -0x18(%rbp),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
jg     1448 <pinv+0x308>
movq   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
jg     1432 <pinv+0x2f2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movq   $0x0,-0x80(%rbp)
mov    -0x80(%rbp),%rax
cmp    -0x78(%rbp),%rax
jg     1403 <pinv+0x2c3>
mov    -0x48(%rbp),%rax
mov    -0x80(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x70(%rbp),%rcx
imul   -0x10(%rbp),%rcx
add    -0x80(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x80(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x80(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x80(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x80(%rbp)
jmp    1392 <pinv+0x252>
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%rcx
imul   -0x10(%rbp),%rcx
add    -0x70(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x70(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
jmp    136a <pinv+0x22a>
jmp    1437 <pinv+0x2f7>
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    134d <pinv+0x20d>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pinv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8,0x28(%rsp)
mov    %rcx,%r14
mov    %rdx,%r12
mov    %rsi,%r15
mov    %rdi,0x30(%rsp)
mov    %rsi,%rbp
imul   %rsi,%rbp
mov    %rdx,%rax
imul   %rsi,%rax
mov    %rax,0x20(%rsp)
add    %rbp,%rax
mov    %rdx,%rcx
imul   %rdx,%rcx
mov    %rcx,0x8(%rsp)
add    %rcx,%rax
cmp    %rdx,%rsi
mov    %rdx,%rcx
cmovg  %rsi,%rcx
mov    %rcx,(%rsp)
lea    0x3(%rcx,%rcx,2),%rdi
mov    %rdx,%rbx
cmovl  %rsi,%rbx
add    %rax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,0x10(%rsp)
lea    (%rax,%rbp,8),%r13
mov    0x20(%rsp),%rdx
lea    0x0(,%rdx,8),%rbp
add    %r13,%rbp
mov    0x8(%rsp),%rax
lea    0x0(%rbp,%rax,8),%rax
mov    %rax,0x18(%rsp)
mov    (%rsp),%rcx
lea    0x8(%rax,%rcx,8),%rax
mov    %rax,0x8(%rsp)
lea    0x8(%rax,%rcx,8),%rax
mov    %rax,(%rsp)
mov    0x28(%rsp),%rax
mov    %r12,(%rax)
mov    %r15,0x8(%rax)
shl    $0x3,%rdx
mov    %r13,%rdi
mov    0x30(%rsp),%rsi
call   1040 <memcpy@plt>
sub    $0x8,%rsp
mov    %r15,%rdi
mov    %r12,%rsi
mov    0x18(%rsp),%rdx
mov    %r13,%rcx
mov    %rbp,%r8
mov    0x20(%rsp),%r9
push   $0x64
push   0x10(%rsp)
push   0x20(%rsp)
call   1060 <uav@plt>
add    $0x20,%rsp
test   %rbx,%rbx
jle    1275 <pinv+0x135>
lea    0x8(,%r12,8),%rcx
xor    %eax,%eax
xorpd  %xmm0,%xmm0
mov    %r13,%rdx
nopl   0x0(%rax,%rax,1)
movsd  (%rdx),%xmm1
ucomisd %xmm0,%xmm1
jne    125c <pinv+0x11c>
jnp    126b <pinv+0x12b>
add    $0x1,%rax
add    %rcx,%rdx
cmp    %rax,%rbx
jne    1250 <pinv+0x110>
mov    %rbx,%rax
test   %r12,%r12
jg     1280 <pinv+0x140>
jmp    1324 <pinv+0x1e4>
xor    %eax,%eax
test   %r12,%r12
jle    1324 <pinv+0x1e4>
lea    0x8(,%r12,8),%rcx
lea    0x0(,%r15,8),%r9
xor    %r8d,%r8d
jmp    12ad <pinv+0x16d>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r8
add    $0x8,%rbp
cmp    %r12,%r8
je     1324 <pinv+0x1e4>
test   %r15,%r15
jle    12a0 <pinv+0x160>
mov    %r8,%r10
imul   %r15,%r10
mov    0x10(%rsp),%rsi
xor    %r11d,%r11d
jmp    12e6 <pinv+0x1a6>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    (%r11,%r10,1),%rdx
movsd  %xmm0,(%r14,%rdx,8)
add    $0x1,%r11
add    %r9,%rsi
cmp    %r15,%r11
je     12a0 <pinv+0x160>
xorpd  %xmm0,%xmm0
test   %rax,%rax
je     12d0 <pinv+0x190>
xor    %edi,%edi
mov    %r13,%rbx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rdi,8),%xmm1
mulsd  (%rsi,%rdx,8),%xmm1
divsd  (%rbx),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
add    %rcx,%rbx
add    %r12,%rdi
cmp    %rdx,%rax
jne    1300 <pinv+0x1c0>
jmp    12d0 <pinv+0x190>
mov    0x10(%rsp),%rdi
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

LLVM IR (-O2):
define void @pinv(double* %0, i64 %1, i64 %2, double* nocapture writeonly %3, i64* nocapture writeonly %4) {
%6 = icmp sgt i64 %1, %2
%7 = select i1 %6, i64 %1, i64 %2
%8 = add nsw i64 %7, 1
%9 = mul nsw i64 %1, %1
%10 = mul nsw i64 %2, %1
%11 = add nsw i64 %10, %9
%12 = mul nsw i64 %2, %2
%13 = add nsw i64 %11, %12
%14 = mul nsw i64 %8, 3
%15 = add nsw i64 %13, %14
%16 = shl i64 %15, 3
%17 = tail call i64 @malloc(i64 %16)
%18 = inttoptr i64 %17 to double*
%19 = getelementptr inbounds double, double* %18, i64 %9
%20 = getelementptr inbounds double, double* %19, i64 %10
%21 = getelementptr inbounds double, double* %20, i64 %12
%22 = getelementptr inbounds double, double* %21, i64 %8
%23 = getelementptr inbounds double, double* %22, i64 %8
store i64 %2, i64* %4
%24 = getelementptr inbounds i64, i64* %4, i64 1
store i64 %1, i64* %24
%25 = shl i64 %10, 3
%26 = tail call i32 @memcpy(double* %19, double* %0, i64 %25)
%27 = tail call i32 @uav(i64 %1, i64 %2, double* %18, double* %19, double* %20, double* %21, double* %22, double* %23, i32 100)
%28 = icmp slt i64 %1, %2
%29 = select i1 %28, i64 %1, i64 %2
%30 = icmp sgt i64 %29, 0
br i1 %30, label %31, label %41
31: ; preds = %5, %38
%32 = phi i64 [ %39, %38 ], [ 0, %5 ]
%33 = mul nsw i64 %32, %2
%34 = add nsw i64 %33, %32
%35 = getelementptr inbounds double, double* %19, i64 %34
%36 = load double, double* %35
%37 = fcmp une double %36, 0.000000e+00
br i1 %37, label %38, label %41
38: ; preds = %31
%39 = add nuw nsw i64 %32, 1
%40 = icmp eq i64 %39, %29
br i1 %40, label %41, label %31
41: ; preds = %31, %38, %5
%42 = phi i64 [ 0, %5 ], [ %29, %38 ], [ %32, %31 ]
%43 = icmp sgt i64 %2, 0
br i1 %43, label %44, label %82
44: ; preds = %41
%45 = icmp sgt i64 %1, 0
%46 = icmp eq i64 %42, 0
br label %47
47: ; preds = %44, %79
%48 = phi i64 [ 0, %44 ], [ %80, %79 ]
br i1 %45, label %49, label %79
49: ; preds = %47
%50 = mul nsw i64 %48, %1
br label %51
51: ; preds = %49, %73
%52 = phi i64 [ 0, %49 ], [ %77, %73 ]
br i1 %46, label %73, label %53
53: ; preds = %51
%54 = mul nsw i64 %52, %1
br label %55
55: ; preds = %53, %55
%56 = phi double [ 0.000000e+00, %53 ], [ %70, %55 ]
%57 = phi i64 [ 0, %53 ], [ %71, %55 ]
%58 = mul nsw i64 %57, %2
%59 = add nsw i64 %58, %48
%60 = getelementptr inbounds double, double* %20, i64 %59
%61 = load double, double* %60
%62 = add nsw i64 %57, %54
%63 = getelementptr inbounds double, double* %18, i64 %62
%64 = load double, double* %63
%65 = fmul double %61, %64
%66 = add nsw i64 %58, %57
%67 = getelementptr inbounds double, double* %19, i64 %66
%68 = load double, double* %67
%69 = fdiv double %65, %68
%70 = fadd double %56, %69
%71 = add nuw nsw i64 %57, 1
%72 = icmp eq i64 %71, %42
br i1 %72, label %73, label %55
73: ; preds = %55, %51
%74 = phi double [ 0.000000e+00, %51 ], [ %70, %55 ]
%75 = add nsw i64 %52, %50
%76 = getelementptr inbounds double, double* %3, i64 %75
store double %74, double* %76
%77 = add nuw nsw i64 %52, 1
%78 = icmp eq i64 %77, %1
br i1 %78, label %79, label %51
79: ; preds = %73, %47
%80 = add nuw nsw i64 %48, 1
%81 = icmp eq i64 %80, %2
br i1 %81, label %82, label %47
82: ; preds = %79, %41
%83 = tail call i32 @free(double* %18)
ret void
}