void cbspl_Bk(double *x, double *G, long *kk, int n, long fix) {
  long i, i4;
  double pp, qq, tmp;
  for (i = 0; i < n; ++i) {
    i4 = 4 * i;
    kk[i] = (long)x[i];
    if (fix && kk[i] > fix)
      kk[i] = fix;
    pp = x[i] - kk[i];
    qq = 1 - pp;
    tmp = pp * pp;
    pp *= tmp;
    G[i4 + 3] = 1. / 6 * pp;
    G[i4 + 1] = 3. / 6 * pp - tmp + 4. / 6;
    tmp = qq * qq;
    qq *= tmp;
    G[i4] = 1. / 6 * qq;
    G[i4 + 2] = 3. / 6 * qq - tmp + 4. / 6;
  }
}


x86 assembly (clang, -O0):
<cbspl_Bk>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cmp    %rcx,%rax
jge    12d5 <cbspl_Bk+0x1d5>
mov    -0x30(%rbp),%rax
shl    $0x2,%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%rdx
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
cmpq   $0x0,-0x28(%rbp)
je     1187 <cbspl_Bk+0x87>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
cmp    -0x28(%rbp),%rax
jle    1187 <cbspl_Bk+0x87>
mov    -0x28(%rbp),%rdx
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
cvtsi2sdq (%rax,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe65(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe2d(%rip),%xmm0
mulsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,0x18(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xde9(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xdd1(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xda9(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xd66(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xd4e(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,0x10(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    111f <cbspl_Bk+0x1f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<cbspl_Bk>:
push   %rbp
push   %r15
push   %r14
push   %rbx
test   %ecx,%ecx
jle    13af <cbspl_Bk+0x2af>
movslq %ecx,%r9
cmp    $0x1,%ecx
jne    111d <cbspl_Bk+0x1d>
xor    %ecx,%ecx
jmp    1300 <cbspl_Bk+0x200>
mov    %r9,%rcx
shl    $0x5,%rcx
add    %rsi,%rcx
cmp    %rdx,%rcx
seta   %r14b
lea    (%rdx,%r9,8),%rbp
cmp    %rsi,%rbp
seta   %r15b
lea    (%rdi,%r9,8),%rbx
cmp    %rdx,%rbx
seta   %al
cmp    %rdi,%rbp
seta   %r11b
cmp    %rsi,%rbx
seta   %r10b
cmp    %rdi,%rcx
seta   %bpl
xor    %ecx,%ecx
test   %r15b,%r14b
jne    1300 <cbspl_Bk+0x200>
and    %r11b,%al
jne    1300 <cbspl_Bk+0x200>
and    %bpl,%r10b
jne    1300 <cbspl_Bk+0x200>
mov    %r9,%rcx
and    $0xfffffffffffffffe,%rcx
movq   %r8,%xmm0
pshufd $0x44,%xmm0,%xmm11
xor    %r10d,%r10d
movdqa 0xe6d(%rip),%xmm8
movdqa %xmm11,%xmm2
pxor   %xmm8,%xmm2
movapd 0xe6a(%rip),%xmm9
movapd 0xe71(%rip),%xmm10
movapd 0xe79(%rip),%xmm5
movapd 0xe81(%rip),%xmm6
mov    %rsi,%rax
jmp    129b <cbspl_Bk+0x19b>
nopw   0x0(%rax,%rax,1)
pxor   %xmm1,%xmm1
movdqa %xmm11,%xmm0
pand   %xmm1,%xmm0
pandn  %xmm3,%xmm1
por    %xmm0,%xmm1
movdqu %xmm1,(%rdx,%r10,8)
movupd (%rdi,%r10,8),%xmm7
movq   %xmm1,%rbp
xorps  %xmm0,%xmm0
cvtsi2sd %rbp,%xmm0
pshufd $0xee,%xmm1,%xmm1
movq   %xmm1,%rbp
xorps  %xmm1,%xmm1
cvtsi2sd %rbp,%xmm1
unpcklpd %xmm1,%xmm0
subpd  %xmm0,%xmm7
movapd %xmm9,%xmm3
subpd  %xmm7,%xmm3
movapd %xmm7,%xmm0
mulpd  %xmm7,%xmm0
mulpd  %xmm0,%xmm7
movapd %xmm7,%xmm1
mulpd  %xmm10,%xmm1
mulpd  %xmm5,%xmm7
subpd  %xmm0,%xmm7
addpd  %xmm6,%xmm7
movapd %xmm3,%xmm0
mulpd  %xmm3,%xmm0
mulpd  %xmm0,%xmm3
movapd %xmm3,%xmm4
mulpd  %xmm10,%xmm4
mulpd  %xmm5,%xmm3
subpd  %xmm0,%xmm3
addpd  %xmm6,%xmm3
movapd %xmm3,%xmm0
unpcklpd %xmm1,%xmm0
unpckhpd %xmm1,%xmm3
movapd %xmm4,%xmm1
unpcklpd %xmm7,%xmm1
unpckhpd %xmm7,%xmm4
movupd %xmm4,0x20(%rax)
movupd %xmm1,(%rax)
movupd %xmm3,0x30(%rax)
movupd %xmm0,0x10(%rax)
add    $0x2,%r10
add    $0x40,%rax
cmp    %r10,%rcx
je     12f7 <cbspl_Bk+0x1f7>
movupd (%rdi,%r10,8),%xmm1
cvttsd2si %xmm1,%rbp
unpckhpd %xmm1,%xmm1
cvttsd2si %xmm1,%rbx
movq   %rbp,%xmm3
movq   %rbx,%xmm1
punpcklqdq %xmm1,%xmm3
test   %r8,%r8
je     11d0 <cbspl_Bk+0xd0>
movdqa %xmm3,%xmm1
pxor   %xmm8,%xmm1
movdqa %xmm1,%xmm7
pcmpgtd %xmm2,%xmm7
pshufd $0xa0,%xmm7,%xmm4
pcmpeqd %xmm2,%xmm1
pshufd $0xf5,%xmm1,%xmm0
pand   %xmm4,%xmm0
pshufd $0xf5,%xmm7,%xmm1
por    %xmm0,%xmm1
jmp    11d4 <cbspl_Bk+0xd4>
cmp    %r9,%rcx
je     13af <cbspl_Bk+0x2af>
mov    %rcx,%rax
shl    $0x5,%rax
add    %rax,%rsi
movsd  0xd3e(%rip),%xmm0
movsd  0xd3e(%rip),%xmm1
movapd 0xd0e(%rip),%xmm2
movapd 0xd16(%rip),%xmm3
nopw   0x0(%rax,%rax,1)
cvttsd2si (%rdi,%rcx,8),%rbp
cmp    %r8,%rbp
mov    %rbp,%rax
cmovg  %r8,%rax
test   %r8,%r8
cmove  %rbp,%rax
mov    %rax,(%rdx,%rcx,8)
movsd  (%rdi,%rcx,8),%xmm4
xorps  %xmm5,%xmm5
cvtsi2sd %rax,%xmm5
subsd  %xmm5,%xmm4
movapd %xmm0,%xmm5
subsd  %xmm4,%xmm5
unpcklpd %xmm5,%xmm4
movapd %xmm4,%xmm5
mulpd  %xmm4,%xmm5
mulpd  %xmm5,%xmm4
movapd %xmm4,%xmm6
mulsd  %xmm1,%xmm6
movsd  %xmm6,0x18(%rsi)
movapd %xmm4,%xmm6
unpckhpd %xmm4,%xmm6
mulsd  %xmm1,%xmm6
movsd  %xmm6,(%rsi)
mulpd  %xmm2,%xmm4
subpd  %xmm5,%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,0x8(%rsi)
add    $0x1,%rcx
add    $0x20,%rsi
cmp    %rcx,%r9
jne    1330 <cbspl_Bk+0x230>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @cbspl_Bk(double* %0, double* %1, i64* %2, i32 %3, i64 %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca i64*
%9 = alloca i32
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
%14 = alloca double
%15 = alloca double
store double* %0, double** %6
store double* %1, double** %7
store i64* %2, i64** %8
store i32 %3, i32* %9
store i64 %4, i64* %10
store i64 0, i64* %11
br label %16
16: ; preds = %100, %5
%17 = load i64, i64* %11
%18 = load i32, i32* %9
%19 = sext i32 %18 to i64
%20 = icmp slt i64 %17, %19
br i1 %20, label %21, label %103
21: ; preds = %16
%22 = load i64, i64* %11
%23 = mul nsw i64 4, %22
store i64 %23, i64* %12
%24 = load double*, double** %6
%25 = load i64, i64* %11
%26 = getelementptr inbounds double, double* %24, i64 %25
%27 = load double, double* %26
%28 = fptosi double %27 to i64
%29 = load i64*, i64** %8
%30 = load i64, i64* %11
%31 = getelementptr inbounds i64, i64* %29, i64 %30
store i64 %28, i64* %31
%32 = load i64, i64* %10
%33 = icmp ne i64 %32, 0
br i1 %33, label %34, label %46
34: ; preds = %21
%35 = load i64*, i64** %8
%36 = load i64, i64* %11
%37 = getelementptr inbounds i64, i64* %35, i64 %36
%38 = load i64, i64* %37
%39 = load i64, i64* %10
%40 = icmp sgt i64 %38, %39
br i1 %40, label %41, label %46
41: ; preds = %34
%42 = load i64, i64* %10
%43 = load i64*, i64** %8
%44 = load i64, i64* %11
%45 = getelementptr inbounds i64, i64* %43, i64 %44
store i64 %42, i64* %45
br label %46
46: ; preds = %41, %34, %21
%47 = load double*, double** %6
%48 = load i64, i64* %11
%49 = getelementptr inbounds double, double* %47, i64 %48
%50 = load double, double* %49
%51 = load i64*, i64** %8
%52 = load i64, i64* %11
%53 = getelementptr inbounds i64, i64* %51, i64 %52
%54 = load i64, i64* %53
%55 = sitofp i64 %54 to double
%56 = fsub double %50, %55
store double %56, double* %13
%57 = load double, double* %13
%58 = fsub double 1.000000e+00, %57
store double %58, double* %14
%59 = load double, double* %13
%60 = load double, double* %13
%61 = fmul double %59, %60
store double %61, double* %15
%62 = load double, double* %15
%63 = load double, double* %13
%64 = fmul double %63, %62
store double %64, double* %13
%65 = load double, double* %13
%66 = fmul double 0x3FC5555555555555, %65
%67 = load double*, double** %7
%68 = load i64, i64* %12
%69 = add nsw i64 %68, 3
%70 = getelementptr inbounds double, double* %67, i64 %69
store double %66, double* %70
%71 = load double, double* %13
%72 = load double, double* %15
%73 = fneg double %72
%74 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %71, double %73)
%75 = fadd double %74, 0x3FE5555555555555
%76 = load double*, double** %7
%77 = load i64, i64* %12
%78 = add nsw i64 %77, 1
%79 = getelementptr inbounds double, double* %76, i64 %78
store double %75, double* %79
%80 = load double, double* %14
%81 = load double, double* %14
%82 = fmul double %80, %81
store double %82, double* %15
%83 = load double, double* %15
%84 = load double, double* %14
%85 = fmul double %84, %83
store double %85, double* %14
%86 = load double, double* %14
%87 = fmul double 0x3FC5555555555555, %86
%88 = load double*, double** %7
%89 = load i64, i64* %12
%90 = getelementptr inbounds double, double* %88, i64 %89
store double %87, double* %90
%91 = load double, double* %14
%92 = load double, double* %15
%93 = fneg double %92
%94 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %91, double %93)
%95 = fadd double %94, 0x3FE5555555555555
%96 = load double*, double** %7
%97 = load i64, i64* %12
%98 = add nsw i64 %97, 2
%99 = getelementptr inbounds double, double* %96, i64 %98
store double %95, double* %99
br label %100
100: ; preds = %46
%101 = load i64, i64* %11
%102 = add nsw i64 %101, 1
store i64 %102, i64* %11
br label %16
103: ; preds = %16
ret void
}

LLVM IR (-O2):
define void @cbspl_Bk(double* nocapture %0, double* nocapture writeonly %1, i64* nocapture %2, i32 %3, i64 %4) {
%6 = sext i32 %3 to i64
%7 = icmp sgt i32 %3, 0
br i1 %7, label %8, label %110
8: ; preds = %5
%9 = icmp ne i64 %4, 0
%10 = icmp eq i32 %3, 1
br i1 %10, label %75, label %11
11: ; preds = %8
%12 = getelementptr i64, i64* %2, i64 %6
%13 = shl nsw i64 %6, 2
%14 = getelementptr double, double* %1, i64 %13
%15 = getelementptr double, double* %0, i64 %6
%16 = bitcast double* %14 to i64*
%17 = icmp ugt i64* %16, %2
%18 = bitcast i64* %12 to double*
%19 = icmp ugt double* %18, %1
%20 = and i1 %17, %19
%21 = bitcast double* %15 to i64*
%22 = icmp ugt i64* %21, %2
%23 = bitcast i64* %12 to double*
%24 = icmp ugt double* %23, %0
%25 = and i1 %22, %24
%26 = or i1 %20, %25
%27 = icmp ugt double* %15, %1
%28 = icmp ugt double* %14, %0
%29 = and i1 %27, %28
%30 = or i1 %26, %29
br i1 %30, label %75, label %31
31: ; preds = %11
%32 = and i64 %6, -2
%33 = insertelement <2 x i64> poison, i64 %4, i64 0
%34 = shufflevector <2 x i64> %33, <2 x i64> poison, <2 x i32> zeroinitializer
%35 = getelementptr inbounds double, double* %1, i64 -2
br label %36
36: ; preds = %36, %31
%37 = phi i64 [ 0, %31 ], [ %71, %36 ]
%38 = shl nsw i64 %37, 2
%39 = getelementptr inbounds double, double* %0, i64 %37
%40 = bitcast double* %39 to <2 x double>*
%41 = load <2 x double>, <2 x double>* %40
%42 = fptosi <2 x double> %41 to <2 x i64>
%43 = getelementptr inbounds i64, i64* %2, i64 %37
%44 = icmp slt <2 x i64> %34, %42
%45 = select i1 %9, <2 x i1> %44, <2 x i1> zeroinitializer
%46 = select <2 x i1> %45, <2 x i64> %34, <2 x i64> %42
%47 = bitcast i64* %43 to <2 x i64>*
store <2 x i64> %46, <2 x i64>* %47
%48 = bitcast double* %39 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = sitofp <2 x i64> %46 to <2 x double>
%51 = fsub <2 x double> %49, %50
%52 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %51
%53 = fmul <2 x double> %51, %51
%54 = fmul <2 x double> %51, %53
%55 = fmul <2 x double> %54, <double 0x3FC5555555555555, double 0x3FC5555555555555>
%56 = fneg <2 x double> %53
%57 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %54, <2 x double> <double 5.000000e-01, double 5.000000e-01>, <2 x double> %56)
%58 = fadd <2 x double> %57, <double 0x3FE5555555555555, double 0x3FE5555555555555>
%59 = fmul <2 x double> %52, %52
%60 = fmul <2 x double> %52, %59
%61 = fmul <2 x double> %60, <double 0x3FC5555555555555, double 0x3FC5555555555555>
%62 = fneg <2 x double> %59
%63 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %60, <2 x double> <double 5.000000e-01, double 5.000000e-01>, <2 x double> %62)
%64 = fadd <2 x double> %63, <double 0x3FE5555555555555, double 0x3FE5555555555555>
%65 = or i64 %38, 2
%66 = getelementptr inbounds double, double* %35, i64 %65
%67 = bitcast double* %66 to <8 x double>*
%68 = shufflevector <2 x double> %61, <2 x double> %58, <4 x i32> <i32 0, i32 1, i32 2, i32 3>
%69 = shufflevector <2 x double> %64, <2 x double> %55, <4 x i32> <i32 0, i32 1, i32 2, i32 3>
%70 = shufflevector <4 x double> %68, <4 x double> %69, <8 x i32> <i32 0, i32 2, i32 4, i32 6, i32 1, i32 3, i32 5, i32 7>
store <8 x double> %70, <8 x double>* %67
%71 = add nuw i64 %37, 2
%72 = icmp eq i64 %71, %32
br i1 %72, label %73, label %36
73: ; preds = %36
%74 = icmp eq i64 %32, %6
br i1 %74, label %110, label %75
75: ; preds = %11, %8, %73
%76 = phi i64 [ 0, %11 ], [ 0, %8 ], [ %32, %73 ]
br label %77
77: ; preds = %75, %77
%78 = phi i64 [ %108, %77 ], [ %76, %75 ]
%79 = shl nsw i64 %78, 2
%80 = getelementptr inbounds double, double* %0, i64 %78
%81 = load double, double* %80
%82 = fptosi double %81 to i64
%83 = getelementptr inbounds i64, i64* %2, i64 %78
%84 = icmp sgt i64 %82, %4
%85 = select i1 %9, i1 %84, i1 false
%86 = select i1 %85, i64 %4, i64 %82
store i64 %86, i64* %83
%87 = load double, double* %80
%88 = sitofp i64 %86 to double
%89 = fsub double %87, %88
%90 = fsub double 1.000000e+00, %89
%91 = insertelement <2 x double> poison, double %89, i64 0
%92 = insertelement <2 x double> %91, double %90, i64 1
%93 = fmul <2 x double> %92, %92
%94 = fmul <2 x double> %92, %93
%95 = extractelement <2 x double> %94, i64 0
%96 = fmul double %95, 0x3FC5555555555555
%97 = or i64 %79, 3
%98 = getelementptr inbounds double, double* %1, i64 %97
store double %96, double* %98
%99 = or i64 %79, 1
%100 = getelementptr inbounds double, double* %1, i64 %99
%101 = extractelement <2 x double> %94, i64 1
%102 = fmul double %101, 0x3FC5555555555555
%103 = getelementptr inbounds double, double* %1, i64 %79
store double %102, double* %103
%104 = fneg <2 x double> %93
%105 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %94, <2 x double> <double 5.000000e-01, double 5.000000e-01>, <2 x double> %104)
%106 = fadd <2 x double> %105, <double 0x3FE5555555555555, double 0x3FE5555555555555>
%107 = bitcast double* %100 to <2 x double>*
store <2 x double> %106, <2 x double>* %107
%108 = add nuw nsw i64 %78, 1
%109 = icmp eq i64 %108, %6
br i1 %109, label %110, label %77
110: ; preds = %77, %73, %5
ret void
}