void init2d(double **array, int dimX, int dimY) {
  int i, j;
  for (i = 0; i < dimX; i++)
    for (j = 0; j < dimY; j++)
      array[i][j] = (i == 0 || i == dimX - 1 || j == 0 || j == dimY - 1)
                        ? 0.01 * (i + 1) + 0.001 * (j + 1)
                        : 0.0;
}


x86 assembly (clang, -O0):
<init2d>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11ec <init2d+0xec>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11d9 <init2d+0xd9>
cmpl   $0x0,-0x14(%rbp)
je     116a <init2d+0x6a>
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
je     116a <init2d+0x6a>
cmpl   $0x0,-0x18(%rbp)
je     116a <init2d+0x6a>
mov    -0x18(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    11a4 <init2d+0xa4>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
mov    -0x18(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
movsd  0xe82(%rip),%xmm1
mulsd  %xmm2,%xmm1
movsd  0xe6e(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11b1 <init2d+0xb1>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11b1 <init2d+0xb1>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1128 <init2d+0x28>
jmp    11de <init2d+0xde>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1115 <init2d+0x15>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init2d>:
push   %r14
push   %rbx
test   %esi,%esi
jle    11ac <init2d+0xac>
lea    -0x1(%rsi),%r11d
lea    -0x1(%rdx),%r9d
mov    %esi,%r8d
mov    %edx,%r14d
xor    %ecx,%ecx
movsd  0xedd(%rip),%xmm0
movsd  0xedd(%rip),%xmm1
jmp    1138 <init2d+0x38>
nopl   (%rax)
mov    %r10,%rcx
cmp    %r8,%r10
je     11ac <init2d+0xac>
test   %edx,%edx
jle    11a0 <init2d+0xa0>
lea    0x1(%rcx),%r10
xorps  %xmm2,%xmm2
cvtsi2sd %r10d,%xmm2
mov    (%rdi,%rcx,8),%rax
xor    %ebx,%ebx
jmp    1178 <init2d+0x78>
lea    0x1(%rbx),%rsi
xorps  %xmm4,%xmm4
cvtsi2sd %esi,%xmm4
mulsd  %xmm0,%xmm4
movapd %xmm2,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%rax,%rbx,8)
mov    %rsi,%rbx
cmp    %r14,%rsi
je     1130 <init2d+0x30>
test   %rcx,%rcx
je     1150 <init2d+0x50>
cmp    %r11,%rcx
je     1150 <init2d+0x50>
test   %rbx,%rbx
je     1150 <init2d+0x50>
cmp    %r9,%rbx
je     1150 <init2d+0x50>
lea    0x1(%rbx),%rsi
xorpd  %xmm3,%xmm3
jmp    116b <init2d+0x6b>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rcx
mov    %rcx,%r10
cmp    %r8,%r10
jne    1138 <init2d+0x38>
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define void @init2d(double** nocapture %0, i32 %1, i32 %2) {
%4 = icmp sgt i32 %1, 0
br i1 %4, label %5, label %48
5: ; preds = %3
%6 = icmp sgt i32 %2, 0
%7 = add nsw i32 %1, -1
%8 = add nsw i32 %2, -1
%9 = zext i32 %8 to i64
%10 = zext i32 %7 to i64
%11 = zext i32 %1 to i64
%12 = zext i32 %2 to i64
br label %13
13: ; preds = %5, %45
%14 = phi i64 [ 0, %5 ], [ %46, %45 ]
br i1 %6, label %17, label %15
15: ; preds = %13
%16 = add nuw nsw i64 %14, 1
br label %45
17: ; preds = %13
%18 = icmp eq i64 %14, 0
%19 = icmp eq i64 %14, %10
%20 = add nuw nsw i64 %14, 1
%21 = trunc i64 %20 to i32
%22 = sitofp i32 %21 to double
%23 = getelementptr inbounds double*, double** %0, i64 %14
%24 = load double*, double** %23
br label %25
25: ; preds = %17, %40
%26 = phi i64 [ 0, %17 ], [ %41, %40 ]
br i1 %18, label %34, label %27
27: ; preds = %25
%28 = icmp eq i64 %26, 0
%29 = select i1 %19, i1 true, i1 %28
%30 = icmp eq i64 %26, %9
%31 = select i1 %29, i1 true, i1 %30
br i1 %31, label %34, label %32
32: ; preds = %27
%33 = add nuw nsw i64 %26, 1
br label %40
34: ; preds = %27, %25
%35 = add nuw nsw i64 %26, 1
%36 = trunc i64 %35 to i32
%37 = sitofp i32 %36 to double
%38 = fmul double %37, 1.000000e-03
%39 = tail call double @llvm.fmuladd.f64(double %22, double 1.000000e-02, double %38)
br label %40
40: ; preds = %32, %34
%41 = phi i64 [ %33, %32 ], [ %35, %34 ]
%42 = phi double [ 0.000000e+00, %32 ], [ %39, %34 ]
%43 = getelementptr inbounds double, double* %24, i64 %26
store double %42, double* %43
%44 = icmp eq i64 %41, %12
br i1 %44, label %45, label %25
45: ; preds = %40, %15
%46 = phi i64 [ %16, %15 ], [ %20, %40 ]
%47 = icmp eq i64 %46, %11
br i1 %47, label %48, label %13
48: ; preds = %45, %3
ret void
}