void init(double a[128][128], double b[128][128], int n) {
  int i, j, s = -1;
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      a[i][j] = i * j * s;
      b[i][j] = i + j + s;
      s = -s;
    }
  }
}


x86 assembly (clang, -O0):
<init>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0xffffffff,-0x20(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11b1 <init+0xb1>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    119e <init+0x9e>
mov    -0x18(%rbp),%eax
imul   -0x1c(%rbp),%eax
imul   -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
shl    $0xa,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    -0x1c(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
shl    $0xa,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
xor    %eax,%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1130 <init+0x30>
jmp    11a3 <init+0xa3>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    111d <init+0x1d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init>:
test   %edx,%edx
jle    116b <init+0x6b>
mov    %edx,%r8d
mov    $0xffffffff,%r10d
xor    %r9d,%r9d
xor    %ecx,%ecx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %ecx,%eax
imul   %r10d,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rdi,%rdx,8)
lea    (%r9,%rdx,1),%eax
add    %r10d,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
neg    %r10d
add    $0x1,%rdx
add    %r9d,%ecx
cmp    %rdx,%r8
jne    1120 <init+0x20>
add    $0x1,%r9
add    $0x400,%rsi
add    $0x400,%rdi
cmp    %r8,%r9
jne    1110 <init+0x10>
ret

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

LLVM IR (-O2):
define void @init([128 x double]* nocapture writeonly %0, [128 x double]* nocapture writeonly %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %30
5: ; preds = %3
%6 = zext i32 %2 to i64
br label %7
7: ; preds = %27, %5
%8 = phi i64 [ 0, %5 ], [ %28, %27 ]
%9 = phi i32 [ -1, %5 ], [ %24, %27 ]
%10 = trunc i64 %8 to i32
br label %11
11: ; preds = %7, %11
%12 = phi i64 [ 0, %7 ], [ %25, %11 ]
%13 = phi i32 [ %9, %7 ], [ %24, %11 ]
%14 = trunc i64 %12 to i32
%15 = mul nsw i32 %14, %10
%16 = mul nsw i32 %15, %13
%17 = sitofp i32 %16 to double
%18 = getelementptr inbounds [128 x double], [128 x double]* %0, i64 %8, i64 %12
store double %17, double* %18
%19 = add nuw nsw i64 %12, %8
%20 = trunc i64 %19 to i32
%21 = add nsw i32 %13, %20
%22 = sitofp i32 %21 to double
%23 = getelementptr inbounds [128 x double], [128 x double]* %1, i64 %8, i64 %12
store double %22, double* %23
%24 = sub nsw i32 0, %13
%25 = add nuw nsw i64 %12, 1
%26 = icmp eq i64 %25, %6
br i1 %26, label %27, label %11
27: ; preds = %11
%28 = add nuw nsw i64 %8, 1
%29 = icmp eq i64 %28, %6
br i1 %29, label %30, label %7
30: ; preds = %27, %3
ret void
}