void initGrid(float u0[], float u1[], int w) {
  int i, j;
  for (i = 0; i < w; i++) {
    for (j = 0; j < w; j++) {
      u0[(i) * (w) + j] = 0.;
      u1[(i) * (w) + j] = 0.;
    }
  }
}


x86 assembly (clang, -O0):
<initGrid>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1188 <initGrid+0x88>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1175 <initGrid+0x75>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1129 <initGrid+0x29>
jmp    117a <initGrid+0x7a>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1116 <initGrid+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<initGrid>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    1394 <initGrid+0x294>
mov    %rsi,%rbx
mov    %edx,%r14d
mov    %r14d,%eax
and    $0xfffffff8,%eax
mov    %rax,-0x18(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x20(%rsp)
shr    $0x3,%rax
add    $0x1,%rax
mov    %eax,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%rax
mov    %rax,-0x30(%rsp)
mov    %r14d,%eax
and    $0x3,%eax
mov    %rax,-0x40(%rsp)
lea    0x70(%rbx),%r9
lea    0x0(,%r14,4),%r15
lea    0x70(%rdi),%rcx
lea    0x10(%rbx),%r12
lea    0x10(%rdi),%r13
mov    %rsi,-0x28(%rsp)
mov    %rsi,%r11
shl    $0x5,%r11
xor    %r8d,%r8d
xorps  %xmm0,%xmm0
mov    %rdi,%rsi
mov    %rbx,%r10
mov    %rdi,-0x8(%rsp)
mov    %r11,-0x10(%rsp)
mov    %edx,-0x44(%rsp)
mov    %rbx,-0x38(%rsp)
jmp    11c9 <initGrid+0xc9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r8
add    %r15,%r9
add    %r15,%rcx
add    %r15,%r12
add    %r15,%r13
add    %r15,%r10
add    %r15,%rsi
cmp    %r14,%r8
mov    -0x8(%rsp),%rdi
mov    -0x10(%rsp),%r11
je     1394 <initGrid+0x294>
cmp    $0x8,%edx
jae    11e0 <initGrid+0xe0>
xor    %r11d,%r11d
mov    -0x40(%rsp),%rdi
jmp    12f0 <initGrid+0x1f0>
nopl   0x0(%rax,%rax,1)
mov    %r8,%rax
imul   %r14,%rax
lea    (%rax,%r14,1),%rbp
mov    -0x38(%rsp),%rdx
lea    (%rdx,%rbp,4),%rbx
lea    (%rdi,%rax,4),%rdx
cmp    %rbx,%rdx
mov    -0x38(%rsp),%rbx
jae    1220 <initGrid+0x120>
lea    (%rbx,%rax,4),%rax
lea    (%rdi,%rbp,4),%rdx
cmp    %rdx,%rax
jae    1220 <initGrid+0x120>
xor    %r11d,%r11d
mov    -0x44(%rsp),%edx
mov    -0x40(%rsp),%rdi
jmp    12f0 <initGrid+0x1f0>
cmpq   $0x18,-0x20(%rsp)
jae    1230 <initGrid+0x130>
xor    %ebp,%ebp
mov    -0x44(%rsp),%edx
jmp    12a0 <initGrid+0x1a0>
mov    -0x30(%rsp),%rax
xor    %ebp,%ebp
mov    -0x44(%rsp),%edx
nopl   0x0(%rax,%rax,1)
movups %xmm0,-0x70(%rcx,%rbp,4)
movups %xmm0,-0x60(%rcx,%rbp,4)
movups %xmm0,-0x70(%r9,%rbp,4)
movups %xmm0,-0x60(%r9,%rbp,4)
movups %xmm0,-0x50(%rcx,%rbp,4)
movups %xmm0,-0x40(%rcx,%rbp,4)
movups %xmm0,-0x50(%r9,%rbp,4)
movups %xmm0,-0x40(%r9,%rbp,4)
movups %xmm0,-0x30(%rcx,%rbp,4)
movups %xmm0,-0x20(%rcx,%rbp,4)
movups %xmm0,-0x30(%r9,%rbp,4)
movups %xmm0,-0x20(%r9,%rbp,4)
movups %xmm0,-0x10(%rcx,%rbp,4)
movups %xmm0,(%rcx,%rbp,4)
movups %xmm0,-0x10(%r9,%rbp,4)
movups %xmm0,(%r9,%rbp,4)
add    $0x20,%rbp
add    $0xfffffffffffffffc,%rax
jne    1240 <initGrid+0x140>
cmpq   $0x0,-0x28(%rsp)
mov    -0x40(%rsp),%rdi
je     12dc <initGrid+0x1dc>
lea    (%r12,%rbp,4),%rax
lea    0x0(,%rbp,4),%rbp
add    %r13,%rbp
xor    %ebx,%ebx
xchg   %ax,%ax
movups %xmm0,-0x10(%rbp,%rbx,1)
movups %xmm0,0x0(%rbp,%rbx,1)
movups %xmm0,-0x10(%rax,%rbx,1)
movups %xmm0,(%rax,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%r11
jne    12c0 <initGrid+0x1c0>
mov    -0x18(%rsp),%rax
mov    %rax,%r11
cmp    %r14,%rax
je     11a0 <initGrid+0xa0>
nopl   (%rax)
mov    %r11,%rax
not    %rax
add    %r14,%rax
test   %rdi,%rdi
je     132a <initGrid+0x22a>
mov    %rdi,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movl   $0x0,(%rsi,%r11,4)
movl   $0x0,(%r10,%r11,4)
add    $0x1,%r11
add    $0xffffffffffffffff,%rbp
jne    1310 <initGrid+0x210>
cmp    $0x3,%rax
jb     11a0 <initGrid+0xa0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movl   $0x0,(%rsi,%r11,4)
movl   $0x0,(%r10,%r11,4)
movl   $0x0,0x4(%rsi,%r11,4)
movl   $0x0,0x4(%r10,%r11,4)
movl   $0x0,0x8(%rsi,%r11,4)
movl   $0x0,0x8(%r10,%r11,4)
movl   $0x0,0xc(%rsi,%r11,4)
movl   $0x0,0xc(%r10,%r11,4)
add    $0x4,%r11
cmp    %r11,%r14
jne    1340 <initGrid+0x240>
jmp    11a0 <initGrid+0xa0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @initGrid(float* %0, float* %1, i32 %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
store float* %0, float** %4
store float* %1, float** %5
store i32 %2, i32* %6
store i32 0, i32* %7
br label %9
9: ; preds = %39, %3
%10 = load i32, i32* %7
%11 = load i32, i32* %6
%12 = icmp slt i32 %10, %11
br i1 %12, label %13, label %42
13: ; preds = %9
store i32 0, i32* %8
br label %14
14: ; preds = %35, %13
%15 = load i32, i32* %8
%16 = load i32, i32* %6
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %38
18: ; preds = %14
%19 = load float*, float** %4
%20 = load i32, i32* %7
%21 = load i32, i32* %6
%22 = mul nsw i32 %20, %21
%23 = load i32, i32* %8
%24 = add nsw i32 %22, %23
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds float, float* %19, i64 %25
store float 0.000000e+00, float* %26
%27 = load float*, float** %5
%28 = load i32, i32* %7
%29 = load i32, i32* %6
%30 = mul nsw i32 %28, %29
%31 = load i32, i32* %8
%32 = add nsw i32 %30, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %27, i64 %33
store float 0.000000e+00, float* %34
br label %35
35: ; preds = %18
%36 = load i32, i32* %8
%37 = add nsw i32 %36, 1
store i32 %37, i32* %8
br label %14
38: ; preds = %14
br label %39
39: ; preds = %38
%40 = load i32, i32* %7
%41 = add nsw i32 %40, 1
store i32 %41, i32* %7
br label %9
42: ; preds = %9
ret void
}

LLVM IR (-O2):
define void @initGrid(float* nocapture writeonly %0, float* nocapture writeonly %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %136
5: ; preds = %3
%6 = zext i32 %2 to i64
%7 = zext i32 %2 to i64
%8 = and i64 %7, 4294967288
%9 = add nsw i64 %8, -8
%10 = lshr exact i64 %9, 3
%11 = add nuw nsw i64 %10, 1
%12 = icmp ult i32 %2, 8
%13 = and i64 %7, 4294967288
%14 = and i64 %11, 3
%15 = icmp ult i64 %9, 24
%16 = and i64 %11, 4611686018427387900
%17 = icmp eq i64 %14, 0
%18 = icmp eq i64 %13, %7
%19 = and i64 %7, 3
%20 = icmp eq i64 %19, 0
br label %21
21: ; preds = %133, %5
%22 = phi i64 [ 0, %5 ], [ %134, %133 ]
%23 = mul i64 %22, %7
%24 = add i64 %23, %7
%25 = mul nsw i64 %22, %6
br i1 %12, label %98, label %26
26: ; preds = %21
%27 = getelementptr float, float* %1, i64 %24
%28 = getelementptr float, float* %1, i64 %23
%29 = getelementptr float, float* %0, i64 %24
%30 = getelementptr float, float* %0, i64 %23
%31 = icmp ult float* %30, %27
%32 = icmp ult float* %28, %29
%33 = and i1 %31, %32
br i1 %33, label %98, label %34
34: ; preds = %26
br i1 %15, label %80, label %35
35: ; preds = %34, %35
%36 = phi i64 [ %77, %35 ], [ 0, %34 ]
%37 = phi i64 [ %78, %35 ], [ 0, %34 ]
%38 = add nuw nsw i64 %36, %25
%39 = getelementptr inbounds float, float* %0, i64 %38
%40 = bitcast float* %39 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %40
%41 = getelementptr inbounds float, float* %39, i64 4
%42 = bitcast float* %41 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %42
%43 = getelementptr inbounds float, float* %1, i64 %38
%44 = bitcast float* %43 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %44
%45 = getelementptr inbounds float, float* %43, i64 4
%46 = bitcast float* %45 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %46
%47 = or i64 %36, 8
%48 = add nuw nsw i64 %47, %25
%49 = getelementptr inbounds float, float* %0, i64 %48
%50 = bitcast float* %49 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %50
%51 = getelementptr inbounds float, float* %49, i64 4
%52 = bitcast float* %51 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %52
%53 = getelementptr inbounds float, float* %1, i64 %48
%54 = bitcast float* %53 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %54
%55 = getelementptr inbounds float, float* %53, i64 4
%56 = bitcast float* %55 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %56
%57 = or i64 %36, 16
%58 = add nuw nsw i64 %57, %25
%59 = getelementptr inbounds float, float* %0, i64 %58
%60 = bitcast float* %59 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %60
%61 = getelementptr inbounds float, float* %59, i64 4
%62 = bitcast float* %61 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %62
%63 = getelementptr inbounds float, float* %1, i64 %58
%64 = bitcast float* %63 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %64
%65 = getelementptr inbounds float, float* %63, i64 4
%66 = bitcast float* %65 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %66
%67 = or i64 %36, 24
%68 = add nuw nsw i64 %67, %25
%69 = getelementptr inbounds float, float* %0, i64 %68
%70 = bitcast float* %69 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %70
%71 = getelementptr inbounds float, float* %69, i64 4
%72 = bitcast float* %71 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %72
%73 = getelementptr inbounds float, float* %1, i64 %68
%74 = bitcast float* %73 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %74
%75 = getelementptr inbounds float, float* %73, i64 4
%76 = bitcast float* %75 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %76
%77 = add nuw i64 %36, 32
%78 = add i64 %37, 4
%79 = icmp eq i64 %78, %16
br i1 %79, label %80, label %35
80: ; preds = %35, %34
%81 = phi i64 [ 0, %34 ], [ %77, %35 ]
br i1 %17, label %97, label %82
82: ; preds = %80, %82
%83 = phi i64 [ %94, %82 ], [ %81, %80 ]
%84 = phi i64 [ %95, %82 ], [ 0, %80 ]
%85 = add nuw nsw i64 %83, %25
%86 = getelementptr inbounds float, float* %0, i64 %85
%87 = bitcast float* %86 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %87
%88 = getelementptr inbounds float, float* %86, i64 4
%89 = bitcast float* %88 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %89
%90 = getelementptr inbounds float, float* %1, i64 %85
%91 = bitcast float* %90 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %91
%92 = getelementptr inbounds float, float* %90, i64 4
%93 = bitcast float* %92 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %93
%94 = add nuw i64 %83, 8
%95 = add i64 %84, 1
%96 = icmp eq i64 %95, %14
br i1 %96, label %97, label %82
97: ; preds = %82, %80
br i1 %18, label %133, label %98
98: ; preds = %26, %21, %97
%99 = phi i64 [ 0, %26 ], [ 0, %21 ], [ %13, %97 ]
%100 = xor i64 %99, -1
%101 = add nsw i64 %100, %7
br i1 %20, label %111, label %102
102: ; preds = %98, %102
%103 = phi i64 [ %108, %102 ], [ %99, %98 ]
%104 = phi i64 [ %109, %102 ], [ 0, %98 ]
%105 = add nuw nsw i64 %103, %25
%106 = getelementptr inbounds float, float* %0, i64 %105
store float 0.000000e+00, float* %106
%107 = getelementptr inbounds float, float* %1, i64 %105
store float 0.000000e+00, float* %107
%108 = add nuw nsw i64 %103, 1
%109 = add i64 %104, 1
%110 = icmp eq i64 %109, %19
br i1 %110, label %111, label %102
111: ; preds = %102, %98
%112 = phi i64 [ %99, %98 ], [ %108, %102 ]
%113 = icmp ult i64 %101, 3
br i1 %113, label %133, label %114
114: ; preds = %111, %114
%115 = phi i64 [ %131, %114 ], [ %112, %111 ]
%116 = add nuw nsw i64 %115, %25
%117 = getelementptr inbounds float, float* %0, i64 %116
store float 0.000000e+00, float* %117
%118 = getelementptr inbounds float, float* %1, i64 %116
store float 0.000000e+00, float* %118
%119 = add nuw nsw i64 %115, 1
%120 = add nuw nsw i64 %119, %25
%121 = getelementptr inbounds float, float* %0, i64 %120
store float 0.000000e+00, float* %121
%122 = getelementptr inbounds float, float* %1, i64 %120
store float 0.000000e+00, float* %122
%123 = add nuw nsw i64 %115, 2
%124 = add nuw nsw i64 %123, %25
%125 = getelementptr inbounds float, float* %0, i64 %124
store float 0.000000e+00, float* %125
%126 = getelementptr inbounds float, float* %1, i64 %124
store float 0.000000e+00, float* %126
%127 = add nuw nsw i64 %115, 3
%128 = add nuw nsw i64 %127, %25
%129 = getelementptr inbounds float, float* %0, i64 %128
store float 0.000000e+00, float* %129
%130 = getelementptr inbounds float, float* %1, i64 %128
store float 0.000000e+00, float* %130
%131 = add nuw nsw i64 %115, 4
%132 = icmp eq i64 %131, %7
br i1 %132, label %133, label %114
133: ; preds = %111, %114, %97
%134 = add nuw nsw i64 %22, 1
%135 = icmp eq i64 %134, %7
br i1 %135, label %136, label %21
136: ; preds = %133, %3
ret void
}