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

```asm
<s_update_hessian_box_mpc>:
  endbr64
  push   %r15
  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,-0x70(%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,-0x38(%rsp)
  mov    (%rax),%r9
  mov    0x68(%rsp),%rax
  mov    %r10,-0x30(%rsp)
  mov    %r9,-0x28(%rsp)
  mov    (%rax),%rbp
  mov    0x70(%rsp),%rax
  mov    %rbp,-0x20(%rsp)
  mov    (%rax),%r11
  mov    0x90(%rsp),%rax
  mov    (%rax),%rcx
  mov    %r11,-0x10(%rsp)
  mov    %rcx,-0x8(%rsp)
  cmp    $0x7,%esi
  jle    L151
  lea    0x40(,%r8,8),%eax
  mov    %r14,-0x48(%rsp)
  mov    %rcx,%r8
  movss  D4(%rip),%xmm1
  cltq
  mov    %r12,-0x40(%rsp)
  mov    %r14,%rcx
  mov    %r13,%rdx
  shl    $0x2,%rax
  mov    %rax,-0x58(%rsp)
  lea    -0x8(%rsi),%eax
  mov    -0x58(%rsp),%r14
  mov    %r15,%rsi
  shr    $0x3,%eax
  add    $0x1,%eax
  mov    %eax,-0x50(%rsp)
  shl    $0x6,%rax
  mov    %rax,-0x60(%rsp)
  add    %rdi,%rax
  mov    %rax,-0x68(%rsp)
  mov    %r12,%rax
  mov    -0x68(%rsp),%r12
  nopw   0x0(%rax,%rax,1)
L70:
  movaps %xmm1,%xmm2
  movaps %xmm1,%xmm3
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  add    $0x40,%rdi
  add    $0x40,%rcx
  add    $0x40,%rax
  add    $0x40,%rsi
  add    $0x40,%rdx
  add    $0x20,%rbp
  add    $0x20,%r11
  add    $0x40,%r8
  divss  -0x3c(%rdi),%xmm3
  movss  %xmm3,-0x3c(%rdx)
  mulss  -0x40(%rcx),%xmm2
  add    $0x100,%r10
  add    $0x20,%r9
  movss  %xmm2,-0x40(%rax)
  movss  -0x3c(%rcx),%xmm2
  mulss  -0x3c(%rdx),%xmm2
  movss  %xmm2,-0x3c(%rax)
  movss  -0x40(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x40(%rsi)
  movss  -0x3c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x3c(%rsi)
  movss  -0x20(%rbp),%xmm2
  addss  -0x40(%rax),%xmm2
  addss  -0x3c(%rax),%xmm2
  movss  %xmm2,(%rbx)
  movss  -0x3c(%rax),%xmm3
  mulss  -0x3c(%r8),%xmm3
  movss  -0x20(%r11),%xmm2
  addss  -0x3c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x40(%rax),%xmm3
  mulss  -0x40(%r8),%xmm3
  addss  -0x3c(%rsi),%xmm2
  subss  -0x40(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x40(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0x100(%r10)
  movss  %xmm2,-0x20(%r9)
  movaps %xmm1,%xmm2
  divss  -0x38(%rdi),%xmm2
  movss  %xmm2,-0x38(%rdx)
  divss  -0x34(%rdi),%xmm3
  movss  %xmm3,-0x34(%rdx)
  mulss  -0x38(%rcx),%xmm2
  movss  %xmm2,-0x38(%rax)
  movss  -0x34(%rcx),%xmm2
  mulss  -0x34(%rdx),%xmm2
  movss  %xmm2,-0x34(%rax)
  movss  -0x38(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x38(%rsi)
  movss  -0x34(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x34(%rsi)
  movss  -0x1c(%rbp),%xmm2
  addss  -0x38(%rax),%xmm2
  addss  -0x34(%rax),%xmm2
  movss  %xmm2,0x24(%rbx)
  movss  -0x34(%rax),%xmm3
  mulss  -0x34(%r8),%xmm3
  movss  -0x1c(%r11),%xmm2
  addss  -0x34(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x38(%rax),%xmm3
  mulss  -0x38(%r8),%xmm3
  addss  -0x34(%rsi),%xmm2
  subss  -0x38(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x38(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0xe0(%r10)
  movss  %xmm2,-0x1c(%r9)
  movaps %xmm1,%xmm2
  divss  -0x30(%rdi),%xmm2
  movss  %xmm2,-0x30(%rdx)
  divss  -0x2c(%rdi),%xmm3
  movss  %xmm3,-0x2c(%rdx)
  mulss  -0x30(%rcx),%xmm2
  movss  %xmm2,-0x30(%rax)
  movss  -0x2c(%rcx),%xmm2
  mulss  -0x2c(%rdx),%xmm2
  movss  %xmm2,-0x2c(%rax)
  movss  -0x30(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x30(%rsi)
  movss  -0x2c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x2c(%rsi)
  movss  -0x18(%rbp),%xmm2
  addss  -0x30(%rax),%xmm2
  addss  -0x2c(%rax),%xmm2
  movss  %xmm2,0x48(%rbx)
  movss  -0x2c(%rax),%xmm3
  mulss  -0x2c(%r8),%xmm3
  movss  -0x18(%r11),%xmm2
  addss  -0x2c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x30(%rax),%xmm3
  mulss  -0x30(%r8),%xmm3
  addss  -0x2c(%rsi),%xmm2
  subss  -0x30(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x30(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0xc0(%r10)
  movss  %xmm2,-0x18(%r9)
  movaps %xmm1,%xmm2
  divss  -0x28(%rdi),%xmm2
  movss  %xmm2,-0x28(%rdx)
  divss  -0x24(%rdi),%xmm3
  movss  %xmm3,-0x24(%rdx)
  mulss  -0x28(%rcx),%xmm2
  movss  %xmm2,-0x28(%rax)
  movss  -0x24(%rcx),%xmm2
  mulss  -0x24(%rdx),%xmm2
  movss  %xmm2,-0x24(%rax)
  movss  -0x28(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x28(%rsi)
  movss  -0x24(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x24(%rsi)
  movss  -0x14(%rbp),%xmm2
  addss  -0x28(%rax),%xmm2
  addss  -0x24(%rax),%xmm2
  movss  %xmm2,0x6c(%rbx)
  movss  -0x24(%rax),%xmm3
  mulss  -0x24(%r8),%xmm3
  movss  -0x14(%r11),%xmm2
  addss  -0x24(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x28(%rax),%xmm3
  mulss  -0x28(%r8),%xmm3
  addss  -0x24(%rsi),%xmm2
  subss  -0x28(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x28(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0xa0(%r10)
  movss  %xmm2,-0x14(%r9)
  movaps %xmm1,%xmm2
  divss  -0x20(%rdi),%xmm2
  movss  %xmm2,-0x20(%rdx)
  divss  -0x1c(%rdi),%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  %xmm0,%xmm2
  movss  %xmm2,-0x20(%rsi)
  movss  -0x1c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x1c(%rsi)
  movss  -0x10(%rbp),%xmm2
  addss  -0x20(%rax),%xmm2
  addss  -0x1c(%rax),%xmm2
  movss  %xmm2,0x90(%rbx)
  movss  -0x10(%r11),%xmm2
  movss  -0x1c(%rax),%xmm3
  mulss  -0x1c(%r8),%xmm3
  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 %xmm1,%xmm3
  movss  %xmm2,-0x80(%r10)
  movss  %xmm2,-0x10(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x18(%rsi)
  movss  -0x14(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x14(%rsi)
  movss  -0xc(%rbp),%xmm2
  addss  -0x18(%rax),%xmm2
  addss  -0x14(%rax),%xmm2
  movss  %xmm2,0xb4(%rbx)
  movss  -0x14(%rax),%xmm3
  mulss  -0x14(%r8),%xmm3
  movss  -0xc(%r11),%xmm2
  addss  -0x14(%rcx),%xmm2
  addss  %xmm3,%xmm2
  addss  -0x14(%rsi),%xmm2
  movss  -0x18(%rax),%xmm3
  mulss  -0x18(%r8),%xmm3
  subss  -0x18(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x18(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0x60(%r10)
  movss  %xmm2,-0xc(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x10(%rsi)
  movss  -0xc(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0xc(%rsi)
  movss  -0x8(%rbp),%xmm2
  addss  -0x10(%rax),%xmm2
  addss  -0xc(%rax),%xmm2
  movss  %xmm2,0xd8(%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 %xmm1,%xmm3
  movss  %xmm2,-0x40(%r10)
  movss  %xmm2,-0x8(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x8(%rsi)
  movss  -0x4(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x4(%rsi)
  movss  -0x4(%rbp),%xmm2
  addss  -0x8(%rax),%xmm2
  addss  -0x4(%rax),%xmm2
  movss  %xmm2,0xfc(%rbx)
  movss  -0x4(%rax),%xmm3
  add    %r14,%rbx
  movss  -0x4(%r11),%xmm2
  mulss  -0x4(%r8),%xmm3
  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,-0x20(%r10)
  movss  %xmm2,-0x4(%r9)
  cmp    %r12,%rdi
  jne    L70
  mov    -0x60(%rsp),%rax
  mov    -0x48(%rsp),%r14
  mov    -0x40(%rsp),%r12
  mov    -0x50(%rsp),%ecx
  add    %rax,%r14
  add    %rax,%r15
  add    %rax,%r13
  add    %rax,%r12
  shl    $0x3,%ecx
L92:
  cmp    %ecx,-0x70(%rsp)
  jle    L10
  mov    0x8(%rsp),%eax
  lea    0x0(,%rcx,8),%edx
  mov    -0x70(%rsp),%esi
  movslq %ecx,%rdi
  mov    -0x38(%rsp),%rbx
  movslq %edx,%rdx
  mov    -0x30(%rsp),%r8
  imul   %ecx,%eax
  sub    %ecx,%esi
  mov    -0x28(%rsp),%r9
  mov    -0x20(%rsp),%r10
  lea    -0x1(%rsi),%ecx
  movss  D4(%rip),%xmm1
  mov    -0x10(%rsp),%r11
  lea    0x1(%rdi,%rcx,1),%rcx
  cltq
  shl    $0x2,%rcx
  add    %rdx,%rax
  lea    (%rbx,%rax,4),%rdx
  mov    -0x68(%rsp),%rbx
  lea    0x0(,%rdi,4),%rax
  neg    %rdi
  lea    (%rbx,%rdi,8),%rsi
  mov    -0x8(%rsp),%rdi
  nopl   0x0(%rax)
L170:
  movaps %xmm1,%xmm2
  movaps %xmm1,%xmm3
  divss  (%rsi,%rax,2),%xmm2
  add    $0x8,%r14
  add    $0x8,%r12
  add    $0x8,%r15
  add    $0x8,%r13
  add    $0x24,%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  %xmm0,%xmm2
  movss  %xmm2,-0x8(%r15)
  movss  -0x4(%r13),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x4(%r15)
  movss  (%r10,%rax,1),%xmm2
  addss  -0x8(%r12),%xmm2
  addss  -0x4(%r12),%xmm2
  movss  %xmm2,-0x24(%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,8)
  movss  %xmm2,(%r9,%rax,1)
  add    $0x4,%rax
  cmp    %rax,%rcx
  jne    L170
L10:
  mov    (%rsp),%eax
  sub    $0x7,%eax
  cmpl   $0x1,-0x78(%rsp)
  mov    %eax,-0x20(%rsp)
  movslq -0x78(%rsp),%rax
  mov    %rax,-0x10(%rsp)
  jle    L48
  mov    0x8(%rsp),%eax
  movss  D4(%rip),%xmm1
  movq   $0x1,-0x78(%rsp)
  lea    0x40(,%rax,8),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x8(%rsp)
  mov    (%rsp),%eax
  sub    $0x8,%eax
  shr    $0x3,%eax
  lea    0x1(%rax),%edx
  mov    %rdx,%rax
  shl    $0x6,%rdx
  shl    $0x3,%eax
  mov    %rdx,-0x68(%rsp)
  mov    %eax,0x4(%rsp)
  nopw   0x0(%rax,%rax,1)
L68:
  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,-0x48(%rsp)
  mov    (%rbx,%rax,8),%r12
  mov    0x60(%rsp),%rbx
  mov    0x70(%rsp),%rcx
  mov    0x90(%rsp),%rsi
  mov    %r9,-0x40(%rsp)
  mov    (%rbx,%rax,8),%r15
  mov    0x48(%rsp),%rbx
  mov    %rbp,-0x38(%rsp)
  mov    (%rcx,%rax,8),%r11
  mov    (%rsi,%rax,8),%r8
  mov    (%rbx,%rax,8),%r13
  mov    0x78(%rsp),%rbx
  mov    %r11,-0x30(%rsp)
  mov    (%rbx,%rax,8),%rbx
  mov    -0x20(%rsp),%eax
  mov    %r8,-0x28(%rsp)
  mov    %rbx,-0x50(%rsp)
  test   %eax,%eax
  jle    L66
  mov    -0x68(%rsp),%rax
  mov    %r14,-0x60(%rsp)
  mov    %r14,%rcx
  mov    %r13,%rdx
  mov    %r12,-0x58(%rsp)
  mov    -0x8(%rsp),%r14
  mov    %r15,%rsi
  add    %rdi,%rax
  mov    %rax,-0x70(%rsp)
  mov    %r12,%rax
  mov    -0x70(%rsp),%r12
  nopl   0x0(%rax,%rax,1)
L45:
  movaps %xmm1,%xmm2
  movaps %xmm1,%xmm3
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  add    $0x40,%rdi
  add    $0x40,%rcx
  add    $0x40,%rax
  add    $0x40,%rsi
  add    $0x40,%rdx
  add    $0x20,%rbp
  add    $0x20,%r11
  add    $0x40,%r8
  divss  -0x3c(%rdi),%xmm3
  movss  %xmm3,-0x3c(%rdx)
  mulss  -0x40(%rcx),%xmm2
  add    $0x100,%r10
  add    $0x20,%r9
  movss  %xmm2,-0x40(%rax)
  movss  -0x3c(%rcx),%xmm2
  mulss  -0x3c(%rdx),%xmm2
  movss  %xmm2,-0x3c(%rax)
  movss  -0x40(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x40(%rsi)
  movss  -0x3c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x3c(%rsi)
  movss  -0x20(%rbp),%xmm2
  addss  -0x40(%rax),%xmm2
  addss  -0x3c(%rax),%xmm2
  movss  %xmm2,(%rbx)
  movss  -0x3c(%rax),%xmm3
  mulss  -0x3c(%r8),%xmm3
  movss  -0x20(%r11),%xmm2
  addss  -0x3c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x40(%rax),%xmm3
  mulss  -0x40(%r8),%xmm3
  addss  -0x3c(%rsi),%xmm2
  subss  -0x40(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x40(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0x100(%r10)
  movss  %xmm2,-0x20(%r9)
  movaps %xmm1,%xmm2
  divss  -0x38(%rdi),%xmm2
  movss  %xmm2,-0x38(%rdx)
  divss  -0x34(%rdi),%xmm3
  movss  %xmm3,-0x34(%rdx)
  mulss  -0x38(%rcx),%xmm2
  movss  %xmm2,-0x38(%rax)
  movss  -0x34(%rcx),%xmm2
  mulss  -0x34(%rdx),%xmm2
  movss  %xmm2,-0x34(%rax)
  movss  -0x38(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x38(%rsi)
  movss  -0x34(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x34(%rsi)
  movss  -0x1c(%rbp),%xmm2
  addss  -0x38(%rax),%xmm2
  addss  -0x34(%rax),%xmm2
  movss  %xmm2,0x24(%rbx)
  movss  -0x34(%rax),%xmm3
  mulss  -0x34(%r8),%xmm3
  movss  -0x1c(%r11),%xmm2
  addss  -0x34(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x38(%rax),%xmm3
  mulss  -0x38(%r8),%xmm3
  addss  -0x34(%rsi),%xmm2
  subss  -0x38(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x38(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0xe0(%r10)
  movss  %xmm2,-0x1c(%r9)
  movaps %xmm1,%xmm2
  divss  -0x30(%rdi),%xmm2
  movss  %xmm2,-0x30(%rdx)
  divss  -0x2c(%rdi),%xmm3
  movss  %xmm3,-0x2c(%rdx)
  mulss  -0x30(%rcx),%xmm2
  movss  %xmm2,-0x30(%rax)
  movss  -0x2c(%rcx),%xmm2
  mulss  -0x2c(%rdx),%xmm2
  movss  %xmm2,-0x2c(%rax)
  movss  -0x30(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x30(%rsi)
  movss  -0x2c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x2c(%rsi)
  movss  -0x18(%rbp),%xmm2
  addss  -0x30(%rax),%xmm2
  addss  -0x2c(%rax),%xmm2
  movss  %xmm2,0x48(%rbx)
  movss  -0x2c(%rax),%xmm3
  mulss  -0x2c(%r8),%xmm3
  movss  -0x18(%r11),%xmm2
  addss  -0x2c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x30(%rax),%xmm3
  mulss  -0x30(%r8),%xmm3
  addss  -0x2c(%rsi),%xmm2
  subss  -0x30(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x30(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0xc0(%r10)
  movss  %xmm2,-0x18(%r9)
  movaps %xmm1,%xmm2
  divss  -0x28(%rdi),%xmm2
  movss  %xmm2,-0x28(%rdx)
  divss  -0x24(%rdi),%xmm3
  movss  %xmm3,-0x24(%rdx)
  mulss  -0x28(%rcx),%xmm2
  movss  %xmm2,-0x28(%rax)
  movss  -0x24(%rcx),%xmm2
  mulss  -0x24(%rdx),%xmm2
  movss  %xmm2,-0x24(%rax)
  movss  -0x28(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x28(%rsi)
  movss  -0x24(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x24(%rsi)
  movss  -0x14(%rbp),%xmm2
  addss  -0x28(%rax),%xmm2
  addss  -0x24(%rax),%xmm2
  movss  %xmm2,0x6c(%rbx)
  movss  -0x24(%rax),%xmm3
  mulss  -0x24(%r8),%xmm3
  movss  -0x14(%r11),%xmm2
  addss  -0x24(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x28(%rax),%xmm3
  mulss  -0x28(%r8),%xmm3
  addss  -0x24(%rsi),%xmm2
  subss  -0x28(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x28(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0xa0(%r10)
  movss  %xmm2,-0x14(%r9)
  movaps %xmm1,%xmm2
  divss  -0x20(%rdi),%xmm2
  movss  %xmm2,-0x20(%rdx)
  divss  -0x1c(%rdi),%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  %xmm0,%xmm2
  movss  %xmm2,-0x20(%rsi)
  movss  -0x1c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x1c(%rsi)
  movss  -0x10(%rbp),%xmm2
  addss  -0x20(%rax),%xmm2
  addss  -0x1c(%rax),%xmm2
  movss  %xmm2,0x90(%rbx)
  movss  -0x10(%r11),%xmm2
  movss  -0x1c(%rax),%xmm3
  mulss  -0x1c(%r8),%xmm3
  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 %xmm1,%xmm3
  movss  %xmm2,-0x80(%r10)
  movss  %xmm2,-0x10(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x18(%rsi)
  movss  -0x14(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x14(%rsi)
  movss  -0xc(%rbp),%xmm2
  addss  -0x18(%rax),%xmm2
  addss  -0x14(%rax),%xmm2
  movss  %xmm2,0xb4(%rbx)
  movss  -0x14(%rax),%xmm3
  mulss  -0x14(%r8),%xmm3
  movss  -0xc(%r11),%xmm2
  addss  -0x14(%rcx),%xmm2
  addss  %xmm3,%xmm2
  addss  -0x14(%rsi),%xmm2
  movss  -0x18(%rax),%xmm3
  mulss  -0x18(%r8),%xmm3
  subss  -0x18(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x18(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,-0x60(%r10)
  movss  %xmm2,-0xc(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x10(%rsi)
  movss  -0xc(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0xc(%rsi)
  movss  -0x8(%rbp),%xmm2
  addss  -0x10(%rax),%xmm2
  addss  -0xc(%rax),%xmm2
  movss  %xmm2,0xd8(%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 %xmm1,%xmm3
  movss  %xmm2,-0x40(%r10)
  movss  %xmm2,-0x8(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x8(%rsi)
  movss  -0x4(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x4(%rsi)
  movss  -0x4(%rbp),%xmm2
  addss  -0x8(%rax),%xmm2
  addss  -0x4(%rax),%xmm2
  movss  %xmm2,0xfc(%rbx)
  movss  -0x4(%rax),%xmm3
  add    %r14,%rbx
  movss  -0x4(%r11),%xmm2
  mulss  -0x4(%r8),%xmm3
  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,-0x20(%r10)
  movss  %xmm2,-0x4(%r9)
  cmp    %r12,%rdi
  jne    L45
  mov    -0x68(%rsp),%rax
  mov    -0x60(%rsp),%r14
  mov    -0x58(%rsp),%r12
  mov    0x4(%rsp),%ecx
  add    %rax,%r14
  add    %rax,%r15
  add    %rax,%r13
  add    %rax,%r12
L38:
  cmp    %ecx,(%rsp)
  jg     L120
L9:
  addq   $0x1,-0x78(%rsp)
  mov    -0x78(%rsp),%rax
  cmp    %rax,-0x10(%rsp)
  jne    L68
L48:
  mov    0xc(%rsp),%ecx
  mov    -0x10(%rsp),%rbx
  mov    -0x18(%rsp),%rsi
  mov    0x50(%rsp),%rdx
  lea    (%rcx,%rcx,1),%eax
  mov    (%rsi,%rbx,8),%rdi
  cltq
  mov    (%rdx,%rbx,8),%rsi
  shl    $0x2,%rax
  mov    0x60(%rsp),%rdx
  add    %rax,%rsi
  add    %rax,%rdi
  mov    %rsi,-0x70(%rsp)
  mov    0x58(%rsp),%rsi
  mov    (%rsi,%rbx,8),%r14
  mov    (%rdx,%rbx,8),%rsi
  add    %rax,%rsi
  add    %rax,%r14
  mov    %rsi,-0x78(%rsp)
  mov    0x48(%rsp),%rsi
  mov    (%rsi,%rbx,8),%r15
  mov    0x78(%rsp),%rsi
  mov    (%rsi,%rbx,8),%rsi
  add    %rax,%r15
  mov    %rsi,-0x60(%rsp)
  mov    0x80(%rsp),%rsi
  mov    (%rsi,%rbx,8),%rsi
  mov    %rsi,-0x58(%rsp)
  mov    0x88(%rsp),%rsi
  mov    (%rsi,%rbx,8),%rsi
  mov    %rsi,-0x50(%rsp)
  mov    0x68(%rsp),%rsi
  mov    (%rsi,%rbx,8),%rsi
  mov    %rsi,-0x48(%rsp)
  mov    0x70(%rsp),%rsi
  mov    (%rsi,%rbx,8),%rsi
  mov    %rsi,-0x40(%rsp)
  mov    0x90(%rsp),%rsi
  mov    (%rsi,%rbx,8),%rbx
  mov    %rbx,-0x38(%rsp)
  cmp    -0x20(%rsp),%ecx
  jge    L150
  mov    %rbx,%rsi
  mov    (%rsp),%r10d
  movslq %ecx,%r12
  mov    -0x48(%rsp),%rbx
  lea    (%rsi,%rax,1),%r8
  mov    0x8(%rsp),%esi
  mov    %r14,-0x18(%rsp)
  xor    %r11d,%r11d
  lea    0x0(,%r12,4),%r9
  mov    %r15,-0x10(%rsp)
  shl    $0x5,%r12
  movss  D4(%rip),%xmm1
  lea    0x40(,%rsi,8),%edx
  lea    (%rbx,%r9,1),%rbp
  mov    -0x40(%rsp),%rbx
  add    -0x58(%rsp),%r12
  movslq %edx,%rdx
  imul   %ecx,%esi
  shl    $0x2,%rdx
  add    %r9,%rbx
  add    -0x50(%rsp),%r9
  mov    %rdx,-0x20(%rsp)
  lea    -0x8(%r10),%edx
  xor    %r10d,%r10d
  sub    %ecx,%edx
  mov    -0x70(%rsp),%rcx
  shr    $0x3,%edx
  mov    %edx,-0x28(%rsp)
  add    $0x1,%edx
  shl    $0x6,%rdx
  mov    %rdx,-0x30(%rsp)
  add    %rdi,%rdx
  mov    %rdx,-0x68(%rsp)
  movslq %esi,%rdx
  mov    -0x78(%rsp),%rsi
  add    %rax,%rdx
  mov    -0x60(%rsp),%rax
  lea    (%rax,%rdx,4),%r13
  mov    %r15,%rdx
  mov    %r14,%rax
  mov    -0x68(%rsp),%r15
  mov    -0x20(%rsp),%r14
  nopw   0x0(%rax,%rax,1)
L136:
  movaps %xmm1,%xmm2
  movaps %xmm1,%xmm3
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  add    $0x40,%rdi
  add    $0x40,%rcx
  add    $0x40,%rax
  add    $0x40,%rsi
  add    $0x40,%rdx
  add    $0x20,%rbp
  add    $0x20,%rbx
  add    $0x40,%r8
  divss  -0x3c(%rdi),%xmm3
  movss  %xmm3,-0x3c(%rdx)
  mulss  -0x40(%rcx),%xmm2
  add    $0x20,%r9
  movss  %xmm2,-0x40(%rax)
  movss  -0x3c(%rcx),%xmm2
  mulss  -0x3c(%rdx),%xmm2
  movss  %xmm2,-0x3c(%rax)
  movss  -0x40(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x40(%rsi)
  movss  -0x3c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x3c(%rsi)
  movss  -0x20(%rbp),%xmm2
  addss  -0x40(%rax),%xmm2
  addss  -0x3c(%rax),%xmm2
  movss  %xmm2,0x0(%r13,%r10,1)
  movss  -0x3c(%rax),%xmm3
  mulss  -0x3c(%r8),%xmm3
  movss  -0x20(%rbx),%xmm2
  addss  -0x3c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x40(%rax),%xmm3
  mulss  -0x40(%r8),%xmm3
  addss  -0x3c(%rsi),%xmm2
  subss  -0x40(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x40(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,(%r12,%r11,1)
  movss  %xmm2,-0x20(%r9)
  movaps %xmm1,%xmm2
  divss  -0x38(%rdi),%xmm2
  movss  %xmm2,-0x38(%rdx)
  divss  -0x34(%rdi),%xmm3
  movss  %xmm3,-0x34(%rdx)
  mulss  -0x38(%rcx),%xmm2
  movss  %xmm2,-0x38(%rax)
  movss  -0x34(%rcx),%xmm2
  mulss  -0x34(%rdx),%xmm2
  movss  %xmm2,-0x34(%rax)
  movss  -0x38(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x38(%rsi)
  movss  -0x34(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x34(%rsi)
  movss  -0x1c(%rbp),%xmm2
  addss  -0x38(%rax),%xmm2
  addss  -0x34(%rax),%xmm2
  movss  %xmm2,0x24(%r13,%r10,1)
  movss  -0x34(%rax),%xmm3
  mulss  -0x34(%r8),%xmm3
  movss  -0x1c(%rbx),%xmm2
  addss  -0x34(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x38(%rax),%xmm3
  mulss  -0x38(%r8),%xmm3
  addss  -0x34(%rsi),%xmm2
  subss  -0x38(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x38(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,0x20(%r12,%r11,1)
  movss  %xmm2,-0x1c(%r9)
  movaps %xmm1,%xmm2
  divss  -0x30(%rdi),%xmm2
  movss  %xmm2,-0x30(%rdx)
  divss  -0x2c(%rdi),%xmm3
  movss  %xmm3,-0x2c(%rdx)
  mulss  -0x30(%rcx),%xmm2
  movss  %xmm2,-0x30(%rax)
  movss  -0x2c(%rcx),%xmm2
  mulss  -0x2c(%rdx),%xmm2
  movss  %xmm2,-0x2c(%rax)
  movss  -0x30(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x30(%rsi)
  movss  -0x2c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x2c(%rsi)
  movss  -0x18(%rbp),%xmm2
  addss  -0x30(%rax),%xmm2
  addss  -0x2c(%rax),%xmm2
  movss  %xmm2,0x48(%r13,%r10,1)
  movss  -0x2c(%rax),%xmm3
  mulss  -0x2c(%r8),%xmm3
  movss  -0x18(%rbx),%xmm2
  addss  -0x2c(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x30(%rax),%xmm3
  mulss  -0x30(%r8),%xmm3
  addss  -0x2c(%rsi),%xmm2
  subss  -0x30(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x30(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,0x40(%r12,%r11,1)
  movss  %xmm2,-0x18(%r9)
  movaps %xmm1,%xmm2
  divss  -0x28(%rdi),%xmm2
  movss  %xmm2,-0x28(%rdx)
  divss  -0x24(%rdi),%xmm3
  movss  %xmm3,-0x24(%rdx)
  mulss  -0x28(%rcx),%xmm2
  movss  %xmm2,-0x28(%rax)
  movss  -0x24(%rcx),%xmm2
  mulss  -0x24(%rdx),%xmm2
  movss  %xmm2,-0x24(%rax)
  movss  -0x28(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x28(%rsi)
  movss  -0x24(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x24(%rsi)
  movss  -0x14(%rbp),%xmm2
  addss  -0x28(%rax),%xmm2
  addss  -0x24(%rax),%xmm2
  movss  %xmm2,0x6c(%r13,%r10,1)
  movss  -0x24(%rax),%xmm3
  mulss  -0x24(%r8),%xmm3
  movss  -0x14(%rbx),%xmm2
  addss  -0x24(%rcx),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x28(%rax),%xmm3
  mulss  -0x28(%r8),%xmm3
  addss  -0x24(%rsi),%xmm2
  subss  -0x28(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x28(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,0x60(%r12,%r11,1)
  movss  %xmm2,-0x14(%r9)
  movaps %xmm1,%xmm2
  divss  -0x20(%rdi),%xmm2
  movss  %xmm2,-0x20(%rdx)
  divss  -0x1c(%rdi),%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  %xmm0,%xmm2
  movss  %xmm2,-0x20(%rsi)
  movss  -0x1c(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x1c(%rsi)
  movss  -0x10(%rbp),%xmm2
  addss  -0x20(%rax),%xmm2
  addss  -0x1c(%rax),%xmm2
  movss  %xmm2,0x90(%r13,%r10,1)
  movss  -0x10(%rbx),%xmm2
  movss  -0x1c(%rax),%xmm3
  mulss  -0x1c(%r8),%xmm3
  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 %xmm1,%xmm3
  movss  %xmm2,0x80(%r12,%r11,1)
  movss  %xmm2,-0x10(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x18(%rsi)
  movss  -0x14(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x14(%rsi)
  movss  -0xc(%rbp),%xmm2
  addss  -0x18(%rax),%xmm2
  addss  -0x14(%rax),%xmm2
  movss  %xmm2,0xb4(%r13,%r10,1)
  movss  -0x14(%rax),%xmm3
  mulss  -0x14(%r8),%xmm3
  movss  -0xc(%rbx),%xmm2
  addss  -0x14(%rcx),%xmm2
  addss  %xmm3,%xmm2
  addss  -0x14(%rsi),%xmm2
  movss  -0x18(%rax),%xmm3
  mulss  -0x18(%r8),%xmm3
  subss  -0x18(%rcx),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x18(%rsi),%xmm2
  movaps %xmm1,%xmm3
  movss  %xmm2,0xa0(%r12,%r11,1)
  movss  %xmm2,-0xc(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x10(%rsi)
  movss  -0xc(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0xc(%rsi)
  movss  -0x8(%rbp),%xmm2
  addss  -0x10(%rax),%xmm2
  addss  -0xc(%rax),%xmm2
  movss  %xmm2,0xd8(%r13,%r10,1)
  movss  -0xc(%rax),%xmm3
  mulss  -0xc(%r8),%xmm3
  movss  -0x8(%rbx),%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 %xmm1,%xmm3
  movss  %xmm2,0xc0(%r12,%r11,1)
  movss  %xmm2,-0x8(%r9)
  movaps %xmm1,%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  %xmm0,%xmm2
  movss  %xmm2,-0x8(%rsi)
  movss  -0x4(%rdx),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x4(%rsi)
  movss  -0x4(%rbp),%xmm2
  addss  -0x8(%rax),%xmm2
  addss  -0x4(%rax),%xmm2
  movss  %xmm2,0xfc(%r13,%r10,1)
  movss  -0x4(%rax),%xmm3
  add    %r14,%r10
  movss  -0x4(%rbx),%xmm2
  mulss  -0x4(%r8),%xmm3
  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,0xe0(%r12,%r11,1)
  add    $0x100,%r11
  movss  %xmm2,-0x4(%r9)
  cmp    %r15,%rdi
  jne    L136
  mov    -0x30(%rsp),%rax
  mov    -0x18(%rsp),%r14
  mov    -0x10(%rsp),%r15
  mov    -0x28(%rsp),%ebx
  add    %rax,-0x70(%rsp)
  add    %rax,%r14
  add    %rax,-0x78(%rsp)
  add    %rax,%r15
  mov    0xc(%rsp),%eax
  lea    0x8(%rax,%rbx,8),%eax
  mov    %eax,0xc(%rsp)
  mov    %eax,%ebx
L137:
  mov    (%rsp),%edi
  cmp    %ebx,%edi
  jle    L44
  mov    0x8(%rsp),%edx
  lea    0x0(,%rbx,8),%eax
  sub    %ebx,%edi
  mov    -0x60(%rsp),%rcx
  cltq
  mov    %edi,%esi
  mov    -0x70(%rsp),%r8
  mov    -0x78(%rsp),%rdi
  imul   %ebx,%edx
  sub    $0x1,%esi
  mov    -0x58(%rsp),%r10
  mov    -0x50(%rsp),%r11
  movss  D4(%rip),%xmm1
  mov    -0x40(%rsp),%rbp
  mov    -0x38(%rsp),%r9
  movslq %edx,%rdx
  add    %rdx,%rax
  lea    (%rcx,%rax,4),%rdx
  movslq %ebx,%rcx
  mov    -0x68(%rsp),%rbx
  lea    0x1(%rcx,%rsi,1),%rsi
  lea    0x0(,%rcx,4),%rax
  neg    %rcx
  lea    (%rbx,%rcx,8),%rcx
  mov    -0x48(%rsp),%rbx
  shl    $0x2,%rsi
  nopl   0x0(%rax)
L81:
  movaps %xmm1,%xmm2
  movaps %xmm1,%xmm3
  divss  (%rcx,%rax,2),%xmm2
  movss  %xmm2,(%r15)
  add    $0x8,%r8
  add    $0x8,%r14
  add    $0x8,%rdi
  add    $0x8,%r15
  divss  0x4(%rcx,%rax,2),%xmm3
  movss  %xmm3,-0x4(%r15)
  add    $0x24,%rdx
  mulss  -0x8(%r8),%xmm2
  movss  %xmm2,-0x8(%r14)
  movss  -0x4(%r8),%xmm2
  mulss  -0x4(%r15),%xmm2
  movss  %xmm2,-0x4(%r14)
  movss  -0x8(%r15),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x8(%rdi)
  movss  -0x4(%r15),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x4(%rdi)
  movss  (%rbx,%rax,1),%xmm2
  addss  -0x8(%r14),%xmm2
  addss  -0x4(%r14),%xmm2
  movss  %xmm2,-0x24(%rdx)
  movss  -0x4(%r14),%xmm3
  mulss  0x4(%r9,%rax,2),%xmm3
  movss  0x0(%rbp,%rax,1),%xmm2
  addss  -0x4(%r8),%xmm2
  addss  %xmm3,%xmm2
  movss  -0x8(%r14),%xmm3
  mulss  (%r9,%rax,2),%xmm3
  addss  -0x4(%rdi),%xmm2
  subss  -0x8(%r8),%xmm2
  subss  %xmm3,%xmm2
  subss  -0x8(%rdi),%xmm2
  movss  %xmm2,(%r10,%rax,8)
  movss  %xmm2,(%r11,%rax,1)
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L81
L44:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L120:
  mov    0x8(%rsp),%eax
  lea    0x0(,%rcx,8),%edx
  mov    (%rsp),%edi
  movslq %ecx,%rsi
  mov    -0x50(%rsp),%rbx
  movslq %edx,%rdx
  mov    -0x48(%rsp),%r8
  imul   %ecx,%eax
  sub    %ecx,%edi
  mov    -0x40(%rsp),%r9
  mov    -0x38(%rsp),%r10
  lea    -0x1(%rdi),%ecx
  mov    -0x30(%rsp),%r11
  mov    -0x28(%rsp),%rdi
  lea    0x1(%rsi,%rcx,1),%rcx
  cltq
  shl    $0x2,%rcx
  add    %rdx,%rax
  lea    (%rbx,%rax,4),%rdx
  mov    -0x70(%rsp),%rbx
  lea    0x0(,%rsi,4),%rax
  neg    %rsi
  lea    (%rbx,%rsi,8),%rsi
L58:
  movaps %xmm1,%xmm2
  movaps %xmm1,%xmm3
  divss  (%rsi,%rax,2),%xmm2
  add    $0x8,%r14
  add    $0x8,%r12
  add    $0x8,%r15
  add    $0x8,%r13
  add    $0x24,%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  %xmm0,%xmm2
  movss  %xmm2,-0x8(%r15)
  movss  -0x4(%r13),%xmm2
  mulss  %xmm0,%xmm2
  movss  %xmm2,-0x4(%r15)
  movss  (%r10,%rax,1),%xmm2
  addss  -0x8(%r12),%xmm2
  addss  -0x4(%r12),%xmm2
  movss  %xmm2,-0x24(%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,8)
  movss  %xmm2,(%r9,%rax,1)
  add    $0x4,%rax
  cmp    %rax,%rcx
  jne    L58
  jmp    L9
  nopw   0x0(%rax,%rax,1)
L66:
  mov    %rdi,-0x70(%rsp)
  xor    %ecx,%ecx
  jmp    L38
L151:
  mov    %rdi,-0x68(%rsp)
  xor    %ecx,%ecx
  jmp    L92
L150:
  mov    %rdi,-0x68(%rsp)
  mov    %ecx,%ebx
  jmp    L137
```

tool_calls: name: extract_symbolic_data
arguments: tag: D4
data_type: f32

tool: D4: 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 = 8;
  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 - 7; ii += 8) {
    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_tinv[8] = 1.0 / ptr_t[8];
    ptr_tinv[9] = 1.0 / ptr_t[9];
    ptr_lamt[8] = ptr_lam[8] * ptr_tinv[8];
    ptr_lamt[9] = ptr_lam[9] * ptr_tinv[9];
    ptr_dlam[8] = ptr_tinv[8] * sigma_mu;
    ptr_dlam[9] = ptr_tinv[9] * sigma_mu;
    ptr_pd[4 + (ii + 4) * bs + ii * cnz] = ptr_bd[ii + 4] + ptr_lamt[8] + ptr_lamt[9];
    ptr_pl[(ii + 4) * bs] = ptr_bl[ii + 4] + ptr_lam[9] + ptr_lamt[9] * ptr_db[2 * ii + 9] + ptr_dlam[9] - ptr_lam[8] - ptr_lamt[8] * ptr_db[2 * ii + 8] - ptr_dlam[8];
    ptr_pl2[ii + 4] = ptr_pl[(ii + 4) * bs];
    ptr_tinv[10] = 1.0 / ptr_t[10];
    ptr_tinv[11] = 1.0 / ptr_t[11];
    ptr_lamt[10] = ptr_lam[10] * ptr_tinv[10];
    ptr_lamt[11] = ptr_lam[11] * ptr_tinv[11];
    ptr_dlam[10] = ptr_tinv[10] * sigma_mu;
    ptr_dlam[11] = ptr_tinv[11] * sigma_mu;
    ptr_pd[5 + (ii + 5) * bs + ii * cnz] = ptr_bd[ii + 5] + ptr_lamt[10] + ptr_lamt[11];
    ptr_pl[(ii + 5) * bs] = ptr_bl[ii + 5] + ptr_lam[11] + ptr_lamt[11] * ptr_db[2 * ii + 11] + ptr_dlam[11] - ptr_lam[10] - ptr_lamt[10] * ptr_db[2 * ii + 10] - ptr_dlam[10];
    ptr_pl2[ii + 5] = ptr_pl[(ii + 5) * bs];
    ptr_tinv[12] = 1.0 / ptr_t[12];
    ptr_tinv[13] = 1.0 / ptr_t[13];
    ptr_lamt[12] = ptr_lam[12] * ptr_tinv[12];
    ptr_lamt[13] = ptr_lam[13] * ptr_tinv[13];
    ptr_dlam[12] = ptr_tinv[12] * sigma_mu;
    ptr_dlam[13] = ptr_tinv[13] * sigma_mu;
    ptr_pd[6 + (ii + 6) * bs + ii * cnz] = ptr_bd[ii + 6] + ptr_lamt[12] + ptr_lamt[13];
    ptr_pl[(ii + 6) * bs] = ptr_bl[ii + 6] + ptr_lam[13] + ptr_lamt[13] * ptr_db[2 * ii + 13] + ptr_dlam[13] - ptr_lam[12] - ptr_lamt[12] * ptr_db[2 * ii + 12] - ptr_dlam[12];
    ptr_pl2[ii + 6] = ptr_pl[(ii + 6) * bs];
    ptr_tinv[14] = 1.0 / ptr_t[14];
    ptr_tinv[15] = 1.0 / ptr_t[15];
    ptr_lamt[14] = ptr_lam[14] * ptr_tinv[14];
    ptr_lamt[15] = ptr_lam[15] * ptr_tinv[15];
    ptr_dlam[14] = ptr_tinv[14] * sigma_mu;
    ptr_dlam[15] = ptr_tinv[15] * sigma_mu;
    ptr_pd[7 + (ii + 7) * bs + ii * cnz] = ptr_bd[ii + 7] + ptr_lamt[14] + ptr_lamt[15];
    ptr_pl[(ii + 7) * bs] = ptr_bl[ii + 7] + ptr_lam[15] + ptr_lamt[15] * ptr_db[2 * ii + 15] + ptr_dlam[15] - ptr_lam[14] - ptr_lamt[14] * ptr_db[2 * ii + 14] - ptr_dlam[14];
    ptr_pl2[ii + 7] = ptr_pl[(ii + 7) * bs];
    ptr_t += 16;
    ptr_lam += 16;
    ptr_lamt += 16;
    ptr_dlam += 16;
    ptr_tinv += 16;
  }
  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 - 7; ii += 8) {
      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_tinv[8] = 1.0 / ptr_t[8];
      ptr_tinv[9] = 1.0 / ptr_t[9];
      ptr_lamt[8] = ptr_lam[8] * ptr_tinv[8];
      ptr_lamt[9] = ptr_lam[9] * ptr_tinv[9];
      ptr_dlam[8] = ptr_tinv[8] * sigma_mu;
      ptr_dlam[9] = ptr_tinv[9] * sigma_mu;
      ptr_pd[4 + (ii + 4) * bs + ii * cnz] = ptr_bd[ii + 4] + ptr_lamt[8] + ptr_lamt[9];
      ptr_pl[(ii + 4) * bs] = ptr_bl[ii + 4] + ptr_lam[9] + ptr_lamt[9] * ptr_db[2 * ii + 9] + ptr_dlam[9] - ptr_lam[8] - ptr_lamt[8] * ptr_db[2 * ii + 8] - ptr_dlam[8];
      ptr_pl2[ii + 4] = ptr_pl[(ii + 4) * bs];
      ptr_tinv[10] = 1.0 / ptr_t[10];
      ptr_tinv[11] = 1.0 / ptr_t[11];
      ptr_lamt[10] = ptr_lam[10] * ptr_tinv[10];
      ptr_lamt[11] = ptr_lam[11] * ptr_tinv[11];
      ptr_dlam[10] = ptr_tinv[10] * sigma_mu;
      ptr_dlam[11] = ptr_tinv[11] * sigma_mu;
      ptr_pd[5 + (ii + 5) * bs + ii * cnz] = ptr_bd[ii + 5] + ptr_lamt[10] + ptr_lamt[11];
      ptr_pl[(ii + 5) * bs] = ptr_bl[ii + 5] + ptr_lam[11] + ptr_lamt[11] * ptr_db[2 * ii + 11] + ptr_dlam[11] - ptr_lam[10] - ptr_lamt[10] * ptr_db[2 * ii + 10] - ptr_dlam[10];
      ptr_pl2[ii + 5] = ptr_pl[(ii + 5) * bs];
      ptr_tinv[12] = 1.0 / ptr_t[12];
      ptr_tinv[13] = 1.0 / ptr_t[13];
      ptr_lamt[12] = ptr_lam[12] * ptr_tinv[12];
      ptr_lamt[13] = ptr_lam[13] * ptr_tinv[13];
      ptr_dlam[12] = ptr_tinv[12] * sigma_mu;
      ptr_dlam[13] = ptr_tinv[13] * sigma_mu;
      ptr_pd[6 + (ii + 6) * bs + ii * cnz] = ptr_bd[ii + 6] + ptr_lamt[12] + ptr_lamt[13];
      ptr_pl[(ii + 6) * bs] = ptr_bl[ii + 6] + ptr_lam[13] + ptr_lamt[13] * ptr_db[2 * ii + 13] + ptr_dlam[13] - ptr_lam[12] - ptr_lamt[12] * ptr_db[2 * ii + 12] - ptr_dlam[12];
      ptr_pl2[ii + 6] = ptr_pl[(ii + 6) * bs];
      ptr_tinv[14] = 1.0 / ptr_t[14];
      ptr_tinv[15] = 1.0 / ptr_t[15];
      ptr_lamt[14] = ptr_lam[14] * ptr_tinv[14];
      ptr_lamt[15] = ptr_lam[15] * ptr_tinv[15];
      ptr_dlam[14] = ptr_tinv[14] * sigma_mu;
      ptr_dlam[15] = ptr_tinv[15] * sigma_mu;
      ptr_pd[7 + (ii + 7) * bs + ii * cnz] = ptr_bd[ii + 7] + ptr_lamt[14] + ptr_lamt[15];
      ptr_pl[(ii + 7) * bs] = ptr_bl[ii + 7] + ptr_lam[15] + ptr_lamt[15] * ptr_db[2 * ii + 15] + ptr_dlam[15] - ptr_lam[14] - ptr_lamt[14] * ptr_db[2 * ii + 14] - ptr_dlam[14];
      ptr_pl2[ii + 7] = ptr_pl[(ii + 7) * bs];
      ptr_t += 16;
      ptr_lam += 16;
      ptr_lamt += 16;
      ptr_dlam += 16;
      ptr_tinv += 16;
    }
    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 - 7; ii += 8) {
    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_tinv[8] = 1.0 / ptr_t[8];
    ptr_tinv[9] = 1.0 / ptr_t[9];
    ptr_lamt[8] = ptr_lam[8] * ptr_tinv[8];
    ptr_lamt[9] = ptr_lam[9] * ptr_tinv[9];
    ptr_dlam[8] = ptr_tinv[8] * sigma_mu;
    ptr_dlam[9] = ptr_tinv[9] * sigma_mu;
    ptr_pd[4 + (ii + 4) * bs + ii * cnz] = ptr_bd[ii + 4] + ptr_lamt[8] + ptr_lamt[9];
    ptr_pl[(ii + 4) * bs] = ptr_bl[ii + 4] + ptr_lam[9] + ptr_lamt[9] * ptr_db[2 * ii + 9] + ptr_dlam[9] - ptr_lam[8] - ptr_lamt[8] * ptr_db[2 * ii + 8] - ptr_dlam[8];
    ptr_pl2[ii + 4] = ptr_pl[(ii + 4) * bs];
    ptr_tinv[10] = 1.0 / ptr_t[10];
    ptr_tinv[11] = 1.0 / ptr_t[11];
    ptr_lamt[10] = ptr_lam[10] * ptr_tinv[10];
    ptr_lamt[11] = ptr_lam[11] * ptr_tinv[11];
    ptr_dlam[10] = ptr_tinv[10] * sigma_mu;
    ptr_dlam[11] = ptr_tinv[11] * sigma_mu;
    ptr_pd[5 + (ii + 5) * bs + ii * cnz] = ptr_bd[ii + 5] + ptr_lamt[10] + ptr_lamt[11];
    ptr_pl[(ii + 5) * bs] = ptr_bl[ii + 5] + ptr_lam[11] + ptr_lamt[11] * ptr_db[2 * ii + 11] + ptr_dlam[11] - ptr_lam[10] - ptr_lamt[10] * ptr_db[2 * ii + 10] - ptr_dlam[10];
    ptr_pl2[ii + 5] = ptr_pl[(ii + 5) * bs];
    ptr_tinv[12] = 1.0 / ptr_t[12];
    ptr_tinv[13] = 1.0 / ptr_t[13];
    ptr_lamt[12] = ptr_lam[12] * ptr_tinv[12];
    ptr_lamt[13] = ptr_lam[13] * ptr_tinv[13];
    ptr_dlam[12] = ptr_tinv[12] * sigma_mu;
    ptr_dlam[13] = ptr_tinv[13] * sigma_mu;
    ptr_pd[6 + (ii + 6) * bs + ii * cnz] = ptr_bd[ii + 6] + ptr_lamt[12] + ptr_lamt[13];
    ptr_pl[(ii + 6) * bs] = ptr_bl[ii + 6] + ptr_lam[13] + ptr_lamt[13] * ptr_db[2 * ii + 13] + ptr_dlam[13] - ptr_lam[12] - ptr_lamt[12] * ptr_db[2 * ii + 12] - ptr_dlam[12];
    ptr_pl2[ii + 6] = ptr_pl[(ii + 6) * bs];
    ptr_tinv[14] = 1.0 / ptr_t[14];
    ptr_tinv[15] = 1.0 / ptr_t[15];
    ptr_lamt[14] = ptr_lam[14] * ptr_tinv[14];
    ptr_lamt[15] = ptr_lam[15] * ptr_tinv[15];
    ptr_dlam[14] = ptr_tinv[14] * sigma_mu;
    ptr_dlam[15] = ptr_tinv[15] * sigma_mu;
    ptr_pd[7 + (ii + 7) * bs + ii * cnz] = ptr_bd[ii + 7] + ptr_lamt[14] + ptr_lamt[15];
    ptr_pl[(ii + 7) * bs] = ptr_bl[ii + 7] + ptr_lam[15] + ptr_lamt[15] * ptr_db[2 * ii + 15] + ptr_dlam[15] - ptr_lam[14] - ptr_lamt[14] * ptr_db[2 * ii + 14] - ptr_dlam[14];
    ptr_pl2[ii + 7] = ptr_pl[(ii + 7) * bs];
    ptr_t += 16;
    ptr_lam += 16;
    ptr_lamt += 16;
    ptr_dlam += 16;
    ptr_tinv += 16;
  }
  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;
    }
  }
}
```