void dt(double *m, int rows, int cols) {
  const int n = ((rows) > (cols) ? (rows) : (cols));
  double *f = (double *)malloc(sizeof(f[0]) * n);
  double *d = (double *)malloc(sizeof(d[0]) * n);
  double *z = (double *)malloc(sizeof(z[0]) * (n + 1));
  int *v = (int *)malloc(sizeof(v[0]) * n);
  int x, y;
  for (x = 0; x < cols; ++x) {
    for (y = 0; y < rows; ++y) {
      f[y] = m[y * cols + x];
    }
    dt_row(f, rows, d, z, v);
    for (y = 0; y < rows; ++y) {
      m[y * cols + x] = d[y];
    }
  }
  for (y = 0; y < rows; ++y) {
    for (x = 0; x < cols; ++x) {
      f[x] = m[y * cols + x];
    }
    dt_row(f, cols, d, z, v);
    for (x = 0; x < cols; ++x) {
      m[y * cols + x] = d[x];
    }
  }
  free(f);
  free(d);
  free(z);
  free(v);
}


x86 assembly (clang, -O0):
<dt>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jle    1159 <dt+0x29>
mov    -0xc(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    115f <dt+0x2f>
mov    -0x10(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x14(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x20(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x30(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x2,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1276 <dt+0x146>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1205 <dt+0xd5>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11c8 <dt+0x98>
mov    -0x20(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%r8
mov    $0x0,%al
call   1050 <dt_row@plt>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1263 <dt+0x133>
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1226 <dt+0xf6>
jmp    1268 <dt+0x138>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11b5 <dt+0x85>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    133e <dt+0x20e>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12cd <dt+0x19d>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1290 <dt+0x160>
mov    -0x20(%rbp),%rdi
mov    -0x10(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%r8
mov    $0x0,%al
call   1050 <dt_row@plt>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    132b <dt+0x1fb>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12ee <dt+0x1be>
jmp    1330 <dt+0x200>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    127d <dt+0x14d>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dt>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %edx,%r15d
mov    %esi,%r13d
mov    %rdi,0x28(%rsp)
cmp    %edx,%esi
mov    %edx,%eax
cmovg  %esi,%eax
movslq %eax,%r14
lea    0x0(,%r14,8),%rbp
mov    %rbp,%rdi
call   1050 <malloc@plt>
mov    %rax,%rbx
mov    %rbp,%rdi
call   1050 <malloc@plt>
mov    %rax,%r12
lea    0x8(,%r14,8),%rdi
call   1050 <malloc@plt>
mov    %rax,0x20(%rsp)
shl    $0x2,%r14
mov    %r14,%rdi
call   1050 <malloc@plt>
mov    %r13d,%r9d
mov    %rax,0x18(%rsp)
mov    %r15d,%r8d
mov    %r13d,%eax
mov    %rax,0x30(%rsp)
mov    %r15d,0x3c(%rsp)
test   %r15d,%r15d
jle    13ea <dt+0x2aa>
mov    0x30(%rsp),%rax
lea    -0x1(%rax),%rcx
mov    %rcx,0x50(%rsp)
mov    %eax,%r13d
and    $0x3,%r13d
mov    %eax,%edi
and    $0xfffffffc,%edi
mov    %r8,%rsi
shl    $0x4,%rsi
mov    %r8,%rbp
shl    $0x5,%rbp
lea    0x0(,%r8,8),%r15
lea    (%r15,%r15,2),%r14
mov    0x28(%rsp),%r10
xor    %r11d,%r11d
mov    %r8,0x10(%rsp)
mov    %r9d,0xc(%rsp)
mov    %rdi,0x48(%rsp)
mov    %rsi,0x40(%rsp)
jmp    126a <dt+0x12a>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rbx,%rdi
mov    0xc(%rsp),%esi
mov    %r12,%rdx
mov    0x20(%rsp),%rcx
mov    0x18(%rsp),%r8
xor    %eax,%eax
call   1060 <dt_row@plt>
mov    0x58(%rsp),%r11
mov    0x60(%rsp),%r10
mov    0x40(%rsp),%rsi
mov    0x48(%rsp),%rdi
mov    0xc(%rsp),%r9d
mov    0x10(%rsp),%r8
add    $0x1,%r11
add    $0x8,%r10
cmp    %r8,%r11
je     13ea <dt+0x2aa>
test   %r9d,%r9d
mov    %r10,0x60(%rsp)
mov    %r11,0x58(%rsp)
jle    1220 <dt+0xe0>
cmpq   $0x3,0x50(%rsp)
jae    1290 <dt+0x150>
xor    %eax,%eax
jmp    12d8 <dt+0x198>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r10,%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
movsd  (%rcx,%r8,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rax,8)
movsd  (%rcx,%rsi,1),%xmm0
movsd  %xmm0,0x10(%rbx,%rax,8)
movsd  (%rcx,%r14,1),%xmm0
movsd  %xmm0,0x18(%rbx,%rax,8)
add    $0x4,%rax
add    %rbp,%rcx
cmp    %rax,%rdi
jne    12a0 <dt+0x160>
test   %r13,%r13
je     1315 <dt+0x1d5>
lea    (%rbx,%rax,8),%rcx
imul   %r8,%rax
add    %r11,%rax
mov    0x28(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rax),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%rdx
add    %r15,%rax
cmp    %rdx,%r13
jne    1300 <dt+0x1c0>
mov    %rbx,%rdi
mov    0xc(%rsp),%esi
mov    %r12,%rdx
mov    0x20(%rsp),%rcx
mov    0x18(%rsp),%r8
xor    %eax,%eax
call   1060 <dt_row@plt>
mov    0x58(%rsp),%r11
mov    0xc(%rsp),%r9d
mov    0x10(%rsp),%r8
test   %r9d,%r9d
mov    0x48(%rsp),%rdi
mov    0x40(%rsp),%rsi
mov    0x60(%rsp),%r10
jle    1259 <dt+0x119>
cmpq   $0x3,0x50(%rsp)
jae    1363 <dt+0x223>
xor    %eax,%eax
jmp    13ac <dt+0x26c>
mov    %r10,%rcx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movsd  (%r12,%rax,8),%xmm0
movsd  %xmm0,(%rcx)
movsd  0x8(%r12,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%r15,1)
movsd  0x10(%r12,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rsi,1)
movsd  0x18(%r12,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%r14,1)
add    $0x4,%rax
add    %rbp,%rcx
cmp    %rax,%rdi
jne    1370 <dt+0x230>
testb  $0x3,0x30(%rsp)
je     1259 <dt+0x119>
mov    %r8,%rcx
imul   %rax,%rcx
add    %r11,%rcx
mov    0x28(%rsp),%rdx
lea    (%rdx,%rcx,8),%rcx
lea    (%r12,%rax,8),%rax
xor    %edx,%edx
movsd  (%rax,%rdx,8),%xmm0
movsd  %xmm0,(%rcx)
add    $0x1,%rdx
add    %r15,%rcx
cmp    %rdx,%r13
jne    13d0 <dt+0x290>
jmp    1259 <dt+0x119>
test   %r9d,%r9d
mov    0x3c(%rsp),%r13d
mov    0x30(%rsp),%rbp
jle    1483 <dt+0x343>
shl    $0x3,%r8
xor    %r15d,%r15d
mov    %r8,0x10(%rsp)
jmp    1462 <dt+0x322>
nopl   0x0(%rax,%rax,1)
movslq %r15d,%rax
mov    0x28(%rsp),%rcx
lea    (%rcx,%rax,8),%r14
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r8,%rdx
call   1040 <memcpy@plt>
mov    %rbx,%rdi
mov    %r13d,%esi
mov    %r12,%rdx
mov    0x20(%rsp),%rcx
mov    0x18(%rsp),%r8
xor    %eax,%eax
call   1060 <dt_row@plt>
mov    %r14,%rdi
mov    %r12,%rsi
mov    0x10(%rsp),%rdx
call   1040 <memcpy@plt>
mov    0x10(%rsp),%r8
add    %r13d,%r15d
add    $0xffffffffffffffff,%rbp
je     1483 <dt+0x343>
test   %r13d,%r13d
jg     1410 <dt+0x2d0>
mov    %rbx,%rdi
mov    %r13d,%esi
mov    %r12,%rdx
mov    0x20(%rsp),%rcx
mov    0x18(%rsp),%r8
xor    %eax,%eax
call   1060 <dt_row@plt>
jmp    1454 <dt+0x314>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    0x20(%rsp),%rdi
call   1030 <free@plt>
mov    0x18(%rsp),%rdi
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @dt(double* %0, i32 %1, i32 %2) {
%4 = alloca double*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32*
%12 = alloca i32
%13 = alloca i32
store double* %0, double** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
%14 = load i32, i32* %5
%15 = load i32, i32* %6
%16 = icmp sgt i32 %14, %15
br i1 %16, label %17, label %19
17: ; preds = %3
%18 = load i32, i32* %5
br label %21
19: ; preds = %3
%20 = load i32, i32* %6
br label %21
21: ; preds = %19, %17
%22 = phi i32 [ %18, %17 ], [ %20, %19 ]
store i32 %22, i32* %7
%23 = load i32, i32* %7
%24 = sext i32 %23 to i64
%25 = mul i64 8, %24
%26 = call i8* @malloc(i64 %25)
%27 = bitcast i8* %26 to double*
store double* %27, double** %8
%28 = load i32, i32* %7
%29 = sext i32 %28 to i64
%30 = mul i64 8, %29
%31 = call i8* @malloc(i64 %30)
%32 = bitcast i8* %31 to double*
store double* %32, double** %9
%33 = load i32, i32* %7
%34 = add nsw i32 %33, 1
%35 = sext i32 %34 to i64
%36 = mul i64 8, %35
%37 = call i8* @malloc(i64 %36)
%38 = bitcast i8* %37 to double*
store double* %38, double** %10
%39 = load i32, i32* %7
%40 = sext i32 %39 to i64
%41 = mul i64 4, %40
%42 = call i8* @malloc(i64 %41)
%43 = bitcast i8* %42 to i32*
store i32* %43, i32** %11
store i32 0, i32* %12
br label %44
44: ; preds = %99, %21
%45 = load i32, i32* %12
%46 = load i32, i32* %6
%47 = icmp slt i32 %45, %46
br i1 %47, label %48, label %102
48: ; preds = %44
store i32 0, i32* %13
br label %49
49: ; preds = %67, %48
%50 = load i32, i32* %13
%51 = load i32, i32* %5
%52 = icmp slt i32 %50, %51
br i1 %52, label %53, label %70
53: ; preds = %49
%54 = load double*, double** %4
%55 = load i32, i32* %13
%56 = load i32, i32* %6
%57 = mul nsw i32 %55, %56
%58 = load i32, i32* %12
%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 = load double*, double** %8
%64 = load i32, i32* %13
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
store double %62, double* %66
br label %67
67: ; preds = %53
%68 = load i32, i32* %13
%69 = add nsw i32 %68, 1
store i32 %69, i32* %13
br label %49
70: ; preds = %49
%71 = load double*, double** %8
%72 = load i32, i32* %5
%73 = load double*, double** %9
%74 = load double*, double** %10
%75 = load i32*, i32** %11
%76 = call i32 (double*, i32, double*, double*, i32*, ...) bitcast (i32 (...)* @dt_row to i32 (double*, i32, double*, double*, i32*, ...)*)(double* %71, i32 %72, double* %73, double* %74, i32* %75)
store i32 0, i32* %13
br label %77
77: ; preds = %95, %70
%78 = load i32, i32* %13
%79 = load i32, i32* %5
%80 = icmp slt i32 %78, %79
br i1 %80, label %81, label %98
81: ; preds = %77
%82 = load double*, double** %9
%83 = load i32, i32* %13
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %82, i64 %84
%86 = load double, double* %85
%87 = load double*, double** %4
%88 = load i32, i32* %13
%89 = load i32, i32* %6
%90 = mul nsw i32 %88, %89
%91 = load i32, i32* %12
%92 = add nsw i32 %90, %91
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %87, i64 %93
store double %86, double* %94
br label %95
95: ; preds = %81
%96 = load i32, i32* %13
%97 = add nsw i32 %96, 1
store i32 %97, i32* %13
br label %77
98: ; preds = %77
br label %99
99: ; preds = %98
%100 = load i32, i32* %12
%101 = add nsw i32 %100, 1
store i32 %101, i32* %12
br label %44
102: ; preds = %44
store i32 0, i32* %13
br label %103
103: ; preds = %158, %102
%104 = load i32, i32* %13
%105 = load i32, i32* %5
%106 = icmp slt i32 %104, %105
br i1 %106, label %107, label %161
107: ; preds = %103
store i32 0, i32* %12
br label %108
108: ; preds = %126, %107
%109 = load i32, i32* %12
%110 = load i32, i32* %6
%111 = icmp slt i32 %109, %110
br i1 %111, label %112, label %129
112: ; preds = %108
%113 = load double*, double** %4
%114 = load i32, i32* %13
%115 = load i32, i32* %6
%116 = mul nsw i32 %114, %115
%117 = load i32, i32* %12
%118 = add nsw i32 %116, %117
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %113, i64 %119
%121 = load double, double* %120
%122 = load double*, double** %8
%123 = load i32, i32* %12
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %122, i64 %124
store double %121, double* %125
br label %126
126: ; preds = %112
%127 = load i32, i32* %12
%128 = add nsw i32 %127, 1
store i32 %128, i32* %12
br label %108
129: ; preds = %108
%130 = load double*, double** %8
%131 = load i32, i32* %6
%132 = load double*, double** %9
%133 = load double*, double** %10
%134 = load i32*, i32** %11
%135 = call i32 (double*, i32, double*, double*, i32*, ...) bitcast (i32 (...)* @dt_row to i32 (double*, i32, double*, double*, i32*, ...)*)(double* %130, i32 %131, double* %132, double* %133, i32* %134)
store i32 0, i32* %12
br label %136
136: ; preds = %154, %129
%137 = load i32, i32* %12
%138 = load i32, i32* %6
%139 = icmp slt i32 %137, %138
br i1 %139, label %140, label %157
140: ; preds = %136
%141 = load double*, double** %9
%142 = load i32, i32* %12
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %141, i64 %143
%145 = load double, double* %144
%146 = load double*, double** %4
%147 = load i32, i32* %13
%148 = load i32, i32* %6
%149 = mul nsw i32 %147, %148
%150 = load i32, i32* %12
%151 = add nsw i32 %149, %150
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %146, i64 %152
store double %145, double* %153
br label %154
154: ; preds = %140
%155 = load i32, i32* %12
%156 = add nsw i32 %155, 1
store i32 %156, i32* %12
br label %136
157: ; preds = %136
br label %158
158: ; preds = %157
%159 = load i32, i32* %13
%160 = add nsw i32 %159, 1
store i32 %160, i32* %13
br label %103
161: ; preds = %103
%162 = load double*, double** %8
%163 = bitcast double* %162 to i8*
%164 = load double*, double** %9
%165 = bitcast double* %164 to i8*
%166 = load double*, double** %10
%167 = bitcast double* %166 to i8*
%168 = load i32*, i32** %11
%169 = bitcast i32* %168 to i8*
ret void
}

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