// File: /workspace/github_cpp_makefile_star2_O2/tomfunkhouser[P]gaps/pkgs/R3Shapes/R3CatmullRomSpline.cpp
// Function: gaps::R3CatmullRomSpline::PointDerivative(double) const
R3Vector R3CatmullRomSpline::
PointDerivative(RNScalar u) const
{
  // Return tangent vector at parametric value u

  // Get index variables
  int i = VertexIndex(u);
  double denom = VertexParameter(i+1) - VertexParameter(i);
  if (RNIsZero(denom)) return R3zero_vector;
  double t = (u - VertexParameter(i)) / denom;

  // Get control vertex_positions
  R3Point P0 = (i > 0) ? VertexPosition(i-1) : VertexPosition(0);
  R3Point P1 = VertexPosition(i);
  R3Point P2 = VertexPosition(i+1);
  R3Point P3 = (i < NVertices()-2) ? VertexPosition(i+2) : VertexPosition(NVertices()-1);

  // Compute derivative vector
  R3Vector v = R3zero_vector;
  v += BlendingWeightDerivative(t, -1) * P0.Vector();
  v += BlendingWeightDerivative(t,  0) * P1.Vector();
  v += BlendingWeightDerivative(t,  1) * P2.Vector();
  v += BlendingWeightDerivative(t,  2) * P3.Vector();

  // Return derivative vector
  return v;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xa8, %rsp
movq %rsi, %r14
movq %rdi, %rbx
movq (%rsi), %rax
movq %rsi, %rdi
movsd %xmm0, (%rsp)
callq *0x268(%rax)
movl %eax, %r15d
leal 0x1(%r15), %ebp
movq 0x10(%r14), %rax
movslq %r15d, %rcx
testq %rax, %rax
je 0x47f5a
movslq %ebp, %rdx
movsd (%rax,%rdx,8), %xmm0
movsd (%rax,%rcx,8), %xmm1
jmp 0x47f63
cvtsi2sd %ebp, %xmm0
cvtsi2sd %r15d, %xmm1
subsd %xmm1, %xmm0
leaq 0x92c3a(%rip), %rdx     # 0xdaba8
movsd (%rdx), %xmm1
movapd 0x4b0b6(%rip), %xmm2    # 0x93030
xorpd %xmm1, %xmm2
ucomisd %xmm2, %xmm0
jb 0x47f9e
ucomisd %xmm0, %xmm1
jb 0x47f9e
leaq 0x516ff(%rip), %rsi     # 0x99690
movq %rbx, %rdi
callq 0x3831c
jmp 0x48200
testq %rax, %rax
je 0x47faa
movsd (%rax,%rcx,8), %xmm1
jmp 0x47faf
cvtsi2sd %r15d, %xmm1
movsd (%rsp), %xmm2
subsd %xmm1, %xmm2
divsd %xmm0, %xmm2
movsd %xmm2, (%rsp)
testl %r15d, %r15d
jle 0x47fd9
leal -0x1(%r15), %edx
leaq 0x48(%rsp), %rdi
movq %r14, %rsi
callq 0x4844c
jmp 0x47fe7
movq 0x8(%r14), %rsi
leaq 0x48(%rsp), %rdi
callq 0x37b5c
leaq 0x90(%rsp), %r13
movq %r13, %rdi
movq %r14, %rsi
movl %r15d, %edx
callq 0x4844c
leaq 0x78(%rsp), %r12
movq %r12, %rdi
movq %r14, %rsi
movl %ebp, %edx
callq 0x4844c
movl 0x20(%r14), %eax
leal -0x2(%rax), %ecx
decl %eax
leal 0x2(%r15), %edx
cmpl %ecx, %r15d
cmovgel %eax, %edx
leaq 0x60(%rsp), %rdi
movq %r14, %rsi
callq 0x4844c
leaq 0x5165a(%rip), %rsi     # 0x99690
movq %rbx, %rdi
callq 0x3831c
movsd (%rsp), %xmm1
movapd %xmm1, %xmm0
mulsd %xmm1, %xmm0
movapd %xmm0, %xmm2
movsd %xmm0, 0x10(%rsp)
movsd 0x60(%r14), %xmm3
movsd 0x4b18d(%rip), %xmm0    # 0x931f0
mulsd %xmm3, %xmm0
mulsd %xmm1, %xmm0
subsd %xmm3, %xmm0
mulsd 0x51f51(%rip), %xmm3    # 0x99fc8
mulsd %xmm2, %xmm3
addsd %xmm0, %xmm3
movsd %xmm3, 0x8(%rsp)
leaq 0x18(%rsp), %rbp
leaq 0x48(%rsp), %rsi
movq %rbp, %rdi
callq 0x212f8
leaq 0x30(%rsp), %r15
movq %r15, %rdi
movq %rbp, %rsi
movsd 0x8(%rsp), %xmm0
callq 0x38a02
movq %rbx, %rdi
movq %r15, %rsi
callq 0x38a24
movsd 0x60(%r14), %xmm0
movsd 0x51f02(%rip), %xmm2    # 0x99fc8
addsd %xmm0, %xmm2
addsd %xmm2, %xmm2
movsd 0x4b132(%rip), %xmm1    # 0x93208
subsd %xmm0, %xmm1
mulsd 0x4b0f6(%rip), %xmm1    # 0x931d8
mulsd 0x10(%rsp), %xmm1
mulsd (%rsp), %xmm2
addsd %xmm1, %xmm2
movsd %xmm2, 0x8(%rsp)
leaq 0x18(%rsp), %r15
movq %r15, %rdi
movq %r13, %rsi
callq 0x212f8
leaq 0x30(%rsp), %r13
movq %r13, %rdi
movq %r15, %rsi
movsd 0x8(%rsp), %xmm0
callq 0x38a02
movq %rbx, %rdi
movq %r13, %rsi
callq 0x38a24
movsd 0x60(%r14), %xmm3
movapd %xmm3, %xmm0
addsd %xmm3, %xmm0
movsd 0x4b09a(%rip), %xmm2    # 0x931d8
movapd %xmm2, %xmm1
subsd %xmm0, %xmm1
addsd %xmm1, %xmm1
mulsd (%rsp), %xmm1
addsd %xmm3, %xmm1
addsd 0x4b0a5(%rip), %xmm3    # 0x93200
mulsd %xmm2, %xmm3
mulsd 0x10(%rsp), %xmm3
addsd %xmm1, %xmm3
movsd %xmm3, 0x8(%rsp)
leaq 0x18(%rsp), %r15
movq %r15, %rdi
movq %r12, %rsi
callq 0x212f8
leaq 0x30(%rsp), %r12
movq %r12, %rdi
movq %r15, %rsi
movsd 0x8(%rsp), %xmm0
callq 0x38a02
movq %rbx, %rdi
movq %r12, %rsi
callq 0x38a24
movsd 0x60(%r14), %xmm0
movsd 0x4b052(%rip), %xmm1    # 0x93200
mulsd %xmm0, %xmm1
mulsd 0x4b01e(%rip), %xmm0    # 0x931d8
mulsd 0x10(%rsp), %xmm0
mulsd (%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
leaq 0x18(%rsp), %r14
movq %r14, %rdi
leaq 0x60(%rsp), %rsi
callq 0x212f8
leaq 0x30(%rsp), %r15
movq %r15, %rdi
movq %r14, %rsi
movsd (%rsp), %xmm0
callq 0x38a02
movq %rbx, %rdi
movq %r15, %rsi
callq 0x38a24
movq %rbx, %rax
addq $0xa8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop