void init_positions(double **coord, int nAtoms, double *box) {
  double cbrt(double x);
  int iBoxD;
  double fBoxD;
  int x, y, z;
  double xPos, yPos, zPos;
  int atomCount;
  iBoxD = (int)cbrt((double)nAtoms);
  if (iBoxD * iBoxD * iBoxD < nAtoms) {
    iBoxD++;
  }
  fBoxD = 3.55;
  *box = iBoxD * fBoxD;
  atomCount = 0;
  for (x = 0; x < iBoxD; x++) {
    xPos = (x + 0.5) * fBoxD;
    for (y = 0; y < iBoxD; y++) {
      yPos = (y + 0.5) * fBoxD;
      for (z = 0; z < iBoxD; z++) {
        if (atomCount < nAtoms) {
          zPos = (z + 0.5) * fBoxD;
          coord[atomCount][0] = xPos;
          coord[atomCount][1] = yPos;
          coord[atomCount][2] = zPos;
          atomCount++;
        } else {
          break;
        }
      }
      if (atomCount >= nAtoms) {
        break;
      }
    }
  }
}


x86 assembly (clang, -O0):
<init_positions>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
call   1030 <cbrt@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
imul   -0x1c(%rbp),%eax
imul   -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1151 <init_positions+0x41>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
movsd  0xea7(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x54(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12ab <init_positions+0x19b>
cvtsi2sdl -0x2c(%rbp),%xmm0
movsd  0xe71(%rip),%xmm1
addsd  %xmm1,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1298 <init_positions+0x188>
cvtsi2sdl -0x30(%rbp),%xmm0
movsd  0xe43(%rip),%xmm1
addsd  %xmm1,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1274 <init_positions+0x164>
mov    -0x54(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    125c <init_positions+0x14c>
cvtsi2sdl -0x34(%rbp),%xmm0
movsd  0xe09(%rip),%xmm1
addsd  %xmm1,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  %xmm0,(%rax)
movsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1261 <init_positions+0x151>
jmp    1274 <init_positions+0x164>
jmp    1266 <init_positions+0x156>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11da <init_positions+0xca>
mov    -0x54(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     1285 <init_positions+0x175>
jmp    1298 <init_positions+0x188>
jmp    128a <init_positions+0x17a>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11ac <init_positions+0x9c>
jmp    129d <init_positions+0x18d>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    117e <init_positions+0x6e>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_positions>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %rdx,%r14
mov    %esi,%r15d
mov    %rdi,%rbx
cvtsi2sd %esi,%xmm0
call   1030 <cbrt@plt>
cvttsd2si %xmm0,%ecx
mov    %ecx,%edx
imul   %ecx,%edx
imul   %ecx,%edx
xor    %eax,%eax
cmp    %r15d,%edx
setl   %al
add    %ecx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xeb2(%rip),%xmm0
movsd  %xmm0,(%r14)
test   %eax,%eax
jle    1230 <init_positions+0x120>
movslq %r15d,%rcx
xorpd  %xmm0,%xmm0
xor    %edx,%edx
movsd  0xe9c(%rip),%xmm1
movsd  0xe8c(%rip),%xmm2
movsd  0xe94(%rip),%xmm3
xor    %r8d,%r8d
jmp    11a1 <init_positions+0x91>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  %xmm3,%xmm0
add    $0x1,%r8d
cmp    %eax,%r8d
je     1230 <init_positions+0x120>
movapd %xmm0,%xmm4
addsd  %xmm1,%xmm4
mulsd  %xmm2,%xmm4
xor    %r9d,%r9d
xorps  %xmm5,%xmm5
cvtsi2sd %r9d,%xmm5
addsd  %xmm1,%xmm5
mulsd  %xmm2,%xmm5
cmp    %r15d,%edx
jge    1190 <init_positions+0x80>
movslq %edx,%rdi
add    $0x1,%rdi
xorpd  %xmm6,%xmm6
mov    $0x1,%esi
cs nopw 0x0(%rax,%rax,1)
nop
mov    %rdi,%rbp
movapd %xmm6,%xmm7
addsd  %xmm1,%xmm7
mulsd  %xmm2,%xmm7
mov    -0x8(%rbx,%rdi,8),%rdi
movsd  %xmm4,(%rdi)
movsd  %xmm5,0x8(%rdi)
movsd  %xmm7,0x10(%rdi)
addsd  %xmm3,%xmm6
add    $0x1,%edx
cmp    %eax,%esi
jge    1219 <init_positions+0x109>
lea    0x1(%rbp),%rdi
add    $0x1,%esi
cmp    %rcx,%rbp
jl     11e0 <init_positions+0xd0>
cmp    %rcx,%rbp
jge    1190 <init_positions+0x80>
add    $0x1,%r9d
cmp    %eax,%r9d
jl     11b0 <init_positions+0xa0>
jmp    1190 <init_positions+0x80>
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @init_positions(double** %0, i32 %1, double* %2) {
%4 = alloca double**
%5 = alloca i32
%6 = alloca double*
%7 = alloca i32
%8 = alloca double
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca i32
store double** %0, double*** %4
store i32 %1, i32* %5
store double* %2, double** %6
%16 = load i32, i32* %5
%17 = sitofp i32 %16 to double
%18 = call double @cbrt(double %17)
%19 = fptosi double %18 to i32
store i32 %19, i32* %7
%20 = load i32, i32* %7
%21 = load i32, i32* %7
%22 = mul nsw i32 %20, %21
%23 = load i32, i32* %7
%24 = mul nsw i32 %22, %23
%25 = load i32, i32* %5
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %30
27: ; preds = %3
%28 = load i32, i32* %7
%29 = add nsw i32 %28, 1
store i32 %29, i32* %7
br label %30
30: ; preds = %27, %3
store double 3.550000e+00, double* %8
%31 = load i32, i32* %7
%32 = sitofp i32 %31 to double
%33 = load double, double* %8
%34 = fmul double %32, %33
%35 = load double*, double** %6
store double %34, double* %35
store i32 0, i32* %15
store i32 0, i32* %9
br label %36
36: ; preds = %108, %30
%37 = load i32, i32* %9
%38 = load i32, i32* %7
%39 = icmp slt i32 %37, %38
br i1 %39, label %40, label %111
40: ; preds = %36
%41 = load i32, i32* %9
%42 = sitofp i32 %41 to double
%43 = fadd double %42, 5.000000e-01
%44 = load double, double* %8
%45 = fmul double %43, %44
store double %45, double* %12
store i32 0, i32* %10
br label %46
46: ; preds = %104, %40
%47 = load i32, i32* %10
%48 = load i32, i32* %7
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %107
50: ; preds = %46
%51 = load i32, i32* %10
%52 = sitofp i32 %51 to double
%53 = fadd double %52, 5.000000e-01
%54 = load double, double* %8
%55 = fmul double %53, %54
store double %55, double* %13
store i32 0, i32* %11
br label %56
56: ; preds = %95, %50
%57 = load i32, i32* %11
%58 = load i32, i32* %7
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %98
60: ; preds = %56
%61 = load i32, i32* %15
%62 = load i32, i32* %5
%63 = icmp slt i32 %61, %62
br i1 %63, label %64, label %93
64: ; preds = %60
%65 = load i32, i32* %11
%66 = sitofp i32 %65 to double
%67 = fadd double %66, 5.000000e-01
%68 = load double, double* %8
%69 = fmul double %67, %68
store double %69, double* %14
%70 = load double, double* %12
%71 = load double**, double*** %4
%72 = load i32, i32* %15
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double*, double** %71, i64 %73
%75 = load double*, double** %74
%76 = getelementptr inbounds double, double* %75, i64 0
store double %70, double* %76
%77 = load double, double* %13
%78 = load double**, double*** %4
%79 = load i32, i32* %15
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double*, double** %78, i64 %80
%82 = load double*, double** %81
%83 = getelementptr inbounds double, double* %82, i64 1
store double %77, double* %83
%84 = load double, double* %14
%85 = load double**, double*** %4
%86 = load i32, i32* %15
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double*, double** %85, i64 %87
%89 = load double*, double** %88
%90 = getelementptr inbounds double, double* %89, i64 2
store double %84, double* %90
%91 = load i32, i32* %15
%92 = add nsw i32 %91, 1
store i32 %92, i32* %15
br label %94
93: ; preds = %60
br label %98
94: ; preds = %64
br label %95
95: ; preds = %94
%96 = load i32, i32* %11
%97 = add nsw i32 %96, 1
store i32 %97, i32* %11
br label %56
98: ; preds = %93, %56
%99 = load i32, i32* %15
%100 = load i32, i32* %5
%101 = icmp sge i32 %99, %100
br i1 %101, label %102, label %103
102: ; preds = %98
br label %107
103: ; preds = %98
br label %104
104: ; preds = %103
%105 = load i32, i32* %10
%106 = add nsw i32 %105, 1
store i32 %106, i32* %10
br label %46
107: ; preds = %102, %46
br label %108
108: ; preds = %107
%109 = load i32, i32* %9
%110 = add nsw i32 %109, 1
store i32 %110, i32* %9
br label %36
111: ; preds = %36
ret void
}

LLVM IR (-O2):
define void @init_positions(double** nocapture %0, i32 %1, double* nocapture writeonly %2) {
%4 = sitofp i32 %1 to double
%5 = tail call double @cbrt(double %4)
%6 = fptosi double %5 to i32
%7 = mul nsw i32 %6, %6
%8 = mul nsw i32 %7, %6
%9 = icmp slt i32 %8, %1
%10 = zext i1 %9 to i32
%11 = add i32 %10, %6
%12 = sitofp i32 %11 to double
%13 = fmul double %12, 3.550000e+00
store double %13, double* %2
%14 = icmp sgt i32 %11, 0
br i1 %14, label %15, label %58
15: ; preds = %3
%16 = sext i32 %1 to i64
br label %17
17: ; preds = %54, %15
%18 = phi i32 [ 0, %15 ], [ %55, %54 ]
%19 = phi i32 [ 0, %15 ], [ %56, %54 ]
%20 = sitofp i32 %19 to double
%21 = fadd double %20, 5.000000e-01
%22 = fmul double %21, 3.550000e+00
br label %23
23: ; preds = %17, %47
%24 = phi i32 [ %48, %47 ], [ %18, %17 ]
%25 = phi i32 [ %49, %47 ], [ 0, %17 ]
%26 = sitofp i32 %25 to double
%27 = fadd double %26, 5.000000e-01
%28 = fmul double %27, 3.550000e+00
%29 = icmp slt i32 %24, %1
br i1 %29, label %30, label %54
30: ; preds = %23
%31 = sext i32 %24 to i64
br label %32
32: ; preds = %30, %32
%33 = phi i64 [ %31, %30 ], [ %42, %32 ]
%34 = phi i32 [ 0, %30 ], [ %43, %32 ]
%35 = sitofp i32 %34 to double
%36 = fadd double %35, 5.000000e-01
%37 = fmul double %36, 3.550000e+00
%38 = getelementptr inbounds double*, double** %0, i64 %33
%39 = load double*, double** %38
store double %22, double* %39
%40 = getelementptr inbounds double, double* %39, i64 1
store double %28, double* %40
%41 = getelementptr inbounds double, double* %39, i64 2
store double %37, double* %41
%42 = add nsw i64 %33, 1
%43 = add nuw nsw i32 %34, 1
%44 = icmp slt i32 %43, %11
%45 = icmp slt i64 %42, %16
%46 = select i1 %44, i1 %45, i1 false
br i1 %46, label %32, label %47
47: ; preds = %32
%48 = trunc i64 %42 to i32
%49 = add nuw nsw i32 %25, 1
%50 = icmp slt i32 %49, %11
%51 = select i1 %45, i1 %50, i1 false
br i1 %51, label %23, label %52
52: ; preds = %47
%53 = trunc i64 %42 to i32
br label %54
54: ; preds = %23, %52
%55 = phi i32 [ %53, %52 ], [ %24, %23 ]
%56 = add nuw nsw i32 %19, 1
%57 = icmp eq i32 %56, %11
br i1 %57, label %58, label %17
58: ; preds = %54, %3
ret void
}