#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 NonlinearADR_Decay_DiracIC(int *, double *, int, double, double *,
                               double *);
double pow(double, double);
int NonlinearADR_Decay_DiracIC_r(int *iwork, double *rwork, int nPoints,
                                 double T, double *x, double *u, double *r) {
  int i, iret;
  double c = rwork[10], d = rwork[11];
  iret = NonlinearADR_Decay_DiracIC(iwork, rwork, nPoints, T, x, u);
  for (i = 0; i < nPoints; i++) {
    r[i] = c * pow(u[i], d);
  }
  return iret;
}


x86 assembly (clang, -O0):
<NonlinearADR_Decay_DiracIC_r>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x50(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x58(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%edx
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rcx
mov    -0x30(%rbp),%r8
call   1030 <NonlinearADR_Decay_DiracIC@plt>
mov    %eax,-0x40(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11db <NonlinearADR_Decay_DiracIC_r+0xbb>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm1
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1187 <NonlinearADR_Decay_DiracIC_r+0x67>
mov    -0x40(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<NonlinearADR_Decay_DiracIC_r>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,%r13
mov    %r8,%rbx
mov    %edx,%r15d
movsd  0x50(%rsi),%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  0x58(%rsi),%xmm1
movsd  %xmm1,0x8(%rsp)
call   1030 <NonlinearADR_Decay_DiracIC@plt>
mov    %eax,%r14d
test   %r15d,%r15d
jle    1249 <NonlinearADR_Decay_DiracIC_r+0x129>
mov    %r15d,%r12d
lea    -0x1(%r12),%rax
mov    %r12d,%r15d
and    $0x3,%r15d
cmp    $0x3,%rax
jae    117a <NonlinearADR_Decay_DiracIC_r+0x5a>
xor    %ebp,%ebp
jmp    1208 <NonlinearADR_Decay_DiracIC_r+0xe8>
and    $0xfffffffc,%r12d
xor    %ebp,%ebp
movsd  (%rbx,%rbp,8),%xmm0
movsd  0x8(%rsp),%xmm1
call   1040 <pow@plt>
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x0(%r13,%rbp,8)
movsd  0x8(%rbx,%rbp,8),%xmm0
movsd  0x8(%rsp),%xmm1
call   1040 <pow@plt>
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x8(%r13,%rbp,8)
movsd  0x10(%rbx,%rbp,8),%xmm0
movsd  0x8(%rsp),%xmm1
call   1040 <pow@plt>
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%r13,%rbp,8)
movsd  0x18(%rbx,%rbp,8),%xmm0
movsd  0x8(%rsp),%xmm1
call   1040 <pow@plt>
movsd  0x10(%rsp),%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm0,0x18(%r13,%rbp,8)
add    $0x4,%rbp
cmp    %rbp,%r12
jne    1180 <NonlinearADR_Decay_DiracIC_r+0x60>
test   %r15,%r15
je     1249 <NonlinearADR_Decay_DiracIC_r+0x129>
lea    0x0(,%rbp,8),%r12
add    %r13,%r12
lea    (%rbx,%rbp,8),%rbx
xor    %ebp,%ebp
xchg   %ax,%ax
movsd  (%rbx,%rbp,8),%xmm0
movsd  0x8(%rsp),%xmm1
call   1040 <pow@plt>
movsd  0x10(%rsp),%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm0,(%r12,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r15
jne    1220 <NonlinearADR_Decay_DiracIC_r+0x100>
mov    %r14d,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @NonlinearADR_Decay_DiracIC_r(i32* %0, double* %1, i32 %2, double %3, double* %4, double* %5, double* %6) {
%8 = alloca i32*
%9 = alloca double*
%10 = alloca i32
%11 = alloca double
%12 = alloca double*
%13 = alloca double*
%14 = alloca double*
%15 = alloca i32
%16 = alloca i32
%17 = alloca double
%18 = alloca double
store i32* %0, i32** %8
store double* %1, double** %9
store i32 %2, i32* %10
store double %3, double* %11
store double* %4, double** %12
store double* %5, double** %13
store double* %6, double** %14
%19 = load double*, double** %9
%20 = getelementptr inbounds double, double* %19, i64 10
%21 = load double, double* %20
store double %21, double* %17
%22 = load double*, double** %9
%23 = getelementptr inbounds double, double* %22, i64 11
%24 = load double, double* %23
store double %24, double* %18
%25 = load i32*, i32** %8
%26 = load double*, double** %9
%27 = load i32, i32* %10
%28 = load double, double* %11
%29 = load double*, double** %12
%30 = load double*, double** %13
%31 = call i32 @NonlinearADR_Decay_DiracIC(i32* %25, double* %26, i32 %27, double %28, double* %29, double* %30)
store i32 %31, i32* %16
store i32 0, i32* %15
br label %32
32: ; preds = %50, %7
%33 = load i32, i32* %15
%34 = load i32, i32* %10
%35 = icmp slt i32 %33, %34
br i1 %35, label %36, label %53
36: ; preds = %32
%37 = load double, double* %17
%38 = load double*, double** %13
%39 = load i32, i32* %15
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = load double, double* %18
%44 = call double @pow(double %42, double %43)
%45 = fmul double %37, %44
%46 = load double*, double** %14
%47 = load i32, i32* %15
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
store double %45, double* %49
br label %50
50: ; preds = %36
%51 = load i32, i32* %15
%52 = add nsw i32 %51, 1
store i32 %52, i32* %15
br label %32
53: ; preds = %32
%54 = load i32, i32* %16
ret i32 %54
}

LLVM IR (-O2):
define i32 @NonlinearADR_Decay_DiracIC_r(i32* %0, double* %1, i32 %2, double %3, double* %4, double* %5, double* nocapture writeonly %6) {
%8 = getelementptr inbounds double, double* %1, i64 10
%9 = load double, double* %8
%10 = getelementptr inbounds double, double* %1, i64 11
%11 = load double, double* %10
%12 = tail call i32 @NonlinearADR_Decay_DiracIC(i32* %0, double* %1, i32 %2, double %3, double* %4, double* %5)
%13 = icmp sgt i32 %2, 0
br i1 %13, label %14, label %64
14: ; preds = %7
%15 = zext i32 %2 to i64
%16 = add nsw i64 %15, -1
%17 = and i64 %15, 3
%18 = icmp ult i64 %16, 3
br i1 %18, label %50, label %19
19: ; preds = %14
%20 = and i64 %15, 4294967292
br label %21
21: ; preds = %21, %19
%22 = phi i64 [ 0, %19 ], [ %47, %21 ]
%23 = phi i64 [ 0, %19 ], [ %48, %21 ]
%24 = getelementptr inbounds double, double* %5, i64 %22
%25 = load double, double* %24
%26 = tail call double @pow(double %25, double %11)
%27 = fmul double %9, %26
%28 = getelementptr inbounds double, double* %6, i64 %22
store double %27, double* %28
%29 = or i64 %22, 1
%30 = getelementptr inbounds double, double* %5, i64 %29
%31 = load double, double* %30
%32 = tail call double @pow(double %31, double %11)
%33 = fmul double %9, %32
%34 = getelementptr inbounds double, double* %6, i64 %29
store double %33, double* %34
%35 = or i64 %22, 2
%36 = getelementptr inbounds double, double* %5, i64 %35
%37 = load double, double* %36
%38 = tail call double @pow(double %37, double %11)
%39 = fmul double %9, %38
%40 = getelementptr inbounds double, double* %6, i64 %35
store double %39, double* %40
%41 = or i64 %22, 3
%42 = getelementptr inbounds double, double* %5, i64 %41
%43 = load double, double* %42
%44 = tail call double @pow(double %43, double %11)
%45 = fmul double %9, %44
%46 = getelementptr inbounds double, double* %6, i64 %41
store double %45, double* %46
%47 = add nuw nsw i64 %22, 4
%48 = add i64 %23, 4
%49 = icmp eq i64 %48, %20
br i1 %49, label %50, label %21
50: ; preds = %21, %14
%51 = phi i64 [ 0, %14 ], [ %47, %21 ]
%52 = icmp eq i64 %17, 0
br i1 %52, label %64, label %53
53: ; preds = %50, %53
%54 = phi i64 [ %61, %53 ], [ %51, %50 ]
%55 = phi i64 [ %62, %53 ], [ 0, %50 ]
%56 = getelementptr inbounds double, double* %5, i64 %54
%57 = load double, double* %56
%58 = tail call double @pow(double %57, double %11)
%59 = fmul double %9, %58
%60 = getelementptr inbounds double, double* %6, i64 %54
store double %59, double* %60
%61 = add nuw nsw i64 %54, 1
%62 = add i64 %55, 1
%63 = icmp eq i64 %62, %17
br i1 %63, label %64, label %53
64: ; preds = %50, %53, %7
ret i32 %12
}