#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double *malloc(int);
double *hydroloops_linspace(double a, double b, int n) {
  int i;
  double *linspace = malloc(sizeof(double[n]));
  double N = n;
  double interval = (b - a) / (N - 1);
  linspace[0] = a;
  for (i = 1; i < n; i++) {
    linspace[i] = interval + *(linspace + (i - 1));
  }
  return linspace;
}


x86 assembly (clang, -O0):
<hydroloops_linspace>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x3,%rax
mov    %eax,%edi
call   1030 <malloc@plt>
mov    %rax,-0x20(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x10(%rbp),%xmm0
subsd  -0x8(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
movsd  0xea8(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11b7 <hydroloops_linspace+0xa7>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1179 <hydroloops_linspace+0x69>
mov    -0x20(%rbp),%rax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<hydroloops_linspace>:
push   %rbx
sub    $0x10,%rsp
mov    %edi,%ebx
movsd  %xmm1,0x8(%rsp)
movsd  %xmm0,(%rsp)
lea    0x0(,%rbx,8),%edi
call   1030 <malloc@plt>
movsd  (%rsp),%xmm2
movsd  0x8(%rsp),%xmm1
movsd  %xmm2,(%rax)
cmp    $0x2,%ebx
jl     11d2 <hydroloops_linspace+0xc2>
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
addsd  0xeab(%rip),%xmm0
subsd  %xmm2,%xmm1
divsd  %xmm0,%xmm1
mov    %ebx,%edi
lea    -0x1(%rdi),%rsi
add    $0xfffffffffffffffe,%rdi
mov    %esi,%ecx
and    $0x3,%ecx
mov    $0x1,%edx
cmp    $0x3,%rdi
jb     11b5 <hydroloops_linspace+0xa5>
and    $0xfffffffffffffffc,%rsi
xor    %edx,%edx
nopl   (%rax)
addsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%rax,%rdx,8)
addsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rax,%rdx,8)
addsd  %xmm1,%xmm2
movsd  %xmm2,0x18(%rax,%rdx,8)
addsd  %xmm1,%xmm2
movsd  %xmm2,0x20(%rax,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1180 <hydroloops_linspace+0x70>
add    $0x1,%rdx
test   %rcx,%rcx
je     11d2 <hydroloops_linspace+0xc2>
lea    (%rax,%rdx,8),%rdx
xor    %esi,%esi
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rdx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    11c0 <hydroloops_linspace+0xb0>
add    $0x10,%rsp
pop    %rbx
ret

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

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