void linspace(double array[], double xi, double xf, int n) {
  int i;
  double dx, range = xf - xi;
  dx = range / (n - 1);
  array[0] = xi;
  for (i = 1; i < n; i++) {
    array[i] = array[i - 1] + dx;
  }
}


x86 assembly (clang, -O0):
<linspace>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %esi,-0x1c(%rbp)
movsd  -0x18(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x10(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    118e <linspace+0x8e>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1150 <linspace+0x50>
pop    %rbp
ret

x86 assembly (clang, -O2):
<linspace>:
movsd  %xmm0,(%rdi)
cmp    $0x2,%esi
jl     1192 <linspace+0x92>
lea    -0x1(%rsi),%eax
cvtsi2sd %eax,%xmm2
subsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
mov    %esi,%esi
lea    -0x1(%rsi),%rdx
add    $0xfffffffffffffffe,%rsi
mov    %edx,%eax
and    $0x3,%eax
mov    $0x1,%ecx
cmp    $0x3,%rsi
jb     1175 <linspace+0x75>
and    $0xfffffffffffffffc,%rdx
xor    %ecx,%ecx
nopl   0x0(%rax)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rdi,%rcx,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x10(%rdi,%rcx,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rdi,%rcx,8)
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rdi,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1140 <linspace+0x40>
add    $0x1,%rcx
test   %rax,%rax
je     1192 <linspace+0x92>
lea    (%rdi,%rcx,8),%rcx
xor    %edx,%edx
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1180 <linspace+0x80>
ret

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

LLVM IR (-O2):
define void @linspace(double* nocapture %0, double %1, double %2, i32 %3) {
%5 = fsub double %2, %1
%6 = add nsw i32 %3, -1
%7 = sitofp i32 %6 to double
%8 = fdiv double %5, %7
store double %1, double* %0
%9 = icmp sgt i32 %3, 1
br i1 %9, label %10, label %49
10: ; preds = %4
%11 = zext i32 %3 to i64
%12 = add nsw i64 %11, -1
%13 = add nsw i64 %11, -2
%14 = and i64 %12, 3
%15 = icmp ult i64 %13, 3
br i1 %15, label %36, label %16
16: ; preds = %10
%17 = and i64 %12, -4
br label %18
18: ; preds = %18, %16
%19 = phi double [ %1, %16 ], [ %31, %18 ]
%20 = phi i64 [ 1, %16 ], [ %33, %18 ]
%21 = phi i64 [ 0, %16 ], [ %34, %18 ]
%22 = fadd double %8, %19
%23 = getelementptr inbounds double, double* %0, i64 %20
store double %22, double* %23
%24 = add nuw nsw i64 %20, 1
%25 = fadd double %8, %22
%26 = getelementptr inbounds double, double* %0, i64 %24
store double %25, double* %26
%27 = add nuw nsw i64 %20, 2
%28 = fadd double %8, %25
%29 = getelementptr inbounds double, double* %0, i64 %27
store double %28, double* %29
%30 = add nuw nsw i64 %20, 3
%31 = fadd double %8, %28
%32 = getelementptr inbounds double, double* %0, i64 %30
store double %31, double* %32
%33 = add nuw nsw i64 %20, 4
%34 = add i64 %21, 4
%35 = icmp eq i64 %34, %17
br i1 %35, label %36, label %18
36: ; preds = %18, %10
%37 = phi double [ %1, %10 ], [ %31, %18 ]
%38 = phi i64 [ 1, %10 ], [ %33, %18 ]
%39 = icmp eq i64 %14, 0
br i1 %39, label %49, label %40
40: ; preds = %36, %40
%41 = phi double [ %44, %40 ], [ %37, %36 ]
%42 = phi i64 [ %46, %40 ], [ %38, %36 ]
%43 = phi i64 [ %47, %40 ], [ 0, %36 ]
%44 = fadd double %8, %41
%45 = getelementptr inbounds double, double* %0, i64 %42
store double %44, double* %45
%46 = add nuw nsw i64 %42, 1
%47 = add i64 %43, 1
%48 = icmp eq i64 %47, %14
br i1 %48, label %49, label %40
49: ; preds = %36, %40, %4
ret void
}