user: Demonstrate the C code that generates the assembly output below:

```asm
<s_update_hessian_box_mpc>:
  endbr64
  push   %r15
  movaps %xmm0,%xmm1
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x10,%rsp
  mov    0x50(%rsp),%rax
  mov    %edi,-0x78(%rsp)
  mov    %r9,-0x18(%rsp)
  mov    (%r9),%rdi
  mov    (%rax),%r14
  mov    0x58(%rsp),%rax
  mov    %ecx,(%rsp)
  mov    %esi,-0x6c(%rsp)
  mov    (%rax),%r12
  mov    0x60(%rsp),%rax
  mov    %edx,0xc(%rsp)
  mov    %r8d,0x8(%rsp)
  mov    (%rax),%r15
  mov    0x48(%rsp),%rax
  mov    (%rax),%r13
  mov    0x78(%rsp),%rax
  mov    (%rax),%rbx
  mov    0x80(%rsp),%rax
  mov    (%rax),%r10
  mov    0x88(%rsp),%rax
  mov    %rbx,-0x40(%rsp)
  mov    (%rax),%r9
  mov    0x68(%rsp),%rax
  mov    %r10,-0x38(%rsp)
  mov    %r9,-0x30(%rsp)
  mov    (%rax),%rbp
  mov    0x70(%rsp),%rax
  mov    %rbp,-0x28(%rsp)
  mov    (%rax),%r11
  mov    0x90(%rsp),%rax
  mov    (%rax),%rcx
  mov    %r11,-0x20(%rsp)
  mov    %rcx,-0x10(%rsp)
  cmp    $0x3,%esi
  jle    L46
  lea    0x10(,%r8,4),%eax
  mov    %r14,-0x48(%rsp)
  mov    %rcx,%r8
  mov    %r14,%rcx
  cltq
  movss  D21(%rip),%xmm0
  mov    %r13,%rdx
  shl    $0x2,%rax
  mov    %rax,-0x60(%rsp)
  lea    -0x4(%rsi),%eax
  mov    -0x60(%rsp),%r14
  mov    %r15,%rsi
  shr    $0x2,%eax
  mov    %r12,-0x60(%rsp)
  add    $0x1,%eax
  mov    %eax,-0x50(%rsp)
  shl    $0x5,%rax
  mov    %rax,-0x58(%rsp)
  add    %rdi,%rax
  mov    %rax,-0x68(%rsp)
  mov    %r12,%rax
  mov    -0x68(%rsp),%r12
  nopl   (%rax)
L18:
  movaps %xmm0,%xmm2
  movaps %xmm0,%xmm3
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  add    $0x20,%rdi
  add    $0x20,%rcx
  add    $0x20,%rax
  add    $0x20,%rsi
  add    $0x20,%rdx
  add    $0x10,%rbp
  add    $0x10,%r11
  add    $0x20,%r8
  divss  -0x1c(%rdi),%xmm3
  movss  %xmm3,-0x1c(%rdx)
  mulss  -0x20(%rcx),%xmm2
  add    $0x40,%r10
  add    $0x10,%r9
  movss  %xmm2,-0x20(%rax)
  movss  -0x1c(%rcx),%xmm2
  mulss  -0x1c(%rdx),%xmm2
  movss  %xmm2,-0x1c(%rax)
  movss  -0x20(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x20(%rsi)
  movss  -0x1c(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x1c(%rsi)
  movss  -0x10(%rbp),%xmm2
  addss  -0x20(%rax),%xmm2
  addss  -0x1c(%rax),%xmm2
  movss  %xmm2,(%rbx)
  movss  -0x1c(%rax),%xmm3
  mulss  -0x1c(%r8),%xmm3
  movss  -0x10(%r11),%xmm2
  addss  -0x1c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x20(%rax),%xmm3
  mulss  -0x20(%r8),%xmm3
  addss  -0x1c(%rsi),%xmm2
  subss  -0x20(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x20(%rsi),%xmm2
  movaps %xmm0,%xmm3
  movss  %xmm2,-0x40(%r10)
  movss  %xmm2,-0x10(%r9)
  movaps %xmm0,%xmm2
  divss  -0x18(%rdi),%xmm2
  movss  %xmm2,-0x18(%rdx)
  divss  -0x14(%rdi),%xmm3
  movss  %xmm3,-0x14(%rdx)
  mulss  -0x18(%rcx),%xmm2
  movss  %xmm2,-0x18(%rax)
  movss  -0x14(%rcx),%xmm2
  mulss  -0x14(%rdx),%xmm2
  movss  %xmm2,-0x14(%rax)
  movss  -0x18(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x18(%rsi)
  movss  -0x14(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x14(%rsi)
  movss  -0xc(%rbp),%xmm2
  addss  -0x18(%rax),%xmm2
  addss  -0x14(%rax),%xmm2
  movss  %xmm2,0x14(%rbx)
  movss  -0x14(%rax),%xmm3
  mulss  -0x14(%r8),%xmm3
  movss  -0xc(%r11),%xmm2
  addss  -0x14(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x18(%rax),%xmm3
  mulss  -0x18(%r8),%xmm3
  addss  -0x14(%rsi),%xmm2
  subss  -0x18(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x18(%rsi),%xmm2
  movaps %xmm0,%xmm3
  movss  %xmm2,-0x30(%r10)
  movss  %xmm2,-0xc(%r9)
  movaps %xmm0,%xmm2
  divss  -0x10(%rdi),%xmm2
  movss  %xmm2,-0x10(%rdx)
  divss  -0xc(%rdi),%xmm3
  movss  %xmm3,-0xc(%rdx)
  mulss  -0x10(%rcx),%xmm2
  movss  %xmm2,-0x10(%rax)
  movss  -0xc(%rcx),%xmm2
  mulss  -0xc(%rdx),%xmm2
  movss  %xmm2,-0xc(%rax)
  movss  -0x10(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x10(%rsi)
  movss  -0xc(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0xc(%rsi)
  movss  -0x8(%rbp),%xmm2
  addss  -0x10(%rax),%xmm2
  addss  -0xc(%rax),%xmm2
  movss  %xmm2,0x28(%rbx)
  movss  -0xc(%rax),%xmm3
  mulss  -0xc(%r8),%xmm3
  movss  -0x8(%r11),%xmm2
  addss  -0xc(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x10(%rax),%xmm3
  mulss  -0x10(%r8),%xmm3
  addss  -0xc(%rsi),%xmm2
  subss  -0x10(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x10(%rsi),%xmm2
  movaps %xmm0,%xmm3
  movss  %xmm2,-0x20(%r10)
  movss  %xmm2,-0x8(%r9)
  movaps %xmm0,%xmm2
  divss  -0x8(%rdi),%xmm2
  movss  %xmm2,-0x8(%rdx)
  divss  -0x4(%rdi),%xmm3
  movss  %xmm3,-0x4(%rdx)
  mulss  -0x8(%rcx),%xmm2
  movss  %xmm2,-0x8(%rax)
  movss  -0x4(%rcx),%xmm2
  mulss  -0x4(%rdx),%xmm2
  movss  %xmm2,-0x4(%rax)
  movss  -0x8(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x8(%rsi)
  movss  -0x4(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x4(%rsi)
  movss  -0x4(%rbp),%xmm2
  addss  -0x8(%rax),%xmm2
  addss  -0x4(%rax),%xmm2
  movss  %xmm2,0x3c(%rbx)
  movss  -0x4(%rax),%xmm3
  mulss  -0x4(%r8),%xmm3
  add    %r14,%rbx
  movss  -0x4(%r11),%xmm2
  addss  -0x4(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x8(%rax),%xmm3
  mulss  -0x8(%r8),%xmm3
  addss  -0x4(%rsi),%xmm2
  subss  -0x8(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x8(%rsi),%xmm2
  movss  %xmm2,-0x10(%r10)
  movss  %xmm2,-0x4(%r9)
  cmp    %r12,%rdi
  jne    L18
  mov    -0x58(%rsp),%rax
  mov    -0x48(%rsp),%r14
  mov    -0x60(%rsp),%r12
  mov    -0x50(%rsp),%ecx
  add    %rax,%r14
  add    %rax,%r15
  add    %rax,%r13
  add    %rax,%r12
  shl    $0x2,%ecx
L96:
  cmp    %ecx,-0x6c(%rsp)
  jle    L66
  mov    0x8(%rsp),%edx
  lea    0x0(,%rcx,4),%eax
  mov    -0x6c(%rsp),%esi
  movslq %ecx,%rdi
  mov    -0x40(%rsp),%rbx
  cltq
  mov    -0x38(%rsp),%r8
  imul   %ecx,%edx
  sub    %ecx,%esi
  mov    -0x30(%rsp),%r9
  mov    -0x28(%rsp),%r10
  lea    -0x1(%rsi),%ecx
  movss  D21(%rip),%xmm0
  mov    -0x20(%rsp),%r11
  lea    0x1(%rdi,%rcx,1),%rcx
  neg    %rdi
  movslq %edx,%rdx
  shl    $0x2,%rcx
  add    %rax,%rdx
  lea    (%rbx,%rdx,4),%rdx
  mov    -0x68(%rsp),%rbx
  lea    (%rbx,%rdi,8),%rsi
  mov    -0x10(%rsp),%rdi
  nopw   0x0(%rax,%rax,1)
L124:
  movaps %xmm0,%xmm2
  movaps %xmm0,%xmm3
  divss  (%rsi,%rax,2),%xmm2
  add    $0x8,%r14
  add    $0x8,%r12
  add    $0x8,%r15
  add    $0x8,%r13
  add    $0x14,%rdx
  movss  %xmm2,-0x8(%r13)
  divss  0x4(%rsi,%rax,2),%xmm3
  movss  %xmm3,-0x4(%r13)
  mulss  -0x8(%r14),%xmm2
  movss  %xmm2,-0x8(%r12)
  movss  -0x4(%r14),%xmm2
  mulss  -0x4(%r13),%xmm2
  movss  %xmm2,-0x4(%r12)
  movss  -0x8(%r13),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x8(%r15)
  movss  -0x4(%r13),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x4(%r15)
  movss  (%r10,%rax,1),%xmm2
  addss  -0x8(%r12),%xmm2
  addss  -0x4(%r12),%xmm2
  movss  %xmm2,-0x14(%rdx)
  movss  -0x4(%r12),%xmm3
  mulss  0x4(%rdi,%rax,2),%xmm3
  movss  (%r11,%rax,1),%xmm2
  addss  -0x4(%r14),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x8(%r12),%xmm3
  mulss  (%rdi,%rax,2),%xmm3
  addss  -0x4(%r15),%xmm2
  subss  -0x8(%r14),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x8(%r15),%xmm2
  movss  %xmm2,(%r8,%rax,4)
  movss  %xmm2,(%r9,%rax,1)
  add    $0x4,%rax
  cmp    %rax,%rcx
  jne    L124
L66:
  mov    (%rsp),%eax
  sub    $0x3,%eax
  cmpl   $0x1,-0x78(%rsp)
  mov    %eax,-0x6c(%rsp)
  movslq -0x78(%rsp),%rax
  mov    %rax,-0x10(%rsp)
  jle    L62
  mov    0x8(%rsp),%eax
  movss  D21(%rip),%xmm2
  movq   $0x1,-0x78(%rsp)
  lea    0x10(,%rax,4),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x8(%rsp)
  mov    (%rsp),%eax
  sub    $0x4,%eax
  shr    $0x2,%eax
  lea    0x1(%rax),%edx
  mov    %rdx,%rax
  shl    $0x5,%rdx
  shl    $0x2,%eax
  mov    %rdx,-0x60(%rsp)
  mov    %eax,0x4(%rsp)
  nopw   0x0(%rax,%rax,1)
L15:
  mov    -0x78(%rsp),%rax
  mov    -0x18(%rsp),%rbx
  mov    0x80(%rsp),%rsi
  mov    0x88(%rsp),%rcx
  mov    (%rbx,%rax,8),%rdi
  mov    0x50(%rsp),%rbx
  mov    (%rsi,%rax,8),%r10
  mov    0x68(%rsp),%rsi
  mov    (%rbx,%rax,8),%r14
  mov    0x58(%rsp),%rbx
  mov    (%rcx,%rax,8),%r9
  mov    (%rsi,%rax,8),%rbp
  mov    %r10,-0x40(%rsp)
  mov    (%rbx,%rax,8),%r12
  mov    0x60(%rsp),%rbx
  mov    0x70(%rsp),%rcx
  mov    0x90(%rsp),%rsi
  mov    %r9,-0x38(%rsp)
  mov    (%rbx,%rax,8),%r15
  mov    0x48(%rsp),%rbx
  mov    %rbp,-0x30(%rsp)
  mov    (%rcx,%rax,8),%r11
  mov    (%rsi,%rax,8),%r8
  mov    (%rbx,%rax,8),%r13
  mov    0x78(%rsp),%rbx
  mov    %r11,-0x28(%rsp)
  mov    (%rbx,%rax,8),%rbx
  mov    -0x6c(%rsp),%eax
  mov    %r8,-0x20(%rsp)
  mov    %rbx,-0x48(%rsp)
  test   %eax,%eax
  jle    L159
  mov    -0x60(%rsp),%rax
  mov    %r14,-0x58(%rsp)
  mov    %r14,%rcx
  mov    %r13,%rdx
  mov    %r12,-0x50(%rsp)
  mov    -0x8(%rsp),%r14
  mov    %r15,%rsi
  add    %rdi,%rax
  mov    %rax,-0x68(%rsp)
  mov    %r12,%rax
  mov    -0x68(%rsp),%r12
  nopl   0x0(%rax,%rax,1)
L164:
  movaps %xmm2,%xmm0
  movaps %xmm2,%xmm3
  divss  (%rdi),%xmm0
  movss  %xmm0,(%rdx)
  add    $0x20,%rdi
  add    $0x20,%rcx
  add    $0x20,%rax
  add    $0x20,%rsi
  add    $0x20,%rdx
  add    $0x10,%rbp
  add    $0x10,%r11
  add    $0x20,%r8
  divss  -0x1c(%rdi),%xmm3
  movss  %xmm3,-0x1c(%rdx)
  mulss  -0x20(%rcx),%xmm0
  add    $0x40,%r10
  add    $0x10,%r9
  movss  %xmm0,-0x20(%rax)
  movss  -0x1c(%rcx),%xmm0
  mulss  -0x1c(%rdx),%xmm0
  movss  %xmm0,-0x1c(%rax)
  movss  -0x20(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x20(%rsi)
  movss  -0x1c(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x1c(%rsi)
  movss  -0x10(%rbp),%xmm0
  addss  -0x20(%rax),%xmm0
  addss  -0x1c(%rax),%xmm0
  movss  %xmm0,(%rbx)
  movss  -0x1c(%rax),%xmm3
  mulss  -0x1c(%r8),%xmm3
  movss  -0x10(%r11),%xmm0
  addss  -0x1c(%rcx),%xmm0
  addss  %xmm3,%xmm0
  movss  -0x20(%rax),%xmm3
  mulss  -0x20(%r8),%xmm3
  addss  -0x1c(%rsi),%xmm0
  subss  -0x20(%rcx),%xmm0
  subss  %xmm3,%xmm0
  subss  -0x20(%rsi),%xmm0
  movaps %xmm2,%xmm3
  movss  %xmm0,-0x40(%r10)
  movss  %xmm0,-0x10(%r9)
  movaps %xmm2,%xmm0
  divss  -0x18(%rdi),%xmm0
  movss  %xmm0,-0x18(%rdx)
  divss  -0x14(%rdi),%xmm3
  movss  %xmm3,-0x14(%rdx)
  mulss  -0x18(%rcx),%xmm0
  movss  %xmm0,-0x18(%rax)
  movss  -0x14(%rcx),%xmm0
  mulss  -0x14(%rdx),%xmm0
  movss  %xmm0,-0x14(%rax)
  movss  -0x18(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x18(%rsi)
  movss  -0x14(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x14(%rsi)
  movss  -0xc(%rbp),%xmm0
  addss  -0x18(%rax),%xmm0
  addss  -0x14(%rax),%xmm0
  movss  %xmm0,0x14(%rbx)
  movss  -0x14(%rax),%xmm3
  mulss  -0x14(%r8),%xmm3
  movss  -0xc(%r11),%xmm0
  addss  -0x14(%rcx),%xmm0
  addss  %xmm3,%xmm0
  movss  -0x18(%rax),%xmm3
  mulss  -0x18(%r8),%xmm3
  addss  -0x14(%rsi),%xmm0
  subss  -0x18(%rcx),%xmm0
  subss  %xmm3,%xmm0
  subss  -0x18(%rsi),%xmm0
  movaps %xmm2,%xmm3
  movss  %xmm0,-0x30(%r10)
  movss  %xmm0,-0xc(%r9)
  movaps %xmm2,%xmm0
  divss  -0x10(%rdi),%xmm0
  movss  %xmm0,-0x10(%rdx)
  divss  -0xc(%rdi),%xmm3
  movss  %xmm3,-0xc(%rdx)
  mulss  -0x10(%rcx),%xmm0
  movss  %xmm0,-0x10(%rax)
  movss  -0xc(%rcx),%xmm0
  mulss  -0xc(%rdx),%xmm0
  movss  %xmm0,-0xc(%rax)
  movss  -0x10(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x10(%rsi)
  movss  -0xc(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0xc(%rsi)
  movss  -0x8(%rbp),%xmm0
  addss  -0x10(%rax),%xmm0
  addss  -0xc(%rax),%xmm0
  movss  %xmm0,0x28(%rbx)
  movss  -0xc(%rax),%xmm3
  mulss  -0xc(%r8),%xmm3
  movss  -0x8(%r11),%xmm0
  addss  -0xc(%rcx),%xmm0
  addss  %xmm3,%xmm0
  movss  -0x10(%rax),%xmm3
  mulss  -0x10(%r8),%xmm3
  addss  -0xc(%rsi),%xmm0
  subss  -0x10(%rcx),%xmm0
  subss  %xmm3,%xmm0
  subss  -0x10(%rsi),%xmm0
  movaps %xmm2,%xmm3
  movss  %xmm0,-0x20(%r10)
  movss  %xmm0,-0x8(%r9)
  movaps %xmm2,%xmm0
  divss  -0x8(%rdi),%xmm0
  movss  %xmm0,-0x8(%rdx)
  divss  -0x4(%rdi),%xmm3
  movss  %xmm3,-0x4(%rdx)
  mulss  -0x8(%rcx),%xmm0
  movss  %xmm0,-0x8(%rax)
  movss  -0x4(%rcx),%xmm0
  mulss  -0x4(%rdx),%xmm0
  movss  %xmm0,-0x4(%rax)
  movss  -0x8(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x8(%rsi)
  movss  -0x4(%rdx),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rsi)
  movss  -0x4(%rbp),%xmm0
  addss  -0x8(%rax),%xmm0
  addss  -0x4(%rax),%xmm0
  movss  %xmm0,0x3c(%rbx)
  movss  -0x4(%rax),%xmm3
  mulss  -0x4(%r8),%xmm3
  add    %r14,%rbx
  movss  -0x4(%r11),%xmm0
  addss  -0x4(%rcx),%xmm0
  addss  %xmm3,%xmm0
  movss  -0x8(%rax),%xmm3
  mulss  -0x8(%r8),%xmm3
  addss  -0x4(%rsi),%xmm0
  subss  -0x8(%rcx),%xmm0
  subss  %xmm3,%xmm0
  subss  -0x8(%rsi),%xmm0
  movss  %xmm0,-0x10(%r10)
  movss  %xmm0,-0x4(%r9)
  cmp    %rdi,%r12
  jne    L164
  mov    -0x60(%rsp),%rax
  mov    -0x58(%rsp),%r14
  mov    -0x50(%rsp),%r12
  mov    0x4(%rsp),%ecx
  add    %rax,%r14
  add    %rax,%r15
  add    %rax,%r13
  add    %rax,%r12
L165:
  cmp    %ecx,(%rsp)
  jg     L161
L127:
  addq   $0x1,-0x78(%rsp)
  mov    -0x78(%rsp),%rax
  cmp    %rax,-0x10(%rsp)
  jne    L15
L62:
  mov    0xc(%rsp),%edi
  mov    -0x10(%rsp),%rax
  mov    -0x18(%rsp),%rcx
  mov    0x90(%rsp),%rsi
  lea    (%rdi,%rdi,1),%r8d
  mov    (%rcx,%rax,8),%rbx
  movslq %r8d,%r8
  mov    0x50(%rsp),%rcx
  shl    $0x2,%r8
  add    %r8,%rbx
  mov    %rbx,-0x60(%rsp)
  mov    (%rcx,%rax,8),%rbx
  mov    0x60(%rsp),%rcx
  add    %r8,%rbx
  mov    %rbx,-0x68(%rsp)
  mov    0x58(%rsp),%rbx
  mov    (%rbx,%rax,8),%r14
  mov    (%rcx,%rax,8),%rbx
  mov    0x48(%rsp),%rcx
  add    %r8,%rbx
  add    %r8,%r14
  mov    %rbx,-0x78(%rsp)
  mov    (%rcx,%rax,8),%rbx
  mov    0x88(%rsp),%rcx
  add    %r8,%rbx
  mov    %rbx,-0x58(%rsp)
  mov    0x78(%rsp),%rbx
  mov    (%rcx,%rax,8),%rbp
  mov    0x68(%rsp),%rcx
  mov    (%rbx,%rax,8),%rbx
  mov    (%rcx,%rax,8),%r12
  mov    0x70(%rsp),%rcx
  mov    %rbx,-0x50(%rsp)
  mov    0x80(%rsp),%rbx
  mov    (%rcx,%rax,8),%r13
  mov    (%rbx,%rax,8),%rbx
  mov    (%rsi,%rax,8),%rax
  mov    %rax,-0x48(%rsp)
  cmp    -0x6c(%rsp),%edi
  jge    L4
  mov    0x8(%rsp),%ecx
  movslq 0xc(%rsp),%rdi
  mov    %r14,-0x40(%rsp)
  xor    %r10d,%r10d
  mov    -0x58(%rsp),%rdx
  mov    -0x78(%rsp),%rsi
  lea    0x10(,%rcx,4),%eax
  mov    -0x60(%rsp),%r9
  add    -0x48(%rsp),%r8
  cltq
  movss  D21(%rip),%xmm0
  lea    0x0(,%rax,4),%r11
  mov    %ecx,%eax
  mov    -0x50(%rsp),%rcx
  imul   %edi,%eax
  cltq
  lea    (%rax,%rdi,4),%rax
  lea    (%rcx,%rax,4),%r15
  mov    -0x68(%rsp),%rcx
  mov    %r14,%rax
  mov    %r11,%r14
  nop
L32:
  movaps %xmm0,%xmm2
  movaps %xmm0,%xmm3
  divss  (%r9),%xmm2
  movss  %xmm2,(%rdx)
  mov    %rdi,%r11
  add    $0x20,%r9
  add    $0x20,%rcx
  add    $0x20,%rax
  shl    $0x4,%r11
  add    $0x20,%rsi
  add    $0x20,%rdx
  add    $0x20,%r8
  divss  -0x1c(%r9),%xmm3
  movss  %xmm3,-0x1c(%rdx)
  mulss  -0x20(%rcx),%xmm2
  movss  %xmm2,-0x20(%rax)
  movss  -0x1c(%rcx),%xmm2
  mulss  -0x1c(%rdx),%xmm2
  movss  %xmm2,-0x1c(%rax)
  movss  -0x20(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x20(%rsi)
  movss  -0x1c(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x1c(%rsi)
  movss  (%r12,%rdi,4),%xmm2
  addss  -0x20(%rax),%xmm2
  addss  -0x1c(%rax),%xmm2
  movss  %xmm2,(%r15,%r10,1)
  movss  -0x1c(%rax),%xmm3
  mulss  -0x1c(%r8),%xmm3
  movss  0x0(%r13,%rdi,4),%xmm2
  addss  -0x1c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x20(%rax),%xmm3
  mulss  -0x20(%r8),%xmm3
  addss  -0x1c(%rsi),%xmm2
  subss  -0x20(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x20(%rsi),%xmm2
  movaps %xmm0,%xmm3
  movss  %xmm2,(%rbx,%r11,1)
  movss  %xmm2,0x0(%rbp,%rdi,4)
  movaps %xmm0,%xmm2
  divss  -0x18(%r9),%xmm2
  movss  %xmm2,-0x18(%rdx)
  divss  -0x14(%r9),%xmm3
  movss  %xmm3,-0x14(%rdx)
  mulss  -0x18(%rcx),%xmm2
  movss  %xmm2,-0x18(%rax)
  movss  -0x14(%rcx),%xmm2
  mulss  -0x14(%rdx),%xmm2
  movss  %xmm2,-0x14(%rax)
  movss  -0x18(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x18(%rsi)
  movss  -0x14(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x14(%rsi)
  movss  0x4(%r12,%rdi,4),%xmm2
  addss  -0x18(%rax),%xmm2
  addss  -0x14(%rax),%xmm2
  movss  %xmm2,0x14(%r15,%r10,1)
  movss  -0x14(%rax),%xmm3
  mulss  -0x14(%r8),%xmm3
  movss  0x4(%r13,%rdi,4),%xmm2
  addss  -0x14(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x18(%rax),%xmm3
  mulss  -0x18(%r8),%xmm3
  addss  -0x14(%rsi),%xmm2
  subss  -0x18(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x18(%rsi),%xmm2
  movaps %xmm0,%xmm3
  movss  %xmm2,0x10(%rbx,%r11,1)
  movss  %xmm2,0x4(%rbp,%rdi,4)
  movaps %xmm0,%xmm2
  divss  -0x10(%r9),%xmm2
  movss  %xmm2,-0x10(%rdx)
  divss  -0xc(%r9),%xmm3
  movss  %xmm3,-0xc(%rdx)
  mulss  -0x10(%rcx),%xmm2
  movss  %xmm2,-0x10(%rax)
  movss  -0xc(%rcx),%xmm2
  mulss  -0xc(%rdx),%xmm2
  movss  %xmm2,-0xc(%rax)
  movss  -0x10(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x10(%rsi)
  movss  -0xc(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0xc(%rsi)
  movss  0x8(%r12,%rdi,4),%xmm2
  addss  -0x10(%rax),%xmm2
  addss  -0xc(%rax),%xmm2
  movss  %xmm2,0x28(%r15,%r10,1)
  movss  -0xc(%rax),%xmm3
  mulss  -0xc(%r8),%xmm3
  movss  0x8(%r13,%rdi,4),%xmm2
  addss  -0xc(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x10(%rax),%xmm3
  mulss  -0x10(%r8),%xmm3
  addss  -0xc(%rsi),%xmm2
  subss  -0x10(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x10(%rsi),%xmm2
  movaps %xmm0,%xmm3
  movss  %xmm2,0x20(%rbx,%r11,1)
  movss  %xmm2,0x8(%rbp,%rdi,4)
  movaps %xmm0,%xmm2
  divss  -0x8(%r9),%xmm2
  movss  %xmm2,-0x8(%rdx)
  divss  -0x4(%r9),%xmm3
  movss  %xmm3,-0x4(%rdx)
  mulss  -0x8(%rcx),%xmm2
  movss  %xmm2,-0x8(%rax)
  movss  -0x4(%rcx),%xmm2
  mulss  -0x4(%rdx),%xmm2
  movss  %xmm2,-0x4(%rax)
  movss  -0x8(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x8(%rsi)
  movss  -0x4(%rdx),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x4(%rsi)
  movss  0xc(%r12,%rdi,4),%xmm2
  addss  -0x8(%rax),%xmm2
  addss  -0x4(%rax),%xmm2
  movss  %xmm2,0x3c(%r15,%r10,1)
  movss  -0x4(%rax),%xmm3
  mulss  -0x4(%r8),%xmm3
  add    %r14,%r10
  movss  0xc(%r13,%rdi,4),%xmm2
  addss  -0x4(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x8(%rax),%xmm3
  mulss  -0x8(%r8),%xmm3
  addss  -0x4(%rsi),%xmm2
  subss  -0x8(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x8(%rsi),%xmm2
  movss  %xmm2,0x30(%rbx,%r11,1)
  movss  %xmm2,0xc(%rbp,%rdi,4)
  add    $0x4,%rdi
  cmp    %edi,-0x6c(%rsp)
  jg     L32
  mov    (%rsp),%eax
  mov    0xc(%rsp),%ecx
  mov    -0x40(%rsp),%r14
  lea    -0x4(%rax),%edx
  sub    %ecx,%edx
  shr    $0x2,%edx
  lea    0x1(%rdx),%eax
  shl    $0x5,%rax
  add    %rax,-0x60(%rsp)
  add    %rax,-0x68(%rsp)
  add    %rax,%r14
  add    %rax,-0x78(%rsp)
  add    %rax,-0x58(%rsp)
  lea    0x4(%rcx,%rdx,4),%eax
  mov    %eax,0xc(%rsp)
L4:
  mov    (%rsp),%esi
  mov    0xc(%rsp),%edi
  cmp    %edi,%esi
  jle    L28
  mov    0x8(%rsp),%edx
  lea    0x0(,%rdi,4),%eax
  mov    -0x50(%rsp),%rcx
  sub    %edi,%esi
  cltq
  sub    $0x1,%esi
  mov    -0x78(%rsp),%r8
  mov    -0x68(%rsp),%r9
  imul   %edi,%edx
  movss  D21(%rip),%xmm0
  mov    -0x48(%rsp),%r10
  movslq %edx,%rdx
  add    %rax,%rdx
  lea    (%rcx,%rdx,4),%rdx
  movslq %edi,%rcx
  mov    -0x60(%rsp),%rdi
  lea    0x1(%rcx,%rsi,1),%rsi
  neg    %rcx
  lea    (%rdi,%rcx,8),%rcx
  mov    -0x58(%rsp),%rdi
  shl    $0x2,%rsi
  nopl   0x0(%rax,%rax,1)
L71:
  movaps %xmm0,%xmm2
  movaps %xmm0,%xmm3
  divss  (%rcx,%rax,2),%xmm2
  movss  %xmm2,(%rdi)
  add    $0x8,%r9
  add    $0x8,%r14
  add    $0x8,%r8
  add    $0x8,%rdi
  divss  0x4(%rcx,%rax,2),%xmm3
  movss  %xmm3,-0x4(%rdi)
  add    $0x14,%rdx
  mulss  -0x8(%r9),%xmm2
  movss  %xmm2,-0x8(%r14)
  movss  -0x4(%r9),%xmm2
  mulss  -0x4(%rdi),%xmm2
  movss  %xmm2,-0x4(%r14)
  movss  -0x8(%rdi),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x8(%r8)
  movss  -0x4(%rdi),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,-0x4(%r8)
  movss  (%r12,%rax,1),%xmm2
  addss  -0x8(%r14),%xmm2
  addss  -0x4(%r14),%xmm2
  movss  %xmm2,-0x14(%rdx)
  movss  -0x4(%r14),%xmm3
  mulss  0x4(%r10,%rax,2),%xmm3
  movss  0x0(%r13,%rax,1),%xmm2
  addss  -0x4(%r9),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x8(%r14),%xmm3
  mulss  (%r10,%rax,2),%xmm3
  addss  -0x4(%r8),%xmm2
  subss  -0x8(%r9),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x8(%r8),%xmm2
  movss  %xmm2,(%rbx,%rax,4)
  movss  %xmm2,0x0(%rbp,%rax,1)
  add    $0x4,%rax
  cmp    %rax,%rsi
  jne    L71
L28:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L161:
  mov    0x8(%rsp),%edx
  lea    0x0(,%rcx,4),%eax
  mov    (%rsp),%edi
  movslq %ecx,%rsi
  mov    -0x48(%rsp),%rbx
  cltq
  mov    -0x40(%rsp),%r8
  imul   %ecx,%edx
  sub    %ecx,%edi
  mov    -0x38(%rsp),%r9
  mov    -0x30(%rsp),%r10
  lea    -0x1(%rdi),%ecx
  mov    -0x28(%rsp),%r11
  mov    -0x20(%rsp),%rdi
  lea    0x1(%rsi,%rcx,1),%rcx
  neg    %rsi
  movslq %edx,%rdx
  shl    $0x2,%rcx
  add    %rax,%rdx
  lea    (%rbx,%rdx,4),%rdx
  mov    -0x68(%rsp),%rbx
  lea    (%rbx,%rsi,8),%rsi
  nopl   0x0(%rax,%rax,1)
L162:
  movaps %xmm2,%xmm0
  movaps %xmm2,%xmm3
  divss  (%rsi,%rax,2),%xmm0
  add    $0x8,%r14
  add    $0x8,%r12
  add    $0x8,%r15
  add    $0x8,%r13
  add    $0x14,%rdx
  movss  %xmm0,-0x8(%r13)
  divss  0x4(%rsi,%rax,2),%xmm3
  movss  %xmm3,-0x4(%r13)
  mulss  -0x8(%r14),%xmm0
  movss  %xmm0,-0x8(%r12)
  movss  -0x4(%r14),%xmm0
  mulss  -0x4(%r13),%xmm0
  movss  %xmm0,-0x4(%r12)
  movss  -0x8(%r13),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x8(%r15)
  movss  -0x4(%r13),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%r15)
  movss  (%r10,%rax,1),%xmm0
  addss  -0x8(%r12),%xmm0
  addss  -0x4(%r12),%xmm0
  movss  %xmm0,-0x14(%rdx)
  movss  -0x4(%r12),%xmm3
  mulss  0x4(%rdi,%rax,2),%xmm3
  movss  (%r11,%rax,1),%xmm0
  addss  -0x4(%r14),%xmm0
  addss  %xmm3,%xmm0
  movss  -0x8(%r12),%xmm3
  mulss  (%rdi,%rax,2),%xmm3
  addss  -0x4(%r15),%xmm0
  subss  -0x8(%r14),%xmm0
  subss  %xmm3,%xmm0
  subss  -0x8(%r15),%xmm0
  movss  %xmm0,(%r8,%rax,4)
  movss  %xmm0,(%r9,%rax,1)
  add    $0x4,%rax
  cmp    %rax,%rcx
  jne    L162
  jmp    L127
  nopw   0x0(%rax,%rax,1)
L159:
  mov    %rdi,-0x68(%rsp)
  xor    %ecx,%ecx
  jmp    L165
L46:
  mov    %rdi,-0x68(%rsp)
  xor    %ecx,%ecx
  jmp    L96
```

tool_calls: name: parse_data
arguments: label: D21
type_guessed: f32

tool: D21: 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;
    }
  }
}
```