void calculate_linear_interpolation(int length_of_interpolation,
                                    double **interpolation,
                                    const double min_value,
                                    const double max_value) {
  double interval = (max_value - min_value) / ((double)length_of_interpolation);
  int coordinate = 0;
  for (coordinate = 0; coordinate < length_of_interpolation; coordinate++) {
    *interpolation[coordinate] = min_value + interval * ((double)coordinate);
  }
}


x86 assembly (clang, -O0):
<calculate_linear_interpolation>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0x4(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    117c <calculate_linear_interpolation+0x7c>
movsd  -0x18(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    113b <calculate_linear_interpolation+0x3b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculate_linear_interpolation>:
test   %edi,%edi
jle    11d1 <calculate_linear_interpolation+0xd1>
subsd  %xmm0,%xmm1
cvtsi2sd %edi,%xmm2
divsd  %xmm2,%xmm1
mov    %edi,%edx
lea    -0x1(%rdx),%rcx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1129 <calculate_linear_interpolation+0x29>
xor    %ecx,%ecx
jmp    11a1 <calculate_linear_interpolation+0xa1>
and    $0xfffffffc,%edx
xor    %ecx,%ecx
xchg   %ax,%ax
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mov    (%rsi,%rcx,8),%rdi
movsd  %xmm2,(%rdi)
lea    0x1(%rcx),%edi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mov    0x8(%rsi,%rcx,8),%rdi
movsd  %xmm2,(%rdi)
lea    0x2(%rcx),%edi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mov    0x10(%rsi,%rcx,8),%rdi
movsd  %xmm2,(%rdi)
lea    0x3(%rcx),%edi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mov    0x18(%rsi,%rcx,8),%rdi
movsd  %xmm2,(%rdi)
add    $0x4,%rcx
cmp    %rdx,%rcx
jne    1130 <calculate_linear_interpolation+0x30>
test   %rax,%rax
je     11d1 <calculate_linear_interpolation+0xd1>
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mov    (%rsi,%rcx,8),%rdx
movsd  %xmm2,(%rdx)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rax
jne    11b0 <calculate_linear_interpolation+0xb0>
ret

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

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