#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*/ free(double *);
double *malloc(int);
int omp_get_max_threads();
int omp_get_thread_num();
int sqrt(double);
double RB_GS_OmpLoopV3err(double **black, double **red, int m, int n,
                          double eps, int maxit, int iterations_print,
                          int *iterations) {
  int i, j;
  double v;
  double *my_error;
  int tid, nthreads;
  int my_iterations = *iterations;
  double error = 10.0 * eps;
  nthreads = omp_get_max_threads();
  my_error = malloc(sizeof(double) * nthreads);
#pragma omp parallel private(j, v, tid) firstprivate(my_iterations)
  {
    tid = omp_get_thread_num();
    while (error >= eps && my_iterations < maxit) {
      my_error[tid] = 0.0;
#pragma omp for schedule(OMP_SCHED)
      for (i = 1; i < m - 1; i++)
        for (j = i % 2; j < n / 2 - (i + 1) % 2; j++) {
          v = 0.25 * (red[i - 1][j] + red[i + 1][j] + red[i][j - i % 2] +
                      red[i][j + (i + 1) % 2]);
          my_error[tid] += (v - black[i][j]) * (v - black[i][j]);
          black[i][j] = v;
        }
#pragma omp for schedule(OMP_SCHED)
      for (i = 1; i < m - 1; i++)
        for (j = (i + 1) % 2; j < n / 2 - i % 2; j++) {
          v = 0.25 * (black[i - 1][j] + black[i + 1][j] +
                      black[i][j - (i + 1) % 2] + black[i][j + i % 2]);
          my_error[tid] += (v - red[i][j]) * (v - red[i][j]);
          red[i][j] = v;
        }
#pragma omp single
      {
        for (j = 0; j < nthreads; j++)
          error += my_error[j];
        error = sqrt(error) / (n * m);
      }
      my_iterations++;
    }
#pragma omp master
    *iterations = my_iterations;
  }
  free(my_error);
  return error;
}


x86 assembly (clang, -O0):
<RB_GS_OmpLoopV3err>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
mov    0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x4c(%rbp)
movsd  0xe75(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    $0x0,%al
call   1060 <omp_get_max_threads@plt>
mov    %eax,-0x48(%rbp)
movslq -0x48(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
mov    $0x0,%al
call   1040 <omp_get_thread_num@plt>
mov    %eax,-0x44(%rbp)
movsd  -0x58(%rbp),%xmm0
xor    %eax,%eax
ucomisd -0x20(%rbp),%xmm0
mov    %al,-0x59(%rbp)
jb     11dd <RB_GS_OmpLoopV3err+0x8d>
mov    -0x4c(%rbp),%eax
cmp    -0x24(%rbp),%eax
setl   %al
mov    %al,-0x59(%rbp)
mov    -0x59(%rbp),%al
test   $0x1,%al
jne    11ed <RB_GS_OmpLoopV3err+0x9d>
jmp    1560 <RB_GS_OmpLoopV3err+0x410>
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    137a <RB_GS_OmpLoopV3err+0x22a>
mov    -0x2c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %edx,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x60(%rbp),%eax
sub    %edx,%ecx
cmp    %ecx,%eax
jge    1367 <RB_GS_OmpLoopV3err+0x217>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x70(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x70(%rbp),%rax
sub    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x68(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x68(%rbp),%rax
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
movsd  0xd20(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1223 <RB_GS_OmpLoopV3err+0xd3>
jmp    136c <RB_GS_OmpLoopV3err+0x21c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1204 <RB_GS_OmpLoopV3err+0xb4>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    14fd <RB_GS_OmpLoopV3err+0x3ad>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %edx,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x2c(%rbp),%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x74(%rbp),%eax
sub    %edx,%ecx
cmp    %ecx,%eax
jge    14ea <RB_GS_OmpLoopV3err+0x39a>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x88(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x88(%rbp),%rax
sub    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x80(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    -0x80(%rbp),%rax
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
movsd  0xb9d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13a3 <RB_GS_OmpLoopV3err+0x253>
jmp    14ef <RB_GS_OmpLoopV3err+0x39f>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1381 <RB_GS_OmpLoopV3err+0x231>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    1535 <RB_GS_OmpLoopV3err+0x3e5>
mov    -0x40(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1504 <RB_GS_OmpLoopV3err+0x3b4>
movsd  -0x58(%rbp),%xmm0
call   1070 <sqrt@plt>
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11bc <RB_GS_OmpLoopV3err+0x6c>
mov    -0x4c(%rbp),%ecx
mov    0x10(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x58(%rbp),%xmm0
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<RB_GS_OmpLoopV3err>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8d,%ebp
movsd  %xmm0,0x30(%rsp)
mov    %ecx,%r15d
mov    %rdx,0x10(%rsp)
mov    %rsi,%r14
mov    %rdi,0x50(%rsp)
mov    0x90(%rsp),%rax
mov    (%rax),%ebx
movsd  0xe77(%rip),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
xor    %eax,%eax
call   1060 <omp_get_max_threads@plt>
mov    %rax,0x28(%rsp)
lea    0x0(,%rax,8),%edi
call   1050 <malloc@plt>
mov    %rax,%r13
xor    %eax,%eax
call   1040 <omp_get_thread_num@plt>
movsd  0x8(%rsp),%xmm0
ucomisd 0x30(%rsp),%xmm0
jb     14a2 <RB_GS_OmpLoopV3err+0x352>
cmp    %ebp,%ebx
mov    %ebx,%ecx
jge    14a4 <RB_GS_OmpLoopV3err+0x354>
mov    %ebp,0x20(%rsp)
cltq
mov    %rax,0x18(%rsp)
mov    0x10(%rsp),%rdx
lea    -0x1(%rdx),%r12d
mov    %r15d,%eax
shr    $0x1f,%r15d
add    %eax,%r15d
sar    %r15d
imul   %edx,%eax
mov    %eax,0x24(%rsp)
mov    0x28(%rsp),%eax
lea    -0x1(%rax),%rdx
mov    %rdx,0x48(%rsp)
mov    %eax,%edx
and    $0x7,%edx
mov    %rdx,0x40(%rsp)
and    $0xfffffff8,%eax
mov    %rax,0x38(%rsp)
movsd  0xde8(%rip),%xmm4
mov    %ecx,0x8(%rsp)
mov    0x18(%rsp),%rax
movq   $0x0,0x0(%r13,%rax,8)
cmpl   $0x3,0x10(%rsp)
mov    0x50(%rsp),%r11
jl     13d0 <RB_GS_OmpLoopV3err+0x280>
mov    $0x1,%r10d
mov    $0x1,%r8d
xorpd  %xmm3,%xmm3
xor    %r9d,%r9d
jmp    1271 <RB_GS_OmpLoopV3err+0x121>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r10d
add    $0x1,%r9b
cmp    %r12,%r8
je     1310 <RB_GS_OmpLoopV3err+0x1c0>
mov    %r8,%rcx
mov    %ecx,%edx
and    $0x1,%edx
add    $0x1,%r8
mov    %r8d,%esi
and    $0x1,%esi
mov    %r15d,%eax
sub    %esi,%eax
cmp    %eax,%edx
jge    1260 <RB_GS_OmpLoopV3err+0x110>
mov    %r10d,%edx
and    $0x1,%edx
movzbl %r9b,%esi
and    $0x1,%esi
add    %rdx,%rsi
mov    %r10d,%edi
and    $0x1,%edi
mov    -0x8(%r14,%rcx,8),%rbp
mov    (%r14,%rcx,8),%rbx
mov    (%r14,%r8,8),%rdx
mov    (%r11,%rcx,8),%rcx
cltq
cs nopw 0x0(%rax,%rax,1)
nop
movsd  0x0(%rbp,%rdi,8),%xmm1
addsd  (%rdx,%rdi,8),%xmm1
addsd  (%rbx),%xmm1
addsd  (%rbx,%rsi,8),%xmm1
mulsd  %xmm4,%xmm1
movapd %xmm1,%xmm2
subsd  (%rcx,%rdi,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm2,%xmm3
movsd  %xmm1,(%rcx,%rdi,8)
add    $0x1,%rdi
add    $0x8,%rbx
cmp    %rax,%rdi
jl     12c0 <RB_GS_OmpLoopV3err+0x170>
mov    0x18(%rsp),%rax
movsd  %xmm3,0x0(%r13,%rax,8)
jmp    1260 <RB_GS_OmpLoopV3err+0x110>
nopl   0x0(%rax)
cmpl   $0x3,0x10(%rsp)
jl     13d0 <RB_GS_OmpLoopV3err+0x280>
mov    $0x1,%r8d
xor    %r9d,%r9d
jmp    133d <RB_GS_OmpLoopV3err+0x1ed>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r9b
cmp    %r12,%r8
je     13d0 <RB_GS_OmpLoopV3err+0x280>
mov    %r8,%rcx
mov    %ecx,%edi
and    $0x1,%edi
add    $0x1,%r8
mov    %r8d,%edx
and    $0x1,%edx
mov    %r15d,%eax
sub    %edi,%eax
cmp    %eax,%edx
jge    1330 <RB_GS_OmpLoopV3err+0x1e0>
movzbl %r9b,%edx
mov    %edx,%esi
and    $0x1,%esi
add    %rdi,%rsi
and    $0x1,%edx
mov    -0x8(%r11,%rcx,8),%rdi
mov    (%r11,%rcx,8),%rbp
mov    (%r11,%r8,8),%rbx
mov    (%r14,%rcx,8),%rcx
cltq
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rdx,8),%xmm1
addsd  (%rbx,%rdx,8),%xmm1
addsd  0x0(%rbp),%xmm1
addsd  0x0(%rbp,%rsi,8),%xmm1
mulsd  %xmm4,%xmm1
movapd %xmm1,%xmm2
subsd  (%rcx,%rdx,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm2,%xmm3
movsd  %xmm1,(%rcx,%rdx,8)
add    $0x1,%rdx
add    $0x8,%rbp
cmp    %rax,%rdx
jl     1380 <RB_GS_OmpLoopV3err+0x230>
mov    0x18(%rsp),%rax
movsd  %xmm3,0x0(%r13,%rax,8)
jmp    1330 <RB_GS_OmpLoopV3err+0x1e0>
nopl   (%rax)
cmpl   $0x0,0x28(%rsp)
jle    146e <RB_GS_OmpLoopV3err+0x31e>
cmpq   $0x7,0x48(%rsp)
mov    0x40(%rsp),%rdx
jae    13f0 <RB_GS_OmpLoopV3err+0x2a0>
xor    %eax,%eax
jmp    1441 <RB_GS_OmpLoopV3err+0x2f1>
nopl   0x0(%rax)
xor    %eax,%eax
mov    0x38(%rsp),%rcx
nopw   0x0(%rax,%rax,1)
addsd  0x0(%r13,%rax,8),%xmm0
addsd  0x8(%r13,%rax,8),%xmm0
addsd  0x10(%r13,%rax,8),%xmm0
addsd  0x18(%r13,%rax,8),%xmm0
addsd  0x20(%r13,%rax,8),%xmm0
addsd  0x28(%r13,%rax,8),%xmm0
addsd  0x30(%r13,%rax,8),%xmm0
addsd  0x38(%r13,%rax,8),%xmm0
add    $0x8,%rax
cmp    %rax,%rcx
jne    1400 <RB_GS_OmpLoopV3err+0x2b0>
test   %rdx,%rdx
je     146e <RB_GS_OmpLoopV3err+0x31e>
lea    0x0(,%rax,8),%rax
add    %r13,%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rax,%rcx,8),%xmm0
add    $0x1,%rcx
cmp    %rcx,%rdx
jne    1460 <RB_GS_OmpLoopV3err+0x310>
call   1070 <sqrt@plt>
cltd
idivl  0x24(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    0x8(%rsp),%ecx
add    $0x1,%ecx
ucomisd 0x30(%rsp),%xmm0
jb     14a4 <RB_GS_OmpLoopV3err+0x354>
cmp    0x20(%rsp),%ecx
movsd  0xb6e(%rip),%xmm4
jl     1220 <RB_GS_OmpLoopV3err+0xd0>
jmp    14a4 <RB_GS_OmpLoopV3err+0x354>
mov    %ebx,%ecx
movsd  %xmm0,0x8(%rsp)
mov    0x90(%rsp),%rax
mov    %ecx,(%rax)
mov    %r13,%rdi
call   1030 <free@plt>
movsd  0x8(%rsp),%xmm0
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @RB_GS_OmpLoopV3err(double** %0, double** %1, i32 %2, i32 %3, double %4, i32 %5, i32 %6, i32* %7) {
%9 = alloca double**
%10 = alloca double**
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32*
%17 = alloca i32
%18 = alloca i32
%19 = alloca double
%20 = alloca double*
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca double
store double** %0, double*** %9
store double** %1, double*** %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store double %4, double* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store i32* %7, i32** %16
%25 = load i32*, i32** %16
%26 = load i32, i32* %25
store i32 %26, i32* %23
%27 = load double, double* %13
%28 = fmul double 1.000000e+01, %27
store double %28, double* %24
%29 = call i32 (...) @omp_get_max_threads()
store i32 %29, i32* %22
%30 = load i32, i32* %22
%31 = sext i32 %30 to i64
%32 = mul i64 8, %31
%33 = trunc i64 %32 to i32
%34 = call double* @malloc(i32 %33)
store double* %34, double** %20
%35 = call i32 (...) @omp_get_thread_num()
store i32 %35, i32* %21
br label %36
36: ; preds = %290, %8
%37 = load double, double* %24
%38 = load double, double* %13
%39 = fcmp oge double %37, %38
br i1 %39, label %40, label %44
40: ; preds = %36
%41 = load i32, i32* %23
%42 = load i32, i32* %14
%43 = icmp slt i32 %41, %42
br label %44
44: ; preds = %40, %36
%45 = phi i1 [ false, %36 ], [ %43, %40 ]
br i1 %45, label %46, label %300
46: ; preds = %44
%47 = load double*, double** %20
%48 = load i32, i32* %21
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
store double 0.000000e+00, double* %50
store i32 1, i32* %17
br label %51
51: ; preds = %159, %46
%52 = load i32, i32* %17
%53 = load i32, i32* %11
%54 = sub nsw i32 %53, 1
%55 = icmp slt i32 %52, %54
br i1 %55, label %56, label %162
56: ; preds = %51
%57 = load i32, i32* %17
%58 = srem i32 %57, 2
store i32 %58, i32* %18
br label %59
59: ; preds = %155, %56
%60 = load i32, i32* %18
%61 = load i32, i32* %12
%62 = sdiv i32 %61, 2
%63 = load i32, i32* %17
%64 = add nsw i32 %63, 1
%65 = srem i32 %64, 2
%66 = sub nsw i32 %62, %65
%67 = icmp slt i32 %60, %66
br i1 %67, label %68, label %158
68: ; preds = %59
%69 = load double**, double*** %10
%70 = load i32, i32* %17
%71 = sub nsw i32 %70, 1
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double*, double** %69, i64 %72
%74 = load double*, double** %73
%75 = load i32, i32* %18
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
%78 = load double, double* %77
%79 = load double**, double*** %10
%80 = load i32, i32* %17
%81 = add nsw i32 %80, 1
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds double*, double** %79, i64 %82
%84 = load double*, double** %83
%85 = load i32, i32* %18
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
%88 = load double, double* %87
%89 = fadd double %78, %88
%90 = load double**, double*** %10
%91 = load i32, i32* %17
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double*, double** %90, i64 %92
%94 = load double*, double** %93
%95 = load i32, i32* %18
%96 = load i32, i32* %17
%97 = srem i32 %96, 2
%98 = sub nsw i32 %95, %97
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %94, i64 %99
%101 = load double, double* %100
%102 = fadd double %89, %101
%103 = load double**, double*** %10
%104 = load i32, i32* %17
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds double*, double** %103, i64 %105
%107 = load double*, double** %106
%108 = load i32, i32* %18
%109 = load i32, i32* %17
%110 = add nsw i32 %109, 1
%111 = srem i32 %110, 2
%112 = add nsw i32 %108, %111
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %107, i64 %113
%115 = load double, double* %114
%116 = fadd double %102, %115
%117 = fmul double 2.500000e-01, %116
store double %117, double* %19
%118 = load double, double* %19
%119 = load double**, double*** %9
%120 = load i32, i32* %17
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double*, double** %119, i64 %121
%123 = load double*, double** %122
%124 = load i32, i32* %18
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %123, i64 %125
%127 = load double, double* %126
%128 = fsub double %118, %127
%129 = load double, double* %19
%130 = load double**, double*** %9
%131 = load i32, i32* %17
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double*, double** %130, i64 %132
%134 = load double*, double** %133
%135 = load i32, i32* %18
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds double, double* %134, i64 %136
%138 = load double, double* %137
%139 = fsub double %129, %138
%140 = load double*, double** %20
%141 = load i32, i32* %21
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %140, i64 %142
%144 = load double, double* %143
%145 = call double @llvm.fmuladd.f64(double %128, double %139, double %144)
store double %145, double* %143
%146 = load double, double* %19
%147 = load double**, double*** %9
%148 = load i32, i32* %17
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds double*, double** %147, i64 %149
%151 = load double*, double** %150
%152 = load i32, i32* %18
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds double, double* %151, i64 %153
store double %146, double* %154
br label %155
155: ; preds = %68
%156 = load i32, i32* %18
%157 = add nsw i32 %156, 1
store i32 %157, i32* %18
br label %59
158: ; preds = %59
br label %159
159: ; preds = %158
%160 = load i32, i32* %17
%161 = add nsw i32 %160, 1
store i32 %161, i32* %17
br label %51
162: ; preds = %51
store i32 1, i32* %17
br label %163
163: ; preds = %271, %162
%164 = load i32, i32* %17
%165 = load i32, i32* %11
%166 = sub nsw i32 %165, 1
%167 = icmp slt i32 %164, %166
br i1 %167, label %168, label %274
168: ; preds = %163
%169 = load i32, i32* %17
%170 = add nsw i32 %169, 1
%171 = srem i32 %170, 2
store i32 %171, i32* %18
br label %172
172: ; preds = %267, %168
%173 = load i32, i32* %18
%174 = load i32, i32* %12
%175 = sdiv i32 %174, 2
%176 = load i32, i32* %17
%177 = srem i32 %176, 2
%178 = sub nsw i32 %175, %177
%179 = icmp slt i32 %173, %178
br i1 %179, label %180, label %270
180: ; preds = %172
%181 = load double**, double*** %9
%182 = load i32, i32* %17
%183 = sub nsw i32 %182, 1
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds double*, double** %181, i64 %184
%186 = load double*, double** %185
%187 = load i32, i32* %18
%188 = sext i32 %187 to i64
%189 = getelementptr inbounds double, double* %186, i64 %188
%190 = load double, double* %189
%191 = load double**, double*** %9
%192 = load i32, i32* %17
%193 = add nsw i32 %192, 1
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds double*, double** %191, i64 %194
%196 = load double*, double** %195
%197 = load i32, i32* %18
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds double, double* %196, i64 %198
%200 = load double, double* %199
%201 = fadd double %190, %200
%202 = load double**, double*** %9
%203 = load i32, i32* %17
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds double*, double** %202, i64 %204
%206 = load double*, double** %205
%207 = load i32, i32* %18
%208 = load i32, i32* %17
%209 = add nsw i32 %208, 1
%210 = srem i32 %209, 2
%211 = sub nsw i32 %207, %210
%212 = sext i32 %211 to i64
%213 = getelementptr inbounds double, double* %206, i64 %212
%214 = load double, double* %213
%215 = fadd double %201, %214
%216 = load double**, double*** %9
%217 = load i32, i32* %17
%218 = sext i32 %217 to i64
%219 = getelementptr inbounds double*, double** %216, i64 %218
%220 = load double*, double** %219
%221 = load i32, i32* %18
%222 = load i32, i32* %17
%223 = srem i32 %222, 2
%224 = add nsw i32 %221, %223
%225 = sext i32 %224 to i64
%226 = getelementptr inbounds double, double* %220, i64 %225
%227 = load double, double* %226
%228 = fadd double %215, %227
%229 = fmul double 2.500000e-01, %228
store double %229, double* %19
%230 = load double, double* %19
%231 = load double**, double*** %10
%232 = load i32, i32* %17
%233 = sext i32 %232 to i64
%234 = getelementptr inbounds double*, double** %231, i64 %233
%235 = load double*, double** %234
%236 = load i32, i32* %18
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double, double* %235, i64 %237
%239 = load double, double* %238
%240 = fsub double %230, %239
%241 = load double, double* %19
%242 = load double**, double*** %10
%243 = load i32, i32* %17
%244 = sext i32 %243 to i64
%245 = getelementptr inbounds double*, double** %242, i64 %244
%246 = load double*, double** %245
%247 = load i32, i32* %18
%248 = sext i32 %247 to i64
%249 = getelementptr inbounds double, double* %246, i64 %248
%250 = load double, double* %249
%251 = fsub double %241, %250
%252 = load double*, double** %20
%253 = load i32, i32* %21
%254 = sext i32 %253 to i64
%255 = getelementptr inbounds double, double* %252, i64 %254
%256 = load double, double* %255
%257 = call double @llvm.fmuladd.f64(double %240, double %251, double %256)
store double %257, double* %255
%258 = load double, double* %19
%259 = load double**, double*** %10
%260 = load i32, i32* %17
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double*, double** %259, i64 %261
%263 = load double*, double** %262
%264 = load i32, i32* %18
%265 = sext i32 %264 to i64
%266 = getelementptr inbounds double, double* %263, i64 %265
store double %258, double* %266
br label %267
267: ; preds = %180
%268 = load i32, i32* %18
%269 = add nsw i32 %268, 1
store i32 %269, i32* %18
br label %172
270: ; preds = %172
br label %271
271: ; preds = %270
%272 = load i32, i32* %17
%273 = add nsw i32 %272, 1
store i32 %273, i32* %17
br label %163
274: ; preds = %163
store i32 0, i32* %18
br label %275
275: ; preds = %287, %274
%276 = load i32, i32* %18
%277 = load i32, i32* %22
%278 = icmp slt i32 %276, %277
br i1 %278, label %279, label %290
279: ; preds = %275
%280 = load double*, double** %20
%281 = load i32, i32* %18
%282 = sext i32 %281 to i64
%283 = getelementptr inbounds double, double* %280, i64 %282
%284 = load double, double* %283
%285 = load double, double* %24
%286 = fadd double %285, %284
store double %286, double* %24
br label %287
287: ; preds = %279
%288 = load i32, i32* %18
%289 = add nsw i32 %288, 1
store i32 %289, i32* %18
br label %275
290: ; preds = %275
%291 = load double, double* %24
%292 = call i32 @sqrt(double %291)
%293 = load i32, i32* %12
%294 = load i32, i32* %11
%295 = mul nsw i32 %293, %294
%296 = sdiv i32 %292, %295
%297 = sitofp i32 %296 to double
store double %297, double* %24
%298 = load i32, i32* %23
%299 = add nsw i32 %298, 1
store i32 %299, i32* %23
br label %36
300: ; preds = %44
%301 = load i32, i32* %23
%302 = load i32*, i32** %16
store i32 %301, i32* %302
%303 = load double*, double** %20
%304 = call i32 @free(double* %303)
%305 = load double, double* %24
ret double %305
}

LLVM IR (-O2):
define double @RB_GS_OmpLoopV3err(double** nocapture %0, double** nocapture %1, i32 %2, i32 %3, double %4, i32 %5, i32 %6, i32* nocapture %7) {
%9 = load i32, i32* %7
%10 = fmul double %4, 1.000000e+01
%11 = tail call i32 (...) @omp_get_max_threads()
%12 = shl i32 %11, 3
%13 = tail call double* @malloc(i32 %12)
%14 = tail call i32 (...) @omp_get_thread_num()
%15 = fcmp oge double %10, %4
%16 = icmp slt i32 %9, %5
%17 = select i1 %15, i1 %16, i1 false
br i1 %17, label %18, label %208
18: ; preds = %8
%19 = sext i32 %14 to i64
%20 = getelementptr inbounds double, double* %13, i64 %19
%21 = add i32 %2, -1
%22 = icmp sgt i32 %2, 2
%23 = sdiv i32 %3, 2
%24 = icmp sgt i32 %11, 0
%25 = mul nsw i32 %3, %2
%26 = zext i32 %21 to i64
%27 = zext i32 %21 to i64
%28 = zext i32 %11 to i64
%29 = add nsw i64 %28, -1
%30 = and i64 %28, 7
%31 = icmp ult i64 %29, 7
%32 = and i64 %28, 4294967288
%33 = icmp eq i64 %30, 0
br label %34
34: ; preds = %18, %199
%35 = phi double [ %10, %18 ], [ %203, %199 ]
%36 = phi i32 [ %9, %18 ], [ %204, %199 ]
store double 0.000000e+00, double* %20
br i1 %22, label %42, label %97
37: ; preds = %69
store double %89, double* %20
br label %38
38: ; preds = %37, %42
%39 = phi double [ %89, %37 ], [ %43, %42 ]
%40 = icmp eq i64 %48, %26
br i1 %40, label %41, label %42
41: ; preds = %38
br i1 %22, label %99, label %97
42: ; preds = %34, %38
%43 = phi double [ %39, %38 ], [ 0.000000e+00, %34 ]
%44 = phi i64 [ %48, %38 ], [ 1, %34 ]
%45 = phi i32 [ %49, %38 ], [ 1, %34 ]
%46 = trunc i64 %44 to i32
%47 = and i32 %46, 1
%48 = add nuw nsw i64 %44, 1
%49 = add nuw nsw i32 %45, 1
%50 = trunc i64 %48 to i32
%51 = and i32 %50, 1
%52 = sub nsw i32 %23, %51
%53 = icmp slt i32 %47, %52
br i1 %53, label %54, label %38
54: ; preds = %42
%55 = and i32 %45, 1
%56 = zext i32 %55 to i64
%57 = add nsw i64 %44, -1
%58 = getelementptr inbounds double*, double** %1, i64 %57
%59 = load double*, double** %58
%60 = getelementptr inbounds double*, double** %1, i64 %48
%61 = load double*, double** %60
%62 = getelementptr inbounds double*, double** %1, i64 %44
%63 = load double*, double** %62
%64 = getelementptr inbounds double*, double** %0, i64 %44
%65 = load double*, double** %64
%66 = zext i32 %47 to i64
%67 = zext i32 %51 to i64
%68 = sext i32 %52 to i64
br label %69
69: ; preds = %54, %69
%70 = phi i64 [ %56, %54 ], [ %90, %69 ]
%71 = phi double [ %43, %54 ], [ %89, %69 ]
%72 = getelementptr inbounds double, double* %59, i64 %70
%73 = load double, double* %72
%74 = getelementptr inbounds double, double* %61, i64 %70
%75 = load double, double* %74
%76 = fadd double %73, %75
%77 = sub nuw nsw i64 %70, %66
%78 = getelementptr inbounds double, double* %63, i64 %77
%79 = load double, double* %78
%80 = fadd double %76, %79
%81 = add nuw nsw i64 %70, %67
%82 = getelementptr inbounds double, double* %63, i64 %81
%83 = load double, double* %82
%84 = fadd double %80, %83
%85 = fmul double %84, 2.500000e-01
%86 = getelementptr inbounds double, double* %65, i64 %70
%87 = load double, double* %86
%88 = fsub double %85, %87
%89 = tail call double @llvm.fmuladd.f64(double %88, double %88, double %71)
store double %85, double* %86
%90 = add nuw nsw i64 %70, 1
%91 = icmp slt i64 %90, %68
br i1 %91, label %69, label %37
92: ; preds = %124
store double %144, double* %20
br label %93
93: ; preds = %92, %99
%94 = phi double [ %144, %92 ], [ %100, %99 ]
%95 = xor i1 %102, true
%96 = icmp eq i64 %103, %27
br i1 %96, label %97, label %99
97: ; preds = %93, %34, %41
br i1 %24, label %98, label %199
98: ; preds = %97
br i1 %31, label %185, label %147
99: ; preds = %41, %93
%100 = phi double [ %94, %93 ], [ %39, %41 ]
%101 = phi i64 [ %103, %93 ], [ 1, %41 ]
%102 = phi i1 [ %95, %93 ], [ false, %41 ]
%103 = add nuw nsw i64 %101, 1
%104 = trunc i64 %103 to i32
%105 = and i32 %104, 1
%106 = trunc i64 %101 to i32
%107 = and i32 %106, 1
%108 = sub nsw i32 %23, %107
%109 = icmp slt i32 %105, %108
br i1 %109, label %110, label %93
110: ; preds = %99
%111 = zext i1 %102 to i64
%112 = add nsw i64 %101, -1
%113 = getelementptr inbounds double*, double** %0, i64 %112
%114 = load double*, double** %113
%115 = getelementptr inbounds double*, double** %0, i64 %103
%116 = load double*, double** %115
%117 = getelementptr inbounds double*, double** %0, i64 %101
%118 = load double*, double** %117
%119 = getelementptr inbounds double*, double** %1, i64 %101
%120 = load double*, double** %119
%121 = zext i32 %105 to i64
%122 = zext i32 %107 to i64
%123 = sext i32 %108 to i64
br label %124
124: ; preds = %110, %124
%125 = phi i64 [ %111, %110 ], [ %145, %124 ]
%126 = phi double [ %100, %110 ], [ %144, %124 ]
%127 = getelementptr inbounds double, double* %114, i64 %125
%128 = load double, double* %127
%129 = getelementptr inbounds double, double* %116, i64 %125
%130 = load double, double* %129
%131 = fadd double %128, %130
%132 = sub nuw nsw i64 %125, %121
%133 = getelementptr inbounds double, double* %118, i64 %132
%134 = load double, double* %133
%135 = fadd double %131, %134
%136 = add nuw nsw i64 %125, %122
%137 = getelementptr inbounds double, double* %118, i64 %136
%138 = load double, double* %137
%139 = fadd double %135, %138
%140 = fmul double %139, 2.500000e-01
%141 = getelementptr inbounds double, double* %120, i64 %125
%142 = load double, double* %141
%143 = fsub double %140, %142
%144 = tail call double @llvm.fmuladd.f64(double %143, double %143, double %126)
store double %140, double* %141
%145 = add nuw nsw i64 %125, 1
%146 = icmp slt i64 %145, %123
br i1 %146, label %124, label %92
147: ; preds = %98, %147
%148 = phi i64 [ %182, %147 ], [ 0, %98 ]
%149 = phi double [ %181, %147 ], [ %35, %98 ]
%150 = phi i64 [ %183, %147 ], [ 0, %98 ]
%151 = getelementptr inbounds double, double* %13, i64 %148
%152 = load double, double* %151
%153 = fadd double %149, %152
%154 = or i64 %148, 1
%155 = getelementptr inbounds double, double* %13, i64 %154
%156 = load double, double* %155
%157 = fadd double %153, %156
%158 = or i64 %148, 2
%159 = getelementptr inbounds double, double* %13, i64 %158
%160 = load double, double* %159
%161 = fadd double %157, %160
%162 = or i64 %148, 3
%163 = getelementptr inbounds double, double* %13, i64 %162
%164 = load double, double* %163
%165 = fadd double %161, %164
%166 = or i64 %148, 4
%167 = getelementptr inbounds double, double* %13, i64 %166
%168 = load double, double* %167
%169 = fadd double %165, %168
%170 = or i64 %148, 5
%171 = getelementptr inbounds double, double* %13, i64 %170
%172 = load double, double* %171
%173 = fadd double %169, %172
%174 = or i64 %148, 6
%175 = getelementptr inbounds double, double* %13, i64 %174
%176 = load double, double* %175
%177 = fadd double %173, %176
%178 = or i64 %148, 7
%179 = getelementptr inbounds double, double* %13, i64 %178
%180 = load double, double* %179
%181 = fadd double %177, %180
%182 = add nuw nsw i64 %148, 8
%183 = add i64 %150, 8
%184 = icmp eq i64 %183, %32
br i1 %184, label %185, label %147
185: ; preds = %147, %98
%186 = phi double [ undef, %98 ], [ %181, %147 ]
%187 = phi i64 [ 0, %98 ], [ %182, %147 ]
%188 = phi double [ %35, %98 ], [ %181, %147 ]
br i1 %33, label %199, label %189
189: ; preds = %185, %189
%190 = phi i64 [ %196, %189 ], [ %187, %185 ]
%191 = phi double [ %195, %189 ], [ %188, %185 ]
%192 = phi i64 [ %197, %189 ], [ 0, %185 ]
%193 = getelementptr inbounds double, double* %13, i64 %190
%194 = load double, double* %193
%195 = fadd double %191, %194
%196 = add nuw nsw i64 %190, 1
%197 = add i64 %192, 1
%198 = icmp eq i64 %197, %30
br i1 %198, label %199, label %189
199: ; preds = %185, %189, %97
%200 = phi double [ %35, %97 ], [ %186, %185 ], [ %195, %189 ]
%201 = tail call i32 @sqrt(double %200)
%202 = sdiv i32 %201, %25
%203 = sitofp i32 %202 to double
%204 = add nsw i32 %36, 1
%205 = fcmp oge double %203, %4
%206 = icmp slt i32 %204, %5
%207 = select i1 %205, i1 %206, i1 false
br i1 %207, label %34, label %208
208: ; preds = %199, %8
%209 = phi i32 [ %9, %8 ], [ %204, %199 ]
%210 = phi double [ %10, %8 ], [ %203, %199 ]
store i32 %209, i32* %7
%211 = tail call i32 @free(double* %13)
ret double %210
}