#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_potential_helper(double *, int, double, double, double,
                                        double *, double *, int, int, double *);
double scf_value(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_potential_helper(&q[0], 1, G, M, r_s, &pars[5], &pars[5 + num_coeff],
                       nmax, lmax, &val[0]);
  _val = val[0];
  return _val;
}


x86 assembly (clang, -O0):
<scf_value>:
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   1040 <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_value+0xf5>
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1202 <scf_value+0xe2>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    11ef <scf_value+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_value+0xa7>
jmp    11f4 <scf_value+0xd4>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11af <scf_value+0x8f>
jmp    1207 <scf_value+0xe7>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1197 <scf_value+0x77>
mov    -0x18(%rbp),%rdi
movsd  -0x28(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm2
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   1030 <scf_potential_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_value>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%r10
mov    %rdi,%r15
cvttsd2si 0x18(%rdi),%r8d
movsd  (%rdi),%xmm0
cvttsd2si 0x20(%rdi),%r9d
movsd  0x8(%rdi),%xmm1
movsd  0x10(%rdi),%xmm2
movq   $0x0,0x10(%rsp)
test   %r8d,%r8d
js     123c <scf_value+0x12c>
lea    0x1(%r9),%esi
lea    -0x7(%r9),%r11d
mov    %r11d,%r14d
shr    $0x3,%r14d
add    $0x1,%r14d
mov    %esi,%edx
and    $0xfffffff8,%edx
movd   %esi,%xmm4
pshufd $0x0,%xmm4,%xmm3
mov    %r14d,%edi
and    $0x7,%edi
and    $0xfffffff8,%r14d
paddd  %xmm4,%xmm4
pshufd $0x0,%xmm4,%xmm4
xor    %ebp,%ebp
xor    %ebx,%ebx
jmp    119e <scf_value+0x8e>
cs nopw 0x0(%rax,%rax,1)
lea    0x1(%rbx),%eax
cmp    %r8d,%ebx
mov    %eax,%ebx
je     123e <scf_value+0x12e>
test   %r9d,%r9d
js     1190 <scf_value+0x80>
xor    %ecx,%ecx
cmp    $0x8,%esi
jb     1221 <scf_value+0x111>
movd   %ebp,%xmm5
pxor   %xmm6,%xmm6
cmp    $0x38,%r11d
jb     11e5 <scf_value+0xd5>
pxor   %xmm6,%xmm6
mov    %r14d,%ebp
nop
paddd  %xmm4,%xmm5
paddd  %xmm4,%xmm6
paddd  %xmm4,%xmm5
paddd  %xmm4,%xmm6
paddd  %xmm4,%xmm5
paddd  %xmm4,%xmm6
paddd  %xmm4,%xmm5
paddd  %xmm4,%xmm6
add    $0xfffffff8,%ebp
jne    11c0 <scf_value+0xb0>
test   %edi,%edi
je     11fd <scf_value+0xed>
mov    %edi,%ebp
nopl   0x0(%rax,%rax,1)
paddd  %xmm3,%xmm5
paddd  %xmm3,%xmm6
add    $0xffffffff,%ebp
jne    11f0 <scf_value+0xe0>
paddd  %xmm6,%xmm5
pshufd $0xee,%xmm5,%xmm6
paddd  %xmm5,%xmm6
pshufd $0x55,%xmm6,%xmm5
paddd  %xmm6,%xmm5
movd   %xmm5,%ebp
mov    %edx,%ecx
cmp    %edx,%esi
je     1190 <scf_value+0x80>
mov    %esi,%eax
sub    %ecx,%eax
cs nopw 0x0(%rax,%rax,1)
nop
add    %esi,%ebp
add    $0xffffffff,%eax
jne    1230 <scf_value+0x120>
jmp    1190 <scf_value+0x80>
xor    %ebp,%ebp
movslq %ebp,%rax
lea    (%r15,%rax,8),%rcx
add    $0x28,%rcx
add    $0x28,%r15
lea    0x10(%rsp),%rax
mov    %rax,(%rsp)
mov    %r10,%rdi
mov    $0x1,%esi
mov    %r15,%rdx
call   1030 <scf_potential_helper@plt>
movsd  0x10(%rsp),%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @scf_value(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* %9
%73 = load double, double* %10
%74 = load double, double* %11
%75 = load double*, double** %6
%76 = getelementptr inbounds double, double* %75, i64 5
%77 = load double*, double** %6
%78 = load i32, i32* %19
%79 = add nsw i32 5, %78
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %77, i64 %80
%82 = load i32, i32* %12
%83 = load i32, i32* %13
%84 = getelementptr inbounds [1 x double], [1 x double]* %14, i64 0, i64 0
%85 = call i32 @scf_potential_helper(double* %71, i32 1, double %72, double %73, double %74, double* %76, double* %81, i32 %82, i32 %83, double* %84)
%86 = getelementptr inbounds [1 x double], [1 x double]* %14, i64 0, i64 0
%87 = load double, double* %86
store double %87, double* %15
%88 = load double, double* %15
ret double %88
}

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