#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 */
int /*<<< orphan*/ MRot(double **, int, int, int, double);
int /*<<< orphan*/ MulMatrix(double **, double **, double **, int, int, int);
double **dmatrix(int, int, int, int);
int /*<<< orphan*/ free_dmatrix(double **, int, int, int, int);
double uniformrng();
void CreateRotationMatrix(double **MResult, int N) {
  int i;
  int k1;
  int k2;
  double **MTurn;
  double **MWork;
  double alpha;
  double myPI;
  myPI = 3.141592653589793;
  MTurn = dmatrix(1, N, 1, N);
  MWork = dmatrix(1, N, 1, N);
  for (k1 = 1; k1 <= N; k1++)
    for (k2 = 1; k2 <= N; k2++)
      MResult[k1][k2] = ((k1 == k2) ? 1.0 : 0.0);
  for (i = 2; i <= N; i++) {
    alpha = (uniformrng() - 0.5) * myPI * 0.5;
    MRot(MTurn, 1, i, N, alpha);
    for (k1 = 1; k1 <= N; k1++)
      for (k2 = 1; k2 <= N; k2++)
        MWork[k1][k2] = MResult[k1][k2];
    MulMatrix(MWork, MTurn, MResult, N, N, N);
  }
  for (i = 2; i <= N - 1; i++) {
    alpha = (uniformrng() - 0.5) * myPI * 0.5;
    MRot(MTurn, i, N, N, alpha);
    for (k1 = 1; k1 <= N; k1++)
      for (k2 = 1; k2 <= N; k2++)
        MWork[k1][k2] = MResult[k1][k2];
    MulMatrix(MWork, MTurn, MResult, N, N, N);
  }
  free_dmatrix(MTurn, 1, N, 1, N);
  free_dmatrix(MWork, 1, N, 1, N);
}


x86 assembly (clang, -O0):
<CreateRotationMatrix>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  0xe99(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0xc(%rbp),%esi
mov    -0xc(%rbp),%ecx
mov    $0x1,%edx
mov    %edx,%edi
call   1050 <dmatrix@plt>
mov    %rax,-0x20(%rbp)
mov    -0xc(%rbp),%esi
mov    -0xc(%rbp),%ecx
mov    $0x1,%edx
mov    %edx,%edi
call   1050 <dmatrix@plt>
mov    %rax,-0x28(%rbp)
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1226 <CreateRotationMatrix+0xd6>
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1213 <CreateRotationMatrix+0xc3>
mov    -0x14(%rbp),%eax
mov    -0x18(%rbp),%ecx
movsd  0xe44(%rip),%xmm0
xorps  %xmm1,%xmm1
movsd  %xmm1,-0x48(%rbp)
cmp    %ecx,%eax
movsd  %xmm0,-0x40(%rbp)
je     11eb <CreateRotationMatrix+0x9b>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11b2 <CreateRotationMatrix+0x62>
jmp    1218 <CreateRotationMatrix+0xc8>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    119f <CreateRotationMatrix+0x4f>
movl   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1316 <CreateRotationMatrix+0x1c6>
mov    $0x0,%al
call   1030 <uniformrng@plt>
movsd  0xdc0(%rip),%xmm1
subsd  %xmm1,%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  0xdaf(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x10(%rbp),%edx
mov    -0xc(%rbp),%ecx
movsd  -0x30(%rbp),%xmm0
mov    $0x1,%esi
call   1060 <MRot@plt>
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12ec <CreateRotationMatrix+0x19c>
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12d9 <CreateRotationMatrix+0x189>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1295 <CreateRotationMatrix+0x145>
jmp    12de <CreateRotationMatrix+0x18e>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1282 <CreateRotationMatrix+0x132>
mov    -0x28(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x8(%rbp),%rdx
mov    -0xc(%rbp),%ecx
mov    -0xc(%rbp),%r8d
mov    -0xc(%rbp),%r9d
call   1040 <MulMatrix@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    122d <CreateRotationMatrix+0xdd>
movl   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1409 <CreateRotationMatrix+0x2b9>
mov    $0x0,%al
call   1030 <uniformrng@plt>
movsd  0xccb(%rip),%xmm1
subsd  %xmm1,%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  0xcba(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x10(%rbp),%esi
mov    -0xc(%rbp),%edx
mov    -0xc(%rbp),%ecx
movsd  -0x30(%rbp),%xmm0
call   1060 <MRot@plt>
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     13df <CreateRotationMatrix+0x28f>
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     13cc <CreateRotationMatrix+0x27c>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1388 <CreateRotationMatrix+0x238>
jmp    13d1 <CreateRotationMatrix+0x281>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1375 <CreateRotationMatrix+0x225>
mov    -0x28(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x8(%rbp),%rdx
mov    -0xc(%rbp),%ecx
mov    -0xc(%rbp),%r8d
mov    -0xc(%rbp),%r9d
call   1040 <MulMatrix@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    131d <CreateRotationMatrix+0x1cd>
mov    -0x20(%rbp),%rdi
mov    -0xc(%rbp),%edx
mov    -0xc(%rbp),%r8d
mov    $0x1,%ecx
mov    %ecx,%esi
call   1070 <free_dmatrix@plt>
mov    -0x28(%rbp),%rdi
mov    -0xc(%rbp),%edx
mov    -0xc(%rbp),%r8d
mov    $0x1,%ecx
mov    %ecx,%esi
call   1070 <free_dmatrix@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<CreateRotationMatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %esi,%ebx
mov    %rdi,%r13
mov    $0x1,%edi
mov    $0x1,%edx
mov    %esi,%ecx
call   1050 <dmatrix@plt>
mov    %rax,0x8(%rsp)
mov    $0x1,%edi
mov    %ebx,%esi
mov    $0x1,%edx
mov    %ebx,%ecx
call   1050 <dmatrix@plt>
mov    %rbx,%rdx
mov    %rax,%r12
test   %edx,%edx
jle    179c <CreateRotationMatrix+0x64c>
mov    %edx,%r9d
lea    0x1(%r9),%r15
lea    -0x1(%r9),%r8
mov    %r9d,%r10d
and    $0x3,%r10d
and    $0xfffffffc,%r9d
mov    $0x1,%eax
xor    %edi,%edi
movsd  0xe41(%rip),%xmm0
mov    %rdx,0x18(%rsp)
mov    %rdx,%rbx
jmp    11e1 <CreateRotationMatrix+0x91>
nopl   0x0(%rax)
add    $0x1,%rax
add    $0x1,%rdi
cmp    %r15,%rax
je     12c2 <CreateRotationMatrix+0x172>
mov    0x0(%r13,%rax,8),%rbp
mov    $0x1,%ecx
cmp    $0x3,%r8
jb     1279 <CreateRotationMatrix+0x129>
mov    $0x1,%ebx
xor    %ecx,%ecx
jmp    1213 <CreateRotationMatrix+0xc3>
xchg   %ax,%ax
movsd  %xmm1,0x20(%rbp,%rcx,8)
add    $0x4,%rbx
add    $0x4,%rcx
cmp    %r9,%rcx
je     1270 <CreateRotationMatrix+0x120>
movapd %xmm0,%xmm1
cmp    %rcx,%rdi
je     1220 <CreateRotationMatrix+0xd0>
xorpd  %xmm1,%xmm1
movsd  %xmm1,0x8(%rbp,%rcx,8)
lea    0x1(%rbx),%rsi
movapd %xmm0,%xmm1
cmp    %rsi,%rax
je     1237 <CreateRotationMatrix+0xe7>
xorpd  %xmm1,%xmm1
movsd  %xmm1,0x10(%rbp,%rcx,8)
lea    0x2(%rbx),%rsi
movapd %xmm0,%xmm1
cmp    %rsi,%rax
je     124e <CreateRotationMatrix+0xfe>
xorpd  %xmm1,%xmm1
movsd  %xmm1,0x18(%rbp,%rcx,8)
lea    0x3(%rbx),%rsi
movapd %xmm0,%xmm1
cmp    %rsi,%rax
je     1200 <CreateRotationMatrix+0xb0>
xorpd  %xmm1,%xmm1
jmp    1200 <CreateRotationMatrix+0xb0>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rcx
mov    0x18(%rsp),%rbx
test   %r10,%r10
je     11d0 <CreateRotationMatrix+0x80>
lea    0x0(%rbp,%rcx,8),%rbp
mov    %rax,%rsi
sub    %rcx,%rsi
xor    %ecx,%ecx
jmp    12b3 <CreateRotationMatrix+0x163>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  %xmm1,0x0(%rbp,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r10
je     11d0 <CreateRotationMatrix+0x80>
movapd %xmm0,%xmm1
cmp    %rcx,%rsi
je     12a0 <CreateRotationMatrix+0x150>
xorpd  %xmm1,%xmm1
jmp    12a0 <CreateRotationMatrix+0x150>
cmp    $0x2,%ebx
mov    %rbx,%rdx
jl     179c <CreateRotationMatrix+0x64c>
movabs $0x7ffffffffffffffc,%rax
ff ff 7f
mov    %rax,0x38(%rsp)
lea    0x1(%rdx),%r14d
lea    -0x1(%r15),%rbp
lea    -0x5(%r15),%rax
mov    %rax,0x40(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rbp,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x30(%rsp)
add    $0x1,%rax
mov    %rax,0x28(%rsp)
mov    %ecx,%eax
and    $0x3,%eax
mov    %rax,0x50(%rsp)
mov    %rcx,0x48(%rsp)
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x20(%rsp)
mov    %rdx,%rbx
mov    $0x2,%eax
jmp    135d <CreateRotationMatrix+0x20d>
mov    %r12,%rdi
mov    0x8(%rsp),%rsi
mov    %r13,%rdx
mov    0x18(%rsp),%rbx
mov    %ebx,%ecx
mov    %ebx,%r8d
mov    %ebx,%r9d
call   1040 <MulMatrix@plt>
mov    0x10(%rsp),%rcx
lea    0x1(%rcx),%eax
cmp    %ebx,%ecx
je     154c <CreateRotationMatrix+0x3fc>
mov    %rax,0x10(%rsp)
xor    %eax,%eax
call   1030 <uniformrng@plt>
addsd  0xc97(%rip),%xmm0
mulsd  0xc97(%rip),%xmm0
mulsd  0xc97(%rip),%xmm0
mov    0x8(%rsp),%rdi
mov    $0x1,%esi
mov    0x10(%rsp),%rdx
mov    %ebx,%ecx
call   1060 <MRot@plt>
mov    $0x1,%eax
mov    0x40(%rsp),%r8
mov    0x30(%rsp),%r9
mov    0x28(%rsp),%r10
mov    0x50(%rsp),%r11
jmp    13cd <CreateRotationMatrix+0x27d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rax
cmp    %r14,%rax
je     1330 <CreateRotationMatrix+0x1e0>
mov    0x0(%r13,%rax,8),%rcx
mov    (%r12,%rax,8),%rdx
mov    $0x1,%esi
cmp    $0x4,%rbp
jb     14d0 <CreateRotationMatrix+0x380>
lea    0x8(%rdx),%rsi
lea    (%rcx,%r15,8),%rdi
cmp    %rdi,%rsi
jae    1408 <CreateRotationMatrix+0x2b8>
lea    (%rdx,%r15,8),%rdi
lea    0x8(%rcx),%rbx
mov    $0x1,%esi
cmp    %rdi,%rbx
jb     14d0 <CreateRotationMatrix+0x380>
cmp    $0xc,%r8
jae    1412 <CreateRotationMatrix+0x2c2>
xor    %esi,%esi
jmp    147e <CreateRotationMatrix+0x32e>
mov    0x20(%rsp),%rdi
xor    %esi,%esi
nopl   0x0(%rax)
movups 0x8(%rcx,%rsi,8),%xmm0
movups 0x18(%rcx,%rsi,8),%xmm1
movups %xmm0,0x8(%rdx,%rsi,8)
movups %xmm1,0x18(%rdx,%rsi,8)
movups 0x28(%rcx,%rsi,8),%xmm0
movups 0x38(%rcx,%rsi,8),%xmm1
movups %xmm0,0x28(%rdx,%rsi,8)
movups %xmm1,0x38(%rdx,%rsi,8)
movups 0x48(%rcx,%rsi,8),%xmm0
movups 0x58(%rcx,%rsi,8),%xmm1
movups %xmm0,0x48(%rdx,%rsi,8)
movups %xmm1,0x58(%rdx,%rsi,8)
movupd 0x68(%rcx,%rsi,8),%xmm0
movupd 0x78(%rcx,%rsi,8),%xmm1
movupd %xmm0,0x68(%rdx,%rsi,8)
movupd %xmm1,0x78(%rdx,%rsi,8)
add    $0x10,%rsi
add    $0xfffffffffffffffc,%rdi
jne    1420 <CreateRotationMatrix+0x2d0>
test   %r11,%r11
je     14bc <CreateRotationMatrix+0x36c>
mov    %r11,%rdi
cs nopw 0x0(%rax,%rax,1)
lea    0x0(,%rsi,8),%rbx
or     $0x8,%rbx
movupd (%rcx,%rbx,1),%xmm0
movupd 0x10(%rcx,%rbx,1),%xmm1
movupd %xmm0,(%rdx,%rbx,1)
movupd %xmm1,0x10(%rdx,%rbx,1)
add    $0x4,%rsi
add    $0xffffffffffffffff,%rdi
jne    1490 <CreateRotationMatrix+0x340>
mov    %r10,%rsi
cmp    %r9,%rbp
je     13c0 <CreateRotationMatrix+0x270>
nopl   0x0(%rax,%rax,1)
mov    %r15d,%ebx
sub    %esi,%ebx
mov    %rsi,%rdi
not    %rdi
add    %r15,%rdi
and    $0x3,%rbx
je     1504 <CreateRotationMatrix+0x3b4>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x1,%rsi
add    $0xffffffffffffffff,%rbx
jne    14f0 <CreateRotationMatrix+0x3a0>
cmp    $0x3,%rdi
jb     13c0 <CreateRotationMatrix+0x270>
xchg   %ax,%ax
movsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
movsd  0x8(%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rsi,8)
movsd  0x10(%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rsi,8)
movsd  0x18(%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x18(%rdx,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%r15
jne    1510 <CreateRotationMatrix+0x3c0>
jmp    13c0 <CreateRotationMatrix+0x270>
cmp    $0x3,%ebx
mov    %rbx,%rdx
jl     179c <CreateRotationMatrix+0x64c>
mov    0x30(%rsp),%rax
add    $0x1,%rax
mov    %rax,0x28(%rsp)
mov    0x48(%rsp),%rax
mov    %eax,%ecx
and    $0x3,%ecx
mov    %rcx,0x20(%rsp)
and    %rax,0x38(%rsp)
mov    $0x2,%ebx
jmp    15bf <CreateRotationMatrix+0x46f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r12,%rdi
mov    0x8(%rsp),%rsi
mov    %r13,%rdx
mov    0x18(%rsp),%rbx
mov    %ebx,%ecx
mov    %ebx,%r8d
mov    %ebx,%r9d
call   1040 <MulMatrix@plt>
mov    %rbx,%rdx
mov    0x10(%rsp),%ebx
add    $0x1,%ebx
cmp    %edx,%ebx
je     179c <CreateRotationMatrix+0x64c>
xor    %eax,%eax
call   1030 <uniformrng@plt>
addsd  0xa3a(%rip),%xmm0
mulsd  0xa3a(%rip),%xmm0
mulsd  0xa3a(%rip),%xmm0
mov    0x8(%rsp),%rdi
mov    %ebx,0x10(%rsp)
mov    %ebx,%esi
mov    0x18(%rsp),%rcx
mov    %ecx,%edx
call   1060 <MRot@plt>
mov    $0x1,%eax
mov    0x38(%rsp),%r8
mov    0x40(%rsp),%r9
mov    0x30(%rsp),%r10
mov    0x28(%rsp),%r11
jmp    161d <CreateRotationMatrix+0x4cd>
add    $0x1,%rax
cmp    %r14,%rax
je     1590 <CreateRotationMatrix+0x440>
mov    0x0(%r13,%rax,8),%rcx
mov    (%r12,%rax,8),%rdx
mov    $0x1,%esi
cmp    $0x4,%rbp
jb     1720 <CreateRotationMatrix+0x5d0>
lea    0x8(%rdx),%rsi
lea    (%rcx,%r15,8),%rdi
cmp    %rdi,%rsi
jae    1658 <CreateRotationMatrix+0x508>
lea    (%rdx,%r15,8),%rdi
lea    0x8(%rcx),%rbx
mov    $0x1,%esi
cmp    %rdi,%rbx
jb     1720 <CreateRotationMatrix+0x5d0>
cmp    $0xc,%r9
jae    1662 <CreateRotationMatrix+0x512>
xor    %esi,%esi
jmp    16ce <CreateRotationMatrix+0x57e>
mov    %r8,%rdi
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movups 0x8(%rcx,%rsi,8),%xmm0
movups 0x18(%rcx,%rsi,8),%xmm1
movups %xmm0,0x8(%rdx,%rsi,8)
movups %xmm1,0x18(%rdx,%rsi,8)
movups 0x28(%rcx,%rsi,8),%xmm0
movups 0x38(%rcx,%rsi,8),%xmm1
movups %xmm0,0x28(%rdx,%rsi,8)
movups %xmm1,0x38(%rdx,%rsi,8)
movups 0x48(%rcx,%rsi,8),%xmm0
movups 0x58(%rcx,%rsi,8),%xmm1
movups %xmm0,0x48(%rdx,%rsi,8)
movups %xmm1,0x58(%rdx,%rsi,8)
movupd 0x68(%rcx,%rsi,8),%xmm0
movupd 0x78(%rcx,%rsi,8),%xmm1
movupd %xmm0,0x68(%rdx,%rsi,8)
movupd %xmm1,0x78(%rdx,%rsi,8)
add    $0x10,%rsi
add    $0xfffffffffffffffc,%rdi
jne    1670 <CreateRotationMatrix+0x520>
cmpq   $0x0,0x20(%rsp)
je     170c <CreateRotationMatrix+0x5bc>
mov    0x20(%rsp),%rdi
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rsi,8),%rbx
or     $0x8,%rbx
movupd (%rcx,%rbx,1),%xmm0
movupd 0x10(%rcx,%rbx,1),%xmm1
movupd %xmm0,(%rdx,%rbx,1)
movupd %xmm1,0x10(%rdx,%rbx,1)
add    $0x4,%rsi
add    $0xffffffffffffffff,%rdi
jne    16e0 <CreateRotationMatrix+0x590>
mov    %r11,%rsi
cmp    %r10,%rbp
je     1610 <CreateRotationMatrix+0x4c0>
nopl   0x0(%rax,%rax,1)
mov    %r14d,%ebx
sub    %esi,%ebx
mov    %rsi,%rdi
not    %rdi
add    %r14,%rdi
and    $0x3,%rbx
je     1754 <CreateRotationMatrix+0x604>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x1,%rsi
add    $0xffffffffffffffff,%rbx
jne    1740 <CreateRotationMatrix+0x5f0>
cmp    $0x3,%rdi
jb     1610 <CreateRotationMatrix+0x4c0>
xchg   %ax,%ax
movsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
movsd  0x8(%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rsi,8)
movsd  0x10(%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rsi,8)
movsd  0x18(%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x18(%rdx,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%r15
jne    1760 <CreateRotationMatrix+0x610>
jmp    1610 <CreateRotationMatrix+0x4c0>
mov    0x8(%rsp),%rdi
mov    $0x1,%esi
mov    $0x1,%ecx
mov    %edx,%r8d
mov    %rdx,%rbx
call   1070 <free_dmatrix@plt>
mov    %r12,%rdi
mov    $0x1,%esi
mov    %ebx,%edx
mov    $0x1,%ecx
mov    %ebx,%r8d
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1070 <free_dmatrix@plt>

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

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