void init(double *p, int n, int m) {
#pragma omp parallel for
  for (int j = 0; j < n; ++j) {
    for (int i = 0; i < n; i++) {
      if (((i == 0) || (j == 0)) || (i == n - 1) || (j == m - 1))
        p[j * n + i] = 1.;
      else
        p[j * n + i] = 0.;
    }
  }
}


x86 assembly (clang, -O0):
<init>:
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    11cc <init+0xcc>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b9 <init+0xb9>
cmpl   $0x0,-0x18(%rbp)
je     116a <init+0x6a>
cmpl   $0x0,-0x14(%rbp)
je     116a <init+0x6a>
mov    -0x18(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
je     116a <init+0x6a>
mov    -0x14(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    118d <init+0x8d>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xe7d(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11a6 <init+0xa6>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11ab <init+0xab>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1128 <init+0x28>
jmp    11be <init+0xbe>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1115 <init+0x15>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init>:
push   %r15
push   %r14
push   %rbx
test   %esi,%esi
jle    11f9 <init+0xf9>
lea    -0x1(%rsi),%r11d
add    $0xffffffff,%edx
mov    %esi,%r8d
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
lea    0x8(%rdi),%r15
lea    0x0(,%r8,8),%r10
xor    %eax,%eax
movsd  0xecc(%rip),%xmm0
jmp    1156 <init+0x56>
cs nopw 0x0(%rax,%rax,1)
movsd  %xmm1,(%rdi,%r14,8)
add    $0x1,%rax
add    %r10,%r15
cmp    %r8,%rax
je     11f9 <init+0xf9>
xor    %ebx,%ebx
cmp    $0x1,%esi
jne    11b6 <init+0xb6>
test   $0x1,%r8b
je     1146 <init+0x46>
movaps %xmm0,%xmm2
cmp    %rdx,%rax
je     116e <init+0x6e>
xorps  %xmm2,%xmm2
movaps %xmm0,%xmm1
cmp    %r11,%rbx
je     1179 <init+0x79>
movaps %xmm2,%xmm1
mov    %rax,%r14
imul   %r8,%r14
movaps %xmm0,%xmm2
test   %rax,%rax
je     118b <init+0x8b>
movaps %xmm1,%xmm2
add    %rbx,%r14
movaps %xmm0,%xmm1
test   %rbx,%rbx
je     1140 <init+0x40>
movaps %xmm2,%xmm1
jmp    1140 <init+0x40>
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,-0x8(%r15,%rbx,8)
movsd  %xmm1,(%r15,%rbx,8)
add    $0x2,%rbx
cmp    %r9,%rbx
je     115d <init+0x5d>
movaps %xmm0,%xmm1
cmp    %rdx,%rax
je     11c1 <init+0xc1>
xorps  %xmm1,%xmm1
movaps %xmm0,%xmm2
cmp    %r11,%rbx
je     11cc <init+0xcc>
movaps %xmm1,%xmm2
lea    0x1(%rbx),%rcx
movaps %xmm0,%xmm4
cmp    %r11,%rcx
je     11db <init+0xdb>
movaps %xmm1,%xmm4
movaps %xmm0,%xmm3
movaps %xmm0,%xmm1
test   %rax,%rax
je     11ec <init+0xec>
movaps %xmm2,%xmm3
movaps %xmm4,%xmm1
movaps %xmm0,%xmm2
test   %rbx,%rbx
je     11a0 <init+0xa0>
movaps %xmm3,%xmm2
jmp    11a0 <init+0xa0>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @init(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 = %57, %3
%10 = load i32, i32* %7
%11 = load i32, i32* %5
%12 = icmp slt i32 %10, %11
br i1 %12, label %13, label %60
13: ; preds = %9
store i32 0, i32* %8
br label %14
14: ; preds = %53, %13
%15 = load i32, i32* %8
%16 = load i32, i32* %5
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %56
18: ; preds = %14
%19 = load i32, i32* %8
%20 = icmp eq i32 %19, 0
br i1 %20, label %34, label %21
21: ; preds = %18
%22 = load i32, i32* %7
%23 = icmp eq i32 %22, 0
br i1 %23, label %34, label %24
24: ; preds = %21
%25 = load i32, i32* %8
%26 = load i32, i32* %5
%27 = sub nsw i32 %26, 1
%28 = icmp eq i32 %25, %27
br i1 %28, label %34, label %29
29: ; preds = %24
%30 = load i32, i32* %7
%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, %24, %21, %18
%35 = load double*, double** %4
%36 = load i32, i32* %7
%37 = load i32, i32* %5
%38 = mul nsw i32 %36, %37
%39 = load i32, i32* %8
%40 = add nsw i32 %38, %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %35, i64 %41
store double 1.000000e+00, double* %42
br label %52
43: ; preds = %29
%44 = load double*, double** %4
%45 = load i32, i32* %7
%46 = load i32, i32* %5
%47 = mul nsw i32 %45, %46
%48 = load i32, i32* %8
%49 = add nsw i32 %47, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %44, i64 %50
store double 0.000000e+00, double* %51
br label %52
52: ; preds = %43, %34
br label %53
53: ; preds = %52
%54 = load i32, i32* %8
%55 = add nsw i32 %54, 1
store i32 %55, i32* %8
br label %14
56: ; preds = %14
br label %57
57: ; preds = %56
%58 = load i32, i32* %7
%59 = add nsw i32 %58, 1
store i32 %59, i32* %7
br label %9
60: ; preds = %9
ret void
}

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