#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 M_PI;
double pow(double, double);
void init_config_BCC(double phi, int nx, int ny, int nz, double *pos,
                     double *lx, double *ly, double *lz) {
  int j;
  int ix, iy, iz;
  (*lx) = (*ly) = (*lz) = pow(8.0 * M_PI / phi / 3.0, 1.0 / 3.0);
  j = 0;
  for (ix = 0; ix < nx; ix++) {
    for (iy = 0; iy < ny; iy++) {
      for (iz = 0; iz < nz; iz++) {
        pos[j * 3 + 0] = (double)ix * (*lx);
        pos[j * 3 + 1] = (double)iy * (*ly);
        pos[j * 3 + 2] = (double)iz * (*lz);
        j++;
        pos[j * 3 + 0] = ((double)ix + 0.5) * (*lx);
        pos[j * 3 + 1] = ((double)iy + 0.5) * (*ly);
        pos[j * 3 + 2] = ((double)iz + 0.5) * (*lz);
        j++;
      }
    }
  }
  (*lx) *= (double)nx;
  (*ly) *= (double)ny;
  (*lz) *= (double)nz;
}


x86 assembly (clang, -O0):
<init_config_BCC>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x10(%rbp),%rax
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  0xed2(%rip),%xmm0
mov    0x2ea3(%rip),%rax
mulsd  (%rax),%xmm0
divsd  -0x8(%rbp),%xmm0
movsd  0xeb2(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xe9e(%rip),%xmm1
call   1030 <pow@plt>
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12e9 <init_config_BCC+0x1d9>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12d6 <init_config_BCC+0x1c6>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12c3 <init_config_BCC+0x1b3>
cvtsi2sdl -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x3c(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x40(%rbp),%xmm0
mov    0x10(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
cvtsi2sdl -0x38(%rbp),%xmm0
movsd  0xde3(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x3c(%rbp),%xmm0
movsd  0xdb7(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x40(%rbp),%xmm0
movsd  0xd8b(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11b3 <init_config_BCC+0xa3>
jmp    12c8 <init_config_BCC+0x1b8>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11a0 <init_config_BCC+0x90>
jmp    12db <init_config_BCC+0x1cb>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118d <init_config_BCC+0x7d>
cvtsi2sdl -0xc(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
cvtsi2sdl -0x10(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
cvtsi2sdl -0x14(%rbp),%xmm0
mov    0x10(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_config_BCC>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,%r13
mov    %r8,%rbp
mov    %rcx,%rbx
mov    %edx,%r15d
mov    %esi,0xc(%rsp)
mov    %edi,%r14d
mov    0x2eb0(%rip),%rax
movsd  (%rax),%xmm2
mulsd  0xebc(%rip),%xmm2
mov    0x50(%rsp),%r12
divsd  %xmm0,%xmm2
divsd  0xeb3(%rip),%xmm2
movsd  0xeb3(%rip),%xmm1
movapd %xmm2,%xmm0
call   1030 <pow@plt>
movsd  %xmm0,(%r12)
movsd  %xmm0,0x0(%r13)
movsd  %xmm0,0x0(%rbp)
mov    %r14d,0x10(%rsp)
test   %r14d,%r14d
jle    12db <init_config_BCC+0x1cb>
movabs $0x600000000,%rax
lea    0x10(%rbx),%r11
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
movsd  0xe77(%rip),%xmm1
movsd  0xe77(%rip),%xmm2
xor    %edi,%edi
jmp    11c5 <init_config_BCC+0xb5>
nopl   (%rax)
addsd  %xmm2,%xmm0
mov    0x14(%rsp),%ecx
add    $0x1,%ecx
cmp    0x10(%rsp),%ecx
je     12d6 <init_config_BCC+0x1c6>
cmpl   $0x0,0xc(%rsp)
mov    %ecx,0x14(%rsp)
jle    11b0 <init_config_BCC+0xa0>
movapd %xmm0,%xmm3
addsd  %xmm1,%xmm3
xor    %esi,%esi
xorpd  %xmm4,%xmm4
jmp    11ed <init_config_BCC+0xdd>
addsd  %xmm2,%xmm4
add    $0x1,%esi
cmp    0xc(%rsp),%esi
je     11b0 <init_config_BCC+0xa0>
test   %r15d,%r15d
jle    11e0 <init_config_BCC+0xd0>
movapd %xmm4,%xmm5
addsd  %xmm1,%xmm5
movslq %edi,%rcx
lea    (%rcx,%rcx,2),%r10
movabs $0x300000000,%r14
imul   %r14,%rcx
movabs $0x500000000,%rdx
lea    (%rcx,%rdx,1),%r8
movabs $0x400000000,%rdx
lea    (%rcx,%rdx,1),%r9
add    %r14,%rcx
lea    (%r11,%r10,8),%r10
xorpd  %xmm6,%xmm6
mov    %r15d,%r14d
nopl   0x0(%rax)
movsd  0x0(%rbp),%xmm7
mulsd  %xmm0,%xmm7
movsd  %xmm7,-0x10(%r10)
movsd  0x0(%r13),%xmm7
mulsd  %xmm4,%xmm7
movsd  %xmm7,-0x8(%r10)
movsd  (%r12),%xmm7
mulsd  %xmm6,%xmm7
movsd  %xmm7,(%r10)
movsd  0x0(%rbp),%xmm7
mulsd  %xmm3,%xmm7
mov    %rcx,%rdx
sar    $0x1d,%rdx
movsd  %xmm7,(%rbx,%rdx,1)
movsd  0x0(%r13),%xmm7
mulsd  %xmm5,%xmm7
mov    %r9,%rdx
sar    $0x1d,%rdx
movsd  %xmm7,(%rbx,%rdx,1)
movapd %xmm6,%xmm7
addsd  %xmm1,%xmm7
mulsd  (%r12),%xmm7
mov    %r8,%rdx
sar    $0x1d,%rdx
movsd  %xmm7,(%rbx,%rdx,1)
addsd  %xmm2,%xmm6
add    %rax,%r8
add    %rax,%r9
add    %rax,%rcx
add    $0x30,%r10
add    $0x2,%edi
add    $0xffffffff,%r14d
jne    1240 <init_config_BCC+0x130>
jmp    11e0 <init_config_BCC+0xd0>
movsd  0x0(%rbp),%xmm0
xorps  %xmm1,%xmm1
cvtsi2sdl 0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0xc(%rsp),%xmm0
mulsd  0x0(%r13),%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r15d,%xmm1
movsd  %xmm0,0x0(%r13)
mulsd  (%r12),%xmm1
movsd  %xmm1,(%r12)
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define void @init_config_BCC(double %0, i32 %1, i32 %2, i32 %3, double* %4, double* %5, double* %6, double* %7) {
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
%14 = alloca double*
%15 = alloca double*
%16 = alloca double*
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
store double %0, double* %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store double* %4, double** %13
store double* %5, double** %14
store double* %6, double** %15
store double* %7, double** %16
%21 = load double, double* @M_PI
%22 = fmul double 8.000000e+00, %21
%23 = load double, double* %9
%24 = fdiv double %22, %23
%25 = fdiv double %24, 3.000000e+00
%26 = call double @pow(double %25, double 0x3FD5555555555555)
%27 = load double*, double** %16
store double %26, double* %27
%28 = load double*, double** %15
store double %26, double* %28
%29 = load double*, double** %14
store double %26, double* %29
store i32 0, i32* %17
store i32 0, i32* %18
br label %30
30: ; preds = %126, %8
%31 = load i32, i32* %18
%32 = load i32, i32* %10
%33 = icmp slt i32 %31, %32
br i1 %33, label %34, label %129
34: ; preds = %30
store i32 0, i32* %19
br label %35
35: ; preds = %122, %34
%36 = load i32, i32* %19
%37 = load i32, i32* %11
%38 = icmp slt i32 %36, %37
br i1 %38, label %39, label %125
39: ; preds = %35
store i32 0, i32* %20
br label %40
40: ; preds = %118, %39
%41 = load i32, i32* %20
%42 = load i32, i32* %12
%43 = icmp slt i32 %41, %42
br i1 %43, label %44, label %121
44: ; preds = %40
%45 = load i32, i32* %18
%46 = sitofp i32 %45 to double
%47 = load double*, double** %14
%48 = load double, double* %47
%49 = fmul double %46, %48
%50 = load double*, double** %13
%51 = load i32, i32* %17
%52 = mul nsw i32 %51, 3
%53 = add nsw i32 %52, 0
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %50, i64 %54
store double %49, double* %55
%56 = load i32, i32* %19
%57 = sitofp i32 %56 to double
%58 = load double*, double** %15
%59 = load double, double* %58
%60 = fmul double %57, %59
%61 = load double*, double** %13
%62 = load i32, i32* %17
%63 = mul nsw i32 %62, 3
%64 = add nsw i32 %63, 1
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %61, i64 %65
store double %60, double* %66
%67 = load i32, i32* %20
%68 = sitofp i32 %67 to double
%69 = load double*, double** %16
%70 = load double, double* %69
%71 = fmul double %68, %70
%72 = load double*, double** %13
%73 = load i32, i32* %17
%74 = mul nsw i32 %73, 3
%75 = add nsw i32 %74, 2
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %72, i64 %76
store double %71, double* %77
%78 = load i32, i32* %17
%79 = add nsw i32 %78, 1
store i32 %79, i32* %17
%80 = load i32, i32* %18
%81 = sitofp i32 %80 to double
%82 = fadd double %81, 5.000000e-01
%83 = load double*, double** %14
%84 = load double, double* %83
%85 = fmul double %82, %84
%86 = load double*, double** %13
%87 = load i32, i32* %17
%88 = mul nsw i32 %87, 3
%89 = add nsw i32 %88, 0
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %86, i64 %90
store double %85, double* %91
%92 = load i32, i32* %19
%93 = sitofp i32 %92 to double
%94 = fadd double %93, 5.000000e-01
%95 = load double*, double** %15
%96 = load double, double* %95
%97 = fmul double %94, %96
%98 = load double*, double** %13
%99 = load i32, i32* %17
%100 = mul nsw i32 %99, 3
%101 = add nsw i32 %100, 1
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %98, i64 %102
store double %97, double* %103
%104 = load i32, i32* %20
%105 = sitofp i32 %104 to double
%106 = fadd double %105, 5.000000e-01
%107 = load double*, double** %16
%108 = load double, double* %107
%109 = fmul double %106, %108
%110 = load double*, double** %13
%111 = load i32, i32* %17
%112 = mul nsw i32 %111, 3
%113 = add nsw i32 %112, 2
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %110, i64 %114
store double %109, double* %115
%116 = load i32, i32* %17
%117 = add nsw i32 %116, 1
store i32 %117, i32* %17
br label %118
118: ; preds = %44
%119 = load i32, i32* %20
%120 = add nsw i32 %119, 1
store i32 %120, i32* %20
br label %40
121: ; preds = %40
br label %122
122: ; preds = %121
%123 = load i32, i32* %19
%124 = add nsw i32 %123, 1
store i32 %124, i32* %19
br label %35
125: ; preds = %35
br label %126
126: ; preds = %125
%127 = load i32, i32* %18
%128 = add nsw i32 %127, 1
store i32 %128, i32* %18
br label %30
129: ; preds = %30
%130 = load i32, i32* %10
%131 = sitofp i32 %130 to double
%132 = load double*, double** %14
%133 = load double, double* %132
%134 = fmul double %133, %131
store double %134, double* %132
%135 = load i32, i32* %11
%136 = sitofp i32 %135 to double
%137 = load double*, double** %15
%138 = load double, double* %137
%139 = fmul double %138, %136
store double %139, double* %137
%140 = load i32, i32* %12
%141 = sitofp i32 %140 to double
%142 = load double*, double** %16
%143 = load double, double* %142
%144 = fmul double %143, %141
store double %144, double* %142
ret void
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define void @init_config_BCC(double %0, i32 %1, i32 %2, i32 %3, double* nocapture writeonly %4, double* nocapture %5, double* nocapture %6, double* nocapture %7) {
%9 = load double, double* @M_PI
%10 = fmul double %9, 8.000000e+00
%11 = fdiv double %10, %0
%12 = fdiv double %11, 3.000000e+00
%13 = tail call double @pow(double %12, double 0x3FD5555555555555)
store double %13, double* %7
store double %13, double* %6
store double %13, double* %5
%14 = icmp sgt i32 %1, 0
br i1 %14, label %15, label %81
15: ; preds = %8
%16 = icmp sgt i32 %2, 0
%17 = icmp sgt i32 %3, 0
br label %18
18: ; preds = %15, %75
%19 = phi i32 [ 0, %15 ], [ %77, %75 ]
%20 = phi i32 [ 0, %15 ], [ %76, %75 ]
br i1 %16, label %21, label %75
21: ; preds = %18
%22 = sitofp i32 %19 to double
%23 = fadd double %22, 5.000000e-01
br label %24
24: ; preds = %21, %71
%25 = phi i32 [ 0, %21 ], [ %73, %71 ]
%26 = phi i32 [ %20, %21 ], [ %72, %71 ]
br i1 %17, label %27, label %71
27: ; preds = %24
%28 = sitofp i32 %25 to double
%29 = fadd double %28, 5.000000e-01
%30 = sext i32 %26 to i64
br label %31
31: ; preds = %27, %31
%32 = phi i64 [ %30, %27 ], [ %66, %31 ]
%33 = phi i32 [ 0, %27 ], [ %67, %31 ]
%34 = load double, double* %5
%35 = fmul double %34, %22
%36 = mul nsw i64 %32, 3
%37 = getelementptr inbounds double, double* %4, i64 %36
store double %35, double* %37
%38 = load double, double* %6
%39 = fmul double %38, %28
%40 = add nsw i64 %36, 1
%41 = getelementptr inbounds double, double* %4, i64 %40
store double %39, double* %41
%42 = sitofp i32 %33 to double
%43 = load double, double* %7
%44 = fmul double %43, %42
%45 = add nsw i64 %36, 2
%46 = getelementptr inbounds double, double* %4, i64 %45
store double %44, double* %46
%47 = load double, double* %5
%48 = fmul double %23, %47
%49 = mul i64 %32, 12884901888
%50 = add i64 %49, 12884901888
%51 = ashr exact i64 %50, 32
%52 = getelementptr inbounds double, double* %4, i64 %51
store double %48, double* %52
%53 = load double, double* %6
%54 = fmul double %29, %53
%55 = mul i64 %32, 12884901888
%56 = add i64 %55, 17179869184
%57 = ashr exact i64 %56, 32
%58 = getelementptr inbounds double, double* %4, i64 %57
store double %54, double* %58
%59 = fadd double %42, 5.000000e-01
%60 = load double, double* %7
%61 = fmul double %59, %60
%62 = mul i64 %32, 12884901888
%63 = add i64 %62, 21474836480
%64 = ashr exact i64 %63, 32
%65 = getelementptr inbounds double, double* %4, i64 %64
store double %61, double* %65
%66 = add nsw i64 %32, 2
%67 = add nuw nsw i32 %33, 1
%68 = icmp eq i32 %67, %3
br i1 %68, label %69, label %31
69: ; preds = %31
%70 = trunc i64 %66 to i32
br label %71
71: ; preds = %69, %24
%72 = phi i32 [ %26, %24 ], [ %70, %69 ]
%73 = add nuw nsw i32 %25, 1
%74 = icmp eq i32 %73, %2
br i1 %74, label %75, label %24
75: ; preds = %71, %18
%76 = phi i32 [ %20, %18 ], [ %72, %71 ]
%77 = add nuw nsw i32 %19, 1
%78 = icmp eq i32 %77, %1
br i1 %78, label %79, label %18
79: ; preds = %75
%80 = load double, double* %5
br label %81
81: ; preds = %79, %8
%82 = phi double [ %80, %79 ], [ %13, %8 ]
%83 = sitofp i32 %1 to double
%84 = fmul double %82, %83
store double %84, double* %5
%85 = sitofp i32 %2 to double
%86 = load double, double* %6
%87 = fmul double %86, %85
store double %87, double* %6
%88 = sitofp i32 %3 to double
%89 = load double, double* %7
%90 = fmul double %89, %88
store double %90, double* %7
ret void
}