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)
  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    $0x7,%esi
  jle    L60
  lea    0x40(,%r8,8),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x78(%rsp)
  mov    %rcx,%r8
  lea    -0x8(%rsi),%eax
  shr    $0x3,%eax
  mov    %eax,-0x58(%rsp)
  add    $0x1,%eax
  mov    %eax,%eax
  shl    $0x6,%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  D41(%rip),%xmm1
  mov    %r12,-0x48(%rsp)
  mov    -0x70(%rsp),%r12
L126:
  movaps %xmm1,%xmm2
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  movaps %xmm1,%xmm3
  divss  0x4(%rdi),%xmm3
  movss  %xmm3,0x4(%rdx)
  mulss  (%rcx),%xmm2
  movss  %xmm2,(%rax)
  movss  0x4(%rcx),%xmm2
  mulss  0x4(%rdx),%xmm2
  movss  %xmm2,0x4(%rax)
  movaps %xmm0,%xmm2
  mulss  (%rdx),%xmm2
  movss  %xmm2,(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x4(%rdx),%xmm2
  movss  %xmm2,0x4(%rsi)
  movss  0x0(%rbp),%xmm2
  addss  (%rax),%xmm2
  addss  0x4(%rax),%xmm2
  movss  %xmm2,(%rbx)
  movss  (%r11),%xmm2
  addss  0x4(%rcx),%xmm2
  movss  0x4(%rax),%xmm3
  mulss  0x4(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x4(%rsi),%xmm2
  subss  (%rcx),%xmm2
  movss  (%rax),%xmm3
  mulss  (%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  (%rsi),%xmm2
  movss  %xmm2,(%r10)
  movss  %xmm2,(%r9)
  movaps %xmm1,%xmm2
  divss  0x8(%rdi),%xmm2
  movss  %xmm2,0x8(%rdx)
  movaps %xmm1,%xmm3
  divss  0xc(%rdi),%xmm3
  movss  %xmm3,0xc(%rdx)
  mulss  0x8(%rcx),%xmm2
  movss  %xmm2,0x8(%rax)
  movss  0xc(%rcx),%xmm2
  mulss  0xc(%rdx),%xmm2
  movss  %xmm2,0xc(%rax)
  movaps %xmm0,%xmm2
  mulss  0x8(%rdx),%xmm2
  movss  %xmm2,0x8(%rsi)
  movaps %xmm0,%xmm2
  mulss  0xc(%rdx),%xmm2
  movss  %xmm2,0xc(%rsi)
  movss  0x4(%rbp),%xmm2
  addss  0x8(%rax),%xmm2
  addss  0xc(%rax),%xmm2
  movss  %xmm2,0x24(%rbx)
  movss  0x4(%r11),%xmm2
  addss  0xc(%rcx),%xmm2
  movss  0xc(%rax),%xmm3
  mulss  0xc(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0xc(%rsi),%xmm2
  subss  0x8(%rcx),%xmm2
  movss  0x8(%rax),%xmm3
  mulss  0x8(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x8(%rsi),%xmm2
  movss  %xmm2,0x20(%r10)
  movss  %xmm2,0x4(%r9)
  movaps %xmm1,%xmm2
  divss  0x10(%rdi),%xmm2
  movss  %xmm2,0x10(%rdx)
  movaps %xmm1,%xmm3
  divss  0x14(%rdi),%xmm3
  movss  %xmm3,0x14(%rdx)
  mulss  0x10(%rcx),%xmm2
  movss  %xmm2,0x10(%rax)
  movss  0x14(%rcx),%xmm2
  mulss  0x14(%rdx),%xmm2
  movss  %xmm2,0x14(%rax)
  movaps %xmm0,%xmm2
  mulss  0x10(%rdx),%xmm2
  movss  %xmm2,0x10(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x14(%rdx),%xmm2
  movss  %xmm2,0x14(%rsi)
  movss  0x8(%rbp),%xmm2
  addss  0x10(%rax),%xmm2
  addss  0x14(%rax),%xmm2
  movss  %xmm2,0x48(%rbx)
  movss  0x8(%r11),%xmm2
  addss  0x14(%rcx),%xmm2
  movss  0x14(%rax),%xmm3
  mulss  0x14(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x14(%rsi),%xmm2
  subss  0x10(%rcx),%xmm2
  movss  0x10(%rax),%xmm3
  mulss  0x10(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x10(%rsi),%xmm2
  movss  %xmm2,0x40(%r10)
  movss  %xmm2,0x8(%r9)
  movaps %xmm1,%xmm2
  divss  0x18(%rdi),%xmm2
  movss  %xmm2,0x18(%rdx)
  movaps %xmm1,%xmm3
  divss  0x1c(%rdi),%xmm3
  movss  %xmm3,0x1c(%rdx)
  mulss  0x18(%rcx),%xmm2
  movss  %xmm2,0x18(%rax)
  movss  0x1c(%rcx),%xmm2
  mulss  0x1c(%rdx),%xmm2
  movss  %xmm2,0x1c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x18(%rdx),%xmm2
  movss  %xmm2,0x18(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x1c(%rdx),%xmm2
  movss  %xmm2,0x1c(%rsi)
  movss  0xc(%rbp),%xmm2
  addss  0x18(%rax),%xmm2
  addss  0x1c(%rax),%xmm2
  movss  %xmm2,0x6c(%rbx)
  movss  0xc(%r11),%xmm2
  addss  0x1c(%rcx),%xmm2
  movss  0x1c(%rax),%xmm3
  mulss  0x1c(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x1c(%rsi),%xmm2
  subss  0x18(%rcx),%xmm2
  movss  0x18(%rax),%xmm3
  mulss  0x18(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x18(%rsi),%xmm2
  movss  %xmm2,0x60(%r10)
  movss  %xmm2,0xc(%r9)
  movaps %xmm1,%xmm2
  divss  0x20(%rdi),%xmm2
  movss  %xmm2,0x20(%rdx)
  movaps %xmm1,%xmm3
  divss  0x24(%rdi),%xmm3
  movss  %xmm3,0x24(%rdx)
  mulss  0x20(%rcx),%xmm2
  movss  %xmm2,0x20(%rax)
  movss  0x24(%rcx),%xmm2
  mulss  0x24(%rdx),%xmm2
  movss  %xmm2,0x24(%rax)
  movaps %xmm0,%xmm2
  mulss  0x20(%rdx),%xmm2
  movss  %xmm2,0x20(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x24(%rdx),%xmm2
  movss  %xmm2,0x24(%rsi)
  movss  0x10(%rbp),%xmm2
  addss  0x20(%rax),%xmm2
  addss  0x24(%rax),%xmm2
  movss  %xmm2,0x90(%rbx)
  movss  0x10(%r11),%xmm2
  addss  0x24(%rcx),%xmm2
  movss  0x24(%rax),%xmm3
  mulss  0x24(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x24(%rsi),%xmm2
  subss  0x20(%rcx),%xmm2
  movss  0x20(%rax),%xmm3
  mulss  0x20(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x20(%rsi),%xmm2
  movss  %xmm2,0x80(%r10)
  movss  %xmm2,0x10(%r9)
  movaps %xmm1,%xmm2
  divss  0x28(%rdi),%xmm2
  movss  %xmm2,0x28(%rdx)
  movaps %xmm1,%xmm3
  divss  0x2c(%rdi),%xmm3
  movss  %xmm3,0x2c(%rdx)
  mulss  0x28(%rcx),%xmm2
  movss  %xmm2,0x28(%rax)
  movss  0x2c(%rcx),%xmm2
  mulss  0x2c(%rdx),%xmm2
  movss  %xmm2,0x2c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x28(%rdx),%xmm2
  movss  %xmm2,0x28(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x2c(%rdx),%xmm2
  movss  %xmm2,0x2c(%rsi)
  movss  0x14(%rbp),%xmm2
  addss  0x28(%rax),%xmm2
  addss  0x2c(%rax),%xmm2
  movss  %xmm2,0xb4(%rbx)
  movss  0x14(%r11),%xmm2
  addss  0x2c(%rcx),%xmm2
  movss  0x2c(%rax),%xmm3
  mulss  0x2c(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x2c(%rsi),%xmm2
  subss  0x28(%rcx),%xmm2
  movss  0x28(%rax),%xmm3
  mulss  0x28(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x28(%rsi),%xmm2
  movss  %xmm2,0xa0(%r10)
  movss  %xmm2,0x14(%r9)
  movaps %xmm1,%xmm2
  divss  0x30(%rdi),%xmm2
  movss  %xmm2,0x30(%rdx)
  movaps %xmm1,%xmm3
  divss  0x34(%rdi),%xmm3
  movss  %xmm3,0x34(%rdx)
  mulss  0x30(%rcx),%xmm2
  movss  %xmm2,0x30(%rax)
  movss  0x34(%rcx),%xmm2
  mulss  0x34(%rdx),%xmm2
  movss  %xmm2,0x34(%rax)
  movaps %xmm0,%xmm2
  mulss  0x30(%rdx),%xmm2
  movss  %xmm2,0x30(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x34(%rdx),%xmm2
  movss  %xmm2,0x34(%rsi)
  movss  0x18(%rbp),%xmm2
  addss  0x30(%rax),%xmm2
  addss  0x34(%rax),%xmm2
  movss  %xmm2,0xd8(%rbx)
  movss  0x18(%r11),%xmm2
  addss  0x34(%rcx),%xmm2
  movss  0x34(%rax),%xmm3
  mulss  0x34(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x34(%rsi),%xmm2
  subss  0x30(%rcx),%xmm2
  movss  0x30(%rax),%xmm3
  mulss  0x30(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x30(%rsi),%xmm2
  movss  %xmm2,0xc0(%r10)
  movss  %xmm2,0x18(%r9)
  movaps %xmm1,%xmm2
  divss  0x38(%rdi),%xmm2
  movss  %xmm2,0x38(%rdx)
  movaps %xmm1,%xmm3
  divss  0x3c(%rdi),%xmm3
  movss  %xmm3,0x3c(%rdx)
  mulss  0x38(%rcx),%xmm2
  movss  %xmm2,0x38(%rax)
  movss  0x3c(%rcx),%xmm2
  mulss  0x3c(%rdx),%xmm2
  movss  %xmm2,0x3c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x38(%rdx),%xmm2
  movss  %xmm2,0x38(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x3c(%rdx),%xmm2
  movss  %xmm2,0x3c(%rsi)
  movss  0x1c(%rbp),%xmm2
  addss  0x38(%rax),%xmm2
  addss  0x3c(%rax),%xmm2
  movss  %xmm2,0xfc(%rbx)
  movss  0x1c(%r11),%xmm2
  addss  0x3c(%rcx),%xmm2
  movss  0x3c(%rax),%xmm3
  mulss  0x3c(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x3c(%rsi),%xmm2
  subss  0x38(%rcx),%xmm2
  movss  0x38(%rax),%xmm3
  mulss  0x38(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x38(%rsi),%xmm2
  movss  %xmm2,0xe0(%r10)
  movss  %xmm2,0x1c(%r9)
  add    $0x40,%rdi
  add    $0x40,%rcx
  add    $0x40,%rax
  add    $0x40,%rsi
  add    $0x40,%rdx
  add    $0x20,%rbp
  add    -0x78(%rsp),%rbx
  add    $0x20,%r11
  add    $0x40,%r8
  add    $0x100,%r10
  add    $0x20,%r9
  cmp    %r12,%rdi
  jne    L126
  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    0x8(,%rax,8),%ecx
L11:
  mov    -0x60(%rsp),%eax
  cmp    %ecx,%eax
  jle    L8
  mov    %eax,%esi
  sub    %ecx,%esi
  test   %esi,%esi
  jle    L8
  lea    0x0(,%rcx,8),%eax
  cltq
  mov    -0x4(%rsp),%edx
  imul   %ecx,%edx
  movslq %edx,%rdx
  add    %rdx,%rax
  mov    -0x40(%rsp),%rbx
  lea    (%rbx,%rax,4),%rdx
  movslq %ecx,%rcx
  lea    0x0(,%rcx,4),%rax
  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  D41(%rip),%xmm3
  mov    -0x38(%rsp),%r8
  mov    -0x30(%rsp),%r9
  mov    -0x28(%rsp),%r10
  mov    -0x20(%rsp),%r11
  mov    -0x10(%rsp),%rdi
L41:
  movaps %xmm3,%xmm1
  divss  (%rcx,%rax,2),%xmm1
  movss  %xmm1,0x0(%r13)
  movaps %xmm3,%xmm2
  divss  0x4(%rcx,%rax,2),%xmm2
  movss  %xmm2,0x4(%r13)
  mulss  (%r15),%xmm1
  movss  %xmm1,(%r12)
  movss  0x4(%r15),%xmm1
  mulss  0x4(%r13),%xmm1
  movss  %xmm1,0x4(%r12)
  movaps %xmm0,%xmm1
  mulss  0x0(%r13),%xmm1
  movss  %xmm1,(%r14)
  movaps %xmm0,%xmm1
  mulss  0x4(%r13),%xmm1
  movss  %xmm1,0x4(%r14)
  movss  (%r10,%rax,1),%xmm1
  addss  (%r12),%xmm1
  addss  0x4(%r12),%xmm1
  movss  %xmm1,(%rdx)
  movss  (%r11,%rax,1),%xmm1
  addss  0x4(%r15),%xmm1
  movss  0x4(%r12),%xmm2
  mulss  0x4(%rdi,%rax,2),%xmm2
  addss  %xmm2,%xmm1
  addss  0x4(%r14),%xmm1
  subss  (%r15),%xmm1
  movss  (%r12),%xmm2
  mulss  (%rdi,%rax,2),%xmm2
  subss  %xmm2,%xmm1
  subss  (%r14),%xmm1
  movss  %xmm1,(%r8,%rax,8)
  movss  %xmm1,(%r9,%rax,1)
  add    $0x8,%r15
  add    $0x8,%r12
  add    $0x8,%r14
  add    $0x8,%r13
  add    $0x24,%rdx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L41
L8:
  mov    0x4(%rsp),%eax
  cmp    $0x1,%eax
  jle    L80
  mov    %eax,%eax
  mov    %rax,-0x10(%rsp)
  mov    -0x4(%rsp),%eax
  lea    0x40(,%rax,8),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,-0x78(%rsp)
  mov    -0x64(%rsp),%eax
  sub    $0x8,%eax
  shr    $0x3,%eax
  add    $0x1,%eax
  mov    %eax,%edx
  shl    $0x6,%rdx
  mov    %rdx,-0x58(%rsp)
  shl    $0x3,%eax
  mov    %eax,-0x8(%rsp)
  movq   $0x1,-0x70(%rsp)
  movss  D41(%rip),%xmm1
  movaps %xmm1,%xmm3
  jmp    L121
L60:
  mov    %rdi,-0x70(%rsp)
  mov    $0x0,%ecx
  jmp    L11
L26:
  mov    %rdi,-0x60(%rsp)
  mov    $0x0,%eax
  jmp    L185
L52:
  movaps %xmm3,%xmm2
  divss  (%rax,%rdx,2),%xmm2
  movss  %xmm2,0x0(%r13)
  movaps %xmm3,%xmm4
  divss  0x4(%rax,%rdx,2),%xmm4
  movss  %xmm4,0x4(%r13)
  mulss  (%r15),%xmm2
  movss  %xmm2,(%r12)
  movss  0x4(%r15),%xmm2
  mulss  0x4(%r13),%xmm2
  movss  %xmm2,0x4(%r12)
  movaps %xmm0,%xmm2
  mulss  0x0(%r13),%xmm2
  movss  %xmm2,(%r14)
  movaps %xmm0,%xmm2
  mulss  0x4(%r13),%xmm2
  movss  %xmm2,0x4(%r14)
  movss  (%r10,%rdx,1),%xmm2
  addss  (%r12),%xmm2
  addss  0x4(%r12),%xmm2
  movss  %xmm2,(%rcx)
  movss  (%r11,%rdx,1),%xmm2
  addss  0x4(%r15),%xmm2
  movss  0x4(%r12),%xmm4
  mulss  0x4(%rdi,%rdx,2),%xmm4
  addss  %xmm4,%xmm2
  addss  0x4(%r14),%xmm2
  subss  (%r15),%xmm2
  movss  (%r12),%xmm4
  mulss  (%rdi,%rdx,2),%xmm4
  subss  %xmm4,%xmm2
  subss  (%r14),%xmm2
  movss  %xmm2,(%r8,%rdx,8)
  movss  %xmm2,(%r9,%rdx,1)
  add    $0x8,%r15
  add    $0x8,%r12
  add    $0x8,%r14
  add    $0x8,%r13
  add    $0x24,%rcx
  add    $0x4,%rdx
  cmp    %rsi,%rdx
  jne    L52
L105:
  addq   $0x1,-0x70(%rsp)
  mov    -0x70(%rsp),%rax
  cmp    -0x10(%rsp),%rax
  je     L80
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   $0x7,-0x64(%rsp)
  jle    L26
  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
L178:
  movaps %xmm1,%xmm2
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  movaps %xmm1,%xmm4
  divss  0x4(%rdi),%xmm4
  movss  %xmm4,0x4(%rdx)
  mulss  (%rcx),%xmm2
  movss  %xmm2,(%rax)
  movss  0x4(%rcx),%xmm2
  mulss  0x4(%rdx),%xmm2
  movss  %xmm2,0x4(%rax)
  movaps %xmm0,%xmm2
  mulss  (%rdx),%xmm2
  movss  %xmm2,(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x4(%rdx),%xmm2
  movss  %xmm2,0x4(%rsi)
  movss  0x0(%rbp),%xmm2
  addss  (%rax),%xmm2
  addss  0x4(%rax),%xmm2
  movss  %xmm2,(%rbx)
  movss  (%r11),%xmm2
  addss  0x4(%rcx),%xmm2
  movss  0x4(%rax),%xmm4
  mulss  0x4(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x4(%rsi),%xmm2
  subss  (%rcx),%xmm2
  movss  (%rax),%xmm4
  mulss  (%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  (%rsi),%xmm2
  movss  %xmm2,(%r10)
  movss  %xmm2,(%r9)
  movaps %xmm1,%xmm2
  divss  0x8(%rdi),%xmm2
  movss  %xmm2,0x8(%rdx)
  movaps %xmm1,%xmm4
  divss  0xc(%rdi),%xmm4
  movss  %xmm4,0xc(%rdx)
  mulss  0x8(%rcx),%xmm2
  movss  %xmm2,0x8(%rax)
  movss  0xc(%rcx),%xmm2
  mulss  0xc(%rdx),%xmm2
  movss  %xmm2,0xc(%rax)
  movaps %xmm0,%xmm2
  mulss  0x8(%rdx),%xmm2
  movss  %xmm2,0x8(%rsi)
  movaps %xmm0,%xmm2
  mulss  0xc(%rdx),%xmm2
  movss  %xmm2,0xc(%rsi)
  movss  0x4(%rbp),%xmm2
  addss  0x8(%rax),%xmm2
  addss  0xc(%rax),%xmm2
  movss  %xmm2,0x24(%rbx)
  movss  0x4(%r11),%xmm2
  addss  0xc(%rcx),%xmm2
  movss  0xc(%rax),%xmm4
  mulss  0xc(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0xc(%rsi),%xmm2
  subss  0x8(%rcx),%xmm2
  movss  0x8(%rax),%xmm4
  mulss  0x8(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x8(%rsi),%xmm2
  movss  %xmm2,0x20(%r10)
  movss  %xmm2,0x4(%r9)
  movaps %xmm1,%xmm2
  divss  0x10(%rdi),%xmm2
  movss  %xmm2,0x10(%rdx)
  movaps %xmm1,%xmm4
  divss  0x14(%rdi),%xmm4
  movss  %xmm4,0x14(%rdx)
  mulss  0x10(%rcx),%xmm2
  movss  %xmm2,0x10(%rax)
  movss  0x14(%rcx),%xmm2
  mulss  0x14(%rdx),%xmm2
  movss  %xmm2,0x14(%rax)
  movaps %xmm0,%xmm2
  mulss  0x10(%rdx),%xmm2
  movss  %xmm2,0x10(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x14(%rdx),%xmm2
  movss  %xmm2,0x14(%rsi)
  movss  0x8(%rbp),%xmm2
  addss  0x10(%rax),%xmm2
  addss  0x14(%rax),%xmm2
  movss  %xmm2,0x48(%rbx)
  movss  0x8(%r11),%xmm2
  addss  0x14(%rcx),%xmm2
  movss  0x14(%rax),%xmm4
  mulss  0x14(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x14(%rsi),%xmm2
  subss  0x10(%rcx),%xmm2
  movss  0x10(%rax),%xmm4
  mulss  0x10(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x10(%rsi),%xmm2
  movss  %xmm2,0x40(%r10)
  movss  %xmm2,0x8(%r9)
  movaps %xmm1,%xmm2
  divss  0x18(%rdi),%xmm2
  movss  %xmm2,0x18(%rdx)
  movaps %xmm1,%xmm4
  divss  0x1c(%rdi),%xmm4
  movss  %xmm4,0x1c(%rdx)
  mulss  0x18(%rcx),%xmm2
  movss  %xmm2,0x18(%rax)
  movss  0x1c(%rcx),%xmm2
  mulss  0x1c(%rdx),%xmm2
  movss  %xmm2,0x1c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x18(%rdx),%xmm2
  movss  %xmm2,0x18(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x1c(%rdx),%xmm2
  movss  %xmm2,0x1c(%rsi)
  movss  0xc(%rbp),%xmm2
  addss  0x18(%rax),%xmm2
  addss  0x1c(%rax),%xmm2
  movss  %xmm2,0x6c(%rbx)
  movss  0xc(%r11),%xmm2
  addss  0x1c(%rcx),%xmm2
  movss  0x1c(%rax),%xmm4
  mulss  0x1c(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x1c(%rsi),%xmm2
  subss  0x18(%rcx),%xmm2
  movss  0x18(%rax),%xmm4
  mulss  0x18(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x18(%rsi),%xmm2
  movss  %xmm2,0x60(%r10)
  movss  %xmm2,0xc(%r9)
  movaps %xmm1,%xmm2
  divss  0x20(%rdi),%xmm2
  movss  %xmm2,0x20(%rdx)
  movaps %xmm1,%xmm4
  divss  0x24(%rdi),%xmm4
  movss  %xmm4,0x24(%rdx)
  mulss  0x20(%rcx),%xmm2
  movss  %xmm2,0x20(%rax)
  movss  0x24(%rcx),%xmm2
  mulss  0x24(%rdx),%xmm2
  movss  %xmm2,0x24(%rax)
  movaps %xmm0,%xmm2
  mulss  0x20(%rdx),%xmm2
  movss  %xmm2,0x20(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x24(%rdx),%xmm2
  movss  %xmm2,0x24(%rsi)
  movss  0x10(%rbp),%xmm2
  addss  0x20(%rax),%xmm2
  addss  0x24(%rax),%xmm2
  movss  %xmm2,0x90(%rbx)
  movss  0x10(%r11),%xmm2
  addss  0x24(%rcx),%xmm2
  movss  0x24(%rax),%xmm4
  mulss  0x24(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x24(%rsi),%xmm2
  subss  0x20(%rcx),%xmm2
  movss  0x20(%rax),%xmm4
  mulss  0x20(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x20(%rsi),%xmm2
  movss  %xmm2,0x80(%r10)
  movss  %xmm2,0x10(%r9)
  movaps %xmm1,%xmm2
  divss  0x28(%rdi),%xmm2
  movss  %xmm2,0x28(%rdx)
  movaps %xmm1,%xmm4
  divss  0x2c(%rdi),%xmm4
  movss  %xmm4,0x2c(%rdx)
  mulss  0x28(%rcx),%xmm2
  movss  %xmm2,0x28(%rax)
  movss  0x2c(%rcx),%xmm2
  mulss  0x2c(%rdx),%xmm2
  movss  %xmm2,0x2c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x28(%rdx),%xmm2
  movss  %xmm2,0x28(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x2c(%rdx),%xmm2
  movss  %xmm2,0x2c(%rsi)
  movss  0x14(%rbp),%xmm2
  addss  0x28(%rax),%xmm2
  addss  0x2c(%rax),%xmm2
  movss  %xmm2,0xb4(%rbx)
  movss  0x14(%r11),%xmm2
  addss  0x2c(%rcx),%xmm2
  movss  0x2c(%rax),%xmm4
  mulss  0x2c(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x2c(%rsi),%xmm2
  subss  0x28(%rcx),%xmm2
  movss  0x28(%rax),%xmm4
  mulss  0x28(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x28(%rsi),%xmm2
  movss  %xmm2,0xa0(%r10)
  movss  %xmm2,0x14(%r9)
  movaps %xmm1,%xmm2
  divss  0x30(%rdi),%xmm2
  movss  %xmm2,0x30(%rdx)
  movaps %xmm1,%xmm4
  divss  0x34(%rdi),%xmm4
  movss  %xmm4,0x34(%rdx)
  mulss  0x30(%rcx),%xmm2
  movss  %xmm2,0x30(%rax)
  movss  0x34(%rcx),%xmm2
  mulss  0x34(%rdx),%xmm2
  movss  %xmm2,0x34(%rax)
  movaps %xmm0,%xmm2
  mulss  0x30(%rdx),%xmm2
  movss  %xmm2,0x30(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x34(%rdx),%xmm2
  movss  %xmm2,0x34(%rsi)
  movss  0x18(%rbp),%xmm2
  addss  0x30(%rax),%xmm2
  addss  0x34(%rax),%xmm2
  movss  %xmm2,0xd8(%rbx)
  movss  0x18(%r11),%xmm2
  addss  0x34(%rcx),%xmm2
  movss  0x34(%rax),%xmm4
  mulss  0x34(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x34(%rsi),%xmm2
  subss  0x30(%rcx),%xmm2
  movss  0x30(%rax),%xmm4
  mulss  0x30(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x30(%rsi),%xmm2
  movss  %xmm2,0xc0(%r10)
  movss  %xmm2,0x18(%r9)
  movaps %xmm1,%xmm2
  divss  0x38(%rdi),%xmm2
  movss  %xmm2,0x38(%rdx)
  movaps %xmm1,%xmm4
  divss  0x3c(%rdi),%xmm4
  movss  %xmm4,0x3c(%rdx)
  mulss  0x38(%rcx),%xmm2
  movss  %xmm2,0x38(%rax)
  movss  0x3c(%rcx),%xmm2
  mulss  0x3c(%rdx),%xmm2
  movss  %xmm2,0x3c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x38(%rdx),%xmm2
  movss  %xmm2,0x38(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x3c(%rdx),%xmm2
  movss  %xmm2,0x3c(%rsi)
  movss  0x1c(%rbp),%xmm2
  addss  0x38(%rax),%xmm2
  addss  0x3c(%rax),%xmm2
  movss  %xmm2,0xfc(%rbx)
  movss  0x1c(%r11),%xmm2
  addss  0x3c(%rcx),%xmm2
  movss  0x3c(%rax),%xmm4
  mulss  0x3c(%r8),%xmm4
  addss  %xmm4,%xmm2
  addss  0x3c(%rsi),%xmm2
  subss  0x38(%rcx),%xmm2
  movss  0x38(%rax),%xmm4
  mulss  0x38(%r8),%xmm4
  subss  %xmm4,%xmm2
  subss  0x38(%rsi),%xmm2
  movss  %xmm2,0xe0(%r10)
  movss  %xmm2,0x1c(%r9)
  add    $0x40,%rdi
  add    $0x40,%rcx
  add    $0x40,%rax
  add    $0x40,%rsi
  add    $0x40,%rdx
  add    $0x20,%rbp
  add    -0x78(%rsp),%rbx
  add    $0x20,%r11
  add    $0x40,%r8
  add    $0x100,%r10
  add    $0x20,%r9
  cmp    %r12,%rdi
  jne    L178
  mov    -0x50(%rsp),%r12
  mov    -0x58(%rsp),%rax
  add    %rax,%r15
  add    %rax,%r12
  add    %rax,%r14
  add    %rax,%r13
  mov    -0x8(%rsp),%eax
L185:
  cmp    %eax,-0x64(%rsp)
  jle    L105
  mov    -0x64(%rsp),%esi
  sub    %eax,%esi
  test   %esi,%esi
  jle    L105
  lea    0x0(,%rax,8),%edx
  movslq %edx,%rdx
  mov    -0x4(%rsp),%ecx
  imul   %eax,%ecx
  movslq %ecx,%rcx
  add    %rcx,%rdx
  mov    -0x48(%rsp),%rbx
  lea    (%rbx,%rdx,4),%rcx
  cltq
  lea    0x0(,%rax,4),%rdx
  lea    -0x1(%rsi),%esi
  lea    0x1(%rax,%rsi,1),%rsi
  shl    $0x2,%rsi
  neg    %rax
  mov    -0x60(%rsp),%rbx
  lea    (%rbx,%rax,8),%rax
  mov    -0x40(%rsp),%r8
  mov    -0x38(%rsp),%r9
  mov    -0x30(%rsp),%r10
  mov    -0x28(%rsp),%r11
  mov    -0x20(%rsp),%rdi
  jmp    L52
L80:
  movslq 0x4(%rsp),%rdx
  mov    (%rsp),%ebx
  lea    (%rbx,%rbx,1),%eax
  cltq
  shl    $0x2,%rax
  mov    -0x18(%rsp),%rsi
  mov    %rax,%rdi
  add    (%rsi,%rdx,8),%rdi
  mov    0x48(%rsp),%rsi
  mov    %rax,%r15
  add    (%rsi,%rdx,8),%r15
  mov    0x50(%rsp),%rsi
  mov    %rax,%r14
  add    (%rsi,%rdx,8),%r14
  mov    0x58(%rsp),%rcx
  mov    %rax,%rsi
  add    (%rcx,%rdx,8),%rsi
  mov    %rsi,-0x70(%rsp)
  mov    0x40(%rsp),%rcx
  mov    %rax,%r10
  add    (%rcx,%rdx,8),%r10
  mov    %r10,-0x60(%rsp)
  mov    0x70(%rsp),%rcx
  mov    (%rcx,%rdx,8),%r13
  mov    %r13,-0x20(%rsp)
  mov    0x78(%rsp),%rcx
  mov    (%rcx,%rdx,8),%rcx
  mov    %rcx,-0x58(%rsp)
  mov    0x80(%rsp),%r11
  mov    (%r11,%rdx,8),%r11
  mov    %r11,-0x18(%rsp)
  mov    0x60(%rsp),%r9
  mov    (%r9,%rdx,8),%r8
  mov    %r8,-0x10(%rsp)
  mov    0x68(%rsp),%r9
  mov    (%r9,%rdx,8),%r9
  mov    %r9,-0x50(%rsp)
  mov    0x88(%rsp),%rbp
  mov    0x0(%rbp,%rdx,8),%rdx
  mov    %rdx,-0x48(%rsp)
  mov    -0x64(%rsp),%ebp
  lea    -0x7(%rbp),%edx
  mov    %ebx,%ecx
  cmp    %edx,%ebx
  jge    L142
  movslq %ebx,%r12
  lea    0x0(,%r12,4),%r9
  lea    (%r8,%r9,1),%rbp
  mov    -0x50(%rsp),%rbx
  add    %r9,%rbx
  mov    -0x48(%rsp),%rdx
  lea    (%rdx,%rax,1),%r8
  add    %r11,%r9
  mov    -0x4(%rsp),%r11d
  lea    0x40(,%r11,8),%edx
  movslq %edx,%rdx
  shl    $0x2,%rdx
  mov    %rdx,-0x78(%rsp)
  mov    -0x64(%rsp),%edx
  sub    $0x8,%edx
  sub    %ecx,%edx
  shr    $0x3,%edx
  mov    %edx,-0x38(%rsp)
  add    $0x1,%edx
  mov    %edx,%edx
  shl    $0x6,%rdx
  mov    %rdx,-0x30(%rsp)
  add    %rdi,%rdx
  mov    %rdx,-0x40(%rsp)
  mov    %r11d,%edx
  imul   %ecx,%edx
  movslq %edx,%rdx
  add    %rax,%rdx
  lea    0x0(%r13,%rdx,4),%r13
  shl    $0x5,%r12
  add    -0x58(%rsp),%r12
  mov    %r10,%rdx
  mov    %r14,%rax
  mov    %r15,%rcx
  mov    $0x0,%r11d
  mov    $0x0,%r10d
  movss  D41(%rip),%xmm1
  mov    %r14,-0x28(%rsp)
  mov    -0x40(%rsp),%r14
L30:
  movaps %xmm1,%xmm2
  divss  (%rdi),%xmm2
  movss  %xmm2,(%rdx)
  movaps %xmm1,%xmm3
  divss  0x4(%rdi),%xmm3
  movss  %xmm3,0x4(%rdx)
  mulss  (%rcx),%xmm2
  movss  %xmm2,(%rax)
  movss  0x4(%rcx),%xmm2
  mulss  0x4(%rdx),%xmm2
  movss  %xmm2,0x4(%rax)
  movaps %xmm0,%xmm2
  mulss  (%rdx),%xmm2
  movss  %xmm2,(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x4(%rdx),%xmm2
  movss  %xmm2,0x4(%rsi)
  movss  0x0(%rbp),%xmm2
  addss  (%rax),%xmm2
  addss  0x4(%rax),%xmm2
  movss  %xmm2,0x0(%r13,%r10,1)
  movss  (%rbx),%xmm2
  addss  0x4(%rcx),%xmm2
  movss  0x4(%rax),%xmm3
  mulss  0x4(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x4(%rsi),%xmm2
  subss  (%rcx),%xmm2
  movss  (%rax),%xmm3
  mulss  (%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  (%rsi),%xmm2
  movss  %xmm2,(%r12,%r11,1)
  movss  %xmm2,(%r9)
  movaps %xmm1,%xmm2
  divss  0x8(%rdi),%xmm2
  movss  %xmm2,0x8(%rdx)
  movaps %xmm1,%xmm3
  divss  0xc(%rdi),%xmm3
  movss  %xmm3,0xc(%rdx)
  mulss  0x8(%rcx),%xmm2
  movss  %xmm2,0x8(%rax)
  movss  0xc(%rcx),%xmm2
  mulss  0xc(%rdx),%xmm2
  movss  %xmm2,0xc(%rax)
  movaps %xmm0,%xmm2
  mulss  0x8(%rdx),%xmm2
  movss  %xmm2,0x8(%rsi)
  movaps %xmm0,%xmm2
  mulss  0xc(%rdx),%xmm2
  movss  %xmm2,0xc(%rsi)
  movss  0x4(%rbp),%xmm2
  addss  0x8(%rax),%xmm2
  addss  0xc(%rax),%xmm2
  movss  %xmm2,0x24(%r13,%r10,1)
  movss  0x4(%rbx),%xmm2
  addss  0xc(%rcx),%xmm2
  movss  0xc(%rax),%xmm3
  mulss  0xc(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0xc(%rsi),%xmm2
  subss  0x8(%rcx),%xmm2
  movss  0x8(%rax),%xmm3
  mulss  0x8(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x8(%rsi),%xmm2
  movss  %xmm2,0x20(%r12,%r11,1)
  movss  %xmm2,0x4(%r9)
  movaps %xmm1,%xmm2
  divss  0x10(%rdi),%xmm2
  movss  %xmm2,0x10(%rdx)
  movaps %xmm1,%xmm3
  divss  0x14(%rdi),%xmm3
  movss  %xmm3,0x14(%rdx)
  mulss  0x10(%rcx),%xmm2
  movss  %xmm2,0x10(%rax)
  movss  0x14(%rcx),%xmm2
  mulss  0x14(%rdx),%xmm2
  movss  %xmm2,0x14(%rax)
  movaps %xmm0,%xmm2
  mulss  0x10(%rdx),%xmm2
  movss  %xmm2,0x10(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x14(%rdx),%xmm2
  movss  %xmm2,0x14(%rsi)
  movss  0x8(%rbp),%xmm2
  addss  0x10(%rax),%xmm2
  addss  0x14(%rax),%xmm2
  movss  %xmm2,0x48(%r13,%r10,1)
  movss  0x8(%rbx),%xmm2
  addss  0x14(%rcx),%xmm2
  movss  0x14(%rax),%xmm3
  mulss  0x14(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x14(%rsi),%xmm2
  subss  0x10(%rcx),%xmm2
  movss  0x10(%rax),%xmm3
  mulss  0x10(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x10(%rsi),%xmm2
  movss  %xmm2,0x40(%r12,%r11,1)
  movss  %xmm2,0x8(%r9)
  movaps %xmm1,%xmm2
  divss  0x18(%rdi),%xmm2
  movss  %xmm2,0x18(%rdx)
  movaps %xmm1,%xmm3
  divss  0x1c(%rdi),%xmm3
  movss  %xmm3,0x1c(%rdx)
  mulss  0x18(%rcx),%xmm2
  movss  %xmm2,0x18(%rax)
  movss  0x1c(%rcx),%xmm2
  mulss  0x1c(%rdx),%xmm2
  movss  %xmm2,0x1c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x18(%rdx),%xmm2
  movss  %xmm2,0x18(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x1c(%rdx),%xmm2
  movss  %xmm2,0x1c(%rsi)
  movss  0xc(%rbp),%xmm2
  addss  0x18(%rax),%xmm2
  addss  0x1c(%rax),%xmm2
  movss  %xmm2,0x6c(%r13,%r10,1)
  movss  0xc(%rbx),%xmm2
  addss  0x1c(%rcx),%xmm2
  movss  0x1c(%rax),%xmm3
  mulss  0x1c(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x1c(%rsi),%xmm2
  subss  0x18(%rcx),%xmm2
  movss  0x18(%rax),%xmm3
  mulss  0x18(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x18(%rsi),%xmm2
  movss  %xmm2,0x60(%r12,%r11,1)
  movss  %xmm2,0xc(%r9)
  movaps %xmm1,%xmm2
  divss  0x20(%rdi),%xmm2
  movss  %xmm2,0x20(%rdx)
  movaps %xmm1,%xmm3
  divss  0x24(%rdi),%xmm3
  movss  %xmm3,0x24(%rdx)
  mulss  0x20(%rcx),%xmm2
  movss  %xmm2,0x20(%rax)
  movss  0x24(%rcx),%xmm2
  mulss  0x24(%rdx),%xmm2
  movss  %xmm2,0x24(%rax)
  movaps %xmm0,%xmm2
  mulss  0x20(%rdx),%xmm2
  movss  %xmm2,0x20(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x24(%rdx),%xmm2
  movss  %xmm2,0x24(%rsi)
  movss  0x10(%rbp),%xmm2
  addss  0x20(%rax),%xmm2
  addss  0x24(%rax),%xmm2
  movss  %xmm2,0x90(%r13,%r10,1)
  movss  0x10(%rbx),%xmm2
  addss  0x24(%rcx),%xmm2
  movss  0x24(%rax),%xmm3
  mulss  0x24(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x24(%rsi),%xmm2
  subss  0x20(%rcx),%xmm2
  movss  0x20(%rax),%xmm3
  mulss  0x20(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x20(%rsi),%xmm2
  movss  %xmm2,0x80(%r12,%r11,1)
  movss  %xmm2,0x10(%r9)
  movaps %xmm1,%xmm2
  divss  0x28(%rdi),%xmm2
  movss  %xmm2,0x28(%rdx)
  movaps %xmm1,%xmm3
  divss  0x2c(%rdi),%xmm3
  movss  %xmm3,0x2c(%rdx)
  mulss  0x28(%rcx),%xmm2
  movss  %xmm2,0x28(%rax)
  movss  0x2c(%rcx),%xmm2
  mulss  0x2c(%rdx),%xmm2
  movss  %xmm2,0x2c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x28(%rdx),%xmm2
  movss  %xmm2,0x28(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x2c(%rdx),%xmm2
  movss  %xmm2,0x2c(%rsi)
  movss  0x14(%rbp),%xmm2
  addss  0x28(%rax),%xmm2
  addss  0x2c(%rax),%xmm2
  movss  %xmm2,0xb4(%r13,%r10,1)
  movss  0x14(%rbx),%xmm2
  addss  0x2c(%rcx),%xmm2
  movss  0x2c(%rax),%xmm3
  mulss  0x2c(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x2c(%rsi),%xmm2
  subss  0x28(%rcx),%xmm2
  movss  0x28(%rax),%xmm3
  mulss  0x28(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x28(%rsi),%xmm2
  movss  %xmm2,0xa0(%r12,%r11,1)
  movss  %xmm2,0x14(%r9)
  movaps %xmm1,%xmm2
  divss  0x30(%rdi),%xmm2
  movss  %xmm2,0x30(%rdx)
  movaps %xmm1,%xmm3
  divss  0x34(%rdi),%xmm3
  movss  %xmm3,0x34(%rdx)
  mulss  0x30(%rcx),%xmm2
  movss  %xmm2,0x30(%rax)
  movss  0x34(%rcx),%xmm2
  mulss  0x34(%rdx),%xmm2
  movss  %xmm2,0x34(%rax)
  movaps %xmm0,%xmm2
  mulss  0x30(%rdx),%xmm2
  movss  %xmm2,0x30(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x34(%rdx),%xmm2
  movss  %xmm2,0x34(%rsi)
  movss  0x18(%rbp),%xmm2
  addss  0x30(%rax),%xmm2
  addss  0x34(%rax),%xmm2
  movss  %xmm2,0xd8(%r13,%r10,1)
  movss  0x18(%rbx),%xmm2
  addss  0x34(%rcx),%xmm2
  movss  0x34(%rax),%xmm3
  mulss  0x34(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x34(%rsi),%xmm2
  subss  0x30(%rcx),%xmm2
  movss  0x30(%rax),%xmm3
  mulss  0x30(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x30(%rsi),%xmm2
  movss  %xmm2,0xc0(%r12,%r11,1)
  movss  %xmm2,0x18(%r9)
  movaps %xmm1,%xmm2
  divss  0x38(%rdi),%xmm2
  movss  %xmm2,0x38(%rdx)
  movaps %xmm1,%xmm3
  divss  0x3c(%rdi),%xmm3
  movss  %xmm3,0x3c(%rdx)
  mulss  0x38(%rcx),%xmm2
  movss  %xmm2,0x38(%rax)
  movss  0x3c(%rcx),%xmm2
  mulss  0x3c(%rdx),%xmm2
  movss  %xmm2,0x3c(%rax)
  movaps %xmm0,%xmm2
  mulss  0x38(%rdx),%xmm2
  movss  %xmm2,0x38(%rsi)
  movaps %xmm0,%xmm2
  mulss  0x3c(%rdx),%xmm2
  movss  %xmm2,0x3c(%rsi)
  movss  0x1c(%rbp),%xmm2
  addss  0x38(%rax),%xmm2
  addss  0x3c(%rax),%xmm2
  movss  %xmm2,0xfc(%r13,%r10,1)
  movss  0x1c(%rbx),%xmm2
  addss  0x3c(%rcx),%xmm2
  movss  0x3c(%rax),%xmm3
  mulss  0x3c(%r8),%xmm3
  addss  %xmm3,%xmm2
  addss  0x3c(%rsi),%xmm2
  subss  0x38(%rcx),%xmm2
  movss  0x38(%rax),%xmm3
  mulss  0x38(%r8),%xmm3
  subss  %xmm3,%xmm2
  subss  0x38(%rsi),%xmm2
  movss  %xmm2,0xe0(%r12,%r11,1)
  movss  %xmm2,0x1c(%r9)
  add    $0x40,%rdi
  add    $0x40,%rcx
  add    $0x40,%rax
  add    $0x40,%rsi
  add    $0x40,%rdx
  add    $0x20,%rbp
  add    $0x20,%rbx
  add    $0x40,%r8
  add    $0x20,%r9
  add    -0x78(%rsp),%r10
  add    $0x100,%r11
  cmp    %r14,%rdi
  jne    L30
  mov    -0x28(%rsp),%r14
  mov    -0x30(%rsp),%rax
  add    %rax,%r15
  add    %rax,%r14
  add    %rax,-0x70(%rsp)
  add    %rax,-0x60(%rsp)
  mov    (%rsp),%eax
  mov    -0x38(%rsp),%ebx
  lea    0x8(%rax,%rbx,8),%eax
  mov    %eax,(%rsp)
L95:
  mov    -0x64(%rsp),%edi
  mov    (%rsp),%ebx
  cmp    %ebx,%edi
  jle    L135
  lea    0x0(,%rbx,8),%eax
  cltq
  mov    -0x4(%rsp),%edx
  imul   %ebx,%edx
  movslq %edx,%rdx
  add    %rdx,%rax
  mov    -0x20(%rsp),%rsi
  lea    (%rsi,%rax,4),%rdx
  movslq %ebx,%rcx
  lea    0x0(,%rcx,4),%rax
  sub    %ebx,%edi
  lea    -0x1(%rdi),%esi
  lea    0x1(%rcx,%rsi,1),%rsi
  shl    $0x2,%rsi
  neg    %rcx
  mov    -0x40(%rsp),%rbx
  lea    (%rbx,%rcx,8),%rcx
  movss  D41(%rip),%xmm3
  mov    -0x70(%rsp),%r8
  mov    -0x60(%rsp),%rdi
  mov    -0x58(%rsp),%r10
  mov    -0x18(%rsp),%r11
  mov    -0x10(%rsp),%rbx
  mov    -0x50(%rsp),%rbp
  mov    -0x48(%rsp),%r9
L173:
  movaps %xmm3,%xmm1
  divss  (%rcx,%rax,2),%xmm1
  movss  %xmm1,(%rdi)
  movaps %xmm3,%xmm2
  divss  0x4(%rcx,%rax,2),%xmm2
  movss  %xmm2,0x4(%rdi)
  mulss  (%r15),%xmm1
  movss  %xmm1,(%r14)
  movss  0x4(%r15),%xmm1
  mulss  0x4(%rdi),%xmm1
  movss  %xmm1,0x4(%r14)
  movaps %xmm0,%xmm1
  mulss  (%rdi),%xmm1
  movss  %xmm1,(%r8)
  movaps %xmm0,%xmm1
  mulss  0x4(%rdi),%xmm1
  movss  %xmm1,0x4(%r8)
  movss  (%rbx,%rax,1),%xmm1
  addss  (%r14),%xmm1
  addss  0x4(%r14),%xmm1
  movss  %xmm1,(%rdx)
  movss  0x0(%rbp,%rax,1),%xmm1
  addss  0x4(%r15),%xmm1
  movss  0x4(%r14),%xmm2
  mulss  0x4(%r9,%rax,2),%xmm2
  addss  %xmm2,%xmm1
  addss  0x4(%r8),%xmm1
  subss  (%r15),%xmm1
  movss  (%r14),%xmm2
  mulss  (%r9,%rax,2),%xmm2
  subss  %xmm2,%xmm1
  subss  (%r8),%xmm1
  movss  %xmm1,(%r10,%rax,8)
  movss  %xmm1,(%r11,%rax,1)
  add    $0x8,%r15
  add    $0x8,%r14
  add    $0x8,%r8
  add    $0x8,%rdi
  add    $0x24,%rdx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L173
L135:
  add    $0x8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L142:
  mov    %rdi,-0x40(%rsp)
  jmp    L95
```

tool_calls: name: extract_symbolic_data
arguments: label: D41
type_guessed: f32

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