void self_adj_grid_points(double *grid_points, int grid_size, double grid_low,
                          double grid_high) {
  double un;
  int i;
  un = (grid_high - grid_low) / (grid_size - 1.0);
  for (i = 1; i <= grid_size; i++)
    grid_points[i] = grid_low + (i - 1.0) * un;
}


x86 assembly (clang, -O0):
<self_adj_grid_points>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  0xed4(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     118a <self_adj_grid_points+0x8a>
movsd  -0x18(%rbp),%xmm1
cvtsi2sdl -0x2c(%rbp),%xmm0
movsd  0xea2(%rip),%xmm2
subsd  %xmm2,%xmm0
movsd  -0x28(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1140 <self_adj_grid_points+0x40>
pop    %rbp
ret

x86 assembly (clang, -O2):
<self_adj_grid_points>:
test   %esi,%esi
jle    1270 <self_adj_grid_points+0x170>
subsd  %xmm0,%xmm1
cvtsi2sd %esi,%xmm2
addsd  0xee8(%rip),%xmm2
divsd  %xmm2,%xmm1
mov    %esi,%eax
mov    $0x1,%ecx
cmp    $0x4,%esi
jb     124a <self_adj_grid_points+0x14a>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
lea    -0x4(%rcx),%rdx
mov    %rdx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rdx,%rdx
je     1271 <self_adj_grid_points+0x171>
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
movdqa 0xea8(%rip),%xmm6
mov    $0x1,%edx
pcmpeqd %xmm5,%xmm5
movdqa 0xea6(%rip),%xmm8
movdqa 0xead(%rip),%xmm9
movdqa 0xeb4(%rip),%xmm10
nopl   0x0(%rax)
movdqa %xmm6,%xmm4
paddd  %xmm5,%xmm4
movdqa %xmm6,%xmm7
cvtdq2pd %xmm4,%xmm4
psubd  %xmm5,%xmm7
cvtdq2pd %xmm7,%xmm7
mulpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
mulpd  %xmm2,%xmm7
addpd  %xmm3,%xmm7
movupd %xmm4,(%rdi,%rdx,8)
movupd %xmm7,0x10(%rdi,%rdx,8)
movdqa %xmm6,%xmm4
paddd  %xmm8,%xmm4
movdqa %xmm6,%xmm7
cvtdq2pd %xmm4,%xmm4
paddd  %xmm9,%xmm7
cvtdq2pd %xmm7,%xmm7
mulpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
mulpd  %xmm2,%xmm7
addpd  %xmm3,%xmm7
movupd %xmm4,0x20(%rdi,%rdx,8)
movupd %xmm7,0x30(%rdi,%rdx,8)
paddd  %xmm10,%xmm6
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1190 <self_adj_grid_points+0x90>
test   $0x1,%r8b
je     1241 <self_adj_grid_points+0x141>
pcmpeqd %xmm4,%xmm4
movdqa %xmm6,%xmm5
paddd  %xmm4,%xmm5
psubd  %xmm4,%xmm6
cvtdq2pd %xmm5,%xmm4
cvtdq2pd %xmm6,%xmm5
mulpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
mulpd  %xmm2,%xmm5
addpd  %xmm3,%xmm5
movupd %xmm4,(%rdi,%rdx,8)
movupd %xmm5,0x10(%rdi,%rdx,8)
cmp    %rcx,%rax
je     1270 <self_adj_grid_points+0x170>
or     $0x1,%rcx
add    $0x1,%rax
xchg   %ax,%ax
lea    -0x1(%rcx),%edx
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rdi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1250 <self_adj_grid_points+0x150>
ret
movdqa 0xd97(%rip),%xmm6
mov    $0x1,%edx
test   $0x1,%r8b
jne    120e <self_adj_grid_points+0x10e>
jmp    1241 <self_adj_grid_points+0x141>

LLVM IR (-O0):
define void @self_adj_grid_points(double* %0, i32 %1, double %2, double %3) {
%5 = alloca double*
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32
store double* %0, double** %5
store i32 %1, i32* %6
store double %2, double* %7
store double %3, double* %8
%11 = load double, double* %8
%12 = load double, double* %7
%13 = fsub double %11, %12
%14 = load i32, i32* %6
%15 = sitofp i32 %14 to double
%16 = fsub double %15, 1.000000e+00
%17 = fdiv double %13, %16
store double %17, double* %9
store i32 1, i32* %10
br label %18
18: ; preds = %33, %4
%19 = load i32, i32* %10
%20 = load i32, i32* %6
%21 = icmp sle i32 %19, %20
br i1 %21, label %22, label %36
22: ; preds = %18
%23 = load double, double* %7
%24 = load i32, i32* %10
%25 = sitofp i32 %24 to double
%26 = fsub double %25, 1.000000e+00
%27 = load double, double* %9
%28 = call double @llvm.fmuladd.f64(double %26, double %27, double %23)
%29 = load double*, double** %5
%30 = load i32, i32* %10
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %29, i64 %31
store double %28, double* %32
br label %33
33: ; preds = %22
%34 = load i32, i32* %10
%35 = add nsw i32 %34, 1
store i32 %35, i32* %10
br label %18
36: ; preds = %18
ret void
}

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