void grid_cell_center_to_bounds_X2D(const char *xunitstr, long xsize,
                                    long ysize,
                                    const double *restrict grid_center_lon,
                                    double *restrict grid_corner_lon,
                                    double dlon) {
  long i, j, index;
  double minlon, maxlon;
  if (!(dlon > 0))
    dlon = 360. / xsize;
  for (i = 0; i < xsize; ++i) {
    minlon = grid_center_lon[i] - 0.5 * dlon;
    maxlon = grid_center_lon[i] + 0.5 * dlon;
    for (j = 0; j < ysize; ++j) {
      index = (j << 2) * xsize + (i << 2);
      grid_corner_lon[index] = minlon;
      grid_corner_lon[index + 1] = maxlon;
      grid_corner_lon[index + 2] = maxlon;
      grid_corner_lon[index + 3] = minlon;
    }
  }
}


x86 assembly (clang, -O0):
<grid_cell_center_to_bounds_X2D>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
ja     1146 <grid_cell_center_to_bounds_X2D+0x46>
cvtsi2sdq -0x10(%rbp),%xmm1
movsd  0xec3(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    124e <grid_cell_center_to_bounds_X2D+0x14e>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
movsd  0xe9a(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
movsd  0xe6b(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1238 <grid_cell_center_to_bounds_X2D+0x138>
mov    -0x40(%rbp),%rax
shl    $0x2,%rax
imul   -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  %xmm0,0x10(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  %xmm0,0x18(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11b2 <grid_cell_center_to_bounds_X2D+0xb2>
jmp    123d <grid_cell_center_to_bounds_X2D+0x13d>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    114e <grid_cell_center_to_bounds_X2D+0x4e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<grid_cell_center_to_bounds_X2D>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
ja     1127 <grid_cell_center_to_bounds_X2D+0x27>
xorps  %xmm1,%xmm1
cvtsi2sd %rsi,%xmm1
movsd  0xedd(%rip),%xmm0
divsd  %xmm1,%xmm0
test   %rsi,%rsi
jle    123b <grid_cell_center_to_bounds_X2D+0x13b>
unpcklpd %xmm0,%xmm0
lea    -0x1(%rdx),%r9
mov    %edx,%r10d
and    $0x3,%r10d
mov    %rdx,%r11
and    $0xfffffffffffffffc,%r11
mov    %rsi,%r15
shl    $0x7,%r15
mov    %rsi,%r12
shl    $0x5,%r12
xor    %r14d,%r14d
mulpd  0xeb1(%rip),%xmm0
jmp    1181 <grid_cell_center_to_bounds_X2D+0x81>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r14
add    $0x20,%r8
cmp    %rsi,%r14
je     123b <grid_cell_center_to_bounds_X2D+0x13b>
test   %rdx,%rdx
jle    1170 <grid_cell_center_to_bounds_X2D+0x70>
movsd  (%rcx,%r14,8),%xmm1
unpcklpd %xmm1,%xmm1
addpd  %xmm0,%xmm1
movapd %xmm1,%xmm2
shufpd $0x1,%xmm1,%xmm2
cmp    $0x3,%r9
jae    11b0 <grid_cell_center_to_bounds_X2D+0xb0>
xor    %ebx,%ebx
jmp    1204 <grid_cell_center_to_bounds_X2D+0x104>
nopw   0x0(%rax,%rax,1)
mov    %r8,%rdi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movupd %xmm1,(%rdi)
movupd %xmm2,0x10(%rdi)
lea    (%rdi,%r12,1),%rax
movupd %xmm1,(%rdi,%r12,1)
movupd %xmm2,0x10(%rdi,%r12,1)
lea    (%rax,%r12,1),%r13
movupd %xmm1,(%r12,%rax,1)
movupd %xmm2,0x10(%r12,%rax,1)
movupd %xmm1,(%r12,%r13,1)
movupd %xmm2,0x10(%r12,%r13,1)
add    $0x4,%rbx
add    %r15,%rdi
cmp    %rbx,%r11
jne    11c0 <grid_cell_center_to_bounds_X2D+0xc0>
test   %r10,%r10
je     1170 <grid_cell_center_to_bounds_X2D+0x70>
imul   %r12,%rbx
mov    %r10,%rdi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd %xmm1,(%r8,%rbx,1)
movupd %xmm2,0x10(%r8,%rbx,1)
add    %r12,%rbx
add    $0xffffffffffffffff,%rdi
jne    1220 <grid_cell_center_to_bounds_X2D+0x120>
jmp    1170 <grid_cell_center_to_bounds_X2D+0x70>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @grid_cell_center_to_bounds_X2D(i8* %0, i64 %1, i64 %2, double* noalias %3, double* noalias %4, double %5) {
%7 = alloca i8*
%8 = alloca i64
%9 = alloca i64
%10 = alloca double*
%11 = alloca double*
%12 = alloca double
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca double
%17 = alloca double
store i8* %0, i8** %7
store i64 %1, i64* %8
store i64 %2, i64* %9
store double* %3, double** %10
store double* %4, double** %11
store double %5, double* %12
%18 = load double, double* %12
%19 = fcmp ogt double %18, 0.000000e+00
br i1 %19, label %24, label %20
20: ; preds = %6
%21 = load i64, i64* %8
%22 = sitofp i64 %21 to double
%23 = fdiv double 3.600000e+02, %22
store double %23, double* %12
br label %24
24: ; preds = %20, %6
store i64 0, i64* %13
br label %25
25: ; preds = %77, %24
%26 = load i64, i64* %13
%27 = load i64, i64* %8
%28 = icmp slt i64 %26, %27
br i1 %28, label %29, label %80
29: ; preds = %25
%30 = load double*, double** %10
%31 = load i64, i64* %13
%32 = getelementptr inbounds double, double* %30, i64 %31
%33 = load double, double* %32
%34 = load double, double* %12
%35 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %34, double %33)
store double %35, double* %16
%36 = load double*, double** %10
%37 = load i64, i64* %13
%38 = getelementptr inbounds double, double* %36, i64 %37
%39 = load double, double* %38
%40 = load double, double* %12
%41 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %40, double %39)
store double %41, double* %17
store i64 0, i64* %14
br label %42
42: ; preds = %73, %29
%43 = load i64, i64* %14
%44 = load i64, i64* %9
%45 = icmp slt i64 %43, %44
br i1 %45, label %46, label %76
46: ; preds = %42
%47 = load i64, i64* %14
%48 = shl i64 %47, 2
%49 = load i64, i64* %8
%50 = mul nsw i64 %48, %49
%51 = load i64, i64* %13
%52 = shl i64 %51, 2
%53 = add nsw i64 %50, %52
store i64 %53, i64* %15
%54 = load double, double* %16
%55 = load double*, double** %11
%56 = load i64, i64* %15
%57 = getelementptr inbounds double, double* %55, i64 %56
store double %54, double* %57
%58 = load double, double* %17
%59 = load double*, double** %11
%60 = load i64, i64* %15
%61 = add nsw i64 %60, 1
%62 = getelementptr inbounds double, double* %59, i64 %61
store double %58, double* %62
%63 = load double, double* %17
%64 = load double*, double** %11
%65 = load i64, i64* %15
%66 = add nsw i64 %65, 2
%67 = getelementptr inbounds double, double* %64, i64 %66
store double %63, double* %67
%68 = load double, double* %16
%69 = load double*, double** %11
%70 = load i64, i64* %15
%71 = add nsw i64 %70, 3
%72 = getelementptr inbounds double, double* %69, i64 %71
store double %68, double* %72
br label %73
73: ; preds = %46
%74 = load i64, i64* %14
%75 = add nsw i64 %74, 1
store i64 %75, i64* %14
br label %42
76: ; preds = %42
br label %77
77: ; preds = %76
%78 = load i64, i64* %13
%79 = add nsw i64 %78, 1
store i64 %79, i64* %13
br label %25
80: ; preds = %25
ret void
}

LLVM IR (-O2):
define void @grid_cell_center_to_bounds_X2D(i8* nocapture readnone %0, i64 %1, i64 %2, double* noalias nocapture %3, double* noalias nocapture writeonly %4, double %5) {
%7 = fcmp ogt double %5, 0.000000e+00
br i1 %7, label %11, label %8
8: ; preds = %6
%9 = sitofp i64 %1 to double
%10 = fdiv double 3.600000e+02, %9
br label %11
11: ; preds = %8, %6
%12 = phi double [ %5, %6 ], [ %10, %8 ]
%13 = icmp sgt i64 %1, 0
br i1 %13, label %14, label %89
14: ; preds = %11
%15 = icmp sgt i64 %2, 0
%16 = shl i64 %1, 2
%17 = insertelement <2 x double> poison, double %12, i64 0
%18 = shufflevector <2 x double> %17, <2 x double> poison, <2 x i32> zeroinitializer
%19 = add i64 %2, -1
%20 = and i64 %2, 3
%21 = icmp ult i64 %19, 3
%22 = and i64 %2, -4
%23 = icmp eq i64 %20, 0
br label %24
24: ; preds = %14, %86
%25 = phi i64 [ 0, %14 ], [ %87, %86 ]
%26 = getelementptr inbounds double, double* %3, i64 %25
%27 = load double, double* %26
%28 = insertelement <2 x double> poison, double %27, i64 0
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> zeroinitializer
%30 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %18, <2 x double> <double -5.000000e-01, double 5.000000e-01>, <2 x double> %29)
br i1 %15, label %31, label %86
31: ; preds = %24
%32 = shl i64 %25, 2
%33 = shufflevector <2 x double> %30, <2 x double> poison, <2 x i32> <i32 1, i32 0>
br i1 %21, label %71, label %34
34: ; preds = %31, %34
%35 = phi i64 [ %68, %34 ], [ 0, %31 ]
%36 = phi i64 [ %69, %34 ], [ 0, %31 ]
%37 = mul i64 %16, %35
%38 = add nsw i64 %37, %32
%39 = getelementptr inbounds double, double* %4, i64 %38
%40 = bitcast double* %39 to <2 x double>*
store <2 x double> %30, <2 x double>* %40
%41 = or i64 %38, 2
%42 = getelementptr inbounds double, double* %4, i64 %41
%43 = bitcast double* %42 to <2 x double>*
store <2 x double> %33, <2 x double>* %43
%44 = or i64 %35, 1
%45 = mul i64 %16, %44
%46 = add nsw i64 %45, %32
%47 = getelementptr inbounds double, double* %4, i64 %46
%48 = bitcast double* %47 to <2 x double>*
store <2 x double> %30, <2 x double>* %48
%49 = or i64 %46, 2
%50 = getelementptr inbounds double, double* %4, i64 %49
%51 = bitcast double* %50 to <2 x double>*
store <2 x double> %33, <2 x double>* %51
%52 = or i64 %35, 2
%53 = mul i64 %16, %52
%54 = add nsw i64 %53, %32
%55 = getelementptr inbounds double, double* %4, i64 %54
%56 = bitcast double* %55 to <2 x double>*
store <2 x double> %30, <2 x double>* %56
%57 = or i64 %54, 2
%58 = getelementptr inbounds double, double* %4, i64 %57
%59 = bitcast double* %58 to <2 x double>*
store <2 x double> %33, <2 x double>* %59
%60 = or i64 %35, 3
%61 = mul i64 %16, %60
%62 = add nsw i64 %61, %32
%63 = getelementptr inbounds double, double* %4, i64 %62
%64 = bitcast double* %63 to <2 x double>*
store <2 x double> %30, <2 x double>* %64
%65 = or i64 %62, 2
%66 = getelementptr inbounds double, double* %4, i64 %65
%67 = bitcast double* %66 to <2 x double>*
store <2 x double> %33, <2 x double>* %67
%68 = add nuw nsw i64 %35, 4
%69 = add i64 %36, 4
%70 = icmp eq i64 %69, %22
br i1 %70, label %71, label %34
71: ; preds = %34, %31
%72 = phi i64 [ 0, %31 ], [ %68, %34 ]
br i1 %23, label %86, label %73
73: ; preds = %71, %73
%74 = phi i64 [ %83, %73 ], [ %72, %71 ]
%75 = phi i64 [ %84, %73 ], [ 0, %71 ]
%76 = mul i64 %16, %74
%77 = add nsw i64 %76, %32
%78 = getelementptr inbounds double, double* %4, i64 %77
%79 = bitcast double* %78 to <2 x double>*
store <2 x double> %30, <2 x double>* %79
%80 = or i64 %77, 2
%81 = getelementptr inbounds double, double* %4, i64 %80
%82 = bitcast double* %81 to <2 x double>*
store <2 x double> %33, <2 x double>* %82
%83 = add nuw nsw i64 %74, 1
%84 = add i64 %75, 1
%85 = icmp eq i64 %84, %20
br i1 %85, label %86, label %73
86: ; preds = %71, %73, %24
%87 = add nuw nsw i64 %25, 1
%88 = icmp eq i64 %87, %1
br i1 %88, label %89, label %24
89: ; preds = %86, %11
ret void
}