#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 fabs(double);
double rdfabsorb(double *, double *, int);
int /*<<< orphan*/ rdfdiffuse(double *, double *, double *, int, double);
int /*<<< orphan*/ rdfreverserxn(double *, double *, int, double, double,
                                 double);
double rdfsteadystate(double *r, double *rdfa, double *rdfd, int n, double step,
                      double b, double eps) {
  const int maxit = 100000;
  const double maxflux = 1e7;
  int i, it;
  double flux, fluxold;
  rdfdiffuse(r, rdfa, rdfd, n, step);
  flux = fluxold = rdfabsorb(r, rdfd, n);
  for (it = 0; it < 30 || (it < maxit && flux < maxflux &&
                           fabs((flux - fluxold) / (fluxold + 1e-20)) > eps);
       it++) {
    fluxold = flux;
    rdfdiffuse(r, rdfa, rdfd, n, step);
    if (b >= 0)
      rdfreverserxn(r, rdfd, n, step, b, flux);
    for (i = 0; i < n; i++)
      rdfa[i] = rdfd[i];
    flux = rdfabsorb(r, rdfa, n);
  }
  if (it >= maxit || flux >= maxflux)
    flux = -1;
  return flux;
}


x86 assembly (clang, -O0):
<rdfsteadystate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  %xmm1,-0x30(%rbp)
movsd  %xmm2,-0x38(%rbp)
movl   $0x186a0,-0x3c(%rbp)
movsd  0xe9b(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
movsd  -0x28(%rbp),%xmm0
call   1040 <rdfdiffuse@plt>
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x1c(%rbp),%edx
call   1050 <rdfabsorb@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x50(%rbp)
mov    $0x1,%al
cmpl   $0x1e,-0x50(%rbp)
mov    %al,-0x61(%rbp)
jl     1214 <rdfsteadystate+0xe4>
xor    %eax,%eax
cmpl   $0x186a0,-0x50(%rbp)
mov    %al,-0x62(%rbp)
jge    120e <rdfsteadystate+0xde>
xor    %eax,%eax
movsd  0xe31(%rip),%xmm0
ucomisd -0x58(%rbp),%xmm0
mov    %al,-0x62(%rbp)
jbe    120e <rdfsteadystate+0xde>
movsd  -0x58(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  0xe19(%rip),%xmm1
addsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movaps 0xe21(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x38(%rbp),%xmm0
seta   %al
mov    %al,-0x62(%rbp)
mov    -0x62(%rbp),%al
mov    %al,-0x61(%rbp)
mov    -0x61(%rbp),%al
test   $0x1,%al
jne    1224 <rdfsteadystate+0xf4>
jmp    12d6 <rdfsteadystate+0x1a6>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
movsd  -0x28(%rbp),%xmm0
call   1040 <rdfdiffuse@plt>
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1278 <rdfsteadystate+0x148>
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x1c(%rbp),%edx
movsd  -0x28(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm2
call   1030 <rdfreverserxn@plt>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12b3 <rdfsteadystate+0x183>
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    127f <rdfsteadystate+0x14f>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x1c(%rbp),%edx
call   1050 <rdfabsorb@plt>
movsd  %xmm0,-0x58(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    11a4 <rdfsteadystate+0x74>
cmpl   $0x186a0,-0x50(%rbp)
jge    12fa <rdfsteadystate+0x1ca>
movsd  -0x58(%rbp),%xmm0
movsd  0xd10(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1307 <rdfsteadystate+0x1d7>
movsd  0xd0e(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rdfsteadystate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
movsd  %xmm2,0x58(%rsp)
movsd  %xmm1,0x18(%rsp)
mov    %ecx,%r14d
mov    %rdx,%rbx
mov    %rsi,%rbp
mov    %rdi,%r15
movsd  %xmm0,0x10(%rsp)
call   1040 <rdfdiffuse@plt>
mov    %r15,0x8(%rsp)
mov    %r15,%rdi
mov    %rbx,%rsi
mov    %r14d,%edx
call   1050 <rdfabsorb@plt>
mov    %r14d,%r12d
lea    (%rbx,%r12,8),%rax
mov    %r12d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0x48(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x40(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
cmp    %rbp,%rax
seta   %al
lea    0x0(,%r12,8),%rcx
add    %rbp,%rcx
cmp    %rbx,%rcx
seta   %cl
and    %al,%cl
mov    %cl,0x7(%rsp)
mov    %edx,%r15d
and    $0x3,%r15d
and    $0xfffffffffffffffc,%rdx
mov    %rdx,0x20(%rsp)
mov    %r12d,%eax
and    $0x3,%eax
lea    0x10(%rbp),%rcx
mov    %rcx,0x30(%rsp)
lea    0x10(%rbx),%rcx
mov    %rcx,0x28(%rsp)
mov    %r15,0x38(%rsp)
shl    $0x5,%r15
xor    %r13d,%r13d
movapd %xmm0,%xmm1
mov    %rax,0x50(%rsp)
jmp    121a <rdfsteadystate+0xea>
cs nopw 0x0(%rax,%rax,1)
mov    0x8(%rsp),%rdi
mov    %rbp,%rsi
mov    %r14d,%edx
call   1050 <rdfabsorb@plt>
add    $0x1,%r13d
movapd 0x60(%rsp),%xmm1
movapd %xmm0,%xmm2
cmp    $0x1e,%r13d
jb     126b <rdfsteadystate+0x13b>
cmp    $0x1869f,%r13d
ja     142d <rdfsteadystate+0x2fd>
movsd  0xdc7(%rip),%xmm0
ucomisd %xmm2,%xmm0
jbe    142d <rdfsteadystate+0x2fd>
movapd %xmm2,%xmm0
subsd  %xmm1,%xmm0
addsd  0xdb5(%rip),%xmm1
divsd  %xmm1,%xmm0
andpd  0xdc1(%rip),%xmm0
ucomisd 0x58(%rsp),%xmm0
jbe    142d <rdfsteadystate+0x2fd>
movapd %xmm2,0x60(%rsp)
mov    0x8(%rsp),%rdi
mov    %rbp,%rsi
mov    %rbx,%rdx
mov    %r14d,%ecx
movsd  0x10(%rsp),%xmm0
call   1040 <rdfdiffuse@plt>
movsd  0x18(%rsp),%xmm0
ucomisd 0xd80(%rip),%xmm0
jb     12bc <rdfsteadystate+0x18c>
mov    0x8(%rsp),%rdi
mov    %rbx,%rsi
mov    %r14d,%edx
movsd  0x10(%rsp),%xmm0
movsd  0x18(%rsp),%xmm1
movapd 0x60(%rsp),%xmm2
call   1030 <rdfreverserxn@plt>
test   %r14d,%r14d
jle    1200 <rdfsteadystate+0xd0>
cmp    $0x4,%r14d
setb   %al
or     0x7(%rsp),%al
je     12e0 <rdfsteadystate+0x1b0>
xor    %eax,%eax
jmp    13b0 <rdfsteadystate+0x280>
nopl   0x0(%rax)
cmpq   $0xc,0x40(%rsp)
jae    12ec <rdfsteadystate+0x1bc>
xor    %eax,%eax
jmp    135d <rdfsteadystate+0x22d>
mov    0x20(%rsp),%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups (%rbx,%rax,8),%xmm0
movups 0x10(%rbx,%rax,8),%xmm1
movups %xmm0,0x0(%rbp,%rax,8)
movups %xmm1,0x10(%rbp,%rax,8)
movups 0x20(%rbx,%rax,8),%xmm0
movups 0x30(%rbx,%rax,8),%xmm1
movups %xmm0,0x20(%rbp,%rax,8)
movups %xmm1,0x30(%rbp,%rax,8)
movups 0x40(%rbx,%rax,8),%xmm0
movups 0x50(%rbx,%rax,8),%xmm1
movups %xmm0,0x40(%rbp,%rax,8)
movups %xmm1,0x50(%rbp,%rax,8)
movupd 0x60(%rbx,%rax,8),%xmm0
movupd 0x70(%rbx,%rax,8),%xmm1
movupd %xmm0,0x60(%rbp,%rax,8)
movupd %xmm1,0x70(%rbp,%rax,8)
add    $0x10,%rax
add    $0xfffffffffffffffc,%rcx
jne    1300 <rdfsteadystate+0x1d0>
cmpq   $0x0,0x38(%rsp)
je     139f <rdfsteadystate+0x26f>
mov    0x30(%rsp),%rcx
lea    (%rcx,%rax,8),%rcx
mov    0x28(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
xor    %edx,%edx
nopl   0x0(%rax)
movupd -0x10(%rax,%rdx,1),%xmm0
movupd (%rax,%rdx,1),%xmm1
movupd %xmm0,-0x10(%rcx,%rdx,1)
movupd %xmm1,(%rcx,%rdx,1)
add    $0x20,%rdx
cmp    %rdx,%r15
jne    1380 <rdfsteadystate+0x250>
mov    0x48(%rsp),%rcx
mov    %rcx,%rax
cmp    %r12,%rcx
je     1200 <rdfsteadystate+0xd0>
mov    %rax,%rcx
not    %rcx
add    %r12,%rcx
mov    0x50(%rsp),%rdx
test   %rdx,%rdx
je     13e5 <rdfsteadystate+0x2b5>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    13d0 <rdfsteadystate+0x2a0>
cmp    $0x3,%rcx
jb     1200 <rdfsteadystate+0xd0>
nop
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
movsd  0x8(%rbx,%rax,8),%xmm0
movsd  %xmm0,0x8(%rbp,%rax,8)
movsd  0x10(%rbx,%rax,8),%xmm0
movsd  %xmm0,0x10(%rbp,%rax,8)
movsd  0x18(%rbx,%rax,8),%xmm0
movsd  %xmm0,0x18(%rbp,%rax,8)
add    $0x4,%rax
cmp    %rax,%r12
jne    13f0 <rdfsteadystate+0x2c0>
jmp    1200 <rdfsteadystate+0xd0>
movsd  0xbdb(%rip),%xmm0
cmp    $0x186a0,%r13d
jae    1457 <rdfsteadystate+0x327>
movsd  0xbba(%rip),%xmm1
cmplesd %xmm2,%xmm1
andpd  %xmm1,%xmm0
andnpd %xmm2,%xmm1
orpd   %xmm1,%xmm0
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @rdfsteadystate(double* %0, double* %1, double* %2, i32 %3, double %4, double %5, double %6) {
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca i32
%16 = alloca double
%17 = alloca i32
%18 = alloca i32
%19 = alloca double
%20 = alloca double
store double* %0, double** %8
store double* %1, double** %9
store double* %2, double** %10
store i32 %3, i32* %11
store double %4, double* %12
store double %5, double* %13
store double %6, double* %14
store i32 100000, i32* %15
store double 1.000000e+07, double* %16
%21 = load double*, double** %8
%22 = load double*, double** %9
%23 = load double*, double** %10
%24 = load i32, i32* %11
%25 = load double, double* %12
%26 = call i32 @rdfdiffuse(double* %21, double* %22, double* %23, i32 %24, double %25)
%27 = load double*, double** %8
%28 = load double*, double** %10
%29 = load i32, i32* %11
%30 = call double @rdfabsorb(double* %27, double* %28, i32 %29)
store double %30, double* %20
store double %30, double* %19
store i32 0, i32* %18
br label %31
31: ; preds = %95, %7
%32 = load i32, i32* %18
%33 = icmp slt i32 %32, 30
br i1 %33, label %52, label %34
34: ; preds = %31
%35 = load i32, i32* %18
%36 = icmp slt i32 %35, 100000
br i1 %36, label %37, label %50
37: ; preds = %34
%38 = load double, double* %19
%39 = fcmp olt double %38, 1.000000e+07
br i1 %39, label %40, label %50
40: ; preds = %37
%41 = load double, double* %19
%42 = load double, double* %20
%43 = fsub double %41, %42
%44 = load double, double* %20
%45 = fadd double %44, 0x3BC79CA10C924223
%46 = fdiv double %43, %45
%47 = call double @llvm.fabs.f64(double %46)
%48 = load double, double* %14
%49 = fcmp ogt double %47, %48
br label %50
50: ; preds = %40, %37, %34
%51 = phi i1 [ false, %37 ], [ false, %34 ], [ %49, %40 ]
br label %52
52: ; preds = %50, %31
%53 = phi i1 [ true, %31 ], [ %51, %50 ]
br i1 %53, label %54, label %98
54: ; preds = %52
%55 = load double, double* %19
store double %55, double* %20
%56 = load double*, double** %8
%57 = load double*, double** %9
%58 = load double*, double** %10
%59 = load i32, i32* %11
%60 = load double, double* %12
%61 = call i32 @rdfdiffuse(double* %56, double* %57, double* %58, i32 %59, double %60)
%62 = load double, double* %13
%63 = fcmp oge double %62, 0.000000e+00
br i1 %63, label %64, label %72
64: ; preds = %54
%65 = load double*, double** %8
%66 = load double*, double** %10
%67 = load i32, i32* %11
%68 = load double, double* %12
%69 = load double, double* %13
%70 = load double, double* %19
%71 = call i32 @rdfreverserxn(double* %65, double* %66, i32 %67, double %68, double %69, double %70)
br label %72
72: ; preds = %64, %54
store i32 0, i32* %17
br label %73
73: ; preds = %87, %72
%74 = load i32, i32* %17
%75 = load i32, i32* %11
%76 = icmp slt i32 %74, %75
br i1 %76, label %77, label %90
77: ; preds = %73
%78 = load double*, double** %10
%79 = load i32, i32* %17
%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* %17
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
store double %82, double* %86
br label %87
87: ; preds = %77
%88 = load i32, i32* %17
%89 = add nsw i32 %88, 1
store i32 %89, i32* %17
br label %73
90: ; preds = %73
%91 = load double*, double** %8
%92 = load double*, double** %9
%93 = load i32, i32* %11
%94 = call double @rdfabsorb(double* %91, double* %92, i32 %93)
store double %94, double* %19
br label %95
95: ; preds = %90
%96 = load i32, i32* %18
%97 = add nsw i32 %96, 1
store i32 %97, i32* %18
br label %31
98: ; preds = %52
%99 = load i32, i32* %18
%100 = icmp sge i32 %99, 100000
br i1 %100, label %104, label %101
101: ; preds = %98
%102 = load double, double* %19
%103 = fcmp oge double %102, 1.000000e+07
br i1 %103, label %104, label %105
104: ; preds = %101, %98
store double -1.000000e+00, double* %19
br label %105
105: ; preds = %104, %101
%106 = load double, double* %19
ret double %106
}

LLVM IR (-O2):
define double @rdfsteadystate(double* %0, double* %1, double* %2, i32 %3, double %4, double %5, double %6) {
%8 = tail call i32 @rdfdiffuse(double* %0, double* %1, double* %2, i32 %3, double %4)
%9 = tail call double @rdfabsorb(double* %0, double* %2, i32 %3)
%10 = fcmp ult double %5, 0.000000e+00
%11 = icmp sgt i32 %3, 0
%12 = zext i32 %3 to i64
%13 = getelementptr double, double* %1, i64 %12
%14 = getelementptr double, double* %2, i64 %12
%15 = and i64 %12, 4294967292
%16 = add nsw i64 %15, -4
%17 = lshr exact i64 %16, 2
%18 = add nuw nsw i64 %17, 1
%19 = icmp ult i32 %3, 4
%20 = icmp ugt double* %14, %1
%21 = icmp ugt double* %13, %2
%22 = and i1 %20, %21
%23 = and i64 %12, 4294967292
%24 = and i64 %18, 3
%25 = icmp ult i64 %16, 12
%26 = and i64 %18, 9223372036854775804
%27 = icmp eq i64 %24, 0
%28 = icmp eq i64 %23, %12
%29 = and i64 %12, 3
%30 = icmp eq i64 %29, 0
br label %31
31: ; preds = %157, %7
%32 = phi i32 [ 0, %7 ], [ %159, %157 ]
%33 = phi double [ %9, %7 ], [ %158, %157 ]
%34 = phi double [ %9, %7 ], [ %33, %157 ]
%35 = icmp ult i32 %32, 30
br i1 %35, label %46, label %36
36: ; preds = %31
%37 = icmp ult i32 %32, 100000
%38 = fcmp olt double %33, 1.000000e+07
%39 = select i1 %37, i1 %38, i1 false
br i1 %39, label %40, label %160
40: ; preds = %36
%41 = fsub double %33, %34
%42 = fadd double %34, 0x3BC79CA10C924223
%43 = fdiv double %41, %42
%44 = tail call double @llvm.fabs.f64(double %43)
%45 = fcmp ogt double %44, %6
br i1 %45, label %46, label %160
46: ; preds = %31, %40
%47 = tail call i32 @rdfdiffuse(double* %0, double* %1, double* %2, i32 %3, double %4)
br i1 %10, label %50, label %48
48: ; preds = %46
%49 = tail call i32 @rdfreverserxn(double* %0, double* %2, i32 %3, double %4, double %5, double %33)
br label %50
50: ; preds = %48, %46
br i1 %11, label %51, label %157
51: ; preds = %50
%52 = or i1 %19, %22
br i1 %52, label %122, label %53
53: ; preds = %51
br i1 %25, label %103, label %54
54: ; preds = %53, %54
%55 = phi i64 [ %100, %54 ], [ 0, %53 ]
%56 = phi i64 [ %101, %54 ], [ 0, %53 ]
%57 = getelementptr inbounds double, double* %2, i64 %55
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = getelementptr inbounds double, double* %57, i64 2
%61 = bitcast double* %60 to <2 x double>*
%62 = load <2 x double>, <2 x double>* %61
%63 = getelementptr inbounds double, double* %1, i64 %55
%64 = bitcast double* %63 to <2 x double>*
store <2 x double> %59, <2 x double>* %64
%65 = getelementptr inbounds double, double* %63, i64 2
%66 = bitcast double* %65 to <2 x double>*
store <2 x double> %62, <2 x double>* %66
%67 = or i64 %55, 4
%68 = getelementptr inbounds double, double* %2, i64 %67
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = getelementptr inbounds double, double* %68, i64 2
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = getelementptr inbounds double, double* %1, i64 %67
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> %70, <2 x double>* %75
%76 = getelementptr inbounds double, double* %74, i64 2
%77 = bitcast double* %76 to <2 x double>*
store <2 x double> %73, <2 x double>* %77
%78 = or i64 %55, 8
%79 = getelementptr inbounds double, double* %2, i64 %78
%80 = bitcast double* %79 to <2 x double>*
%81 = load <2 x double>, <2 x double>* %80
%82 = getelementptr inbounds double, double* %79, i64 2
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = getelementptr inbounds double, double* %1, i64 %78
%86 = bitcast double* %85 to <2 x double>*
store <2 x double> %81, <2 x double>* %86
%87 = getelementptr inbounds double, double* %85, i64 2
%88 = bitcast double* %87 to <2 x double>*
store <2 x double> %84, <2 x double>* %88
%89 = or i64 %55, 12
%90 = getelementptr inbounds double, double* %2, i64 %89
%91 = bitcast double* %90 to <2 x double>*
%92 = load <2 x double>, <2 x double>* %91
%93 = getelementptr inbounds double, double* %90, i64 2
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = getelementptr inbounds double, double* %1, i64 %89
%97 = bitcast double* %96 to <2 x double>*
store <2 x double> %92, <2 x double>* %97
%98 = getelementptr inbounds double, double* %96, i64 2
%99 = bitcast double* %98 to <2 x double>*
store <2 x double> %95, <2 x double>* %99
%100 = add nuw i64 %55, 16
%101 = add nuw i64 %56, 4
%102 = icmp eq i64 %101, %26
br i1 %102, label %103, label %54
103: ; preds = %54, %53
%104 = phi i64 [ 0, %53 ], [ %100, %54 ]
br i1 %27, label %121, label %105
105: ; preds = %103, %105
%106 = phi i64 [ %118, %105 ], [ %104, %103 ]
%107 = phi i64 [ %119, %105 ], [ 0, %103 ]
%108 = getelementptr inbounds double, double* %2, i64 %106
%109 = bitcast double* %108 to <2 x double>*
%110 = load <2 x double>, <2 x double>* %109
%111 = getelementptr inbounds double, double* %108, i64 2
%112 = bitcast double* %111 to <2 x double>*
%113 = load <2 x double>, <2 x double>* %112
%114 = getelementptr inbounds double, double* %1, i64 %106
%115 = bitcast double* %114 to <2 x double>*
store <2 x double> %110, <2 x double>* %115
%116 = getelementptr inbounds double, double* %114, i64 2
%117 = bitcast double* %116 to <2 x double>*
store <2 x double> %113, <2 x double>* %117
%118 = add nuw i64 %106, 4
%119 = add i64 %107, 1
%120 = icmp eq i64 %119, %24
br i1 %120, label %121, label %105
121: ; preds = %105, %103
br i1 %28, label %157, label %122
122: ; preds = %51, %121
%123 = phi i64 [ 0, %51 ], [ %23, %121 ]
%124 = xor i64 %123, -1
%125 = add nsw i64 %124, %12
br i1 %30, label %135, label %126
126: ; preds = %122, %126
%127 = phi i64 [ %132, %126 ], [ %123, %122 ]
%128 = phi i64 [ %133, %126 ], [ 0, %122 ]
%129 = getelementptr inbounds double, double* %2, i64 %127
%130 = load double, double* %129
%131 = getelementptr inbounds double, double* %1, i64 %127
store double %130, double* %131
%132 = add nuw nsw i64 %127, 1
%133 = add i64 %128, 1
%134 = icmp eq i64 %133, %29
br i1 %134, label %135, label %126
135: ; preds = %126, %122
%136 = phi i64 [ %123, %122 ], [ %132, %126 ]
%137 = icmp ult i64 %125, 3
br i1 %137, label %157, label %138
138: ; preds = %135, %138
%139 = phi i64 [ %155, %138 ], [ %136, %135 ]
%140 = getelementptr inbounds double, double* %2, i64 %139
%141 = load double, double* %140
%142 = getelementptr inbounds double, double* %1, i64 %139
store double %141, double* %142
%143 = add nuw nsw i64 %139, 1
%144 = getelementptr inbounds double, double* %2, i64 %143
%145 = load double, double* %144
%146 = getelementptr inbounds double, double* %1, i64 %143
store double %145, double* %146
%147 = add nuw nsw i64 %139, 2
%148 = getelementptr inbounds double, double* %2, i64 %147
%149 = load double, double* %148
%150 = getelementptr inbounds double, double* %1, i64 %147
store double %149, double* %150
%151 = add nuw nsw i64 %139, 3
%152 = getelementptr inbounds double, double* %2, i64 %151
%153 = load double, double* %152
%154 = getelementptr inbounds double, double* %1, i64 %151
store double %153, double* %154
%155 = add nuw nsw i64 %139, 4
%156 = icmp eq i64 %155, %12
br i1 %156, label %157, label %138
157: ; preds = %135, %138, %121, %50
%158 = tail call double @rdfabsorb(double* %0, double* %1, i32 %3)
%159 = add nuw nsw i32 %32, 1
br label %31
160: ; preds = %36, %40
%161 = icmp ugt i32 %32, 99999
%162 = fcmp oge double %33, 1.000000e+07
%163 = select i1 %161, i1 true, i1 %162
%164 = select i1 %163, double -1.000000e+00, double %33
ret double %164
}