double nurbs(unsigned int n, unsigned int p, const double tau[n + 1],
             const double coord[n - p], const double w[n - p], double x) {
  double coordw[n - p];
  for (unsigned int i = 0; i < n + 0 - p; i++)
    coordw[i] = w[i] * coord[i];
  double sum = bspline_curve(n, p, tau, coordw, x);
  double nrm = bspline_curve(n, p, tau, w, x);
  return sum / nrm;
}


x86 assembly (clang, -O0):
<nurbs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    %ecx,%eax
mov    %rsp,%rcx
mov    %rcx,-0x30(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x58(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x4(%rbp),%ecx
add    $0x0,%ecx
sub    -0x8(%rbp),%ecx
cmp    %ecx,%eax
jae    11a8 <nurbs+0x98>
mov    -0x58(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x3c(%rbp),%edx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x18(%rbp),%rcx
mov    -0x3c(%rbp),%edx
mulsd  (%rcx,%rdx,8),%xmm0
mov    -0x3c(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1162 <nurbs+0x52>
mov    -0x58(%rbp),%rcx
mov    -0x4(%rbp),%edi
mov    -0x8(%rbp),%esi
mov    -0x10(%rbp),%rdx
movsd  -0x28(%rbp),%xmm0
mov    $0x1,%al
call   1030 <bspline_curve@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x8(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x20(%rbp),%rcx
movsd  -0x28(%rbp),%xmm0
mov    $0x1,%al
call   1030 <bspline_curve@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
divsd  -0x50(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nurbs>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %r8,%rbx
mov    %rdx,%r14
mov    %esi,%r15d
mov    %edi,%r12d
mov    %edi,%r9d
sub    %esi,%r9d
lea    0xf(,%r9,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rax
sub    %rdx,%rax
mov    %rax,%rsp
test   %r9d,%r9d
je     1248 <nurbs+0x138>
cmp    $0x4,%r9d
jae    115c <nurbs+0x4c>
xor    %esi,%esi
jmp    1230 <nurbs+0x120>
mov    %r9d,%esi
and    $0xfffffffc,%esi
lea    -0x4(%rsi),%rdx
mov    %rdx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rdx,%rdx
je     12a1 <nurbs+0x191>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd (%rbx,%rdi,8),%xmm3
movupd 0x10(%rbx,%rdi,8),%xmm1
movupd (%rcx,%rdi,8),%xmm2
mulpd  %xmm3,%xmm2
movupd 0x10(%rcx,%rdi,8),%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm2,(%rax,%rdi,8)
movapd %xmm3,0x10(%rax,%rdi,8)
movupd 0x20(%rbx,%rdi,8),%xmm3
movupd 0x30(%rbx,%rdi,8),%xmm1
movupd 0x20(%rcx,%rdi,8),%xmm2
mulpd  %xmm3,%xmm2
movupd 0x30(%rcx,%rdi,8),%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm2,0x20(%rax,%rdi,8)
movapd %xmm3,0x30(%rax,%rdi,8)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rdx
jne    1190 <nurbs+0x80>
test   $0x1,%r8b
je     121e <nurbs+0x10e>
movupd (%rbx,%rdi,8),%xmm3
movupd 0x10(%rbx,%rdi,8),%xmm1
movupd (%rcx,%rdi,8),%xmm2
mulpd  %xmm3,%xmm2
movupd 0x10(%rcx,%rdi,8),%xmm3
mulpd  %xmm1,%xmm3
movapd %xmm2,(%rax,%rdi,8)
movapd %xmm3,0x10(%rax,%rdi,8)
cmp    %r9,%rsi
je     1248 <nurbs+0x138>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rbx,%rsi,8),%xmm1
mulsd  (%rcx,%rsi,8),%xmm1
movsd  %xmm1,(%rax,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%r9
jne    1230 <nurbs+0x120>
mov    %r12d,%edi
mov    %r15d,%esi
mov    %r14,%rdx
mov    %rax,%rcx
mov    $0x1,%al
movsd  %xmm0,-0x28(%rbp)
call   1030 <bspline_curve@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    %r12d,%edi
mov    %r15d,%esi
mov    %r14,%rdx
mov    %rbx,%rcx
movsd  -0x28(%rbp),%xmm0
mov    $0x1,%al
call   1030 <bspline_curve@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  -0x30(%rbp),%xmm0
divsd  %xmm1,%xmm0
lea    -0x20(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edi,%edi
test   $0x1,%r8b
jne    11f5 <nurbs+0xe5>
jmp    121e <nurbs+0x10e>

LLVM IR (-O0):
define double @nurbs(i32 %0, i32 %1, double* %2, double* %3, double* %4, double %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca double
%13 = alloca i8*
%14 = alloca i64
%15 = alloca i32
%16 = alloca double
%17 = alloca double
store i32 %0, i32* %7
store i32 %1, i32* %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
store double %5, double* %12
%18 = load i32, i32* %7
%19 = add i32 %18, 1
%20 = zext i32 %19 to i64
%21 = load i32, i32* %7
%22 = load i32, i32* %8
%23 = sub i32 %21, %22
%24 = zext i32 %23 to i64
%25 = load i32, i32* %7
%26 = load i32, i32* %8
%27 = sub i32 %25, %26
%28 = zext i32 %27 to i64
%29 = load i32, i32* %7
%30 = load i32, i32* %8
%31 = sub i32 %29, %30
%32 = zext i32 %31 to i64
%33 = call i8* @llvm.stacksave()
store i8* %33, i8** %13
%34 = alloca double, i64 %32
store i64 %32, i64* %14
store i32 0, i32* %15
br label %35
35: ; preds = %57, %6
%36 = load i32, i32* %15
%37 = load i32, i32* %7
%38 = add i32 %37, 0
%39 = load i32, i32* %8
%40 = sub i32 %38, %39
%41 = icmp ult i32 %36, %40
br i1 %41, label %42, label %60
42: ; preds = %35
%43 = load double*, double** %11
%44 = load i32, i32* %15
%45 = zext i32 %44 to i64
%46 = getelementptr inbounds double, double* %43, i64 %45
%47 = load double, double* %46
%48 = load double*, double** %10
%49 = load i32, i32* %15
%50 = zext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = fmul double %47, %52
%54 = load i32, i32* %15
%55 = zext i32 %54 to i64
%56 = getelementptr inbounds double, double* %34, i64 %55
store double %53, double* %56
br label %57
57: ; preds = %42
%58 = load i32, i32* %15
%59 = add i32 %58, 1
store i32 %59, i32* %15
br label %35
60: ; preds = %35
%61 = load i32, i32* %7
%62 = load i32, i32* %8
%63 = load double*, double** %9
%64 = load double, double* %12
%65 = call i32 (i32, i32, double*, double*, double, ...) bitcast (i32 (...)* @bspline_curve to i32 (i32, i32, double*, double*, double, ...)*)(i32 %61, i32 %62, double* %63, double* %34, double %64)
%66 = sitofp i32 %65 to double
store double %66, double* %16
%67 = load i32, i32* %7
%68 = load i32, i32* %8
%69 = load double*, double** %9
%70 = load double*, double** %11
%71 = load double, double* %12
%72 = call i32 (i32, i32, double*, double*, double, ...) bitcast (i32 (...)* @bspline_curve to i32 (i32, i32, double*, double*, double, ...)*)(i32 %67, i32 %68, double* %69, double* %70, double %71)
%73 = sitofp i32 %72 to double
store double %73, double* %17
%74 = load double, double* %16
%75 = load double, double* %17
%76 = fdiv double %74, %75
%77 = load i8*, i8** %13
ret double %76
}

LLVM IR (-O2):
define double @nurbs(i32 %0, i32 %1, double* %2, double* nocapture %3, double* %4, double %5) {
%7 = sub i32 %0, %1
%8 = zext i32 %7 to i64
%9 = alloca double, i64 %8
%10 = icmp eq i32 %7, 0
br i1 %10, label %91, label %11
11: ; preds = %6
%12 = icmp ult i32 %7, 4
br i1 %12, label %89, label %13
13: ; preds = %11
%14 = and i64 %8, 4294967292
%15 = add nsw i64 %14, -4
%16 = lshr exact i64 %15, 2
%17 = add nuw nsw i64 %16, 1
%18 = and i64 %17, 1
%19 = icmp eq i64 %15, 0
br i1 %19, label %65, label %20
20: ; preds = %13
%21 = and i64 %17, 9223372036854775806
br label %22
22: ; preds = %22, %20
%23 = phi i64 [ 0, %20 ], [ %62, %22 ]
%24 = phi i64 [ 0, %20 ], [ %63, %22 ]
%25 = getelementptr inbounds double, double* %4, i64 %23
%26 = bitcast double* %25 to <2 x double>*
%27 = load <2 x double>, <2 x double>* %26
%28 = getelementptr inbounds double, double* %25, i64 2
%29 = bitcast double* %28 to <2 x double>*
%30 = load <2 x double>, <2 x double>* %29
%31 = getelementptr inbounds double, double* %3, i64 %23
%32 = bitcast double* %31 to <2 x double>*
%33 = load <2 x double>, <2 x double>* %32
%34 = getelementptr inbounds double, double* %31, i64 2
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = fmul <2 x double> %27, %33
%38 = fmul <2 x double> %30, %36
%39 = getelementptr inbounds double, double* %9, i64 %23
%40 = bitcast double* %39 to <2 x double>*
store <2 x double> %37, <2 x double>* %40
%41 = getelementptr inbounds double, double* %39, i64 2
%42 = bitcast double* %41 to <2 x double>*
store <2 x double> %38, <2 x double>* %42
%43 = or i64 %23, 4
%44 = getelementptr inbounds double, double* %4, i64 %43
%45 = bitcast double* %44 to <2 x double>*
%46 = load <2 x double>, <2 x double>* %45
%47 = getelementptr inbounds double, double* %44, i64 2
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = getelementptr inbounds double, double* %3, i64 %43
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = getelementptr inbounds double, double* %50, i64 2
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = fmul <2 x double> %46, %52
%57 = fmul <2 x double> %49, %55
%58 = getelementptr inbounds double, double* %9, i64 %43
%59 = bitcast double* %58 to <2 x double>*
store <2 x double> %56, <2 x double>* %59
%60 = getelementptr inbounds double, double* %58, i64 2
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %57, <2 x double>* %61
%62 = add nuw i64 %23, 8
%63 = add i64 %24, 2
%64 = icmp eq i64 %63, %21
br i1 %64, label %65, label %22
65: ; preds = %22, %13
%66 = phi i64 [ 0, %13 ], [ %62, %22 ]
%67 = icmp eq i64 %18, 0
br i1 %67, label %87, label %68
68: ; preds = %65
%69 = getelementptr inbounds double, double* %4, i64 %66
%70 = bitcast double* %69 to <2 x double>*
%71 = load <2 x double>, <2 x double>* %70
%72 = getelementptr inbounds double, double* %69, i64 2
%73 = bitcast double* %72 to <2 x double>*
%74 = load <2 x double>, <2 x double>* %73
%75 = getelementptr inbounds double, double* %3, i64 %66
%76 = bitcast double* %75 to <2 x double>*
%77 = load <2 x double>, <2 x double>* %76
%78 = getelementptr inbounds double, double* %75, i64 2
%79 = bitcast double* %78 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = fmul <2 x double> %71, %77
%82 = fmul <2 x double> %74, %80
%83 = getelementptr inbounds double, double* %9, i64 %66
%84 = bitcast double* %83 to <2 x double>*
store <2 x double> %81, <2 x double>* %84
%85 = getelementptr inbounds double, double* %83, i64 2
%86 = bitcast double* %85 to <2 x double>*
store <2 x double> %82, <2 x double>* %86
br label %87
87: ; preds = %65, %68
%88 = icmp eq i64 %14, %8
br i1 %88, label %91, label %89
89: ; preds = %11, %87
%90 = phi i64 [ 0, %11 ], [ %14, %87 ]
br label %97
91: ; preds = %97, %87, %6
%92 = call i32 (i32, i32, double*, double*, double, ...) bitcast (i32 (...)* @bspline_curve to i32 (i32, i32, double*, double*, double, ...)*)(i32 %0, i32 %1, double* %2, double* %9, double %5)
%93 = sitofp i32 %92 to double
%94 = call i32 (i32, i32, double*, double*, double, ...) bitcast (i32 (...)* @bspline_curve to i32 (i32, i32, double*, double*, double, ...)*)(i32 %0, i32 %1, double* %2, double* %4, double %5)
%95 = sitofp i32 %94 to double
%96 = fdiv double %93, %95
ret double %96
97: ; preds = %89, %97
%98 = phi i64 [ %105, %97 ], [ %90, %89 ]
%99 = getelementptr inbounds double, double* %4, i64 %98
%100 = load double, double* %99
%101 = getelementptr inbounds double, double* %3, i64 %98
%102 = load double, double* %101
%103 = fmul double %100, %102
%104 = getelementptr inbounds double, double* %9, i64 %98
store double %103, double* %104
%105 = add nuw nsw i64 %98, 1
%106 = icmp eq i64 %105, %8
br i1 %106, label %91, label %97
}