#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 **dmatrix(int, int, int, int);
int /*<<< orphan*/ free_dmatrix(double **, int, int, int, int);
int /*<<< orphan*/ mglin(double **, int, int);
int /*<<< orphan*/ printf(char *, ...);
int main(void) {
  int i, j, midl = 33 / 2 + 1;
  double **f, **u;
  f = dmatrix(1, 33, 1, 33);
  u = dmatrix(1, 33, 1, 33);
  for (i = 1; i <= 33; i++)
    for (j = 1; j <= 33; j++)
      u[i][j] = 0.0;
  u[midl][midl] = 2.0;
  mglin(u, 33, 2);
  printf("MGLIN solution:\n");
  for (i = 1; i <= 33; i += 4) {
    for (j = 1; j <= 33; j += 4)
      printf("%8.4f", u[i][j]);
    printf("\n");
  }
  printf("\n Test that solution satisfies difference equations:\n");
  for (i = 4 + 1; i < 33; i += 4) {
    for (j = 4 + 1; j < 33; j += 4)
      f[i][j] =
          u[i + 1][j] + u[i - 1][j] + u[i][j + 1] + u[i][j - 1] - 4.0 * u[i][j];
    printf("%7s", " ");
    for (j = 4 + 1; j < 33; j += 4)
      printf("%8.4f", f[i][j] * (33 - 1) * (33 - 1));
    printf("\n");
  }
  free_dmatrix(u, 1, 33, 1, 33);
  free_dmatrix(f, 1, 33, 1, 33);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movl   $0x0,-0x4(%rbp)
movl   $0x11,-0x10(%rbp)
mov    $0x1,%edx
mov    $0x21,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   1040 <dmatrix@plt>
mov    %rax,-0x18(%rbp)
mov    $0x1,%edx
mov    $0x21,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   1040 <dmatrix@plt>
mov    %rax,-0x20(%rbp)
movl   $0x1,-0x8(%rbp)
cmpl   $0x21,-0x8(%rbp)
jg     11df <main+0x9f>
movl   $0x1,-0xc(%rbp)
cmpl   $0x21,-0xc(%rbp)
jg     11cc <main+0x8c>
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    119c <main+0x5c>
jmp    11d1 <main+0x91>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    118b <main+0x4b>
mov    -0x20(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  0xe09(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rdi
mov    $0x21,%esi
mov    $0x2,%edx
call   1060 <mglin@plt>
lea    0xe02(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x1,-0x8(%rbp)
cmpl   $0x21,-0x8(%rbp)
jg     128c <main+0x14c>
movl   $0x1,-0xc(%rbp)
cmpl   $0x21,-0xc(%rbp)
jg     1270 <main+0x130>
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xdce(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0xc(%rbp)
jmp    1235 <main+0xf5>
lea    0xdb0(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    -0x8(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x8(%rbp)
jmp    1224 <main+0xe4>
lea    0xd9c(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x5,-0x8(%rbp)
cmpl   $0x21,-0x8(%rbp)
jge    13f3 <main+0x2b3>
movl   $0x5,-0xc(%rbp)
cmpl   $0x21,-0xc(%rbp)
jge    136c <main+0x22c>
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xccf(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0xc(%rbp)
jmp    12b2 <main+0x172>
lea    0xcf2(%rip),%rdi
lea    0xcef(%rip),%rsi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x5,-0xc(%rbp)
cmpl   $0x21,-0xc(%rbp)
jge    13d7 <main+0x297>
mov    -0x18(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  0xc5e(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xc51(%rip),%xmm1
mulsd  %xmm1,%xmm0
lea    0xc67(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0xc(%rbp)
jmp    1388 <main+0x248>
lea    0xc49(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    -0x8(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x8(%rbp)
jmp    12a1 <main+0x161>
mov    -0x20(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x21,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1050 <free_dmatrix@plt>
mov    -0x18(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x21,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1050 <free_dmatrix@plt>
xor    %eax,%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    $0x1,%edi
mov    $0x21,%esi
mov    $0x1,%edx
mov    $0x21,%ecx
call   1070 <dmatrix@plt>
mov    %rax,%r14
mov    $0x1,%edi
mov    $0x21,%esi
mov    $0x1,%edx
mov    $0x21,%ecx
call   1070 <dmatrix@plt>
mov    %rax,%r15
mov    0x8(%rax),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x10(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x18(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x20(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x28(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x30(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x38(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x40(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x48(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x50(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x58(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x60(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x68(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x70(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x78(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x80(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x88(%r15),%rbx
lea    0x8(%rbx),%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x90(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x98(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xa0(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xa8(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xb0(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xb8(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xc0(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xc8(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xd0(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xd8(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xe0(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xe8(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xf0(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xf8(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x100(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x108(%r15),%rdi
add    $0x8,%rdi
mov    $0x108,%edx
xor    %esi,%esi
call   1060 <memset@plt>
movabs $0x4000000000000000,%rax
mov    %rax,0x88(%rbx)
mov    %r15,%rdi
mov    $0x21,%esi
mov    $0x2,%edx
call   1090 <mglin@plt>
lea    0xb75(%rip),%rdi
call   1040 <puts@plt>
mov    $0xfffffffffffffffd,%rbp
lea    0xb56(%rip),%rbx
nopw   0x0(%rax,%rax,1)
mov    0x20(%r15,%rbp,8),%rax
movsd  0x8(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0x28(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0x48(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0x68(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0x88(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0xa8(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0xc8(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0xe8(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    0x20(%r15,%rbp,8),%rax
movsd  0x108(%rax),%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    $0xa,%edi
call   1030 <putchar@plt>
add    $0x4,%rbp
cmp    $0x1e,%rbp
jb     14c0 <main+0x350>
lea    0xa8a(%rip),%rdi
call   1040 <puts@plt>
mov    $0x6,%ebx
lea    0xa63(%rip),%rbp
lea    0xa60(%rip),%r12
lea    0xa4f(%rip),%r13
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%r15,%rbx,8),%rsi
mov    -0x10(%r15,%rbx,8),%rdx
mov    -0x8(%r15,%rbx,8),%rax
movsd  0x28(%rsi),%xmm0
addsd  0x28(%rdx),%xmm0
addsd  0x30(%rax),%xmm0
addsd  0x20(%rax),%xmm0
mov    -0x8(%r14,%rbx,8),%rcx
movsd  0x28(%rax),%xmm1
movsd  0x9fc(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x28(%rcx)
movsd  0x48(%rsi),%xmm0
addsd  0x48(%rdx),%xmm0
addsd  0x50(%rax),%xmm0
addsd  0x40(%rax),%xmm0
movsd  0x48(%rax),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x48(%rcx)
movsd  0x68(%rsi),%xmm0
addsd  0x68(%rdx),%xmm0
addsd  0x70(%rax),%xmm0
addsd  0x60(%rax),%xmm0
movsd  0x68(%rax),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x68(%rcx)
movsd  0x88(%rsi),%xmm0
addsd  0x88(%rdx),%xmm0
addsd  0x90(%rax),%xmm0
addsd  0x80(%rax),%xmm0
movsd  0x88(%rax),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x88(%rcx)
movsd  0xa8(%rsi),%xmm0
addsd  0xa8(%rdx),%xmm0
addsd  0xb0(%rax),%xmm0
addsd  0xa0(%rax),%xmm0
movsd  0xa8(%rax),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0xa8(%rcx)
movsd  0xc8(%rsi),%xmm0
addsd  0xc8(%rdx),%xmm0
addsd  0xd0(%rax),%xmm0
addsd  0xc0(%rax),%xmm0
movsd  0xc8(%rax),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0xc8(%rcx)
movsd  0xe8(%rsi),%xmm0
addsd  0xe8(%rdx),%xmm0
addsd  0xf0(%rax),%xmm0
addsd  0xe0(%rax),%xmm0
movsd  0xe8(%rax),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0xe8(%rcx)
mov    %rbp,%rdi
mov    %r12,%rsi
xor    %eax,%eax
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0x28(%rax),%xmm0
movsd  0x8ac(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0x48(%rax),%xmm0
movsd  0x888(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0x68(%rax),%xmm0
movsd  0x864(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0x88(%rax),%xmm0
movsd  0x83d(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0xa8(%rax),%xmm0
movsd  0x816(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0xc8(%rax),%xmm0
movsd  0x7ef(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0xe8(%rax),%xmm0
movsd  0x7c8(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0x4(%rbx),%rax
add    $0xffffffffffffffff,%rbx
cmp    $0x1d,%rbx
mov    %rax,%rbx
jb     15d0 <main+0x460>
mov    %r15,%rdi
mov    $0x1,%esi
mov    $0x21,%edx
mov    $0x1,%ecx
mov    $0x21,%r8d
call   1080 <free_dmatrix@plt>
mov    %r14,%rdi
mov    $0x1,%esi
mov    $0x21,%edx
mov    $0x1,%ecx
mov    $0x21,%r8d
call   1080 <free_dmatrix@plt>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str.1 = [6 x i8] c"%8.4f\00"
@.str.4 = [4 x i8] c"%7s\00"
@.str.5 = [2 x i8] c" \00"
@str = [16 x i8] c"MGLIN solution:\00"
@str.6 = [53 x i8] c"\0A Test that solution satisfies difference equations:\00"
define i32 @main() {
%1 = tail call double** @dmatrix(i32 1, i32 33, i32 1, i32 33)
%2 = tail call double** @dmatrix(i32 1, i32 33, i32 1, i32 33)
%3 = getelementptr inbounds double*, double** %2, i64 1
%4 = load double*, double** %3
%5 = getelementptr double, double* %4, i64 1
%6 = bitcast double* %5 to i8*
%7 = getelementptr inbounds double*, double** %2, i64 2
%8 = load double*, double** %7
%9 = getelementptr double, double* %8, i64 1
%10 = bitcast double* %9 to i8*
%11 = getelementptr inbounds double*, double** %2, i64 3
%12 = load double*, double** %11
%13 = getelementptr double, double* %12, i64 1
%14 = bitcast double* %13 to i8*
%15 = getelementptr inbounds double*, double** %2, i64 4
%16 = load double*, double** %15
%17 = getelementptr double, double* %16, i64 1
%18 = bitcast double* %17 to i8*
%19 = getelementptr inbounds double*, double** %2, i64 5
%20 = load double*, double** %19
%21 = getelementptr double, double* %20, i64 1
%22 = bitcast double* %21 to i8*
%23 = getelementptr inbounds double*, double** %2, i64 6
%24 = load double*, double** %23
%25 = getelementptr double, double* %24, i64 1
%26 = bitcast double* %25 to i8*
%27 = getelementptr inbounds double*, double** %2, i64 7
%28 = load double*, double** %27
%29 = getelementptr double, double* %28, i64 1
%30 = bitcast double* %29 to i8*
%31 = getelementptr inbounds double*, double** %2, i64 8
%32 = load double*, double** %31
%33 = getelementptr double, double* %32, i64 1
%34 = bitcast double* %33 to i8*
%35 = getelementptr inbounds double*, double** %2, i64 9
%36 = load double*, double** %35
%37 = getelementptr double, double* %36, i64 1
%38 = bitcast double* %37 to i8*
%39 = getelementptr inbounds double*, double** %2, i64 10
%40 = load double*, double** %39
%41 = getelementptr double, double* %40, i64 1
%42 = bitcast double* %41 to i8*
%43 = getelementptr inbounds double*, double** %2, i64 11
%44 = load double*, double** %43
%45 = getelementptr double, double* %44, i64 1
%46 = bitcast double* %45 to i8*
%47 = getelementptr inbounds double*, double** %2, i64 12
%48 = load double*, double** %47
%49 = getelementptr double, double* %48, i64 1
%50 = bitcast double* %49 to i8*
%51 = getelementptr inbounds double*, double** %2, i64 13
%52 = load double*, double** %51
%53 = getelementptr double, double* %52, i64 1
%54 = bitcast double* %53 to i8*
%55 = getelementptr inbounds double*, double** %2, i64 14
%56 = load double*, double** %55
%57 = getelementptr double, double* %56, i64 1
%58 = bitcast double* %57 to i8*
%59 = getelementptr inbounds double*, double** %2, i64 15
%60 = load double*, double** %59
%61 = getelementptr double, double* %60, i64 1
%62 = bitcast double* %61 to i8*
%63 = getelementptr inbounds double*, double** %2, i64 16
%64 = load double*, double** %63
%65 = getelementptr double, double* %64, i64 1
%66 = bitcast double* %65 to i8*
%67 = getelementptr inbounds double*, double** %2, i64 17
%68 = load double*, double** %67
%69 = getelementptr double, double* %68, i64 1
%70 = bitcast double* %69 to i8*
%71 = getelementptr inbounds double*, double** %2, i64 18
%72 = load double*, double** %71
%73 = getelementptr double, double* %72, i64 1
%74 = bitcast double* %73 to i8*
%75 = getelementptr inbounds double*, double** %2, i64 19
%76 = load double*, double** %75
%77 = getelementptr double, double* %76, i64 1
%78 = bitcast double* %77 to i8*
%79 = getelementptr inbounds double*, double** %2, i64 20
%80 = load double*, double** %79
%81 = getelementptr double, double* %80, i64 1
%82 = bitcast double* %81 to i8*
%83 = getelementptr inbounds double*, double** %2, i64 21
%84 = load double*, double** %83
%85 = getelementptr double, double* %84, i64 1
%86 = bitcast double* %85 to i8*
%87 = getelementptr inbounds double*, double** %2, i64 22
%88 = load double*, double** %87
%89 = getelementptr double, double* %88, i64 1
%90 = bitcast double* %89 to i8*
%91 = getelementptr inbounds double*, double** %2, i64 23
%92 = load double*, double** %91
%93 = getelementptr double, double* %92, i64 1
%94 = bitcast double* %93 to i8*
%95 = getelementptr inbounds double*, double** %2, i64 24
%96 = load double*, double** %95
%97 = getelementptr double, double* %96, i64 1
%98 = bitcast double* %97 to i8*
%99 = getelementptr inbounds double*, double** %2, i64 25
%100 = load double*, double** %99
%101 = getelementptr double, double* %100, i64 1
%102 = bitcast double* %101 to i8*
%103 = getelementptr inbounds double*, double** %2, i64 26
%104 = load double*, double** %103
%105 = getelementptr double, double* %104, i64 1
%106 = bitcast double* %105 to i8*
%107 = getelementptr inbounds double*, double** %2, i64 27
%108 = load double*, double** %107
%109 = getelementptr double, double* %108, i64 1
%110 = bitcast double* %109 to i8*
%111 = getelementptr inbounds double*, double** %2, i64 28
%112 = load double*, double** %111
%113 = getelementptr double, double* %112, i64 1
%114 = bitcast double* %113 to i8*
%115 = getelementptr inbounds double*, double** %2, i64 29
%116 = load double*, double** %115
%117 = getelementptr double, double* %116, i64 1
%118 = bitcast double* %117 to i8*
%119 = getelementptr inbounds double*, double** %2, i64 30
%120 = load double*, double** %119
%121 = getelementptr double, double* %120, i64 1
%122 = bitcast double* %121 to i8*
%123 = getelementptr inbounds double*, double** %2, i64 31
%124 = load double*, double** %123
%125 = getelementptr double, double* %124, i64 1
%126 = bitcast double* %125 to i8*
%127 = getelementptr inbounds double*, double** %2, i64 32
%128 = load double*, double** %127
%129 = getelementptr double, double* %128, i64 1
%130 = bitcast double* %129 to i8*
%131 = getelementptr inbounds double*, double** %2, i64 33
%132 = load double*, double** %131
%133 = getelementptr double, double* %132, i64 1
%134 = bitcast double* %133 to i8*
%135 = getelementptr inbounds double, double* %68, i64 17
store double 2.000000e+00, double* %135
%136 = tail call i32 @mglin(double** %2, i32 33, i32 2)
%137 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([16 x i8], [16 x i8]* @str, i64 0, i64 0))
br label %138
138: ; preds = %0, %138
%139 = phi i64 [ 1, %0 ], [ %178, %138 ]
%140 = getelementptr inbounds double*, double** %2, i64 %139
%141 = load double*, double** %140
%142 = getelementptr inbounds double, double* %141, i64 1
%143 = load double, double* %142
%144 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %143)
%145 = load double*, double** %140
%146 = getelementptr inbounds double, double* %145, i64 5
%147 = load double, double* %146
%148 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %147)
%149 = load double*, double** %140
%150 = getelementptr inbounds double, double* %149, i64 9
%151 = load double, double* %150
%152 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %151)
%153 = load double*, double** %140
%154 = getelementptr inbounds double, double* %153, i64 13
%155 = load double, double* %154
%156 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %155)
%157 = load double*, double** %140
%158 = getelementptr inbounds double, double* %157, i64 17
%159 = load double, double* %158
%160 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %159)
%161 = load double*, double** %140
%162 = getelementptr inbounds double, double* %161, i64 21
%163 = load double, double* %162
%164 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %163)
%165 = load double*, double** %140
%166 = getelementptr inbounds double, double* %165, i64 25
%167 = load double, double* %166
%168 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %167)
%169 = load double*, double** %140
%170 = getelementptr inbounds double, double* %169, i64 29
%171 = load double, double* %170
%172 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %171)
%173 = load double*, double** %140
%174 = getelementptr inbounds double, double* %173, i64 33
%175 = load double, double* %174
%176 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %175)
%177 = tail call i32 @putchar(i32 10)
%178 = add nuw nsw i64 %139, 4
%179 = icmp ult i64 %139, 30
br i1 %179, label %138, label %180
180: ; preds = %138
%181 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([53 x i8], [53 x i8]* @str.6, i64 0, i64 0))
br label %182
182: ; preds = %180, %182
%183 = phi i64 [ 5, %180 ], [ %343, %182 ]
%184 = add nuw nsw i64 %183, 1
%185 = getelementptr inbounds double*, double** %2, i64 %184
%186 = load double*, double** %185
%187 = add nsw i64 %183, -1
%188 = getelementptr inbounds double*, double** %2, i64 %187
%189 = load double*, double** %188
%190 = getelementptr inbounds double*, double** %2, i64 %183
%191 = load double*, double** %190
%192 = getelementptr inbounds double*, double** %1, i64 %183
%193 = load double*, double** %192
%194 = getelementptr inbounds double, double* %186, i64 5
%195 = load double, double* %194
%196 = getelementptr inbounds double, double* %189, i64 5
%197 = load double, double* %196
%198 = fadd double %195, %197
%199 = getelementptr inbounds double, double* %191, i64 6
%200 = load double, double* %199
%201 = fadd double %198, %200
%202 = getelementptr inbounds double, double* %191, i64 4
%203 = load double, double* %202
%204 = fadd double %201, %203
%205 = getelementptr inbounds double, double* %191, i64 5
%206 = load double, double* %205
%207 = tail call double @llvm.fmuladd.f64(double %206, double -4.000000e+00, double %204)
%208 = getelementptr inbounds double, double* %193, i64 5
store double %207, double* %208
%209 = getelementptr inbounds double, double* %186, i64 9
%210 = load double, double* %209
%211 = getelementptr inbounds double, double* %189, i64 9
%212 = load double, double* %211
%213 = fadd double %210, %212
%214 = getelementptr inbounds double, double* %191, i64 10
%215 = load double, double* %214
%216 = fadd double %213, %215
%217 = getelementptr inbounds double, double* %191, i64 8
%218 = load double, double* %217
%219 = fadd double %216, %218
%220 = getelementptr inbounds double, double* %191, i64 9
%221 = load double, double* %220
%222 = tail call double @llvm.fmuladd.f64(double %221, double -4.000000e+00, double %219)
%223 = getelementptr inbounds double, double* %193, i64 9
store double %222, double* %223
%224 = getelementptr inbounds double, double* %186, i64 13
%225 = load double, double* %224
%226 = getelementptr inbounds double, double* %189, i64 13
%227 = load double, double* %226
%228 = fadd double %225, %227
%229 = getelementptr inbounds double, double* %191, i64 14
%230 = load double, double* %229
%231 = fadd double %228, %230
%232 = getelementptr inbounds double, double* %191, i64 12
%233 = load double, double* %232
%234 = fadd double %231, %233
%235 = getelementptr inbounds double, double* %191, i64 13
%236 = load double, double* %235
%237 = tail call double @llvm.fmuladd.f64(double %236, double -4.000000e+00, double %234)
%238 = getelementptr inbounds double, double* %193, i64 13
store double %237, double* %238
%239 = getelementptr inbounds double, double* %186, i64 17
%240 = load double, double* %239
%241 = getelementptr inbounds double, double* %189, i64 17
%242 = load double, double* %241
%243 = fadd double %240, %242
%244 = getelementptr inbounds double, double* %191, i64 18
%245 = load double, double* %244
%246 = fadd double %243, %245
%247 = getelementptr inbounds double, double* %191, i64 16
%248 = load double, double* %247
%249 = fadd double %246, %248
%250 = getelementptr inbounds double, double* %191, i64 17
%251 = load double, double* %250
%252 = tail call double @llvm.fmuladd.f64(double %251, double -4.000000e+00, double %249)
%253 = getelementptr inbounds double, double* %193, i64 17
store double %252, double* %253
%254 = getelementptr inbounds double, double* %186, i64 21
%255 = load double, double* %254
%256 = getelementptr inbounds double, double* %189, i64 21
%257 = load double, double* %256
%258 = fadd double %255, %257
%259 = getelementptr inbounds double, double* %191, i64 22
%260 = load double, double* %259
%261 = fadd double %258, %260
%262 = getelementptr inbounds double, double* %191, i64 20
%263 = load double, double* %262
%264 = fadd double %261, %263
%265 = getelementptr inbounds double, double* %191, i64 21
%266 = load double, double* %265
%267 = tail call double @llvm.fmuladd.f64(double %266, double -4.000000e+00, double %264)
%268 = getelementptr inbounds double, double* %193, i64 21
store double %267, double* %268
%269 = getelementptr inbounds double, double* %186, i64 25
%270 = load double, double* %269
%271 = getelementptr inbounds double, double* %189, i64 25
%272 = load double, double* %271
%273 = fadd double %270, %272
%274 = getelementptr inbounds double, double* %191, i64 26
%275 = load double, double* %274
%276 = fadd double %273, %275
%277 = getelementptr inbounds double, double* %191, i64 24
%278 = load double, double* %277
%279 = fadd double %276, %278
%280 = getelementptr inbounds double, double* %191, i64 25
%281 = load double, double* %280
%282 = tail call double @llvm.fmuladd.f64(double %281, double -4.000000e+00, double %279)
%283 = getelementptr inbounds double, double* %193, i64 25
store double %282, double* %283
%284 = getelementptr inbounds double, double* %186, i64 29
%285 = load double, double* %284
%286 = getelementptr inbounds double, double* %189, i64 29
%287 = load double, double* %286
%288 = fadd double %285, %287
%289 = getelementptr inbounds double, double* %191, i64 30
%290 = load double, double* %289
%291 = fadd double %288, %290
%292 = getelementptr inbounds double, double* %191, i64 28
%293 = load double, double* %292
%294 = fadd double %291, %293
%295 = getelementptr inbounds double, double* %191, i64 29
%296 = load double, double* %295
%297 = tail call double @llvm.fmuladd.f64(double %296, double -4.000000e+00, double %294)
%298 = getelementptr inbounds double, double* %193, i64 29
store double %297, double* %298
%299 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.5, i64 0, i64 0))
%300 = load double*, double** %192
%301 = getelementptr inbounds double, double* %300, i64 5
%302 = load double, double* %301
%303 = fmul double %302, 3.200000e+01
%304 = fmul double %303, 3.200000e+01
%305 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %304)
%306 = load double*, double** %192
%307 = getelementptr inbounds double, double* %306, i64 9
%308 = load double, double* %307
%309 = fmul double %308, 3.200000e+01
%310 = fmul double %309, 3.200000e+01
%311 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %310)
%312 = load double*, double** %192
%313 = getelementptr inbounds double, double* %312, i64 13
%314 = load double, double* %313
%315 = fmul double %314, 3.200000e+01
%316 = fmul double %315, 3.200000e+01
%317 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %316)
%318 = load double*, double** %192
%319 = getelementptr inbounds double, double* %318, i64 17
%320 = load double, double* %319
%321 = fmul double %320, 3.200000e+01
%322 = fmul double %321, 3.200000e+01
%323 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %322)
%324 = load double*, double** %192
%325 = getelementptr inbounds double, double* %324, i64 21
%326 = load double, double* %325
%327 = fmul double %326, 3.200000e+01
%328 = fmul double %327, 3.200000e+01
%329 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %328)
%330 = load double*, double** %192
%331 = getelementptr inbounds double, double* %330, i64 25
%332 = load double, double* %331
%333 = fmul double %332, 3.200000e+01
%334 = fmul double %333, 3.200000e+01
%335 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %334)
%336 = load double*, double** %192
%337 = getelementptr inbounds double, double* %336, i64 29
%338 = load double, double* %337
%339 = fmul double %338, 3.200000e+01
%340 = fmul double %339, 3.200000e+01
%341 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %340)
%342 = tail call i32 @putchar(i32 10)
%343 = add nuw nsw i64 %183, 4
%344 = icmp ult i64 %183, 29
br i1 %344, label %182, label %345
345: ; preds = %182
%346 = tail call i32 @free_dmatrix(double** %2, i32 1, i32 33, i32 1, i32 33)
%347 = tail call i32 @free_dmatrix(double** %1, i32 1, i32 33, i32 1, i32 33)
ret i32 0
}