void initialize(double *phi, double *mu, double *V) {
  long i, j, index;
  for (i = 0; i < 100; i++) {
    for (j = 0; j < 50; j++) {
      index = i * 50 + j;
      if ((i - (100 / 2)) * (i - (100 / 2)) + (j - (50 / 2)) * (j - (50 / 2)) <=
          (15.0) * (15.0)) {
        phi[index] = 0;
      } else {
        phi[index] = 1;
      }
      V[index] = 0;
      mu[index] = (-(1.0) / (15.0));
    }
  }
}


x86 assembly (clang, -O0):
<initialize>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
cmpq   $0x64,-0x20(%rbp)
jge    11fb <initialize+0xfb>
movq   $0x0,-0x28(%rbp)
cmpq   $0x32,-0x28(%rbp)
jge    11e5 <initialize+0xe5>
imul   $0x32,-0x20(%rbp),%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
sub    $0x32,%rax
mov    -0x20(%rbp),%rcx
sub    $0x32,%rcx
imul   %rcx,%rax
mov    -0x28(%rbp),%rcx
sub    $0x19,%rcx
mov    -0x28(%rbp),%rdx
sub    $0x19,%rdx
imul   %rdx,%rcx
add    %rcx,%rax
cvtsi2sd %rax,%xmm1
movsd  0xe85(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     119a <initialize+0x9a>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11af <initialize+0xaf>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  0xe5e(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  0xe41(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    112b <initialize+0x2b>
jmp    11ea <initialize+0xea>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1118 <initialize+0x18>
pop    %rbp
ret

x86 assembly (clang, -O2):
<initialize>:
add    $0x8,%rsi
add    $0x8,%rdx
add    $0x8,%rdi
xor    %r8d,%r8d
movsd  0xee9(%rip),%xmm0
movabs $0xbfb1111111111111,%r9
11 b1 bf
jmp    1153 <initialize+0x53>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r8
add    $0x190,%rsi
add    $0x190,%rdx
add    $0x190,%rdi
cmp    $0x64,%r8
je     11e7 <initialize+0xe7>
lea    -0x32(%r8),%r10
imul   %r10,%r10
mov    $0xffffffffffffffe8,%rcx
jmp    1197 <initialize+0x97>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm1,0xc0(%rdi,%rcx,8)
movq   $0x0,0xc0(%rdx,%rcx,8)
mov    %r9,0xc0(%rsi,%rcx,8)
add    $0x2,%rcx
cmp    $0x1a,%rcx
je     1130 <initialize+0x30>
lea    -0x1(%rcx),%rax
imul   %rax,%rax
add    %r10,%rax
xorps  %xmm1,%xmm1
xorps  %xmm2,%xmm2
cmp    $0xe2,%rax
jb     11b3 <initialize+0xb3>
movaps %xmm0,%xmm2
movsd  %xmm2,0xb8(%rdi,%rcx,8)
movq   $0x0,0xb8(%rdx,%rcx,8)
mov    %r9,0xb8(%rsi,%rcx,8)
mov    %rcx,%rax
imul   %rcx,%rax
add    %r10,%rax
cmp    $0xe2,%rax
jb     1170 <initialize+0x70>
movaps %xmm0,%xmm1
jmp    1170 <initialize+0x70>
ret

LLVM IR (-O0):
define void @initialize(double* %0, double* %1, double* %2) {
%4 = alloca double*
%5 = alloca double*
%6 = alloca double*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
store double* %0, double** %4
store double* %1, double** %5
store double* %2, double** %6
store i64 0, i64* %7
br label %10
10: ; preds = %54, %3
%11 = load i64, i64* %7
%12 = icmp slt i64 %11, 100
br i1 %12, label %13, label %57
13: ; preds = %10
store i64 0, i64* %8
br label %14
14: ; preds = %50, %13
%15 = load i64, i64* %8
%16 = icmp slt i64 %15, 50
br i1 %16, label %17, label %53
17: ; preds = %14
%18 = load i64, i64* %7
%19 = mul nsw i64 %18, 50
%20 = load i64, i64* %8
%21 = add nsw i64 %19, %20
store i64 %21, i64* %9
%22 = load i64, i64* %7
%23 = sub nsw i64 %22, 50
%24 = load i64, i64* %7
%25 = sub nsw i64 %24, 50
%26 = mul nsw i64 %23, %25
%27 = load i64, i64* %8
%28 = sub nsw i64 %27, 25
%29 = load i64, i64* %8
%30 = sub nsw i64 %29, 25
%31 = mul nsw i64 %28, %30
%32 = add nsw i64 %26, %31
%33 = sitofp i64 %32 to double
%34 = fcmp ole double %33, 2.250000e+02
br i1 %34, label %35, label %39
35: ; preds = %17
%36 = load double*, double** %4
%37 = load i64, i64* %9
%38 = getelementptr inbounds double, double* %36, i64 %37
store double 0.000000e+00, double* %38
br label %43
39: ; preds = %17
%40 = load double*, double** %4
%41 = load i64, i64* %9
%42 = getelementptr inbounds double, double* %40, i64 %41
store double 1.000000e+00, double* %42
br label %43
43: ; preds = %39, %35
%44 = load double*, double** %6
%45 = load i64, i64* %9
%46 = getelementptr inbounds double, double* %44, i64 %45
store double 0.000000e+00, double* %46
%47 = load double*, double** %5
%48 = load i64, i64* %9
%49 = getelementptr inbounds double, double* %47, i64 %48
store double 0xBFB1111111111111, double* %49
br label %50
50: ; preds = %43
%51 = load i64, i64* %8
%52 = add nsw i64 %51, 1
store i64 %52, i64* %8
br label %14
53: ; preds = %14
br label %54
54: ; preds = %53
%55 = load i64, i64* %7
%56 = add nsw i64 %55, 1
store i64 %56, i64* %7
br label %10
57: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @initialize(double* nocapture writeonly %0, double* nocapture writeonly %1, double* nocapture writeonly %2) {
br label %4
4: ; preds = %3, %32
%5 = phi i64 [ 0, %3 ], [ %33, %32 ]
%6 = mul nuw nsw i64 %5, 50
%7 = add nsw i64 %5, -50
%8 = mul nsw i64 %7, %7
br label %9
9: ; preds = %9, %4
%10 = phi i64 [ 0, %4 ], [ %30, %9 ]
%11 = add nuw nsw i64 %10, %6
%12 = add nsw i64 %10, -25
%13 = mul nsw i64 %12, %12
%14 = add nuw nsw i64 %13, %8
%15 = icmp ult i64 %14, 226
%16 = select i1 %15, double 0.000000e+00, double 1.000000e+00
%17 = getelementptr inbounds double, double* %0, i64 %11
store double %16, double* %17
%18 = getelementptr inbounds double, double* %2, i64 %11
store double 0.000000e+00, double* %18
%19 = getelementptr inbounds double, double* %1, i64 %11
store double 0xBFB1111111111111, double* %19
%20 = or i64 %10, 1
%21 = add nuw nsw i64 %20, %6
%22 = add nsw i64 %10, -24
%23 = mul nsw i64 %22, %22
%24 = add nuw nsw i64 %23, %8
%25 = icmp ult i64 %24, 226
%26 = select i1 %25, double 0.000000e+00, double 1.000000e+00
%27 = getelementptr inbounds double, double* %0, i64 %21
store double %26, double* %27
%28 = getelementptr inbounds double, double* %2, i64 %21
store double 0.000000e+00, double* %28
%29 = getelementptr inbounds double, double* %1, i64 %21
store double 0xBFB1111111111111, double* %29
%30 = add nuw nsw i64 %10, 2
%31 = icmp eq i64 %30, 50
br i1 %31, label %32, label %9
32: ; preds = %9
%33 = add nuw nsw i64 %5, 1
%34 = icmp eq i64 %33, 100
br i1 %34, label %35, label %4
35: ; preds = %32
ret void
}