#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*/ mgfas(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;
  mgfas(u, 33, 2);
  printf("MGFAS 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] + u[i][j] * u[i][j] / ((33 - 1) * (33 - 1));
    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 <mgfas@plt>
lea    0xe0a(%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    0xdd6(%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    0xdb8(%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    0xda4(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x5,-0x8(%rbp)
cmpl   $0x21,-0x8(%rbp)
jge    142d <main+0x2ed>
movl   $0x5,-0xc(%rbp)
cmpl   $0x21,-0xc(%rbp)
jge    13a6 <main+0x266>
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  0xcd7(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  0xc95(%rip),%xmm2
divsd  %xmm2,%xmm1
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    0xcc0(%rip),%rdi
lea    0xcbd(%rip),%rsi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x5,-0xc(%rbp)
cmpl   $0x21,-0xc(%rbp)
jge    1411 <main+0x2d1>
mov    -0x18(%rbp),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  0xc24(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xc17(%rip),%xmm1
mulsd  %xmm1,%xmm0
lea    0xc35(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0xc(%rbp)
jmp    13c2 <main+0x282>
lea    0xc17(%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 <mgfas@plt>
lea    0xb7d(%rip),%rdi
call   1040 <puts@plt>
mov    $0xfffffffffffffffd,%rbp
lea    0xb5e(%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    0xa92(%rip),%rdi
call   1040 <puts@plt>
mov    $0x6,%ebx
lea    0xa6b(%rip),%rbp
lea    0xa68(%rip),%r12
lea    0xa57(%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),%rcx
mov    -0x8(%r14,%rbx,8),%rax
movsd  0x28(%rsi),%xmm0
addsd  0x28(%rdx),%xmm0
addsd  0x30(%rcx),%xmm0
addsd  0x20(%rcx),%xmm0
movsd  0x28(%rcx),%xmm1
movapd %xmm1,%xmm2
movsd  0x9f8(%rip),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
movsd  0x9ec(%rip),%xmm4
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x28(%rax)
movsd  0x48(%rsi),%xmm0
addsd  0x48(%rdx),%xmm0
addsd  0x50(%rcx),%xmm0
addsd  0x40(%rcx),%xmm0
movsd  0x48(%rcx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x48(%rax)
movsd  0x68(%rsi),%xmm0
addsd  0x68(%rdx),%xmm0
addsd  0x70(%rcx),%xmm0
addsd  0x60(%rcx),%xmm0
movsd  0x68(%rcx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x68(%rax)
movsd  0x88(%rsi),%xmm0
addsd  0x88(%rdx),%xmm0
addsd  0x90(%rcx),%xmm0
addsd  0x80(%rcx),%xmm0
movsd  0x88(%rcx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x88(%rax)
movsd  0xa8(%rsi),%xmm0
addsd  0xa8(%rdx),%xmm0
addsd  0xb0(%rcx),%xmm0
addsd  0xa0(%rcx),%xmm0
movsd  0xa8(%rcx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0xa8(%rax)
movsd  0xc8(%rsi),%xmm0
addsd  0xc8(%rdx),%xmm0
addsd  0xd0(%rcx),%xmm0
addsd  0xc0(%rcx),%xmm0
movsd  0xc8(%rcx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0xc8(%rax)
movsd  0xe8(%rsi),%xmm0
addsd  0xe8(%rdx),%xmm0
addsd  0xf0(%rcx),%xmm0
addsd  0xe0(%rcx),%xmm0
movsd  0xe8(%rcx),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0xe8(%rax)
mov    %rbp,%rdi
mov    %r12,%rsi
xor    %eax,%eax
call   1050 <printf@plt>
mov    -0x8(%r14,%rbx,8),%rax
movsd  0x28(%rax),%xmm0
movsd  0x83c(%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  0x818(%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  0x7f4(%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  0x7cd(%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  0x7a6(%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  0x77f(%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  0x758(%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"MGFAS 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 @mgfas(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 = %188, %70
%73 = load i32, i32* %2
%74 = icmp slt i32 %73, 33
br i1 %74, label %75, label %191
75: ; preds = %72
store i32 5, i32* %3
br label %76
76: ; preds = %162, %75
%77 = load i32, i32* %3
%78 = icmp slt i32 %77, 33
br i1 %78, label %79, label %165
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*** %6
%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
%141 = load double, double* %140
%142 = load double**, double*** %6
%143 = load i32, i32* %2
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double*, double** %142, i64 %144
%146 = load double*, double** %145
%147 = load i32, i32* %3
%148 = sext i32 %147 to i64
%149 = getelementptr inbounds double, double* %146, i64 %148
%150 = load double, double* %149
%151 = fmul double %141, %150
%152 = fdiv double %151, 1.024000e+03
%153 = fadd double %132, %152
%154 = load double**, double*** %5
%155 = load i32, i32* %2
%156 = sext i32 %155 to i64
%157 = getelementptr inbounds double*, double** %154, i64 %156
%158 = load double*, double** %157
%159 = load i32, i32* %3
%160 = sext i32 %159 to i64
%161 = getelementptr inbounds double, double* %158, i64 %160
store double %153, double* %161
br label %162
162: ; preds = %79
%163 = load i32, i32* %3
%164 = add nsw i32 %163, 4
store i32 %164, i32* %3
br label %76
165: ; preds = %76
%166 = 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 %167
167: ; preds = %183, %165
%168 = load i32, i32* %3
%169 = icmp slt i32 %168, 33
br i1 %169, label %170, label %186
170: ; preds = %167
%171 = load double**, double*** %5
%172 = load i32, i32* %2
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double*, double** %171, i64 %173
%175 = load double*, double** %174
%176 = load i32, i32* %3
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds double, double* %175, i64 %177
%179 = load double, double* %178
%180 = fmul double %179, 3.200000e+01
%181 = fmul double %180, 3.200000e+01
%182 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %181)
br label %183
183: ; preds = %170
%184 = load i32, i32* %3
%185 = add nsw i32 %184, 4
store i32 %185, i32* %3
br label %167
186: ; preds = %167
%187 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.2, i64 0, i64 0))
br label %188
188: ; preds = %186
%189 = load i32, i32* %2
%190 = add nsw i32 %189, 4
store i32 %190, i32* %2
br label %72
191: ; preds = %72
%192 = load double**, double*** %6
%193 = call i32 @free_dmatrix(double** %192, i32 1, i32 33, i32 1, i32 33)
%194 = load double**, double*** %5
%195 = call i32 @free_dmatrix(double** %194, 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"MGFAS 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 @mgfas(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 ], [ %364, %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 = fmul double %206, %206
%209 = fmul double %208, 0x3F50000000000000
%210 = fadd double %207, %209
%211 = getelementptr inbounds double, double* %193, i64 5
store double %210, double* %211
%212 = getelementptr inbounds double, double* %186, i64 9
%213 = load double, double* %212
%214 = getelementptr inbounds double, double* %189, i64 9
%215 = load double, double* %214
%216 = fadd double %213, %215
%217 = getelementptr inbounds double, double* %191, i64 10
%218 = load double, double* %217
%219 = fadd double %216, %218
%220 = getelementptr inbounds double, double* %191, i64 8
%221 = load double, double* %220
%222 = fadd double %219, %221
%223 = getelementptr inbounds double, double* %191, i64 9
%224 = load double, double* %223
%225 = tail call double @llvm.fmuladd.f64(double %224, double -4.000000e+00, double %222)
%226 = fmul double %224, %224
%227 = fmul double %226, 0x3F50000000000000
%228 = fadd double %225, %227
%229 = getelementptr inbounds double, double* %193, i64 9
store double %228, double* %229
%230 = getelementptr inbounds double, double* %186, i64 13
%231 = load double, double* %230
%232 = getelementptr inbounds double, double* %189, i64 13
%233 = load double, double* %232
%234 = fadd double %231, %233
%235 = getelementptr inbounds double, double* %191, i64 14
%236 = load double, double* %235
%237 = fadd double %234, %236
%238 = getelementptr inbounds double, double* %191, i64 12
%239 = load double, double* %238
%240 = fadd double %237, %239
%241 = getelementptr inbounds double, double* %191, i64 13
%242 = load double, double* %241
%243 = tail call double @llvm.fmuladd.f64(double %242, double -4.000000e+00, double %240)
%244 = fmul double %242, %242
%245 = fmul double %244, 0x3F50000000000000
%246 = fadd double %243, %245
%247 = getelementptr inbounds double, double* %193, i64 13
store double %246, double* %247
%248 = getelementptr inbounds double, double* %186, i64 17
%249 = load double, double* %248
%250 = getelementptr inbounds double, double* %189, i64 17
%251 = load double, double* %250
%252 = fadd double %249, %251
%253 = getelementptr inbounds double, double* %191, i64 18
%254 = load double, double* %253
%255 = fadd double %252, %254
%256 = getelementptr inbounds double, double* %191, i64 16
%257 = load double, double* %256
%258 = fadd double %255, %257
%259 = getelementptr inbounds double, double* %191, i64 17
%260 = load double, double* %259
%261 = tail call double @llvm.fmuladd.f64(double %260, double -4.000000e+00, double %258)
%262 = fmul double %260, %260
%263 = fmul double %262, 0x3F50000000000000
%264 = fadd double %261, %263
%265 = getelementptr inbounds double, double* %193, i64 17
store double %264, double* %265
%266 = getelementptr inbounds double, double* %186, i64 21
%267 = load double, double* %266
%268 = getelementptr inbounds double, double* %189, i64 21
%269 = load double, double* %268
%270 = fadd double %267, %269
%271 = getelementptr inbounds double, double* %191, i64 22
%272 = load double, double* %271
%273 = fadd double %270, %272
%274 = getelementptr inbounds double, double* %191, i64 20
%275 = load double, double* %274
%276 = fadd double %273, %275
%277 = getelementptr inbounds double, double* %191, i64 21
%278 = load double, double* %277
%279 = tail call double @llvm.fmuladd.f64(double %278, double -4.000000e+00, double %276)
%280 = fmul double %278, %278
%281 = fmul double %280, 0x3F50000000000000
%282 = fadd double %279, %281
%283 = getelementptr inbounds double, double* %193, i64 21
store double %282, double* %283
%284 = getelementptr inbounds double, double* %186, i64 25
%285 = load double, double* %284
%286 = getelementptr inbounds double, double* %189, i64 25
%287 = load double, double* %286
%288 = fadd double %285, %287
%289 = getelementptr inbounds double, double* %191, i64 26
%290 = load double, double* %289
%291 = fadd double %288, %290
%292 = getelementptr inbounds double, double* %191, i64 24
%293 = load double, double* %292
%294 = fadd double %291, %293
%295 = getelementptr inbounds double, double* %191, i64 25
%296 = load double, double* %295
%297 = tail call double @llvm.fmuladd.f64(double %296, double -4.000000e+00, double %294)
%298 = fmul double %296, %296
%299 = fmul double %298, 0x3F50000000000000
%300 = fadd double %297, %299
%301 = getelementptr inbounds double, double* %193, i64 25
store double %300, double* %301
%302 = getelementptr inbounds double, double* %186, i64 29
%303 = load double, double* %302
%304 = getelementptr inbounds double, double* %189, i64 29
%305 = load double, double* %304
%306 = fadd double %303, %305
%307 = getelementptr inbounds double, double* %191, i64 30
%308 = load double, double* %307
%309 = fadd double %306, %308
%310 = getelementptr inbounds double, double* %191, i64 28
%311 = load double, double* %310
%312 = fadd double %309, %311
%313 = getelementptr inbounds double, double* %191, i64 29
%314 = load double, double* %313
%315 = tail call double @llvm.fmuladd.f64(double %314, double -4.000000e+00, double %312)
%316 = fmul double %314, %314
%317 = fmul double %316, 0x3F50000000000000
%318 = fadd double %315, %317
%319 = getelementptr inbounds double, double* %193, i64 29
store double %318, double* %319
%320 = 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))
%321 = load double*, double** %192
%322 = getelementptr inbounds double, double* %321, i64 5
%323 = load double, double* %322
%324 = fmul double %323, 3.200000e+01
%325 = fmul double %324, 3.200000e+01
%326 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %325)
%327 = load double*, double** %192
%328 = getelementptr inbounds double, double* %327, i64 9
%329 = load double, double* %328
%330 = fmul double %329, 3.200000e+01
%331 = fmul double %330, 3.200000e+01
%332 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %331)
%333 = load double*, double** %192
%334 = getelementptr inbounds double, double* %333, i64 13
%335 = load double, double* %334
%336 = fmul double %335, 3.200000e+01
%337 = fmul double %336, 3.200000e+01
%338 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %337)
%339 = load double*, double** %192
%340 = getelementptr inbounds double, double* %339, i64 17
%341 = load double, double* %340
%342 = fmul double %341, 3.200000e+01
%343 = fmul double %342, 3.200000e+01
%344 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %343)
%345 = load double*, double** %192
%346 = getelementptr inbounds double, double* %345, i64 21
%347 = load double, double* %346
%348 = fmul double %347, 3.200000e+01
%349 = fmul double %348, 3.200000e+01
%350 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %349)
%351 = load double*, double** %192
%352 = getelementptr inbounds double, double* %351, i64 25
%353 = load double, double* %352
%354 = fmul double %353, 3.200000e+01
%355 = fmul double %354, 3.200000e+01
%356 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %355)
%357 = load double*, double** %192
%358 = getelementptr inbounds double, double* %357, i64 29
%359 = load double, double* %358
%360 = fmul double %359, 3.200000e+01
%361 = fmul double %360, 3.200000e+01
%362 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), double %361)
%363 = tail call i32 @putchar(i32 10)
%364 = add nuw nsw i64 %183, 4
%365 = icmp ult i64 %183, 29
br i1 %365, label %182, label %366
366: ; preds = %182
%367 = tail call i32 @free_dmatrix(double** %2, i32 1, i32 33, i32 1, i32 33)
%368 = tail call i32 @free_dmatrix(double** %1, i32 1, i32 33, i32 1, i32 33)
ret i32 0
}