user: 请从给定的汇编指令推导出C源代码:

```asm
<computeLinGNmatPnt2Pnt>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x140,%rsp
  mov    %edi,-0xa4(%rbp)
  mov    %esi,-0xa8(%rbp)
  mov    %rdx,-0xb0(%rbp)
  mov    %ecx,-0xb4(%rbp)
  mov    %r8d,-0xb8(%rbp)
  mov    %r9,-0xc0(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xc8(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xd0(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0xd8(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0xe0(%rbp)
  mov    0x38(%rbp),%rax
  mov    %rax,-0xe8(%rbp)
  mov    0x40(%rbp),%rax
  mov    %rax,-0xf0(%rbp)
  mov    0x48(%rbp),%rax
  mov    %rax,-0xf8(%rbp)
  mov    0x50(%rbp),%rax
  mov    %rax,-0x100(%rbp)
  mov    0x58(%rbp),%rax
  mov    %rax,-0x108(%rbp)
  mov    0x60(%rbp),%rax
  mov    %rax,-0x110(%rbp)
  mov    0x68(%rbp),%rax
  mov    %rax,-0x118(%rbp)
  mov    0x70(%rbp),%rax
  mov    %rax,-0x120(%rbp)
  mov    0x78(%rbp),%rax
  mov    %rax,-0x128(%rbp)
  mov    0x80(%rbp),%rax
  mov    %rax,-0x130(%rbp)
  mov    0x88(%rbp),%rax
  mov    %rax,-0x138(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0xb4(%rbp),%eax
  shl    $0x2,%eax
  mov    %eax,-0x84(%rbp)
  mov    -0xa4(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x80(%rbp)
  mov    -0xa8(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x7c(%rbp)
  mov    -0x80(%rbp),%eax
  imul   -0x7c(%rbp),%eax
  mov    %eax,-0x78(%rbp)
  mov    -0xb4(%rbp),%eax
  imul   -0xb8(%rbp),%eax
  mov    %eax,-0x74(%rbp)
  movl   $0x0,-0x8c(%rbp)
  jmp    L172
L69:
  mov    -0x8c(%rbp),%eax
  add    %eax,%eax
  mov    %eax,-0x70(%rbp)
  mov    -0x8c(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  mov    %eax,-0x6c(%rbp)
  mov    -0x70(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rsi
  mov    -0xc0(%rbp),%rdx
  mov    -0xa4(%rbp),%ecx
  mov    -0xb4(%rbp),%eax
  movq   %rsi,%xmm0
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <FindSpan@plt>
  mov    %eax,-0x68(%rbp)
  mov    -0x70(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdi
  mov    -0x118(%rbp),%r9
  mov    -0x110(%rbp),%r8
  mov    -0x120(%rbp),%rcx
  mov    -0xc0(%rbp),%rdx
  mov    -0xa4(%rbp),%esi
  mov    -0x68(%rbp),%eax
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <BasisFuns@plt>
  mov    -0x68(%rbp),%eax
  sub    -0xa4(%rbp),%eax
  mov    %eax,-0x64(%rbp)
  mov    -0x70(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rsi
  mov    -0xc8(%rbp),%rdx
  mov    -0xa8(%rbp),%ecx
  mov    -0xb8(%rbp),%eax
  movq   %rsi,%xmm0
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <FindSpan@plt>
  mov    %eax,-0x60(%rbp)
  mov    -0x70(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdi
  mov    -0x130(%rbp),%r9
  mov    -0x128(%rbp),%r8
  mov    -0x138(%rbp),%rcx
  mov    -0xc8(%rbp),%rdx
  mov    -0xa8(%rbp),%esi
  mov    -0x60(%rbp),%eax
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <BasisFuns@plt>
  mov    -0x60(%rbp),%eax
  sub    -0xa8(%rbp),%eax
  mov    %eax,-0x5c(%rbp)
  mov    -0x5c(%rbp),%eax
  imul   -0x84(%rbp),%eax
  mov    %eax,%edx
  mov    -0x64(%rbp),%eax
  shl    $0x2,%eax
  add    %edx,%eax
  mov    %eax,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x0,-0x94(%rbp)
  jmp    L68
L47:
  movl   $0x0,-0x90(%rbp)
  jmp    L137
L175:
  movsd  -0x30(%rbp),%xmm1
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x138(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm2
  mov    -0x94(%rbp),%eax
  imul   -0x84(%rbp),%eax
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  shl    $0x2,%eax
  add    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x28(%rbp),%xmm1
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x138(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm2
  mov    -0x94(%rbp),%eax
  imul   -0x84(%rbp),%eax
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  shl    $0x2,%eax
  add    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x20(%rbp),%xmm1
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x138(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm2
  mov    -0x94(%rbp),%eax
  imul   -0x84(%rbp),%eax
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  shl    $0x2,%eax
  add    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x18(%rbp),%xmm1
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x138(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm2
  mov    -0x94(%rbp),%eax
  imul   -0x84(%rbp),%eax
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  shl    $0x2,%eax
  add    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    %edx,%eax
  cltq
  add    $0x3,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xb0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  addl   $0x1,-0x90(%rbp)
L137:
  mov    -0x90(%rbp),%eax
  cmp    -0xa4(%rbp),%eax
  jle    L175
  addl   $0x1,-0x94(%rbp)
L68:
  mov    -0x94(%rbp),%eax
  cmp    -0xa8(%rbp),%eax
  jle    L47
  movl   $0x0,-0x88(%rbp)
  movl   $0x0,-0x94(%rbp)
  jmp    L138
L98:
  movl   $0x0,-0x90(%rbp)
  jmp    L49
L135:
  mov    -0x5c(%rbp),%edx
  mov    -0x94(%rbp),%eax
  add    %edx,%eax
  imul   -0xb4(%rbp),%eax
  mov    %eax,%edx
  mov    -0x64(%rbp),%eax
  add    %eax,%edx
  mov    -0x90(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x54(%rbp)
  mov    -0x88(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rax,%rdx
  mov    -0x54(%rbp),%eax
  mov    %eax,(%rdx)
  mov    -0x54(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  mov    %eax,-0x54(%rbp)
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x138(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm1
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x88(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm1
  mov    -0x54(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x88(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm1
  mov    -0x54(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x88(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x88(%rbp)
  addl   $0x1,-0x90(%rbp)
L49:
  mov    -0x90(%rbp),%eax
  cmp    -0xa4(%rbp),%eax
  jle    L135
  addl   $0x1,-0x94(%rbp)
L138:
  mov    -0x94(%rbp),%eax
  cmp    -0xa8(%rbp),%eax
  jle    L98
  mov    -0x8c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x6c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x8c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x6c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x8c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x6c(%rbp),%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movl   $0x0,-0x94(%rbp)
  jmp    L32
L141:
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x50(%rbp),%xmm2
  mov    -0x94(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm2
  movsd  -0x48(%rbp),%xmm3
  mov    -0x94(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm3
  mov    -0x94(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x0,-0x90(%rbp)
  jmp    L173
L17:
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x100(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  cmp    %eax,%edx
  jl     L178
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  imul   -0x74(%rbp),%eax
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x100(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x94(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x90(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm2
  mov    -0x94(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm3
  mov    -0x90(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm0,%xmm2
  mov    -0x94(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm3
  mov    -0x90(%rbp),%edx
  mov    %edx,%eax
  add    %eax,%eax
  add    %edx,%eax
  cltq
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  mov    -0x94(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  imul   -0x74(%rbp),%eax
  mov    %eax,%edx
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x100(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L178:
  addl   $0x1,-0x90(%rbp)
L173:
  mov    -0x90(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jl     L17
  addl   $0x1,-0x94(%rbp)
L32:
  mov    -0x94(%rbp),%eax
  cmp    -0x78(%rbp),%eax
  jl     L141
  addl   $0x1,-0x8c(%rbp)
L172:
  mov    -0x8c(%rbp),%eax
  cmp    0x20(%rbp),%eax
  jl     L69
  movl   $0x1,-0x94(%rbp)
  jmp    L192
L107:
  movl   $0x0,-0x90(%rbp)
  jmp    L117
L54:
  mov    -0x90(%rbp),%eax
  imul   -0x74(%rbp),%eax
  mov    %eax,%edx
  mov    -0x94(%rbp),%eax
  add    %edx,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xf0(%rbp),%rax
  add    %rax,%rdx
  mov    -0x94(%rbp),%eax
  imul   -0x74(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x90(%rbp),%eax
  add    %ecx,%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0xf0(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x90(%rbp)
L117:
  mov    -0x90(%rbp),%eax
  cmp    -0x94(%rbp),%eax
  jl     L54
  addl   $0x1,-0x94(%rbp)
L192:
  mov    -0x94(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  jl     L107
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L199
  call   <__stack_chk_fail@plt>
L199:
  leave
  ret
```

assistant: ```c
void computeLinGNmatPnt2Pnt(int degU, int degV, double *cp, int ncp, int kcp, double *knotU, double *knotV, double *us, int nus, double *dirOffset, double *cpDir, double *sgnDistDiff, double *HessApprx, double *gradIent, int *derivInds, double *derivVals, double *leftU, double *rightU, double *NU, double *leftV, double *rightV, double *NV) {
  int i, j, ii, twoii, threeii, spanU, spanV, ind, ind2, dind, fourncp, myint, numCP, orderU, orderV, nnz, cpInd;
  double bspnt[4], doubleVals[3];
  fourncp = 4 * ncp;
  orderU = degU + 1;
  orderV = degV + 1;
  nnz = orderU * orderV;
  numCP = ncp * kcp;
  for (ii = 0; ii < nus; ii++) {
    twoii = 2 * ii;
    threeii = 3 * ii;
    spanU = FindSpan(ncp, degU, us[twoii], knotU);
    BasisFuns(spanU, us[twoii], degU, knotU, NU, leftU, rightU);
    ind = spanU - degU;
    spanV = FindSpan(kcp, degV, us[twoii + 1], knotV);
    BasisFuns(spanV, us[twoii + 1], degV, knotV, NV, leftV, rightV);
    ind2 = spanV - degV;
    myint = ind2 * fourncp + ind * 4;
    bspnt[0] = 0.0;
    bspnt[1] = 0.0;
    bspnt[2] = 0.0;
    bspnt[3] = 0.0;
    for (i = 0; i <= degV; i++) {
      for (j = 0; j <= degU; j++) {
        bspnt[0] += NV[i] * NU[j] * cp[i * fourncp + j * 4 + myint];
        bspnt[1] += NV[i] * NU[j] * cp[i * fourncp + j * 4 + myint + 1];
        bspnt[2] += NV[i] * NU[j] * cp[i * fourncp + j * 4 + myint + 2];
        bspnt[3] += NV[i] * NU[j] * cp[i * fourncp + j * 4 + myint + 3];
      }
    }
    dind = 0;
    for (i = 0; i <= degV; i++) {
      for (j = 0; j <= degU; j++) {
        cpInd = (ind2 + i) * ncp + ind + j;
        derivInds[dind] = cpInd;
        cpInd = 3 * cpInd;
        doubleVals[0] = NV[i] * NU[j] / bspnt[3];
        derivVals[3 * dind] = doubleVals[0] * cpDir[cpInd];
        derivVals[3 * dind + 1] = doubleVals[0] * cpDir[cpInd + 1];
        derivVals[3 * dind + 2] = doubleVals[0] * cpDir[cpInd + 2];
        dind++;
      }
    }
    doubleVals[0] = sgnDistDiff[ii] * dirOffset[threeii];
    doubleVals[1] = sgnDistDiff[ii] * dirOffset[threeii + 1];
    doubleVals[2] = sgnDistDiff[ii] * dirOffset[threeii + 2];
    for (i = 0; i < nnz; i++) {
      gradIent[derivInds[i]] += doubleVals[0] * derivVals[3 * i] + doubleVals[1] * derivVals[3 * i + 1] + doubleVals[2] * derivVals[3 * i + 2];
      for (j = 0; j < nnz; j++) {
        if (derivInds[j] >= derivInds[i]) {
          HessApprx[derivInds[i] * numCP + derivInds[j]] += derivVals[3 * i] * derivVals[3 * j] + derivVals[3 * i + 1] * derivVals[3 * j + 1] + derivVals[3 * i + 2] * derivVals[3 * j + 2];
        }
      }
    }
  }
  for (i = 1; i < numCP; i++) {
    for (j = 0; j < i; j++) {
      HessApprx[i * numCP + j] = HessApprx[j * numCP + i];
    }
  }
  return;
}
```