void creategrid(double *gridx, double *gridy, int xsize, int ysize) {
  int ix, iy;
  int i;
  for (iy = 0, i = 0; iy <= ysize; iy++) {
    for (ix = 0; ix <= xsize; ix++) {
      gridx[i] = ix;
      gridy[i] = iy;
      i++;
    }
  }
}


x86 assembly (clang, -O0):
<creategrid>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     118d <creategrid+0x8d>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     117a <creategrid+0x7a>
cvtsi2sdl -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1133 <creategrid+0x33>
jmp    117f <creategrid+0x7f>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1120 <creategrid+0x20>
pop    %rbp
ret

x86 assembly (clang, -O2):
<creategrid>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %ecx,%ecx
js     134e <creategrid+0x24e>
mov    %edx,%eax
lea    0x1(%rax),%r8
lea    (%rsi,%rax,8),%rbp
add    $0x8,%rbp
mov    %r8,-0x8(%rsp)
and    $0xfffffffffffffffc,%r8
lea    -0x4(%r8),%r10
mov    %r10,-0x10(%rsp)
shr    $0x2,%r10
add    $0x1,%r10
cmp    %rdi,%rbp
seta   %bl
lea    (%rdi,%rax,8),%rax
add    $0x8,%rax
cmp    %rsi,%rax
seta   %r11b
lea    0x1(%rdx),%r13d
and    %bl,%r11b
mov    %r10,%rax
and    $0xfffffffffffffffe,%rax
mov    %rax,-0x18(%rsp)
lea    0x30(%rsi),%rax
mov    %rax,-0x20(%rsp)
lea    0x30(%rdi),%rax
mov    %rax,-0x28(%rsp)
xor    %ebx,%ebx
movdqa 0xe81(%rip),%xmm8
movdqa 0xe89(%rip),%xmm1
movdqa 0xe90(%rip),%xmm9
movdqa 0xe98(%rip),%xmm3
movdqa 0xea0(%rip),%xmm4
xor    %r9d,%r9d
jmp    11c0 <creategrid+0xc0>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0x1(%r9),%eax
cmp    %ecx,%r9d
mov    %eax,%r9d
je     134e <creategrid+0x24e>
test   %edx,%edx
js     11b0 <creategrid+0xb0>
cmp    $0x3,%edx
xorps  %xmm5,%xmm5
cvtsi2sd %r9d,%xmm5
setb   %al
movslq %ebx,%rbx
xor    %ebp,%ebp
or     %r11b,%al
jne    12c8 <creategrid+0x1c8>
movapd %xmm5,%xmm6
unpcklpd %xmm5,%xmm6
cmpq   $0x0,-0x10(%rsp)
je     1337 <creategrid+0x237>
mov    -0x20(%rsp),%rax
lea    (%rax,%rbx,8),%r12
mov    -0x28(%rsp),%rax
lea    (%rax,%rbx,8),%r14
mov    -0x18(%rsp),%rax
xor    %r15d,%r15d
movdqa %xmm8,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movdqa %xmm0,%xmm7
cvtdq2pd %xmm0,%xmm2
paddd  %xmm1,%xmm7
cvtdq2pd %xmm7,%xmm7
movups %xmm2,-0x30(%r14,%r15,8)
movups %xmm7,-0x20(%r14,%r15,8)
movupd %xmm6,-0x30(%r12,%r15,8)
movdqa %xmm0,%xmm2
paddd  %xmm9,%xmm2
movdqa %xmm0,%xmm7
paddd  %xmm3,%xmm7
cvtdq2pd %xmm2,%xmm2
cvtdq2pd %xmm7,%xmm7
movupd %xmm6,-0x20(%r12,%r15,8)
movups %xmm2,-0x10(%r14,%r15,8)
movups %xmm7,(%r14,%r15,8)
movupd %xmm6,-0x10(%r12,%r15,8)
movupd %xmm6,(%r12,%r15,8)
add    $0x8,%r15
paddd  %xmm4,%xmm0
add    $0xfffffffffffffffe,%rax
jne    1220 <creategrid+0x120>
test   $0x1,%r10b
je     12b7 <creategrid+0x1b7>
cvtdq2pd %xmm0,%xmm2
add    %rbx,%r15
paddd  %xmm1,%xmm0
cvtdq2pd %xmm0,%xmm0
movups %xmm2,(%rdi,%r15,8)
movdqu %xmm0,0x10(%rdi,%r15,8)
movupd %xmm6,(%rsi,%r15,8)
movupd %xmm6,0x10(%rsi,%r15,8)
add    %r8,%rbx
mov    %r8d,%ebp
cmp    %r8,-0x8(%rsp)
je     11b0 <creategrid+0xb0>
mov    %r13d,%eax
sub    %ebp,%eax
test   $0x1,%al
jne    12e0 <creategrid+0x1e0>
mov    %ebp,%eax
cmp    %edx,%ebp
jne    1300 <creategrid+0x200>
jmp    11b0 <creategrid+0xb0>
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
movsd  %xmm0,(%rdi,%rbx,8)
movsd  %xmm5,(%rsi,%rbx,8)
add    $0x1,%rbx
lea    0x1(%rbp),%eax
cmp    %edx,%ebp
je     11b0 <creategrid+0xb0>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rdi,%rbx,8)
movsd  %xmm5,(%rsi,%rbx,8)
lea    0x1(%rax),%ebp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
movsd  %xmm0,0x8(%rdi,%rbx,8)
movsd  %xmm5,0x8(%rsi,%rbx,8)
add    $0x2,%rbx
add    $0x2,%eax
cmp    %edx,%ebp
jne    1300 <creategrid+0x200>
jmp    11b0 <creategrid+0xb0>
xor    %r15d,%r15d
movdqa %xmm8,%xmm0
test   $0x1,%r10b
jne    128f <creategrid+0x18f>
jmp    12b7 <creategrid+0x1b7>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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