user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<s_update_hessian_box_mpc>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x8,%rsp
  mov    %edi,0x4(%rsp)
  mov    %esi,-0x60(%rsp)
  mov    %edx,(%rsp)
  mov    %ecx,-0x64(%rsp)
  mov    %r8d,-0x4(%rsp)
  movaps %xmm0,%xmm1
  mov    %r9,-0x18(%rsp)
  mov    (%r9),%rdi
  mov    0x48(%rsp),%rax
  mov    (%rax),%r15
  mov    0x50(%rsp),%rax
  mov    (%rax),%r12
  mov    0x58(%rsp),%rax
  mov    (%rax),%r14
  mov    0x40(%rsp),%rax
  mov    (%rax),%r13
  mov    0x70(%rsp),%rax
  mov    (%rax),%rbx
  mov    %rbx,-0x40(%rsp)
  mov    0x78(%rsp),%rax
  mov    (%rax),%r10
  mov    %r10,-0x38(%rsp)
  mov    0x80(%rsp),%rax
  mov    (%rax),%r9
  mov    %r9,-0x30(%rsp)
  mov    0x60(%rsp),%rax
  mov    (%rax),%rbp
  mov    %rbp,-0x28(%rsp)
  mov    0x68(%rsp),%rax
  mov    (%rax),%r11
  mov    %r11,-0x20(%rsp)
  mov    0x88(%rsp),%rax
  mov    (%rax),%rcx
  mov    %rcx,-0x10(%rsp)
  cmp    $0x3,%esi
  jle    L167
  lea    0x10(,%r8,4),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x78(%rsp)
  mov    %rcx,%r8
  lea    -0x4(%rsi),%eax
  shr    $0x2,%eax
  mov    %eax,-0x58(%rsp)
  add    $0x1,%eax
  mov    %eax,%eax
  shl    $0x5,%rax
  mov    %rax,-0x50(%rsp)
  add    %rdi,%rax
  mov    %rax,-0x70(%rsp)
  mov    %r13,%rdx
  mov    %r14,%rsi
  mov    %r12,%rax
  mov    %r15,%rcx
  movss  D20(%rip),%xmm2
  mov    %r12,-0x48(%rsp)
  mov    -0x70(%rsp),%r12
L183:
  movaps %xmm2,%xmm0
  divss  (%rdi),%xmm0
  movss  %xmm0,(%rdx)
  movaps %xmm2,%xmm3
  divss  0x4(%rdi),%xmm3
  movss  %xmm3,0x4(%rdx)
  mulss  (%rcx),%xmm0
  movss  %xmm0,(%rax)
  movss  0x4(%rcx),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  movaps %xmm1,%xmm0
  mulss  (%rdx),%xmm0
  movss  %xmm0,(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rsi)
  movss  0x0(%rbp),%xmm0
  addss  (%rax),%xmm0
  addss  0x4(%rax),%xmm0
  movss  %xmm0,(%rbx)
  movss  (%r11),%xmm0
  addss  0x4(%rcx),%xmm0
  movss  0x4(%rax),%xmm3
  mulss  0x4(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0x4(%rsi),%xmm0
  subss  (%rcx),%xmm0
  movss  (%rax),%xmm3
  mulss  (%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  (%rsi),%xmm0
  movss  %xmm0,(%r10)
  movss  %xmm0,(%r9)
  movaps %xmm2,%xmm0
  divss  0x8(%rdi),%xmm0
  movss  %xmm0,0x8(%rdx)
  movaps %xmm2,%xmm3
  divss  0xc(%rdi),%xmm3
  movss  %xmm3,0xc(%rdx)
  mulss  0x8(%rcx),%xmm0
  movss  %xmm0,0x8(%rax)
  movss  0xc(%rcx),%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  %xmm0,0xc(%rax)
  movaps %xmm1,%xmm0
  mulss  0x8(%rdx),%xmm0
  movss  %xmm0,0x8(%rsi)
  movaps %xmm1,%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  %xmm0,0xc(%rsi)
  movss  0x4(%rbp),%xmm0
  addss  0x8(%rax),%xmm0
  addss  0xc(%rax),%xmm0
  movss  %xmm0,0x14(%rbx)
  movss  0x4(%r11),%xmm0
  addss  0xc(%rcx),%xmm0
  movss  0xc(%rax),%xmm3
  mulss  0xc(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0xc(%rsi),%xmm0
  subss  0x8(%rcx),%xmm0
  movss  0x8(%rax),%xmm3
  mulss  0x8(%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  0x8(%rsi),%xmm0
  movss  %xmm0,0x10(%r10)
  movss  %xmm0,0x4(%r9)
  movaps %xmm2,%xmm0
  divss  0x10(%rdi),%xmm0
  movss  %xmm0,0x10(%rdx)
  movaps %xmm2,%xmm3
  divss  0x14(%rdi),%xmm3
  movss  %xmm3,0x14(%rdx)
  mulss  0x10(%rcx),%xmm0
  movss  %xmm0,0x10(%rax)
  movss  0x14(%rcx),%xmm0
  mulss  0x14(%rdx),%xmm0
  movss  %xmm0,0x14(%rax)
  movaps %xmm1,%xmm0
  mulss  0x10(%rdx),%xmm0
  movss  %xmm0,0x10(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x14(%rdx),%xmm0
  movss  %xmm0,0x14(%rsi)
  movss  0x8(%rbp),%xmm0
  addss  0x10(%rax),%xmm0
  addss  0x14(%rax),%xmm0
  movss  %xmm0,0x28(%rbx)
  movss  0x8(%r11),%xmm0
  addss  0x14(%rcx),%xmm0
  movss  0x14(%rax),%xmm3
  mulss  0x14(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0x14(%rsi),%xmm0
  subss  0x10(%rcx),%xmm0
  movss  0x10(%rax),%xmm3
  mulss  0x10(%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  0x10(%rsi),%xmm0
  movss  %xmm0,0x20(%r10)
  movss  %xmm0,0x8(%r9)
  movaps %xmm2,%xmm0
  divss  0x18(%rdi),%xmm0
  movss  %xmm0,0x18(%rdx)
  movaps %xmm2,%xmm3
  divss  0x1c(%rdi),%xmm3
  movss  %xmm3,0x1c(%rdx)
  mulss  0x18(%rcx),%xmm0
  movss  %xmm0,0x18(%rax)
  movss  0x1c(%rcx),%xmm0
  mulss  0x1c(%rdx),%xmm0
  movss  %xmm0,0x1c(%rax)
  movaps %xmm1,%xmm0
  mulss  0x18(%rdx),%xmm0
  movss  %xmm0,0x18(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x1c(%rdx),%xmm0
  movss  %xmm0,0x1c(%rsi)
  movss  0xc(%rbp),%xmm0
  addss  0x18(%rax),%xmm0
  addss  0x1c(%rax),%xmm0
  movss  %xmm0,0x3c(%rbx)
  movss  0xc(%r11),%xmm0
  addss  0x1c(%rcx),%xmm0
  movss  0x1c(%rax),%xmm3
  mulss  0x1c(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0x1c(%rsi),%xmm0
  subss  0x18(%rcx),%xmm0
  movss  0x18(%rax),%xmm3
  mulss  0x18(%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  0x18(%rsi),%xmm0
  movss  %xmm0,0x30(%r10)
  movss  %xmm0,0xc(%r9)
  add    $0x20,%rdi
  add    $0x20,%rcx
  add    $0x20,%rax
  add    $0x20,%rsi
  add    $0x20,%rdx
  add    $0x10,%rbp
  add    -0x78(%rsp),%rbx
  add    $0x10,%r11
  add    $0x20,%r8
  add    $0x40,%r10
  add    $0x10,%r9
  cmp    %r12,%rdi
  jne    L183
  mov    -0x48(%rsp),%r12
  mov    -0x50(%rsp),%rax
  add    %rax,%r15
  add    %rax,%r12
  add    %rax,%r14
  add    %rax,%r13
  mov    -0x58(%rsp),%eax
  lea    0x4(,%rax,4),%ecx
L56:
  mov    -0x60(%rsp),%eax
  cmp    %ecx,%eax
  jle    L119
  mov    %eax,%esi
  sub    %ecx,%esi
  test   %esi,%esi
  jle    L119
  lea    0x0(,%rcx,4),%eax
  cltq
  mov    -0x4(%rsp),%edx
  imul   %ecx,%edx
  movslq %edx,%rdx
  add    %rax,%rdx
  mov    -0x40(%rsp),%rbx
  lea    (%rbx,%rdx,4),%rdx
  movslq %ecx,%rcx
  lea    -0x1(%rsi),%esi
  lea    0x1(%rcx,%rsi,1),%rsi
  shl    $0x2,%rsi
  neg    %rcx
  mov    -0x70(%rsp),%rbx
  lea    (%rbx,%rcx,8),%rcx
  movss  D20(%rip),%xmm3
  mov    -0x38(%rsp),%r8
  mov    -0x30(%rsp),%r9
  mov    -0x28(%rsp),%r10
  mov    -0x20(%rsp),%r11
  mov    -0x10(%rsp),%rdi
L164:
  movaps %xmm3,%xmm0
  divss  (%rcx,%rax,2),%xmm0
  movss  %xmm0,0x0(%r13)
  movaps %xmm3,%xmm2
  divss  0x4(%rcx,%rax,2),%xmm2
  movss  %xmm2,0x4(%r13)
  mulss  (%r15),%xmm0
  movss  %xmm0,(%r12)
  movss  0x4(%r15),%xmm0
  mulss  0x4(%r13),%xmm0
  movss  %xmm0,0x4(%r12)
  movaps %xmm1,%xmm0
  mulss  0x0(%r13),%xmm0
  movss  %xmm0,(%r14)
  movaps %xmm1,%xmm0
  mulss  0x4(%r13),%xmm0
  movss  %xmm0,0x4(%r14)
  movss  (%r10,%rax,1),%xmm0
  addss  (%r12),%xmm0
  addss  0x4(%r12),%xmm0
  movss  %xmm0,(%rdx)
  movss  (%r11,%rax,1),%xmm0
  addss  0x4(%r15),%xmm0
  movss  0x4(%r12),%xmm2
  mulss  0x4(%rdi,%rax,2),%xmm2
  addss  %xmm2,%xmm0
  addss  0x4(%r14),%xmm0
  subss  (%r15),%xmm0
  movss  (%r12),%xmm2
  mulss  (%rdi,%rax,2),%xmm2
  subss  %xmm2,%xmm0
  subss  (%r14),%xmm0
  movss  %xmm0,(%r8,%rax,4)
  movss  %xmm0,(%r9,%rax,1)
  add    $0x8,%r15
  add    $0x8,%r12
  add    $0x8,%r14
  add    $0x8,%r13
  add    $0x14,%rdx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L164
L119:
  mov    0x4(%rsp),%eax
  cmp    $0x1,%eax
  jle    L23
  mov    %eax,%eax
  mov    %rax,-0x10(%rsp)
  mov    -0x4(%rsp),%eax
  lea    0x10(,%rax,4),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x78(%rsp)
  mov    -0x64(%rsp),%eax
  sub    $0x4,%eax
  shr    $0x2,%eax
  add    $0x1,%eax
  mov    %eax,%edx
  shl    $0x5,%rdx
  mov    %rdx,-0x58(%rsp)
  shl    $0x2,%eax
  mov    %eax,-0x8(%rsp)
  movq   $0x1,-0x70(%rsp)
  movss  D20(%rip),%xmm2
  movaps %xmm2,%xmm3
  jmp    L121
L167:
  mov    %rdi,-0x70(%rsp)
  mov    $0x0,%ecx
  jmp    L56
L35:
  mov    %rdi,-0x60(%rsp)
  mov    $0x0,%edx
  jmp    L9
L140:
  movaps %xmm3,%xmm0
  divss  (%rdx,%rax,2),%xmm0
  movss  %xmm0,0x0(%r13)
  movaps %xmm3,%xmm4
  divss  0x4(%rdx,%rax,2),%xmm4
  movss  %xmm4,0x4(%r13)
  mulss  (%r15),%xmm0
  movss  %xmm0,(%r12)
  movss  0x4(%r15),%xmm0
  mulss  0x4(%r13),%xmm0
  movss  %xmm0,0x4(%r12)
  movaps %xmm1,%xmm0
  mulss  0x0(%r13),%xmm0
  movss  %xmm0,(%r14)
  movaps %xmm1,%xmm0
  mulss  0x4(%r13),%xmm0
  movss  %xmm0,0x4(%r14)
  movss  (%r10,%rax,1),%xmm0
  addss  (%r12),%xmm0
  addss  0x4(%r12),%xmm0
  movss  %xmm0,(%rcx)
  movss  (%r11,%rax,1),%xmm0
  addss  0x4(%r15),%xmm0
  movss  0x4(%r12),%xmm4
  mulss  0x4(%rdi,%rax,2),%xmm4
  addss  %xmm4,%xmm0
  addss  0x4(%r14),%xmm0
  subss  (%r15),%xmm0
  movss  (%r12),%xmm4
  mulss  (%rdi,%rax,2),%xmm4
  subss  %xmm4,%xmm0
  subss  (%r14),%xmm0
  movss  %xmm0,(%r8,%rax,4)
  movss  %xmm0,(%r9,%rax,1)
  add    $0x8,%r15
  add    $0x8,%r12
  add    $0x8,%r14
  add    $0x8,%r13
  add    $0x14,%rcx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L140
L196:
  addq   $0x1,-0x70(%rsp)
  mov    -0x70(%rsp),%rax
  cmp    -0x10(%rsp),%rax
  je     L23
L121:
  mov    -0x18(%rsp),%rbx
  mov    -0x70(%rsp),%rax
  mov    (%rbx,%rax,8),%rdi
  mov    0x48(%rsp),%rbx
  mov    (%rbx,%rax,8),%r15
  mov    0x50(%rsp),%rbx
  mov    (%rbx,%rax,8),%r12
  mov    0x58(%rsp),%rbx
  mov    (%rbx,%rax,8),%r14
  mov    0x40(%rsp),%rbx
  mov    (%rbx,%rax,8),%r13
  mov    0x70(%rsp),%rbx
  mov    (%rbx,%rax,8),%rbx
  mov    %rbx,-0x48(%rsp)
  mov    0x78(%rsp),%rsi
  mov    (%rsi,%rax,8),%r10
  mov    %r10,-0x40(%rsp)
  mov    0x80(%rsp),%rsi
  mov    (%rsi,%rax,8),%r9
  mov    %r9,-0x38(%rsp)
  mov    0x60(%rsp),%rcx
  mov    (%rcx,%rax,8),%rbp
  mov    %rbp,-0x30(%rsp)
  mov    0x68(%rsp),%rsi
  mov    (%rsi,%rax,8),%r11
  mov    %r11,-0x28(%rsp)
  mov    0x88(%rsp),%rcx
  mov    (%rcx,%rax,8),%r8
  mov    %r8,-0x20(%rsp)
  cmpl   $0x3,-0x64(%rsp)
  jle    L35
  mov    -0x58(%rsp),%rax
  add    %rdi,%rax
  mov    %rax,-0x60(%rsp)
  mov    %r13,%rdx
  mov    %r14,%rsi
  mov    %r12,%rax
  mov    %r15,%rcx
  mov    %r12,-0x50(%rsp)
  mov    -0x60(%rsp),%r12
L65:
  movaps %xmm2,%xmm0
  divss  (%rdi),%xmm0
  movss  %xmm0,(%rdx)
  movaps %xmm2,%xmm4
  divss  0x4(%rdi),%xmm4
  movss  %xmm4,0x4(%rdx)
  mulss  (%rcx),%xmm0
  movss  %xmm0,(%rax)
  movss  0x4(%rcx),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  movaps %xmm1,%xmm0
  mulss  (%rdx),%xmm0
  movss  %xmm0,(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rsi)
  movss  0x0(%rbp),%xmm0
  addss  (%rax),%xmm0
  addss  0x4(%rax),%xmm0
  movss  %xmm0,(%rbx)
  movss  (%r11),%xmm0
  addss  0x4(%rcx),%xmm0
  movss  0x4(%rax),%xmm4
  mulss  0x4(%r8),%xmm4
  addss  %xmm4,%xmm0
  addss  0x4(%rsi),%xmm0
  subss  (%rcx),%xmm0
  movss  (%rax),%xmm4
  mulss  (%r8),%xmm4
  subss  %xmm4,%xmm0
  subss  (%rsi),%xmm0
  movss  %xmm0,(%r10)
  movss  %xmm0,(%r9)
  movaps %xmm2,%xmm0
  divss  0x8(%rdi),%xmm0
  movss  %xmm0,0x8(%rdx)
  movaps %xmm2,%xmm4
  divss  0xc(%rdi),%xmm4
  movss  %xmm4,0xc(%rdx)
  mulss  0x8(%rcx),%xmm0
  movss  %xmm0,0x8(%rax)
  movss  0xc(%rcx),%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  %xmm0,0xc(%rax)
  movaps %xmm1,%xmm0
  mulss  0x8(%rdx),%xmm0
  movss  %xmm0,0x8(%rsi)
  movaps %xmm1,%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  %xmm0,0xc(%rsi)
  movss  0x4(%rbp),%xmm0
  addss  0x8(%rax),%xmm0
  addss  0xc(%rax),%xmm0
  movss  %xmm0,0x14(%rbx)
  movss  0x4(%r11),%xmm0
  addss  0xc(%rcx),%xmm0
  movss  0xc(%rax),%xmm4
  mulss  0xc(%r8),%xmm4
  addss  %xmm4,%xmm0
  addss  0xc(%rsi),%xmm0
  subss  0x8(%rcx),%xmm0
  movss  0x8(%rax),%xmm4
  mulss  0x8(%r8),%xmm4
  subss  %xmm4,%xmm0
  subss  0x8(%rsi),%xmm0
  movss  %xmm0,0x10(%r10)
  movss  %xmm0,0x4(%r9)
  movaps %xmm2,%xmm0
  divss  0x10(%rdi),%xmm0
  movss  %xmm0,0x10(%rdx)
  movaps %xmm2,%xmm4
  divss  0x14(%rdi),%xmm4
  movss  %xmm4,0x14(%rdx)
  mulss  0x10(%rcx),%xmm0
  movss  %xmm0,0x10(%rax)
  movss  0x14(%rcx),%xmm0
  mulss  0x14(%rdx),%xmm0
  movss  %xmm0,0x14(%rax)
  movaps %xmm1,%xmm0
  mulss  0x10(%rdx),%xmm0
  movss  %xmm0,0x10(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x14(%rdx),%xmm0
  movss  %xmm0,0x14(%rsi)
  movss  0x8(%rbp),%xmm0
  addss  0x10(%rax),%xmm0
  addss  0x14(%rax),%xmm0
  movss  %xmm0,0x28(%rbx)
  movss  0x8(%r11),%xmm0
  addss  0x14(%rcx),%xmm0
  movss  0x14(%rax),%xmm4
  mulss  0x14(%r8),%xmm4
  addss  %xmm4,%xmm0
  addss  0x14(%rsi),%xmm0
  subss  0x10(%rcx),%xmm0
  movss  0x10(%rax),%xmm4
  mulss  0x10(%r8),%xmm4
  subss  %xmm4,%xmm0
  subss  0x10(%rsi),%xmm0
  movss  %xmm0,0x20(%r10)
  movss  %xmm0,0x8(%r9)
  movaps %xmm2,%xmm0
  divss  0x18(%rdi),%xmm0
  movss  %xmm0,0x18(%rdx)
  movaps %xmm2,%xmm4
  divss  0x1c(%rdi),%xmm4
  movss  %xmm4,0x1c(%rdx)
  mulss  0x18(%rcx),%xmm0
  movss  %xmm0,0x18(%rax)
  movss  0x1c(%rcx),%xmm0
  mulss  0x1c(%rdx),%xmm0
  movss  %xmm0,0x1c(%rax)
  movaps %xmm1,%xmm0
  mulss  0x18(%rdx),%xmm0
  movss  %xmm0,0x18(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x1c(%rdx),%xmm0
  movss  %xmm0,0x1c(%rsi)
  movss  0xc(%rbp),%xmm0
  addss  0x18(%rax),%xmm0
  addss  0x1c(%rax),%xmm0
  movss  %xmm0,0x3c(%rbx)
  movss  0xc(%r11),%xmm0
  addss  0x1c(%rcx),%xmm0
  movss  0x1c(%rax),%xmm4
  mulss  0x1c(%r8),%xmm4
  addss  %xmm4,%xmm0
  addss  0x1c(%rsi),%xmm0
  subss  0x18(%rcx),%xmm0
  movss  0x18(%rax),%xmm4
  mulss  0x18(%r8),%xmm4
  subss  %xmm4,%xmm0
  subss  0x18(%rsi),%xmm0
  movss  %xmm0,0x30(%r10)
  movss  %xmm0,0xc(%r9)
  add    $0x20,%rdi
  add    $0x20,%rcx
  add    $0x20,%rax
  add    $0x20,%rsi
  add    $0x20,%rdx
  add    $0x10,%rbp
  add    -0x78(%rsp),%rbx
  add    $0x10,%r11
  add    $0x20,%r8
  add    $0x40,%r10
  add    $0x10,%r9
  cmp    %r12,%rdi
  jne    L65
  mov    -0x50(%rsp),%r12
  mov    -0x58(%rsp),%rax
  add    %rax,%r15
  add    %rax,%r12
  add    %rax,%r14
  add    %rax,%r13
  mov    -0x8(%rsp),%edx
L9:
  cmp    %edx,-0x64(%rsp)
  jle    L196
  mov    -0x64(%rsp),%esi
  sub    %edx,%esi
  test   %esi,%esi
  jle    L196
  lea    0x0(,%rdx,4),%eax
  cltq
  mov    -0x4(%rsp),%ecx
  imul   %edx,%ecx
  movslq %ecx,%rcx
  add    %rax,%rcx
  mov    -0x48(%rsp),%rbx
  lea    (%rbx,%rcx,4),%rcx
  movslq %edx,%rdx
  lea    -0x1(%rsi),%esi
  lea    0x1(%rdx,%rsi,1),%rsi
  shl    $0x2,%rsi
  neg    %rdx
  mov    -0x60(%rsp),%rbx
  lea    (%rbx,%rdx,8),%rdx
  mov    -0x40(%rsp),%r8
  mov    -0x38(%rsp),%r9
  mov    -0x30(%rsp),%r10
  mov    -0x28(%rsp),%r11
  mov    -0x20(%rsp),%rdi
  jmp    L140
L23:
  movslq 0x4(%rsp),%rax
  mov    (%rsp),%edi
  lea    (%rdi,%rdi,1),%r8d
  movslq %r8d,%r8
  shl    $0x2,%r8
  mov    -0x18(%rsp),%rbx
  mov    %r8,%r9
  add    (%rbx,%rax,8),%r9
  mov    %r9,-0x70(%rsp)
  mov    0x48(%rsp),%rbx
  mov    %r8,%rcx
  add    (%rbx,%rax,8),%rcx
  mov    %rcx,-0x58(%rsp)
  mov    0x50(%rsp),%rbx
  mov    %r8,%r14
  add    (%rbx,%rax,8),%r14
  mov    0x58(%rsp),%rbx
  mov    %r8,%rsi
  add    (%rbx,%rax,8),%rsi
  mov    %rsi,-0x60(%rsp)
  mov    0x40(%rsp),%rbx
  mov    %r8,%rdx
  add    (%rbx,%rax,8),%rdx
  mov    %rdx,-0x50(%rsp)
  mov    0x70(%rsp),%rbx
  mov    (%rbx,%rax,8),%r11
  mov    %r11,-0x38(%rsp)
  mov    0x78(%rsp),%rbx
  mov    (%rbx,%rax,8),%rbx
  mov    0x80(%rsp),%r15
  mov    (%r15,%rax,8),%rbp
  mov    0x60(%rsp),%r15
  mov    (%r15,%rax,8),%r12
  mov    0x68(%rsp),%r15
  mov    (%r15,%rax,8),%r13
  mov    0x88(%rsp),%r15
  mov    (%r15,%rax,8),%rax
  mov    %rax,-0x30(%rsp)
  mov    -0x64(%rsp),%r15d
  sub    $0x3,%r15d
  mov    %r15d,-0x48(%rsp)
  cmp    %r15d,%edi
  jge    L137
  mov    %edi,%r15d
  movslq %edi,%rdi
  add    %rax,%r8
  mov    -0x4(%rsp),%r10d
  lea    0x10(,%r10,4),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x78(%rsp)
  mov    %r10d,%eax
  imul   %r15d,%eax
  cltq
  lea    (%rax,%rdi,4),%rax
  lea    (%r11,%rax,4),%r15
  mov    %r14,%rax
  mov    $0x0,%r10d
  movss  D20(%rip),%xmm2
  mov    %r14,-0x40(%rsp)
  mov    -0x48(%rsp),%r14d
L5:
  movaps %xmm2,%xmm0
  divss  (%r9),%xmm0
  movss  %xmm0,(%rdx)
  movaps %xmm2,%xmm3
  divss  0x4(%r9),%xmm3
  movss  %xmm3,0x4(%rdx)
  mulss  (%rcx),%xmm0
  movss  %xmm0,(%rax)
  movss  0x4(%rcx),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  movaps %xmm1,%xmm0
  mulss  (%rdx),%xmm0
  movss  %xmm0,(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rsi)
  movss  (%r12,%rdi,4),%xmm0
  addss  (%rax),%xmm0
  addss  0x4(%rax),%xmm0
  movss  %xmm0,(%r15,%r10,1)
  movss  0x0(%r13,%rdi,4),%xmm0
  addss  0x4(%rcx),%xmm0
  movss  0x4(%rax),%xmm3
  mulss  0x4(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0x4(%rsi),%xmm0
  subss  (%rcx),%xmm0
  movss  (%rax),%xmm3
  mulss  (%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  (%rsi),%xmm0
  mov    %rdi,%r11
  shl    $0x4,%r11
  movss  %xmm0,(%rbx,%r11,1)
  movss  %xmm0,0x0(%rbp,%rdi,4)
  movaps %xmm2,%xmm0
  divss  0x8(%r9),%xmm0
  movss  %xmm0,0x8(%rdx)
  movaps %xmm2,%xmm3
  divss  0xc(%r9),%xmm3
  movss  %xmm3,0xc(%rdx)
  mulss  0x8(%rcx),%xmm0
  movss  %xmm0,0x8(%rax)
  movss  0xc(%rcx),%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  %xmm0,0xc(%rax)
  movaps %xmm1,%xmm0
  mulss  0x8(%rdx),%xmm0
  movss  %xmm0,0x8(%rsi)
  movaps %xmm1,%xmm0
  mulss  0xc(%rdx),%xmm0
  movss  %xmm0,0xc(%rsi)
  movss  0x4(%r12,%rdi,4),%xmm0
  addss  0x8(%rax),%xmm0
  addss  0xc(%rax),%xmm0
  movss  %xmm0,0x14(%r15,%r10,1)
  movss  0x4(%r13,%rdi,4),%xmm0
  addss  0xc(%rcx),%xmm0
  movss  0xc(%rax),%xmm3
  mulss  0xc(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0xc(%rsi),%xmm0
  subss  0x8(%rcx),%xmm0
  movss  0x8(%rax),%xmm3
  mulss  0x8(%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  0x8(%rsi),%xmm0
  movss  %xmm0,0x10(%rbx,%r11,1)
  movss  %xmm0,0x4(%rbp,%rdi,4)
  movaps %xmm2,%xmm0
  divss  0x10(%r9),%xmm0
  movss  %xmm0,0x10(%rdx)
  movaps %xmm2,%xmm3
  divss  0x14(%r9),%xmm3
  movss  %xmm3,0x14(%rdx)
  mulss  0x10(%rcx),%xmm0
  movss  %xmm0,0x10(%rax)
  movss  0x14(%rcx),%xmm0
  mulss  0x14(%rdx),%xmm0
  movss  %xmm0,0x14(%rax)
  movaps %xmm1,%xmm0
  mulss  0x10(%rdx),%xmm0
  movss  %xmm0,0x10(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x14(%rdx),%xmm0
  movss  %xmm0,0x14(%rsi)
  movss  0x8(%r12,%rdi,4),%xmm0
  addss  0x10(%rax),%xmm0
  addss  0x14(%rax),%xmm0
  movss  %xmm0,0x28(%r15,%r10,1)
  movss  0x8(%r13,%rdi,4),%xmm0
  addss  0x14(%rcx),%xmm0
  movss  0x14(%rax),%xmm3
  mulss  0x14(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0x14(%rsi),%xmm0
  subss  0x10(%rcx),%xmm0
  movss  0x10(%rax),%xmm3
  mulss  0x10(%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  0x10(%rsi),%xmm0
  movss  %xmm0,0x20(%rbx,%r11,1)
  movss  %xmm0,0x8(%rbp,%rdi,4)
  movaps %xmm2,%xmm0
  divss  0x18(%r9),%xmm0
  movss  %xmm0,0x18(%rdx)
  movaps %xmm2,%xmm3
  divss  0x1c(%r9),%xmm3
  movss  %xmm3,0x1c(%rdx)
  mulss  0x18(%rcx),%xmm0
  movss  %xmm0,0x18(%rax)
  movss  0x1c(%rcx),%xmm0
  mulss  0x1c(%rdx),%xmm0
  movss  %xmm0,0x1c(%rax)
  movaps %xmm1,%xmm0
  mulss  0x18(%rdx),%xmm0
  movss  %xmm0,0x18(%rsi)
  movaps %xmm1,%xmm0
  mulss  0x1c(%rdx),%xmm0
  movss  %xmm0,0x1c(%rsi)
  movss  0xc(%r12,%rdi,4),%xmm0
  addss  0x18(%rax),%xmm0
  addss  0x1c(%rax),%xmm0
  movss  %xmm0,0x3c(%r15,%r10,1)
  movss  0xc(%r13,%rdi,4),%xmm0
  addss  0x1c(%rcx),%xmm0
  movss  0x1c(%rax),%xmm3
  mulss  0x1c(%r8),%xmm3
  addss  %xmm3,%xmm0
  addss  0x1c(%rsi),%xmm0
  subss  0x18(%rcx),%xmm0
  movss  0x18(%rax),%xmm3
  mulss  0x18(%r8),%xmm3
  subss  %xmm3,%xmm0
  subss  0x18(%rsi),%xmm0
  movss  %xmm0,0x30(%rbx,%r11,1)
  movss  %xmm0,0xc(%rbp,%rdi,4)
  add    $0x20,%r9
  add    $0x20,%rcx
  add    $0x20,%rax
  add    $0x20,%rsi
  add    $0x20,%rdx
  add    $0x4,%rdi
  add    $0x20,%r8
  add    -0x78(%rsp),%r10
  cmp    %edi,%r14d
  jg     L5
  mov    -0x40(%rsp),%r14
  mov    -0x64(%rsp),%eax
  lea    -0x4(%rax),%edx
  mov    (%rsp),%ecx
  sub    %ecx,%edx
  shr    $0x2,%edx
  lea    0x1(%rdx),%eax
  mov    %eax,%eax
  shl    $0x5,%rax
  add    %rax,-0x70(%rsp)
  add    %rax,-0x58(%rsp)
  add    %rax,%r14
  add    %rax,-0x60(%rsp)
  add    %rax,-0x50(%rsp)
  lea    0x4(%rcx,%rdx,4),%eax
  mov    %eax,(%rsp)
L137:
  mov    -0x64(%rsp),%esi
  mov    (%rsp),%edi
  cmp    %edi,%esi
  jle    L123
  lea    0x0(,%rdi,4),%eax
  cltq
  mov    -0x4(%rsp),%edx
  imul   %edi,%edx
  movslq %edx,%rdx
  add    %rax,%rdx
  mov    -0x38(%rsp),%rcx
  lea    (%rcx,%rdx,4),%rdx
  movslq %edi,%rcx
  sub    %edi,%esi
  lea    -0x1(%rsi),%esi
  lea    0x1(%rcx,%rsi,1),%rsi
  shl    $0x2,%rsi
  neg    %rcx
  mov    -0x70(%rsp),%rdi
  lea    (%rdi,%rcx,8),%rcx
  movss  D20(%rip),%xmm3
  mov    -0x60(%rsp),%r8
  mov    -0x58(%rsp),%r9
  mov    -0x50(%rsp),%rdi
  mov    -0x30(%rsp),%r10
L129:
  movaps %xmm3,%xmm0
  divss  (%rcx,%rax,2),%xmm0
  movss  %xmm0,(%rdi)
  movaps %xmm3,%xmm2
  divss  0x4(%rcx,%rax,2),%xmm2
  movss  %xmm2,0x4(%rdi)
  mulss  (%r9),%xmm0
  movss  %xmm0,(%r14)
  movss  0x4(%r9),%xmm0
  mulss  0x4(%rdi),%xmm0
  movss  %xmm0,0x4(%r14)
  movaps %xmm1,%xmm0
  mulss  (%rdi),%xmm0
  movss  %xmm0,(%r8)
  movaps %xmm1,%xmm0
  mulss  0x4(%rdi),%xmm0
  movss  %xmm0,0x4(%r8)
  movss  (%r12,%rax,1),%xmm0
  addss  (%r14),%xmm0
  addss  0x4(%r14),%xmm0
  movss  %xmm0,(%rdx)
  movss  0x0(%r13,%rax,1),%xmm0
  addss  0x4(%r9),%xmm0
  movss  0x4(%r14),%xmm2
  mulss  0x4(%r10,%rax,2),%xmm2
  addss  %xmm2,%xmm0
  addss  0x4(%r8),%xmm0
  subss  (%r9),%xmm0
  movss  (%r14),%xmm2
  mulss  (%r10,%rax,2),%xmm2
  subss  %xmm2,%xmm0
  subss  (%r8),%xmm0
  movss  %xmm0,(%rbx,%rax,4)
  movss  %xmm0,0x0(%rbp,%rax,1)
  add    $0x8,%r9
  add    $0x8,%r14
  add    $0x8,%r8
  add    $0x8,%rdi
  add    $0x14,%rdx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L129
L123:
  add    $0x8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D20
datatype: f32

tool: D20: 1

assistant: ```c
void s_update_hessian_box_mpc(int N, int k0, int k1, int kmax, int cnz, float sigma_mu, float **t, float **t_inv, float **lam, float **lamt, float **dlam, float **bd, float **bl, float **pd, float **pl, float **pl2, float **db) {
  const int bs = 4;
  float temp0, temp1;
  float *ptr_t, *ptr_lam, *ptr_lamt, *ptr_dlam, *ptr_tinv, *ptr_pd, *ptr_pl, *ptr_pl2, *ptr_bd, *ptr_bl, *ptr_db;
  int ii, jj, ll, bs0;
  ptr_t = t[0];
  ptr_lam = lam[0];
  ptr_lamt = lamt[0];
  ptr_dlam = dlam[0];
  ptr_tinv = t_inv[0];
  ptr_pd = pd[0];
  ptr_pl = pl[0];
  ptr_pl2 = pl2[0];
  ptr_bd = bd[0];
  ptr_bl = bl[0];
  ptr_db = db[0];
  ii = 0;
  for (; ii < k0 - 3; ii += 4) {
    ptr_tinv[0] = 1.0 / ptr_t[0];
    ptr_tinv[1] = 1.0 / ptr_t[1];
    ptr_lamt[0] = ptr_lam[0] * ptr_tinv[0];
    ptr_lamt[1] = ptr_lam[1] * ptr_tinv[1];
    ptr_dlam[0] = ptr_tinv[0] * sigma_mu;
    ptr_dlam[1] = ptr_tinv[1] * sigma_mu;
    ptr_pd[0 + (ii + 0) * bs + ii * cnz] = ptr_bd[ii + 0] + ptr_lamt[0] + ptr_lamt[1];
    ptr_pl[(ii + 0) * bs] = ptr_bl[ii + 0] + ptr_lam[1] + ptr_lamt[1] * ptr_db[2 * ii + 1] + ptr_dlam[1] - ptr_lam[0] - ptr_lamt[0] * ptr_db[2 * ii + 0] - ptr_dlam[0];
    ptr_pl2[ii + 0] = ptr_pl[(ii + 0) * bs];
    ptr_tinv[2] = 1.0 / ptr_t[2];
    ptr_tinv[3] = 1.0 / ptr_t[3];
    ptr_lamt[2] = ptr_lam[2] * ptr_tinv[2];
    ptr_lamt[3] = ptr_lam[3] * ptr_tinv[3];
    ptr_dlam[2] = ptr_tinv[2] * sigma_mu;
    ptr_dlam[3] = ptr_tinv[3] * sigma_mu;
    ptr_pd[1 + (ii + 1) * bs + ii * cnz] = ptr_bd[ii + 1] + ptr_lamt[2] + ptr_lamt[3];
    ptr_pl[(ii + 1) * bs] = ptr_bl[ii + 1] + ptr_lam[3] + ptr_lamt[3] * ptr_db[2 * ii + 3] + ptr_dlam[3] - ptr_lam[2] - ptr_lamt[2] * ptr_db[2 * ii + 2] - ptr_dlam[2];
    ptr_pl2[ii + 1] = ptr_pl[(ii + 1) * bs];
    ptr_tinv[4] = 1.0 / ptr_t[4];
    ptr_tinv[5] = 1.0 / ptr_t[5];
    ptr_lamt[4] = ptr_lam[4] * ptr_tinv[4];
    ptr_lamt[5] = ptr_lam[5] * ptr_tinv[5];
    ptr_dlam[4] = ptr_tinv[4] * sigma_mu;
    ptr_dlam[5] = ptr_tinv[5] * sigma_mu;
    ptr_pd[2 + (ii + 2) * bs + ii * cnz] = ptr_bd[ii + 2] + ptr_lamt[4] + ptr_lamt[5];
    ptr_pl[(ii + 2) * bs] = ptr_bl[ii + 2] + ptr_lam[5] + ptr_lamt[5] * ptr_db[2 * ii + 5] + ptr_dlam[5] - ptr_lam[4] - ptr_lamt[4] * ptr_db[2 * ii + 4] - ptr_dlam[4];
    ptr_pl2[ii + 2] = ptr_pl[(ii + 2) * bs];
    ptr_tinv[6] = 1.0 / ptr_t[6];
    ptr_tinv[7] = 1.0 / ptr_t[7];
    ptr_lamt[6] = ptr_lam[6] * ptr_tinv[6];
    ptr_lamt[7] = ptr_lam[7] * ptr_tinv[7];
    ptr_dlam[6] = ptr_tinv[6] * sigma_mu;
    ptr_dlam[7] = ptr_tinv[7] * sigma_mu;
    ptr_pd[3 + (ii + 3) * bs + ii * cnz] = ptr_bd[ii + 3] + ptr_lamt[6] + ptr_lamt[7];
    ptr_pl[(ii + 3) * bs] = ptr_bl[ii + 3] + ptr_lam[7] + ptr_lamt[7] * ptr_db[2 * ii + 7] + ptr_dlam[7] - ptr_lam[6] - ptr_lamt[6] * ptr_db[2 * ii + 6] - ptr_dlam[6];
    ptr_pl2[ii + 3] = ptr_pl[(ii + 3) * bs];
    ptr_t += 8;
    ptr_lam += 8;
    ptr_lamt += 8;
    ptr_dlam += 8;
    ptr_tinv += 8;
  }
  if (ii < k0) {
    bs0 = k0 - ii;
    for (ll = 0; ll < bs0; ll++) {
      ptr_tinv[0] = 1.0 / ptr_t[0];
      ptr_tinv[1] = 1.0 / ptr_t[1];
      ptr_lamt[0] = ptr_lam[0] * ptr_tinv[0];
      ptr_lamt[1] = ptr_lam[1] * ptr_tinv[1];
      ptr_dlam[0] = ptr_tinv[0] * sigma_mu;
      ptr_dlam[1] = ptr_tinv[1] * sigma_mu;
      ptr_pd[ll + (ii + ll) * bs + ii * cnz] = ptr_bd[ii + ll] + ptr_lamt[0] + ptr_lamt[1];
      ptr_pl[(ii + ll) * bs] = ptr_bl[ii + ll] + ptr_lam[1] + ptr_lamt[1] * ptr_db[2 * ii + 2 * ll + 1] + ptr_dlam[1] - ptr_lam[0] - ptr_lamt[0] * ptr_db[2 * ii + 2 * ll + 0] - ptr_dlam[0];
      ptr_pl2[ii + ll + 0] = ptr_pl[(ii + ll) * bs];
      ptr_t += 2;
      ptr_lam += 2;
      ptr_lamt += 2;
      ptr_dlam += 2;
      ptr_tinv += 2;
    }
  }
  for (jj = 1; jj < N; jj++) {
    ptr_t = t[jj];
    ptr_lam = lam[jj];
    ptr_lamt = lamt[jj];
    ptr_dlam = dlam[jj];
    ptr_tinv = t_inv[jj];
    ptr_pd = pd[jj];
    ptr_pl = pl[jj];
    ptr_pl2 = pl2[jj];
    ptr_bd = bd[jj];
    ptr_bl = bl[jj];
    ptr_db = db[jj];
    ii = 0;
    for (; ii < kmax - 3; ii += 4) {
      ptr_tinv[0] = 1.0 / ptr_t[0];
      ptr_tinv[1] = 1.0 / ptr_t[1];
      ptr_lamt[0] = ptr_lam[0] * ptr_tinv[0];
      ptr_lamt[1] = ptr_lam[1] * ptr_tinv[1];
      ptr_dlam[0] = ptr_tinv[0] * sigma_mu;
      ptr_dlam[1] = ptr_tinv[1] * sigma_mu;
      ptr_pd[0 + (ii + 0) * bs + ii * cnz] = ptr_bd[ii + 0] + ptr_lamt[0] + ptr_lamt[1];
      ptr_pl[(ii + 0) * bs] = ptr_bl[ii + 0] + ptr_lam[1] + ptr_lamt[1] * ptr_db[2 * ii + 1] + ptr_dlam[1] - ptr_lam[0] - ptr_lamt[0] * ptr_db[2 * ii + 0] - ptr_dlam[0];
      ptr_pl2[ii + 0] = ptr_pl[(ii + 0) * bs];
      ptr_tinv[2] = 1.0 / ptr_t[2];
      ptr_tinv[3] = 1.0 / ptr_t[3];
      ptr_lamt[2] = ptr_lam[2] * ptr_tinv[2];
      ptr_lamt[3] = ptr_lam[3] * ptr_tinv[3];
      ptr_dlam[2] = ptr_tinv[2] * sigma_mu;
      ptr_dlam[3] = ptr_tinv[3] * sigma_mu;
      ptr_pd[1 + (ii + 1) * bs + ii * cnz] = ptr_bd[ii + 1] + ptr_lamt[2] + ptr_lamt[3];
      ptr_pl[(ii + 1) * bs] = ptr_bl[ii + 1] + ptr_lam[3] + ptr_lamt[3] * ptr_db[2 * ii + 3] + ptr_dlam[3] - ptr_lam[2] - ptr_lamt[2] * ptr_db[2 * ii + 2] - ptr_dlam[2];
      ptr_pl2[ii + 1] = ptr_pl[(ii + 1) * bs];
      ptr_tinv[4] = 1.0 / ptr_t[4];
      ptr_tinv[5] = 1.0 / ptr_t[5];
      ptr_lamt[4] = ptr_lam[4] * ptr_tinv[4];
      ptr_lamt[5] = ptr_lam[5] * ptr_tinv[5];
      ptr_dlam[4] = ptr_tinv[4] * sigma_mu;
      ptr_dlam[5] = ptr_tinv[5] * sigma_mu;
      ptr_pd[2 + (ii + 2) * bs + ii * cnz] = ptr_bd[ii + 2] + ptr_lamt[4] + ptr_lamt[5];
      ptr_pl[(ii + 2) * bs] = ptr_bl[ii + 2] + ptr_lam[5] + ptr_lamt[5] * ptr_db[2 * ii + 5] + ptr_dlam[5] - ptr_lam[4] - ptr_lamt[4] * ptr_db[2 * ii + 4] - ptr_dlam[4];
      ptr_pl2[ii + 2] = ptr_pl[(ii + 2) * bs];
      ptr_tinv[6] = 1.0 / ptr_t[6];
      ptr_tinv[7] = 1.0 / ptr_t[7];
      ptr_lamt[6] = ptr_lam[6] * ptr_tinv[6];
      ptr_lamt[7] = ptr_lam[7] * ptr_tinv[7];
      ptr_dlam[6] = ptr_tinv[6] * sigma_mu;
      ptr_dlam[7] = ptr_tinv[7] * sigma_mu;
      ptr_pd[3 + (ii + 3) * bs + ii * cnz] = ptr_bd[ii + 3] + ptr_lamt[6] + ptr_lamt[7];
      ptr_pl[(ii + 3) * bs] = ptr_bl[ii + 3] + ptr_lam[7] + ptr_lamt[7] * ptr_db[2 * ii + 7] + ptr_dlam[7] - ptr_lam[6] - ptr_lamt[6] * ptr_db[2 * ii + 6] - ptr_dlam[6];
      ptr_pl2[ii + 3] = ptr_pl[(ii + 3) * bs];
      ptr_t += 8;
      ptr_lam += 8;
      ptr_lamt += 8;
      ptr_dlam += 8;
      ptr_tinv += 8;
    }
    if (ii < kmax) {
      bs0 = kmax - ii;
      for (ll = 0; ll < bs0; ll++) {
        ptr_tinv[0] = 1.0 / ptr_t[0];
        ptr_tinv[1] = 1.0 / ptr_t[1];
        ptr_lamt[0] = ptr_lam[0] * ptr_tinv[0];
        ptr_lamt[1] = ptr_lam[1] * ptr_tinv[1];
        ptr_dlam[0] = ptr_tinv[0] * sigma_mu;
        ptr_dlam[1] = ptr_tinv[1] * sigma_mu;
        ptr_pd[ll + (ii + ll) * bs + ii * cnz] = ptr_bd[ii + ll] + ptr_lamt[0] + ptr_lamt[1];
        ptr_pl[(ii + ll) * bs] = ptr_bl[ii + ll] + ptr_lam[1] + ptr_lamt[1] * ptr_db[2 * ii + 2 * ll + 1] + ptr_dlam[1] - ptr_lam[0] - ptr_lamt[0] * ptr_db[2 * ii + 2 * ll + 0] - ptr_dlam[0];
        ptr_pl2[ii + ll + 0] = ptr_pl[(ii + ll) * bs];
        ptr_t += 2;
        ptr_lam += 2;
        ptr_lamt += 2;
        ptr_dlam += 2;
        ptr_tinv += 2;
      }
    }
  }
  ptr_t = t[N] + 2 * k1;
  ptr_lam = lam[N] + 2 * k1;
  ptr_lamt = lamt[N] + 2 * k1;
  ptr_dlam = dlam[N] + 2 * k1;
  ptr_tinv = t_inv[N] + 2 * k1;
  ptr_pd = pd[N];
  ptr_pl = pl[N];
  ptr_pl2 = pl2[N];
  ptr_bd = bd[N];
  ptr_bl = bl[N];
  ptr_db = db[N];
  ii = k1;
  for (; ii < kmax - 3; ii += 4) {
    ptr_tinv[0] = 1.0 / ptr_t[0];
    ptr_tinv[1] = 1.0 / ptr_t[1];
    ptr_lamt[0] = ptr_lam[0] * ptr_tinv[0];
    ptr_lamt[1] = ptr_lam[1] * ptr_tinv[1];
    ptr_dlam[0] = ptr_tinv[0] * sigma_mu;
    ptr_dlam[1] = ptr_tinv[1] * sigma_mu;
    ptr_pd[0 + (ii + 0) * bs + ii * cnz] = ptr_bd[ii + 0] + ptr_lamt[0] + ptr_lamt[1];
    ptr_pl[(ii + 0) * bs] = ptr_bl[ii + 0] + ptr_lam[1] + ptr_lamt[1] * ptr_db[2 * ii + 1] + ptr_dlam[1] - ptr_lam[0] - ptr_lamt[0] * ptr_db[2 * ii + 0] - ptr_dlam[0];
    ptr_pl2[ii + 0] = ptr_pl[(ii + 0) * bs];
    ptr_tinv[2] = 1.0 / ptr_t[2];
    ptr_tinv[3] = 1.0 / ptr_t[3];
    ptr_lamt[2] = ptr_lam[2] * ptr_tinv[2];
    ptr_lamt[3] = ptr_lam[3] * ptr_tinv[3];
    ptr_dlam[2] = ptr_tinv[2] * sigma_mu;
    ptr_dlam[3] = ptr_tinv[3] * sigma_mu;
    ptr_pd[1 + (ii + 1) * bs + ii * cnz] = ptr_bd[ii + 1] + ptr_lamt[2] + ptr_lamt[3];
    ptr_pl[(ii + 1) * bs] = ptr_bl[ii + 1] + ptr_lam[3] + ptr_lamt[3] * ptr_db[2 * ii + 3] + ptr_dlam[3] - ptr_lam[2] - ptr_lamt[2] * ptr_db[2 * ii + 2] - ptr_dlam[2];
    ptr_pl2[ii + 1] = ptr_pl[(ii + 1) * bs];
    ptr_tinv[4] = 1.0 / ptr_t[4];
    ptr_tinv[5] = 1.0 / ptr_t[5];
    ptr_lamt[4] = ptr_lam[4] * ptr_tinv[4];
    ptr_lamt[5] = ptr_lam[5] * ptr_tinv[5];
    ptr_dlam[4] = ptr_tinv[4] * sigma_mu;
    ptr_dlam[5] = ptr_tinv[5] * sigma_mu;
    ptr_pd[2 + (ii + 2) * bs + ii * cnz] = ptr_bd[ii + 2] + ptr_lamt[4] + ptr_lamt[5];
    ptr_pl[(ii + 2) * bs] = ptr_bl[ii + 2] + ptr_lam[5] + ptr_lamt[5] * ptr_db[2 * ii + 5] + ptr_dlam[5] - ptr_lam[4] - ptr_lamt[4] * ptr_db[2 * ii + 4] - ptr_dlam[4];
    ptr_pl2[ii + 2] = ptr_pl[(ii + 2) * bs];
    ptr_tinv[6] = 1.0 / ptr_t[6];
    ptr_tinv[7] = 1.0 / ptr_t[7];
    ptr_lamt[6] = ptr_lam[6] * ptr_tinv[6];
    ptr_lamt[7] = ptr_lam[7] * ptr_tinv[7];
    ptr_dlam[6] = ptr_tinv[6] * sigma_mu;
    ptr_dlam[7] = ptr_tinv[7] * sigma_mu;
    ptr_pd[3 + (ii + 3) * bs + ii * cnz] = ptr_bd[ii + 3] + ptr_lamt[6] + ptr_lamt[7];
    ptr_pl[(ii + 3) * bs] = ptr_bl[ii + 3] + ptr_lam[7] + ptr_lamt[7] * ptr_db[2 * ii + 7] + ptr_dlam[7] - ptr_lam[6] - ptr_lamt[6] * ptr_db[2 * ii + 6] - ptr_dlam[6];
    ptr_pl2[ii + 3] = ptr_pl[(ii + 3) * bs];
    ptr_t += 8;
    ptr_lam += 8;
    ptr_lamt += 8;
    ptr_dlam += 8;
    ptr_tinv += 8;
  }
  if (ii < kmax) {
    bs0 = kmax - ii;
    for (ll = 0; ll < bs0; ll++) {
      ptr_tinv[0] = 1.0 / ptr_t[0];
      ptr_tinv[1] = 1.0 / ptr_t[1];
      ptr_lamt[0] = ptr_lam[0] * ptr_tinv[0];
      ptr_lamt[1] = ptr_lam[1] * ptr_tinv[1];
      ptr_dlam[0] = ptr_tinv[0] * sigma_mu;
      ptr_dlam[1] = ptr_tinv[1] * sigma_mu;
      ptr_pd[ll + (ii + ll) * bs + ii * cnz] = ptr_bd[ii + ll] + ptr_lamt[0] + ptr_lamt[1];
      ptr_pl[(ii + ll) * bs] = ptr_bl[ii + ll] + ptr_lam[1] + ptr_lamt[1] * ptr_db[2 * ii + 2 * ll + 1] + ptr_dlam[1] - ptr_lam[0] - ptr_lamt[0] * ptr_db[2 * ii + 2 * ll + 0] - ptr_dlam[0];
      ptr_pl2[ii + ll + 0] = ptr_pl[(ii + ll) * bs];
      ptr_t += 2;
      ptr_lam += 2;
      ptr_lamt += 2;
      ptr_dlam += 2;
      ptr_tinv += 2;
    }
  }
}
```