#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*/ scf_density_helper(double *, int, double, double, double *,
                                      double *, int, int, double *);
double scf_density(double t, double *pars, double *q, int n_dim) {
  double G = pars[0];
  double M = pars[1];
  double r_s = pars[2];
  int nmax = (int)pars[3];
  int lmax = (int)pars[4];
  double val[1] = {0.};
  double _val;
  int n, l, m;
  int num_coeff = 0;
  for (n = 0; n < (nmax + 1); n++) {
    for (l = 0; l < (lmax + 1); l++) {
      for (m = 0; m < (lmax + 1); m++) {
        num_coeff++;
      }
    }
  }
  scf_density_helper(&q[0], 1, M, r_s, &pars[5], &pars[5 + num_coeff], nmax,
                     lmax, &val[0]);
  _val = val[0];
  return _val;
}


x86 assembly (clang, -O0):
<scf_density>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
cvttsd2si 0x18(%rax),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
cvttsd2si 0x20(%rax),%eax
mov    %eax,-0x40(%rbp)
lea    -0x48(%rbp),%rdi
xor    %esi,%esi
mov    $0x8,%edx
call   1030 <memset@plt>
movl   $0x0,-0x60(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1215 <scf_density+0xf5>
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1202 <scf_density+0xe2>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    11ef <scf_density+0xcf>
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    11c7 <scf_density+0xa7>
jmp    11f4 <scf_density+0xd4>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11af <scf_density+0x8f>
jmp    1207 <scf_density+0xe7>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1197 <scf_density+0x77>
mov    -0x18(%rbp),%rdi
movsd  -0x30(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
mov    -0x10(%rbp),%rdx
add    $0x28,%rdx
mov    -0x10(%rbp),%rcx
mov    -0x60(%rbp),%eax
add    $0x5,%eax
cltq
shl    $0x3,%rax
add    %rax,%rcx
mov    -0x3c(%rbp),%r8d
mov    -0x40(%rbp),%r9d
lea    -0x48(%rbp),%rax
mov    $0x1,%esi
mov    %rax,(%rsp)
call   1040 <scf_density_helper@plt>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<scf_density>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x18,%rsp
mov    %rdi,%r11
cvttsd2si 0x18(%rdi),%r8d
mov    %rsi,%rdi
cvttsd2si 0x20(%r11),%r9d
movsd  0x8(%r11),%xmm0
movsd  0x10(%r11),%xmm1
movq   $0x0,0x10(%rsp)
test   %r8d,%r8d
js     123c <scf_density+0x12c>
lea    0x1(%r9),%ecx
lea    -0x7(%r9),%r10d
mov    %r10d,%r14d
shr    $0x3,%r14d
add    $0x1,%r14d
mov    %ecx,%esi
and    $0xfffffff8,%esi
movd   %ecx,%xmm3
pshufd $0x0,%xmm3,%xmm2
mov    %r14d,%r15d
and    $0x7,%r15d
and    $0xfffffff8,%r14d
paddd  %xmm3,%xmm3
pshufd $0x0,%xmm3,%xmm3
xor    %ebp,%ebp
xor    %ebx,%ebx
jmp    119e <scf_density+0x8e>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0x1(%rbx),%eax
cmp    %r8d,%ebx
mov    %eax,%ebx
je     123e <scf_density+0x12e>
test   %r9d,%r9d
js     1190 <scf_density+0x80>
xor    %edx,%edx
cmp    $0x8,%ecx
jb     1221 <scf_density+0x111>
movd   %ebp,%xmm4
pxor   %xmm5,%xmm5
cmp    $0x38,%r10d
jb     11e5 <scf_density+0xd5>
pxor   %xmm5,%xmm5
mov    %r14d,%ebp
nop
paddd  %xmm3,%xmm4
paddd  %xmm3,%xmm5
paddd  %xmm3,%xmm4
paddd  %xmm3,%xmm5
paddd  %xmm3,%xmm4
paddd  %xmm3,%xmm5
paddd  %xmm3,%xmm4
paddd  %xmm3,%xmm5
add    $0xfffffff8,%ebp
jne    11c0 <scf_density+0xb0>
test   %r15d,%r15d
je     11fd <scf_density+0xed>
mov    %r15d,%ebp
nopl   (%rax)
paddd  %xmm2,%xmm4
paddd  %xmm2,%xmm5
add    $0xffffffff,%ebp
jne    11f0 <scf_density+0xe0>
paddd  %xmm5,%xmm4
pshufd $0xee,%xmm4,%xmm5
paddd  %xmm4,%xmm5
pshufd $0x55,%xmm5,%xmm4
paddd  %xmm5,%xmm4
movd   %xmm4,%ebp
mov    %esi,%edx
cmp    %esi,%ecx
je     1190 <scf_density+0x80>
mov    %ecx,%eax
sub    %edx,%eax
cs nopw 0x0(%rax,%rax,1)
nop
add    %ecx,%ebp
add    $0xffffffff,%eax
jne    1230 <scf_density+0x120>
jmp    1190 <scf_density+0x80>
xor    %ebp,%ebp
movslq %ebp,%rax
lea    (%r11,%rax,8),%rcx
add    $0x28,%rcx
add    $0x28,%r11
lea    0x10(%rsp),%rax
mov    %rax,(%rsp)
mov    $0x1,%esi
mov    %r11,%rdx
call   1030 <scf_density_helper@plt>
movsd  0x10(%rsp),%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @scf_density(double %0, double* %1, double* %2, i32 %3) {
%5 = alloca double
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca i32
%13 = alloca i32
%14 = alloca [1 x double]
%15 = alloca double
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
store double %0, double* %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
%20 = load double*, double** %6
%21 = getelementptr inbounds double, double* %20, i64 0
%22 = load double, double* %21
store double %22, double* %9
%23 = load double*, double** %6
%24 = getelementptr inbounds double, double* %23, i64 1
%25 = load double, double* %24
store double %25, double* %10
%26 = load double*, double** %6
%27 = getelementptr inbounds double, double* %26, i64 2
%28 = load double, double* %27
store double %28, double* %11
%29 = load double*, double** %6
%30 = getelementptr inbounds double, double* %29, i64 3
%31 = load double, double* %30
%32 = fptosi double %31 to i32
store i32 %32, i32* %12
%33 = load double*, double** %6
%34 = getelementptr inbounds double, double* %33, i64 4
%35 = load double, double* %34
%36 = fptosi double %35 to i32
store i32 %36, i32* %13
%37 = bitcast [1 x double]* %14 to i8*
store i32 0, i32* %19
store i32 0, i32* %16
br label %38
38: ; preds = %66, %4
%39 = load i32, i32* %16
%40 = load i32, i32* %12
%41 = add nsw i32 %40, 1
%42 = icmp slt i32 %39, %41
br i1 %42, label %43, label %69
43: ; preds = %38
store i32 0, i32* %17
br label %44
44: ; preds = %62, %43
%45 = load i32, i32* %17
%46 = load i32, i32* %13
%47 = add nsw i32 %46, 1
%48 = icmp slt i32 %45, %47
br i1 %48, label %49, label %65
49: ; preds = %44
store i32 0, i32* %18
br label %50
50: ; preds = %58, %49
%51 = load i32, i32* %18
%52 = load i32, i32* %13
%53 = add nsw i32 %52, 1
%54 = icmp slt i32 %51, %53
br i1 %54, label %55, label %61
55: ; preds = %50
%56 = load i32, i32* %19
%57 = add nsw i32 %56, 1
store i32 %57, i32* %19
br label %58
58: ; preds = %55
%59 = load i32, i32* %18
%60 = add nsw i32 %59, 1
store i32 %60, i32* %18
br label %50
61: ; preds = %50
br label %62
62: ; preds = %61
%63 = load i32, i32* %17
%64 = add nsw i32 %63, 1
store i32 %64, i32* %17
br label %44
65: ; preds = %44
br label %66
66: ; preds = %65
%67 = load i32, i32* %16
%68 = add nsw i32 %67, 1
store i32 %68, i32* %16
br label %38
69: ; preds = %38
%70 = load double*, double** %7
%71 = getelementptr inbounds double, double* %70, i64 0
%72 = load double, double* %10
%73 = load double, double* %11
%74 = load double*, double** %6
%75 = getelementptr inbounds double, double* %74, i64 5
%76 = load double*, double** %6
%77 = load i32, i32* %19
%78 = add nsw i32 5, %77
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %76, i64 %79
%81 = load i32, i32* %12
%82 = load i32, i32* %13
%83 = getelementptr inbounds [1 x double], [1 x double]* %14, i64 0, i64 0
%84 = call i32 @scf_density_helper(double* %71, i32 1, double %72, double %73, double* %75, double* %80, i32 %81, i32 %82, double* %83)
%85 = getelementptr inbounds [1 x double], [1 x double]* %14, i64 0, i64 0
%86 = load double, double* %85
store double %86, double* %15
%87 = load double, double* %15
ret double %87
}

LLVM IR (-O2):
define double @scf_density(double %0, double* %1, double* %2, i32 %3) {
%5 = alloca [1 x double]
%6 = getelementptr inbounds double, double* %1, i64 1
%7 = load double, double* %6
%8 = getelementptr inbounds double, double* %1, i64 2
%9 = load double, double* %8
%10 = getelementptr inbounds double, double* %1, i64 3
%11 = load double, double* %10
%12 = fptosi double %11 to i32
%13 = getelementptr inbounds double, double* %1, i64 4
%14 = load double, double* %13
%15 = fptosi double %14 to i32
%16 = bitcast [1 x double]* %5 to i8*
%17 = bitcast [1 x double]* %5 to i64*
store i64 0, i64* %17
%18 = icmp slt i32 %12, 0
br i1 %18, label %95, label %19
19: ; preds = %4
%20 = icmp slt i32 %15, 0
%21 = add i32 %15, 1
%22 = add i32 %15, -7
%23 = lshr i32 %22, 3
%24 = add nuw nsw i32 %23, 1
%25 = icmp ult i32 %21, 8
%26 = and i32 %21, -8
%27 = insertelement <4 x i32> poison, i32 %21, i64 0
%28 = shufflevector <4 x i32> %27, <4 x i32> poison, <4 x i32> zeroinitializer
%29 = insertelement <4 x i32> poison, i32 %21, i64 0
%30 = shufflevector <4 x i32> %29, <4 x i32> poison, <4 x i32> zeroinitializer
%31 = and i32 %24, 7
%32 = icmp ult i32 %22, 56
%33 = and i32 %24, 1073741816
%34 = shl <4 x i32> %28, <i32 1, i32 1, i32 1, i32 1>
%35 = shl <4 x i32> %30, <i32 1, i32 1, i32 1, i32 1>
%36 = shl <4 x i32> %28, <i32 1, i32 1, i32 1, i32 1>
%37 = shl <4 x i32> %30, <i32 1, i32 1, i32 1, i32 1>
%38 = shl <4 x i32> %28, <i32 1, i32 1, i32 1, i32 1>
%39 = shl <4 x i32> %30, <i32 1, i32 1, i32 1, i32 1>
%40 = shl <4 x i32> %28, <i32 1, i32 1, i32 1, i32 1>
%41 = shl <4 x i32> %30, <i32 1, i32 1, i32 1, i32 1>
%42 = icmp eq i32 %31, 0
%43 = icmp eq i32 %21, %26
br label %44
44: ; preds = %19, %91
%45 = phi i32 [ 0, %19 ], [ %92, %91 ]
%46 = phi i32 [ 0, %19 ], [ %93, %91 ]
br i1 %20, label %91, label %47
47: ; preds = %44
br i1 %25, label %82, label %48
48: ; preds = %47
%49 = insertelement <4 x i32> <i32 poison, i32 0, i32 0, i32 0>, i32 %45, i64 0
br i1 %32, label %64, label %50
50: ; preds = %48, %50
%51 = phi <4 x i32> [ %60, %50 ], [ %49, %48 ]
%52 = phi <4 x i32> [ %61, %50 ], [ zeroinitializer, %48 ]
%53 = phi i32 [ %62, %50 ], [ 0, %48 ]
%54 = add <4 x i32> %51, %34
%55 = add <4 x i32> %52, %35
%56 = add <4 x i32> %54, %36
%57 = add <4 x i32> %55, %37
%58 = add <4 x i32> %56, %38
%59 = add <4 x i32> %57, %39
%60 = add <4 x i32> %58, %40
%61 = add <4 x i32> %59, %41
%62 = add nuw i32 %53, 8
%63 = icmp eq i32 %62, %33
br i1 %63, label %64, label %50
64: ; preds = %50, %48
%65 = phi <4 x i32> [ undef, %48 ], [ %60, %50 ]
%66 = phi <4 x i32> [ undef, %48 ], [ %61, %50 ]
%67 = phi <4 x i32> [ %49, %48 ], [ %60, %50 ]
%68 = phi <4 x i32> [ zeroinitializer, %48 ], [ %61, %50 ]
br i1 %42, label %77, label %69
69: ; preds = %64, %69
%70 = phi <4 x i32> [ %73, %69 ], [ %67, %64 ]
%71 = phi <4 x i32> [ %74, %69 ], [ %68, %64 ]
%72 = phi i32 [ %75, %69 ], [ 0, %64 ]
%73 = add <4 x i32> %28, %70
%74 = add <4 x i32> %30, %71
%75 = add i32 %72, 1
%76 = icmp eq i32 %75, %31
br i1 %76, label %77, label %69
77: ; preds = %69, %64
%78 = phi <4 x i32> [ %65, %64 ], [ %73, %69 ]
%79 = phi <4 x i32> [ %66, %64 ], [ %74, %69 ]
%80 = add <4 x i32> %79, %78
%81 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %80)
br i1 %43, label %91, label %82
82: ; preds = %47, %77
%83 = phi i32 [ %45, %47 ], [ %81, %77 ]
%84 = phi i32 [ 0, %47 ], [ %26, %77 ]
br label %85
85: ; preds = %82, %85
%86 = phi i32 [ %88, %85 ], [ %83, %82 ]
%87 = phi i32 [ %89, %85 ], [ %84, %82 ]
%88 = add i32 %21, %86
%89 = add nuw i32 %87, 1
%90 = icmp eq i32 %87, %15
br i1 %90, label %91, label %85
91: ; preds = %85, %77, %44
%92 = phi i32 [ %45, %44 ], [ %81, %77 ], [ %88, %85 ]
%93 = add nuw i32 %46, 1
%94 = icmp eq i32 %46, %12
br i1 %94, label %95, label %44
95: ; preds = %91, %4
%96 = phi i32 [ 0, %4 ], [ %92, %91 ]
%97 = getelementptr inbounds double, double* %1, i64 5
%98 = add nsw i32 %96, 5
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %1, i64 %99
%101 = getelementptr inbounds [1 x double], [1 x double]* %5, i64 0, i64 0
%102 = call i32 @scf_density_helper(double* %2, i32 1, double %7, double %9, double* %97, double* %100, i32 %12, i32 %15, double* %101)
%103 = load double, double* %101
ret double %103
}