double *linspace(double a, double b, int n, double u[]) {
  double c;
  int i;
  if (n < 2 || u == 0)
    return (void *)0;
  c = (b - a) / (n - 1);
  for (i = 0; i < n - 1; ++i)
    u[i] = a + i * c;
  u[n - 1] = b;
  return u;
}


x86 assembly (clang, -O0):
<linspace>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %edi,-0x1c(%rbp)
mov    %rsi,-0x28(%rbp)
cmpl   $0x2,-0x1c(%rbp)
jl     112a <linspace+0x2a>
cmpq   $0x0,-0x28(%rbp)
jne    1137 <linspace+0x37>
movq   $0x0,-0x8(%rbp)
jmp    11bd <linspace+0xbd>
movsd  -0x18(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    119e <linspace+0x9e>
movsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x34(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    115b <linspace+0x5b>
movsd  -0x18(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<linspace>:
xor    %eax,%eax
cmp    $0x2,%edi
jl     1278 <linspace+0x178>
test   %rsi,%rsi
je     1278 <linspace+0x178>
movapd %xmm1,%xmm2
subsd  %xmm0,%xmm2
add    $0xffffffff,%edi
cvtsi2sd %edi,%xmm3
divsd  %xmm3,%xmm2
cmp    $0x4,%edi
jae    1133 <linspace+0x33>
xor    %eax,%eax
jmp    1250 <linspace+0x150>
mov    %edi,%eax
and    $0xfffffffc,%eax
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
movapd %xmm0,%xmm4
unpcklpd %xmm0,%xmm4
lea    -0x4(%rax),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rcx,%rcx
je     1279 <linspace+0x179>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
movdqa 0xe91(%rip),%xmm6
xor    %ecx,%ecx
movdqa 0xe96(%rip),%xmm8
movdqa 0xe9d(%rip),%xmm9
movdqa 0xea4(%rip),%xmm10
movdqa 0xeab(%rip),%xmm11
cs nopw 0x0(%rax,%rax,1)
nop
movdqa %xmm6,%xmm5
cvtdq2pd %xmm6,%xmm7
paddd  %xmm8,%xmm5
cvtdq2pd %xmm5,%xmm5
mulpd  %xmm3,%xmm7
addpd  %xmm4,%xmm7
mulpd  %xmm3,%xmm5
addpd  %xmm4,%xmm5
movupd %xmm7,(%rsi,%rcx,8)
movupd %xmm5,0x10(%rsi,%rcx,8)
movdqa %xmm6,%xmm5
paddd  %xmm9,%xmm5
movdqa %xmm6,%xmm7
cvtdq2pd %xmm5,%xmm5
paddd  %xmm10,%xmm7
cvtdq2pd %xmm7,%xmm7
mulpd  %xmm3,%xmm5
addpd  %xmm4,%xmm5
mulpd  %xmm3,%xmm7
addpd  %xmm4,%xmm7
movupd %xmm5,0x20(%rsi,%rcx,8)
movupd %xmm7,0x30(%rsi,%rcx,8)
add    $0x8,%rcx
paddd  %xmm11,%xmm6
add    $0xfffffffffffffffe,%rdx
jne    11a0 <linspace+0xa0>
test   $0x1,%r8b
je     1242 <linspace+0x142>
cvtdq2pd %xmm6,%xmm5
paddd  0xded(%rip),%xmm6
cvtdq2pd %xmm6,%xmm6
mulpd  %xmm3,%xmm5
addpd  %xmm4,%xmm5
mulpd  %xmm3,%xmm6
addpd  %xmm4,%xmm6
movupd %xmm5,(%rsi,%rcx,8)
movupd %xmm6,0x10(%rsi,%rcx,8)
cmp    %rdi,%rax
je     126d <linspace+0x16d>
nopw   0x0(%rax,%rax,1)
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,(%rsi,%rax,8)
add    $0x1,%rax
cmp    %rax,%rdi
jne    1250 <linspace+0x150>
movslq %edi,%rax
movsd  %xmm1,(%rsi,%rax,8)
mov    %rsi,%rax
ret
movdqa 0xd7f(%rip),%xmm6
xor    %ecx,%ecx
test   $0x1,%r8b
jne    1217 <linspace+0x117>
jmp    1242 <linspace+0x142>

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

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