void init_grid(double **G, int N) {
  int i, j;
  for (i = 0; i < N; i++) {
    for (j = 0; j < N; j++) {
      if (i == 0)
        G[i][j] = 4.56;
      else if (i == N - 1)
        G[i][j] = 9.85;
      else if (j == 0)
        G[i][j] = 7.32;
      else if (j == N - 1)
        G[i][j] = 6.88;
      else
        G[i][j] = 0.0;
    }
  }
}


x86 assembly (clang, -O0):
<init_grid>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    123c <init_grid+0x13c>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1229 <init_grid+0x129>
cmpl   $0x0,-0x10(%rbp)
jne    115d <init_grid+0x5d>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  0xec5(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1216 <init_grid+0x116>
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    1190 <init_grid+0x90>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  0xe8a(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1211 <init_grid+0x111>
cmpl   $0x0,-0x14(%rbp)
jne    11bc <init_grid+0xbc>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  0xe56(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    120c <init_grid+0x10c>
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    11ef <init_grid+0xef>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  0xe1b(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1207 <init_grid+0x107>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    120c <init_grid+0x10c>
jmp    1211 <init_grid+0x111>
jmp    1216 <init_grid+0x116>
jmp    121b <init_grid+0x11b>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1125 <init_grid+0x25>
jmp    122e <init_grid+0x12e>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1112 <init_grid+0x12>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_grid>:
push   %r14
push   %rbx
test   %esi,%esi
jle    11d9 <init_grid+0xd9>
lea    -0x1(%rsi),%eax
mov    %esi,%r14d
xor    %edx,%edx
movabs $0x40123d70a3d70a3d,%r8
3d 12 40
movabs $0x4023b33333333333,%r9
b3 23 40
movabs $0x401b851eb851eb85,%r10
85 1b 40
movabs $0x401d47ae147ae148,%r11
47 1d 40
jmp    114d <init_grid+0x4d>
nopl   (%rax)
add    $0x1,%rdx
cmp    %r14,%rdx
je     11d9 <init_grid+0xd9>
test   %rdx,%rdx
je     1170 <init_grid+0x70>
mov    (%rdi,%rdx,8),%rbx
cmp    %rax,%rdx
jne    1180 <init_grid+0x80>
mov    %r9,(%rbx)
cmp    $0x1,%esi
je     1140 <init_grid+0x40>
jmp    1188 <init_grid+0x88>
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%rdi),%rcx
mov    %r8,(%rcx)
cmp    $0x1,%esi
je     1140 <init_grid+0x40>
jmp    1188 <init_grid+0x88>
nopl   (%rax)
mov    %r11,(%rbx)
cmp    $0x1,%esi
je     1140 <init_grid+0x40>
mov    $0x1,%ebx
jmp    11a0 <init_grid+0xa0>
nop
mov    (%rdi),%rcx
mov    %r8,(%rcx,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%r14
je     1140 <init_grid+0x40>
test   %rdx,%rdx
je     1190 <init_grid+0x90>
cmp    %rax,%rdx
jne    11c0 <init_grid+0xc0>
mov    (%rdi,%rdx,8),%rcx
mov    %r9,(%rcx,%rbx,8)
jmp    1197 <init_grid+0x97>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    (%rdi,%rdx,8),%rcx
cmp    %rbx,%rax
jne    11cf <init_grid+0xcf>
mov    %r10,(%rcx,%rax,8)
jmp    1197 <init_grid+0x97>
movq   $0x0,(%rcx,%rbx,8)
jmp    1197 <init_grid+0x97>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @init_grid(double** %0, i32 %1) {
%3 = alloca double**
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
store double** %0, double*** %3
store i32 %1, i32* %4
store i32 0, i32* %5
br label %7
7: ; preds = %85, %2
%8 = load i32, i32* %5
%9 = load i32, i32* %4
%10 = icmp slt i32 %8, %9
br i1 %10, label %11, label %88
11: ; preds = %7
store i32 0, i32* %6
br label %12
12: ; preds = %81, %11
%13 = load i32, i32* %6
%14 = load i32, i32* %4
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %84
16: ; preds = %12
%17 = load i32, i32* %5
%18 = icmp eq i32 %17, 0
br i1 %18, label %19, label %28
19: ; preds = %16
%20 = load double**, double*** %3
%21 = load i32, i32* %5
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double*, double** %20, i64 %22
%24 = load double*, double** %23
%25 = load i32, i32* %6
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
store double 4.560000e+00, double* %27
br label %80
28: ; preds = %16
%29 = load i32, i32* %5
%30 = load i32, i32* %4
%31 = sub nsw i32 %30, 1
%32 = icmp eq i32 %29, %31
br i1 %32, label %33, label %42
33: ; preds = %28
%34 = load double**, double*** %3
%35 = load i32, i32* %5
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double*, double** %34, i64 %36
%38 = load double*, double** %37
%39 = load i32, i32* %6
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
store double 0x4023B33333333333, double* %41
br label %79
42: ; preds = %28
%43 = load i32, i32* %6
%44 = icmp eq i32 %43, 0
br i1 %44, label %45, label %54
45: ; preds = %42
%46 = load double**, double*** %3
%47 = load i32, i32* %5
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double*, double** %46, i64 %48
%50 = load double*, double** %49
%51 = load i32, i32* %6
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %50, i64 %52
store double 7.320000e+00, double* %53
br label %78
54: ; preds = %42
%55 = load i32, i32* %6
%56 = load i32, i32* %4
%57 = sub nsw i32 %56, 1
%58 = icmp eq i32 %55, %57
br i1 %58, label %59, label %68
59: ; preds = %54
%60 = load double**, double*** %3
%61 = load i32, i32* %5
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double*, double** %60, i64 %62
%64 = load double*, double** %63
%65 = load i32, i32* %6
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %64, i64 %66
store double 6.880000e+00, double* %67
br label %77
68: ; preds = %54
%69 = load double**, double*** %3
%70 = load i32, i32* %5
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double*, double** %69, i64 %71
%73 = load double*, double** %72
%74 = load i32, i32* %6
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %73, i64 %75
store double 0.000000e+00, double* %76
br label %77
77: ; preds = %68, %59
br label %78
78: ; preds = %77, %45
br label %79
79: ; preds = %78, %33
br label %80
80: ; preds = %79, %19
br label %81
81: ; preds = %80
%82 = load i32, i32* %6
%83 = add nsw i32 %82, 1
store i32 %83, i32* %6
br label %12
84: ; preds = %12
br label %85
85: ; preds = %84
%86 = load i32, i32* %5
%87 = add nsw i32 %86, 1
store i32 %87, i32* %5
br label %7
88: ; preds = %7
ret void
}

LLVM IR (-O2):
define void @init_grid(double** nocapture %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %43
4: ; preds = %2
%5 = add nsw i32 %1, -1
%6 = zext i32 %5 to i64
%7 = zext i32 %1 to i64
%8 = icmp eq i32 %1, 1
br label %9
9: ; preds = %40, %4
%10 = phi i64 [ 0, %4 ], [ %41, %40 ]
%11 = icmp eq i64 %10, 0
%12 = icmp eq i64 %10, %6
%13 = getelementptr inbounds double*, double** %0, i64 %10
br i1 %11, label %18, label %14
14: ; preds = %9
%15 = load double*, double** %13
br i1 %12, label %17, label %16
16: ; preds = %14
store double 7.320000e+00, double* %15
br label %20
17: ; preds = %14
store double 0x4023B33333333333, double* %15
br label %20
18: ; preds = %9
%19 = load double*, double** %0
store double 4.560000e+00, double* %19
br label %20
20: ; preds = %18, %17, %16
br i1 %8, label %40, label %21
21: ; preds = %20, %37
%22 = phi i64 [ %38, %37 ], [ 1, %20 ]
br i1 %11, label %23, label %26
23: ; preds = %21
%24 = load double*, double** %0
%25 = getelementptr inbounds double, double* %24, i64 %22
store double 4.560000e+00, double* %25
br label %37
26: ; preds = %21
br i1 %12, label %27, label %30
27: ; preds = %26
%28 = load double*, double** %13
%29 = getelementptr inbounds double, double* %28, i64 %22
store double 0x4023B33333333333, double* %29
br label %37
30: ; preds = %26
%31 = icmp eq i64 %22, %6
%32 = load double*, double** %13
br i1 %31, label %33, label %35
33: ; preds = %30
%34 = getelementptr inbounds double, double* %32, i64 %6
store double 6.880000e+00, double* %34
br label %37
35: ; preds = %30
%36 = getelementptr inbounds double, double* %32, i64 %22
store double 0.000000e+00, double* %36
br label %37
37: ; preds = %23, %35, %33, %27
%38 = add nuw nsw i64 %22, 1
%39 = icmp eq i64 %38, %7
br i1 %39, label %40, label %21
40: ; preds = %37, %20
%41 = add nuw nsw i64 %10, 1
%42 = icmp eq i64 %41, %7
br i1 %42, label %43, label %9
43: ; preds = %40, %2
ret void
}