user: Translate the below assembly code into C programming language:

```asm
<binvcrhs>:
  endbr64
  movsd  0x28(%rdi),%xmm2
  movsd  D1(%rip),%xmm0
  movapd %xmm0,%xmm1
  divsd  (%rdi),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x28(%rdi)
  movsd  0x50(%rdi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x50(%rdi)
  movsd  0x78(%rdi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x78(%rdi)
  movsd  0xa0(%rdi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xa0(%rdi)
  movsd  (%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,(%rsi)
  movsd  0x28(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x28(%rsi)
  movsd  0x50(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x50(%rsi)
  movsd  0x78(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x78(%rsi)
  movsd  0xa0(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xa0(%rsi)
  mulsd  (%rdx),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  0x8(%rdi),%xmm1
  movsd  0x28(%rdi),%xmm3
  movsd  0x30(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x50(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x30(%rdi)
  movsd  0x58(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x58(%rdi)
  movsd  0x80(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0x80(%rdi)
  movsd  0xa0(%rdi),%xmm3
  movsd  0xa8(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa8(%rdi)
  movsd  (%rsi),%xmm3
  movsd  0x8(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x28(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x8(%rsi)
  movsd  0x30(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x50(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x30(%rsi)
  movsd  0x58(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x58(%rsi)
  movsd  0x80(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x80(%rsi)
  movsd  0xa8(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa8(%rsi)
  mulsd  (%rdx),%xmm1
  movsd  0x8(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rdx)
  movsd  0x10(%rdi),%xmm1
  movsd  0x28(%rdi),%xmm3
  movsd  0x38(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x50(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x38(%rdi)
  movsd  0x60(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x60(%rdi)
  movsd  0x88(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x88(%rdi)
  movsd  0xb0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb0(%rdi)
  movsd  (%rsi),%xmm3
  movsd  0x10(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x28(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x10(%rsi)
  movsd  0x38(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x50(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x38(%rsi)
  movsd  0x60(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x60(%rsi)
  movsd  0x88(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x88(%rsi)
  movsd  0xb0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb0(%rsi)
  mulsd  (%rdx),%xmm1
  movsd  0x10(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rdx)
  movsd  0x18(%rdi),%xmm1
  movsd  0x28(%rdi),%xmm3
  movsd  0x40(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x50(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x40(%rdi)
  movsd  0x68(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x68(%rdi)
  movsd  0x90(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x90(%rdi)
  movsd  0xb8(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb8(%rdi)
  movsd  (%rsi),%xmm3
  movsd  0x18(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x28(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x18(%rsi)
  movsd  0x40(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x50(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x40(%rsi)
  movsd  0x68(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x68(%rsi)
  movsd  0x90(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x90(%rsi)
  movsd  0xb8(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb8(%rsi)
  mulsd  (%rdx),%xmm1
  movsd  0x18(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x18(%rdx)
  movsd  0x20(%rdi),%xmm1
  movsd  0x28(%rdi),%xmm3
  movsd  0x48(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x50(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x48(%rdi)
  movsd  0x70(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x70(%rdi)
  movsd  0x98(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x98(%rdi)
  movsd  0xc0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rdi)
  movsd  (%rsi),%xmm3
  movsd  0x20(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x28(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x20(%rsi)
  movsd  0x48(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x50(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x48(%rsi)
  movsd  0x70(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x78(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x70(%rsi)
  movsd  0x98(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x98(%rsi)
  movsd  0xc0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rsi)
  mulsd  (%rdx),%xmm1
  movsd  0x20(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movsd  %xmm2,0x20(%rdx)
  movsd  0x58(%rdi),%xmm2
  divsd  0x30(%rdi),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x58(%rdi)
  movsd  0x80(%rdi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x80(%rdi)
  movsd  0xa8(%rdi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xa8(%rdi)
  movsd  0x8(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rsi)
  movsd  0x30(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x30(%rsi)
  movsd  0x58(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x58(%rsi)
  movsd  0x80(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x80(%rsi)
  movsd  0xa8(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xa8(%rsi)
  mulsd  0x8(%rdx),%xmm1
  movsd  %xmm1,0x8(%rdx)
  movsd  0x28(%rdi),%xmm1
  movsd  0x58(%rdi),%xmm3
  movsd  0x50(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x80(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x50(%rdi)
  movsd  0x78(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x78(%rdi)
  movsd  0xa0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa0(%rdi)
  movsd  0x8(%rsi),%xmm3
  movsd  (%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x30(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,(%rsi)
  movsd  0x28(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x58(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x28(%rsi)
  movsd  0x50(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x80(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x50(%rsi)
  movsd  0x78(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x78(%rsi)
  movsd  0xa0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa0(%rsi)
  mulsd  0x8(%rdx),%xmm1
  movsd  (%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,(%rdx)
  movsd  0x38(%rdi),%xmm1
  movsd  0x58(%rdi),%xmm3
  movsd  0x60(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x80(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x60(%rdi)
  movsd  0x88(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x88(%rdi)
  movsd  0xb0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb0(%rdi)
  movsd  0x8(%rsi),%xmm3
  movsd  0x10(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x30(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x10(%rsi)
  movsd  0x38(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x58(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x38(%rsi)
  movsd  0x60(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x80(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x60(%rsi)
  movsd  0x88(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x88(%rsi)
  movsd  0xb0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb0(%rsi)
  mulsd  0x8(%rdx),%xmm1
  movsd  0x10(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rdx)
  movsd  0x40(%rdi),%xmm1
  movsd  0x58(%rdi),%xmm3
  movsd  0x68(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0x68(%rdi)
  movsd  0x80(%rdi),%xmm3
  movsd  0x90(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x90(%rdi)
  movsd  0xb8(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb8(%rdi)
  movsd  0x8(%rsi),%xmm3
  movsd  0x18(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x30(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x18(%rsi)
  movsd  0x40(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x58(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x40(%rsi)
  movsd  0x68(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x80(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x68(%rsi)
  movsd  0x90(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x90(%rsi)
  movsd  0xb8(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb8(%rsi)
  mulsd  0x8(%rdx),%xmm1
  movsd  0x18(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x18(%rdx)
  movsd  0x48(%rdi),%xmm1
  movsd  0x58(%rdi),%xmm3
  movsd  0x70(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x80(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x70(%rdi)
  movsd  0x98(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x98(%rdi)
  movsd  0xc0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rdi)
  movsd  0x8(%rsi),%xmm3
  movsd  0x20(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x30(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x20(%rsi)
  movsd  0x48(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x58(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x48(%rsi)
  movsd  0x70(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x80(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x70(%rsi)
  movsd  0x98(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xa8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x98(%rsi)
  movsd  0xc0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rsi)
  mulsd  0x8(%rdx),%xmm1
  movsd  0x20(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movsd  %xmm2,0x20(%rdx)
  movsd  0x88(%rdi),%xmm2
  divsd  0x60(%rdi),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x88(%rdi)
  movsd  0xb0(%rdi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xb0(%rdi)
  movsd  0x10(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rsi)
  movsd  0x38(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x38(%rsi)
  movsd  0x60(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x60(%rsi)
  movsd  0x88(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x88(%rsi)
  movsd  0xb0(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xb0(%rsi)
  mulsd  0x10(%rdx),%xmm1
  movsd  %xmm1,0x10(%rdx)
  movsd  0x50(%rdi),%xmm1
  movsd  0x78(%rdi),%xmm2
  movsd  0x88(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0xb0(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x78(%rdi)
  movsd  0xa0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa0(%rdi)
  movsd  0x10(%rsi),%xmm3
  movsd  (%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x38(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,(%rsi)
  movsd  0x28(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x60(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x28(%rsi)
  movsd  0x50(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x88(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x50(%rsi)
  movsd  0x78(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x78(%rsi)
  movsd  0xa0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa0(%rsi)
  mulsd  0x10(%rdx),%xmm1
  movsd  (%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,(%rdx)
  movsd  0x58(%rdi),%xmm1
  movsd  0x88(%rdi),%xmm3
  movsd  0x80(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0x80(%rdi)
  movsd  0xb0(%rdi),%xmm3
  movsd  0xa8(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa8(%rdi)
  movsd  0x10(%rsi),%xmm3
  movsd  0x8(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x38(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x8(%rsi)
  movsd  0x30(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x60(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x30(%rsi)
  movsd  0x58(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x88(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x58(%rsi)
  movsd  0x80(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x80(%rsi)
  movsd  0xa8(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa8(%rsi)
  mulsd  0x10(%rdx),%xmm1
  movsd  0x8(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rdx)
  movsd  0x68(%rdi),%xmm1
  movsd  0x88(%rdi),%xmm3
  movsd  0x90(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0xb0(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x90(%rdi)
  movsd  0xb8(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb8(%rdi)
  movsd  0x10(%rsi),%xmm3
  movsd  0x18(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x38(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x18(%rsi)
  movsd  0x40(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x60(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x40(%rsi)
  movsd  0x68(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x88(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x68(%rsi)
  movsd  0x90(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb0(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x90(%rsi)
  movsd  0xb8(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb8(%rsi)
  mulsd  0x10(%rdx),%xmm1
  movsd  0x18(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x18(%rdx)
  movsd  0x70(%rdi),%xmm1
  movsd  0x88(%rdi),%xmm3
  movsd  0x98(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0xb0(%rdi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x98(%rdi)
  movsd  0xc0(%rdi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rdi)
  movsd  0x10(%rsi),%xmm3
  movsd  0x20(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x38(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x20(%rsi)
  movsd  0x48(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x60(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x48(%rsi)
  movsd  0x70(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x88(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x70(%rsi)
  movsd  0x98(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0x98(%rsi)
  movsd  0xb0(%rsi),%xmm3
  movsd  0xc0(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rsi)
  mulsd  0x10(%rdx),%xmm1
  movsd  0x20(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movsd  %xmm2,0x20(%rdx)
  movsd  0xb8(%rdi),%xmm2
  divsd  0x90(%rdi),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xb8(%rdi)
  movsd  0x18(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x18(%rsi)
  movsd  0x40(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x40(%rsi)
  movsd  0x68(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x68(%rsi)
  movsd  0x90(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x90(%rsi)
  movsd  0xb8(%rsi),%xmm2
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0xb8(%rsi)
  mulsd  0x18(%rdx),%xmm1
  movsd  %xmm1,0x18(%rdx)
  movsd  0x78(%rdi),%xmm1
  movsd  0xb8(%rdi),%xmm3
  movsd  0xa0(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa0(%rdi)
  movsd  0x18(%rsi),%xmm3
  movsd  (%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x40(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,(%rsi)
  movsd  0x28(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x68(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x28(%rsi)
  movsd  0x50(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x90(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x50(%rsi)
  movsd  0x78(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x78(%rsi)
  movsd  0xa0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa0(%rsi)
  mulsd  0x18(%rdx),%xmm1
  movsd  (%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,(%rdx)
  movsd  0x80(%rdi),%xmm1
  movsd  0xb8(%rdi),%xmm3
  movsd  0xa8(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa8(%rdi)
  movsd  0x18(%rsi),%xmm3
  movsd  0x8(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x40(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x8(%rsi)
  movsd  0x30(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x68(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x30(%rsi)
  movsd  0x58(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x90(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x58(%rsi)
  movsd  0x80(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x80(%rsi)
  movsd  0xa8(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xa8(%rsi)
  mulsd  0x18(%rdx),%xmm1
  movsd  0x8(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rdx)
  movsd  0x88(%rdi),%xmm1
  movsd  0xb8(%rdi),%xmm3
  movsd  0xb0(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb0(%rdi)
  movsd  0x18(%rsi),%xmm3
  movsd  0x10(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x40(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x10(%rsi)
  movsd  0x38(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x68(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x38(%rsi)
  movsd  0x60(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0x90(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x60(%rsi)
  movsd  0x88(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x88(%rsi)
  movsd  0xb0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xb0(%rsi)
  mulsd  0x18(%rdx),%xmm1
  movsd  0x10(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rdx)
  movsd  0x98(%rdi),%xmm1
  movsd  0xb8(%rdi),%xmm3
  movsd  0xc0(%rdi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rdi)
  movsd  0x18(%rsi),%xmm3
  movsd  0x20(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x40(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x20(%rsi)
  movsd  0x48(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0x48(%rsi)
  movsd  0x68(%rsi),%xmm3
  movsd  0x70(%rsi),%xmm2
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm2
  movsd  0x90(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x70(%rsi)
  movsd  0x98(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  0xb8(%rsi),%xmm3
  mulsd  %xmm1,%xmm3
  movsd  %xmm2,0x98(%rsi)
  movsd  0xc0(%rsi),%xmm2
  subsd  %xmm3,%xmm2
  movsd  %xmm2,0xc0(%rsi)
  mulsd  0x18(%rdx),%xmm1
  movsd  0x20(%rdx),%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0x20(%rdx)
  movsd  0x20(%rsi),%xmm1
  divsd  0xc0(%rdi),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x20(%rsi)
  movsd  0x48(%rsi),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x48(%rsi)
  movsd  0x70(%rsi),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x70(%rsi)
  movsd  0x98(%rsi),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x98(%rsi)
  movsd  0xc0(%rsi),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0xc0(%rsi)
  mulsd  0x20(%rdx),%xmm0
  movsd  %xmm0,0x20(%rdx)
  movsd  0x20(%rsi),%xmm2
  movsd  (%rsi),%xmm1
  movsd  0xa0(%rdi),%xmm0
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  0x48(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,(%rsi)
  movsd  0x28(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  %xmm1,0x28(%rsi)
  movsd  0x70(%rsi),%xmm2
  movsd  0x50(%rsi),%xmm1
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  0x98(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x50(%rsi)
  movsd  0x78(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0xc0(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x78(%rsi)
  movsd  0xa0(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  %xmm1,0xa0(%rsi)
  mulsd  0x20(%rdx),%xmm0
  movsd  (%rdx),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdx)
  movsd  0x20(%rsi),%xmm2
  movsd  0x8(%rsi),%xmm1
  movsd  0xa8(%rdi),%xmm0
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  0x48(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x8(%rsi)
  movsd  0x30(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0x70(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x30(%rsi)
  movsd  0x58(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0x98(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x58(%rsi)
  movsd  0x80(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0xc0(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x80(%rsi)
  movsd  0xa8(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  %xmm1,0xa8(%rsi)
  mulsd  0x20(%rdx),%xmm0
  movsd  0x8(%rdx),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rdx)
  movsd  0xb0(%rdi),%xmm0
  movsd  0x20(%rsi),%xmm2
  movsd  0x10(%rsi),%xmm1
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  0x48(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x10(%rsi)
  movsd  0x38(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0x70(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x38(%rsi)
  movsd  0x60(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0x98(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x60(%rsi)
  movsd  0x88(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0xc0(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x88(%rsi)
  movsd  0xb0(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  %xmm1,0xb0(%rsi)
  mulsd  0x20(%rdx),%xmm0
  movsd  0x10(%rdx),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rdx)
  movsd  0x20(%rsi),%xmm2
  movsd  0x18(%rsi),%xmm1
  movsd  0xb8(%rdi),%xmm0
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  0x48(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x18(%rsi)
  movsd  0x40(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0x70(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x40(%rsi)
  movsd  0x68(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0x98(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x68(%rsi)
  movsd  0x90(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  0xc0(%rsi),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm1,0x90(%rsi)
  movsd  0xb8(%rsi),%xmm1
  subsd  %xmm2,%xmm1
  movsd  %xmm1,0xb8(%rsi)
  mulsd  0x20(%rdx),%xmm0
  movsd  0x18(%rdx),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,0x18(%rdx)
  ret
```

tool_calls: name: retrieve_data
arguments: label: D1
datatype: f64

tool: D1: 1

assistant: ```c
void binvcrhs(double lhs[], double c[], double r[]) {
  double pivot, coeff;
  pivot = 1.00e0 / lhs[(1 - 1) + 5 * (1 - 1)];
  lhs[(1 - 1) + 5 * (2 - 1)] = lhs[(1 - 1) + 5 * (2 - 1)] * pivot;
  lhs[(1 - 1) + 5 * (3 - 1)] = lhs[(1 - 1) + 5 * (3 - 1)] * pivot;
  lhs[(1 - 1) + 5 * (4 - 1)] = lhs[(1 - 1) + 5 * (4 - 1)] * pivot;
  lhs[(1 - 1) + 5 * (5 - 1)] = lhs[(1 - 1) + 5 * (5 - 1)] * pivot;
  c[(1 - 1) + 5 * (1 - 1)] = c[(1 - 1) + 5 * (1 - 1)] * pivot;
  c[(1 - 1) + 5 * (2 - 1)] = c[(1 - 1) + 5 * (2 - 1)] * pivot;
  c[(1 - 1) + 5 * (3 - 1)] = c[(1 - 1) + 5 * (3 - 1)] * pivot;
  c[(1 - 1) + 5 * (4 - 1)] = c[(1 - 1) + 5 * (4 - 1)] * pivot;
  c[(1 - 1) + 5 * (5 - 1)] = c[(1 - 1) + 5 * (5 - 1)] * pivot;
  r[1 - 1] = r[1 - 1] * pivot;
  coeff = lhs[(2 - 1) + 5 * (1 - 1)];
  lhs[(2 - 1) + 5 * (2 - 1)] = lhs[(2 - 1) + 5 * (2 - 1)] - coeff * lhs[(1 - 1) + 5 * (2 - 1)];
  lhs[(2 - 1) + 5 * (3 - 1)] = lhs[(2 - 1) + 5 * (3 - 1)] - coeff * lhs[(1 - 1) + 5 * (3 - 1)];
  lhs[(2 - 1) + 5 * (4 - 1)] = lhs[(2 - 1) + 5 * (4 - 1)] - coeff * lhs[(1 - 1) + 5 * (4 - 1)];
  lhs[(2 - 1) + 5 * (5 - 1)] = lhs[(2 - 1) + 5 * (5 - 1)] - coeff * lhs[(1 - 1) + 5 * (5 - 1)];
  c[(2 - 1) + 5 * (1 - 1)] = c[(2 - 1) + 5 * (1 - 1)] - coeff * c[(1 - 1) + 5 * (1 - 1)];
  c[(2 - 1) + 5 * (2 - 1)] = c[(2 - 1) + 5 * (2 - 1)] - coeff * c[(1 - 1) + 5 * (2 - 1)];
  c[(2 - 1) + 5 * (3 - 1)] = c[(2 - 1) + 5 * (3 - 1)] - coeff * c[(1 - 1) + 5 * (3 - 1)];
  c[(2 - 1) + 5 * (4 - 1)] = c[(2 - 1) + 5 * (4 - 1)] - coeff * c[(1 - 1) + 5 * (4 - 1)];
  c[(2 - 1) + 5 * (5 - 1)] = c[(2 - 1) + 5 * (5 - 1)] - coeff * c[(1 - 1) + 5 * (5 - 1)];
  r[2 - 1] = r[2 - 1] - coeff * r[1 - 1];
  coeff = lhs[(3 - 1) + 5 * (1 - 1)];
  lhs[(3 - 1) + 5 * (2 - 1)] = lhs[(3 - 1) + 5 * (2 - 1)] - coeff * lhs[(1 - 1) + 5 * (2 - 1)];
  lhs[(3 - 1) + 5 * (3 - 1)] = lhs[(3 - 1) + 5 * (3 - 1)] - coeff * lhs[(1 - 1) + 5 * (3 - 1)];
  lhs[(3 - 1) + 5 * (4 - 1)] = lhs[(3 - 1) + 5 * (4 - 1)] - coeff * lhs[(1 - 1) + 5 * (4 - 1)];
  lhs[(3 - 1) + 5 * (5 - 1)] = lhs[(3 - 1) + 5 * (5 - 1)] - coeff * lhs[(1 - 1) + 5 * (5 - 1)];
  c[(3 - 1) + 5 * (1 - 1)] = c[(3 - 1) + 5 * (1 - 1)] - coeff * c[(1 - 1) + 5 * (1 - 1)];
  c[(3 - 1) + 5 * (2 - 1)] = c[(3 - 1) + 5 * (2 - 1)] - coeff * c[(1 - 1) + 5 * (2 - 1)];
  c[(3 - 1) + 5 * (3 - 1)] = c[(3 - 1) + 5 * (3 - 1)] - coeff * c[(1 - 1) + 5 * (3 - 1)];
  c[(3 - 1) + 5 * (4 - 1)] = c[(3 - 1) + 5 * (4 - 1)] - coeff * c[(1 - 1) + 5 * (4 - 1)];
  c[(3 - 1) + 5 * (5 - 1)] = c[(3 - 1) + 5 * (5 - 1)] - coeff * c[(1 - 1) + 5 * (5 - 1)];
  r[3 - 1] = r[3 - 1] - coeff * r[1 - 1];
  coeff = lhs[(4 - 1) + 5 * (1 - 1)];
  lhs[(4 - 1) + 5 * (2 - 1)] = lhs[(4 - 1) + 5 * (2 - 1)] - coeff * lhs[(1 - 1) + 5 * (2 - 1)];
  lhs[(4 - 1) + 5 * (3 - 1)] = lhs[(4 - 1) + 5 * (3 - 1)] - coeff * lhs[(1 - 1) + 5 * (3 - 1)];
  lhs[(4 - 1) + 5 * (4 - 1)] = lhs[(4 - 1) + 5 * (4 - 1)] - coeff * lhs[(1 - 1) + 5 * (4 - 1)];
  lhs[(4 - 1) + 5 * (5 - 1)] = lhs[(4 - 1) + 5 * (5 - 1)] - coeff * lhs[(1 - 1) + 5 * (5 - 1)];
  c[(4 - 1) + 5 * (1 - 1)] = c[(4 - 1) + 5 * (1 - 1)] - coeff * c[(1 - 1) + 5 * (1 - 1)];
  c[(4 - 1) + 5 * (2 - 1)] = c[(4 - 1) + 5 * (2 - 1)] - coeff * c[(1 - 1) + 5 * (2 - 1)];
  c[(4 - 1) + 5 * (3 - 1)] = c[(4 - 1) + 5 * (3 - 1)] - coeff * c[(1 - 1) + 5 * (3 - 1)];
  c[(4 - 1) + 5 * (4 - 1)] = c[(4 - 1) + 5 * (4 - 1)] - coeff * c[(1 - 1) + 5 * (4 - 1)];
  c[(4 - 1) + 5 * (5 - 1)] = c[(4 - 1) + 5 * (5 - 1)] - coeff * c[(1 - 1) + 5 * (5 - 1)];
  r[4 - 1] = r[4 - 1] - coeff * r[1 - 1];
  coeff = lhs[(5 - 1) + 5 * (1 - 1)];
  lhs[(5 - 1) + 5 * (2 - 1)] = lhs[(5 - 1) + 5 * (2 - 1)] - coeff * lhs[(1 - 1) + 5 * (2 - 1)];
  lhs[(5 - 1) + 5 * (3 - 1)] = lhs[(5 - 1) + 5 * (3 - 1)] - coeff * lhs[(1 - 1) + 5 * (3 - 1)];
  lhs[(5 - 1) + 5 * (4 - 1)] = lhs[(5 - 1) + 5 * (4 - 1)] - coeff * lhs[(1 - 1) + 5 * (4 - 1)];
  lhs[(5 - 1) + 5 * (5 - 1)] = lhs[(5 - 1) + 5 * (5 - 1)] - coeff * lhs[(1 - 1) + 5 * (5 - 1)];
  c[(5 - 1) + 5 * (1 - 1)] = c[(5 - 1) + 5 * (1 - 1)] - coeff * c[(1 - 1) + 5 * (1 - 1)];
  c[(5 - 1) + 5 * (2 - 1)] = c[(5 - 1) + 5 * (2 - 1)] - coeff * c[(1 - 1) + 5 * (2 - 1)];
  c[(5 - 1) + 5 * (3 - 1)] = c[(5 - 1) + 5 * (3 - 1)] - coeff * c[(1 - 1) + 5 * (3 - 1)];
  c[(5 - 1) + 5 * (4 - 1)] = c[(5 - 1) + 5 * (4 - 1)] - coeff * c[(1 - 1) + 5 * (4 - 1)];
  c[(5 - 1) + 5 * (5 - 1)] = c[(5 - 1) + 5 * (5 - 1)] - coeff * c[(1 - 1) + 5 * (5 - 1)];
  r[5 - 1] = r[5 - 1] - coeff * r[1 - 1];
  pivot = 1.00e0 / lhs[(2 - 1) + 5 * (2 - 1)];
  lhs[(2 - 1) + 5 * (3 - 1)] = lhs[(2 - 1) + 5 * (3 - 1)] * pivot;
  lhs[(2 - 1) + 5 * (4 - 1)] = lhs[(2 - 1) + 5 * (4 - 1)] * pivot;
  lhs[(2 - 1) + 5 * (5 - 1)] = lhs[(2 - 1) + 5 * (5 - 1)] * pivot;
  c[(2 - 1) + 5 * (1 - 1)] = c[(2 - 1) + 5 * (1 - 1)] * pivot;
  c[(2 - 1) + 5 * (2 - 1)] = c[(2 - 1) + 5 * (2 - 1)] * pivot;
  c[(2 - 1) + 5 * (3 - 1)] = c[(2 - 1) + 5 * (3 - 1)] * pivot;
  c[(2 - 1) + 5 * (4 - 1)] = c[(2 - 1) + 5 * (4 - 1)] * pivot;
  c[(2 - 1) + 5 * (5 - 1)] = c[(2 - 1) + 5 * (5 - 1)] * pivot;
  r[2 - 1] = r[2 - 1] * pivot;
  coeff = lhs[(1 - 1) + 5 * (2 - 1)];
  lhs[(1 - 1) + 5 * (3 - 1)] = lhs[(1 - 1) + 5 * (3 - 1)] - coeff * lhs[(2 - 1) + 5 * (3 - 1)];
  lhs[(1 - 1) + 5 * (4 - 1)] = lhs[(1 - 1) + 5 * (4 - 1)] - coeff * lhs[(2 - 1) + 5 * (4 - 1)];
  lhs[(1 - 1) + 5 * (5 - 1)] = lhs[(1 - 1) + 5 * (5 - 1)] - coeff * lhs[(2 - 1) + 5 * (5 - 1)];
  c[(1 - 1) + 5 * (1 - 1)] = c[(1 - 1) + 5 * (1 - 1)] - coeff * c[(2 - 1) + 5 * (1 - 1)];
  c[(1 - 1) + 5 * (2 - 1)] = c[(1 - 1) + 5 * (2 - 1)] - coeff * c[(2 - 1) + 5 * (2 - 1)];
  c[(1 - 1) + 5 * (3 - 1)] = c[(1 - 1) + 5 * (3 - 1)] - coeff * c[(2 - 1) + 5 * (3 - 1)];
  c[(1 - 1) + 5 * (4 - 1)] = c[(1 - 1) + 5 * (4 - 1)] - coeff * c[(2 - 1) + 5 * (4 - 1)];
  c[(1 - 1) + 5 * (5 - 1)] = c[(1 - 1) + 5 * (5 - 1)] - coeff * c[(2 - 1) + 5 * (5 - 1)];
  r[1 - 1] = r[1 - 1] - coeff * r[2 - 1];
  coeff = lhs[(3 - 1) + 5 * (2 - 1)];
  lhs[(3 - 1) + 5 * (3 - 1)] = lhs[(3 - 1) + 5 * (3 - 1)] - coeff * lhs[(2 - 1) + 5 * (3 - 1)];
  lhs[(3 - 1) + 5 * (4 - 1)] = lhs[(3 - 1) + 5 * (4 - 1)] - coeff * lhs[(2 - 1) + 5 * (4 - 1)];
  lhs[(3 - 1) + 5 * (5 - 1)] = lhs[(3 - 1) + 5 * (5 - 1)] - coeff * lhs[(2 - 1) + 5 * (5 - 1)];
  c[(3 - 1) + 5 * (1 - 1)] = c[(3 - 1) + 5 * (1 - 1)] - coeff * c[(2 - 1) + 5 * (1 - 1)];
  c[(3 - 1) + 5 * (2 - 1)] = c[(3 - 1) + 5 * (2 - 1)] - coeff * c[(2 - 1) + 5 * (2 - 1)];
  c[(3 - 1) + 5 * (3 - 1)] = c[(3 - 1) + 5 * (3 - 1)] - coeff * c[(2 - 1) + 5 * (3 - 1)];
  c[(3 - 1) + 5 * (4 - 1)] = c[(3 - 1) + 5 * (4 - 1)] - coeff * c[(2 - 1) + 5 * (4 - 1)];
  c[(3 - 1) + 5 * (5 - 1)] = c[(3 - 1) + 5 * (5 - 1)] - coeff * c[(2 - 1) + 5 * (5 - 1)];
  r[3 - 1] = r[3 - 1] - coeff * r[2 - 1];
  coeff = lhs[(4 - 1) + 5 * (2 - 1)];
  lhs[(4 - 1) + 5 * (3 - 1)] = lhs[(4 - 1) + 5 * (3 - 1)] - coeff * lhs[(2 - 1) + 5 * (3 - 1)];
  lhs[(4 - 1) + 5 * (4 - 1)] = lhs[(4 - 1) + 5 * (4 - 1)] - coeff * lhs[(2 - 1) + 5 * (4 - 1)];
  lhs[(4 - 1) + 5 * (5 - 1)] = lhs[(4 - 1) + 5 * (5 - 1)] - coeff * lhs[(2 - 1) + 5 * (5 - 1)];
  c[(4 - 1) + 5 * (1 - 1)] = c[(4 - 1) + 5 * (1 - 1)] - coeff * c[(2 - 1) + 5 * (1 - 1)];
  c[(4 - 1) + 5 * (2 - 1)] = c[(4 - 1) + 5 * (2 - 1)] - coeff * c[(2 - 1) + 5 * (2 - 1)];
  c[(4 - 1) + 5 * (3 - 1)] = c[(4 - 1) + 5 * (3 - 1)] - coeff * c[(2 - 1) + 5 * (3 - 1)];
  c[(4 - 1) + 5 * (4 - 1)] = c[(4 - 1) + 5 * (4 - 1)] - coeff * c[(2 - 1) + 5 * (4 - 1)];
  c[(4 - 1) + 5 * (5 - 1)] = c[(4 - 1) + 5 * (5 - 1)] - coeff * c[(2 - 1) + 5 * (5 - 1)];
  r[4 - 1] = r[4 - 1] - coeff * r[2 - 1];
  coeff = lhs[(5 - 1) + 5 * (2 - 1)];
  lhs[(5 - 1) + 5 * (3 - 1)] = lhs[(5 - 1) + 5 * (3 - 1)] - coeff * lhs[(2 - 1) + 5 * (3 - 1)];
  lhs[(5 - 1) + 5 * (4 - 1)] = lhs[(5 - 1) + 5 * (4 - 1)] - coeff * lhs[(2 - 1) + 5 * (4 - 1)];
  lhs[(5 - 1) + 5 * (5 - 1)] = lhs[(5 - 1) + 5 * (5 - 1)] - coeff * lhs[(2 - 1) + 5 * (5 - 1)];
  c[(5 - 1) + 5 * (1 - 1)] = c[(5 - 1) + 5 * (1 - 1)] - coeff * c[(2 - 1) + 5 * (1 - 1)];
  c[(5 - 1) + 5 * (2 - 1)] = c[(5 - 1) + 5 * (2 - 1)] - coeff * c[(2 - 1) + 5 * (2 - 1)];
  c[(5 - 1) + 5 * (3 - 1)] = c[(5 - 1) + 5 * (3 - 1)] - coeff * c[(2 - 1) + 5 * (3 - 1)];
  c[(5 - 1) + 5 * (4 - 1)] = c[(5 - 1) + 5 * (4 - 1)] - coeff * c[(2 - 1) + 5 * (4 - 1)];
  c[(5 - 1) + 5 * (5 - 1)] = c[(5 - 1) + 5 * (5 - 1)] - coeff * c[(2 - 1) + 5 * (5 - 1)];
  r[5 - 1] = r[5 - 1] - coeff * r[2 - 1];
  pivot = 1.00e0 / lhs[(3 - 1) + 5 * (3 - 1)];
  lhs[(3 - 1) + 5 * (4 - 1)] = lhs[(3 - 1) + 5 * (4 - 1)] * pivot;
  lhs[(3 - 1) + 5 * (5 - 1)] = lhs[(3 - 1) + 5 * (5 - 1)] * pivot;
  c[(3 - 1) + 5 * (1 - 1)] = c[(3 - 1) + 5 * (1 - 1)] * pivot;
  c[(3 - 1) + 5 * (2 - 1)] = c[(3 - 1) + 5 * (2 - 1)] * pivot;
  c[(3 - 1) + 5 * (3 - 1)] = c[(3 - 1) + 5 * (3 - 1)] * pivot;
  c[(3 - 1) + 5 * (4 - 1)] = c[(3 - 1) + 5 * (4 - 1)] * pivot;
  c[(3 - 1) + 5 * (5 - 1)] = c[(3 - 1) + 5 * (5 - 1)] * pivot;
  r[3 - 1] = r[3 - 1] * pivot;
  coeff = lhs[(1 - 1) + 5 * (3 - 1)];
  lhs[(1 - 1) + 5 * (4 - 1)] = lhs[(1 - 1) + 5 * (4 - 1)] - coeff * lhs[(3 - 1) + 5 * (4 - 1)];
  lhs[(1 - 1) + 5 * (5 - 1)] = lhs[(1 - 1) + 5 * (5 - 1)] - coeff * lhs[(3 - 1) + 5 * (5 - 1)];
  c[(1 - 1) + 5 * (1 - 1)] = c[(1 - 1) + 5 * (1 - 1)] - coeff * c[(3 - 1) + 5 * (1 - 1)];
  c[(1 - 1) + 5 * (2 - 1)] = c[(1 - 1) + 5 * (2 - 1)] - coeff * c[(3 - 1) + 5 * (2 - 1)];
  c[(1 - 1) + 5 * (3 - 1)] = c[(1 - 1) + 5 * (3 - 1)] - coeff * c[(3 - 1) + 5 * (3 - 1)];
  c[(1 - 1) + 5 * (4 - 1)] = c[(1 - 1) + 5 * (4 - 1)] - coeff * c[(3 - 1) + 5 * (4 - 1)];
  c[(1 - 1) + 5 * (5 - 1)] = c[(1 - 1) + 5 * (5 - 1)] - coeff * c[(3 - 1) + 5 * (5 - 1)];
  r[1 - 1] = r[1 - 1] - coeff * r[3 - 1];
  coeff = lhs[(2 - 1) + 5 * (3 - 1)];
  lhs[(2 - 1) + 5 * (4 - 1)] = lhs[(2 - 1) + 5 * (4 - 1)] - coeff * lhs[(3 - 1) + 5 * (4 - 1)];
  lhs[(2 - 1) + 5 * (5 - 1)] = lhs[(2 - 1) + 5 * (5 - 1)] - coeff * lhs[(3 - 1) + 5 * (5 - 1)];
  c[(2 - 1) + 5 * (1 - 1)] = c[(2 - 1) + 5 * (1 - 1)] - coeff * c[(3 - 1) + 5 * (1 - 1)];
  c[(2 - 1) + 5 * (2 - 1)] = c[(2 - 1) + 5 * (2 - 1)] - coeff * c[(3 - 1) + 5 * (2 - 1)];
  c[(2 - 1) + 5 * (3 - 1)] = c[(2 - 1) + 5 * (3 - 1)] - coeff * c[(3 - 1) + 5 * (3 - 1)];
  c[(2 - 1) + 5 * (4 - 1)] = c[(2 - 1) + 5 * (4 - 1)] - coeff * c[(3 - 1) + 5 * (4 - 1)];
  c[(2 - 1) + 5 * (5 - 1)] = c[(2 - 1) + 5 * (5 - 1)] - coeff * c[(3 - 1) + 5 * (5 - 1)];
  r[2 - 1] = r[2 - 1] - coeff * r[3 - 1];
  coeff = lhs[(4 - 1) + 5 * (3 - 1)];
  lhs[(4 - 1) + 5 * (4 - 1)] = lhs[(4 - 1) + 5 * (4 - 1)] - coeff * lhs[(3 - 1) + 5 * (4 - 1)];
  lhs[(4 - 1) + 5 * (5 - 1)] = lhs[(4 - 1) + 5 * (5 - 1)] - coeff * lhs[(3 - 1) + 5 * (5 - 1)];
  c[(4 - 1) + 5 * (1 - 1)] = c[(4 - 1) + 5 * (1 - 1)] - coeff * c[(3 - 1) + 5 * (1 - 1)];
  c[(4 - 1) + 5 * (2 - 1)] = c[(4 - 1) + 5 * (2 - 1)] - coeff * c[(3 - 1) + 5 * (2 - 1)];
  c[(4 - 1) + 5 * (3 - 1)] = c[(4 - 1) + 5 * (3 - 1)] - coeff * c[(3 - 1) + 5 * (3 - 1)];
  c[(4 - 1) + 5 * (4 - 1)] = c[(4 - 1) + 5 * (4 - 1)] - coeff * c[(3 - 1) + 5 * (4 - 1)];
  c[(4 - 1) + 5 * (5 - 1)] = c[(4 - 1) + 5 * (5 - 1)] - coeff * c[(3 - 1) + 5 * (5 - 1)];
  r[4 - 1] = r[4 - 1] - coeff * r[3 - 1];
  coeff = lhs[(5 - 1) + 5 * (3 - 1)];
  lhs[(5 - 1) + 5 * (4 - 1)] = lhs[(5 - 1) + 5 * (4 - 1)] - coeff * lhs[(3 - 1) + 5 * (4 - 1)];
  lhs[(5 - 1) + 5 * (5 - 1)] = lhs[(5 - 1) + 5 * (5 - 1)] - coeff * lhs[(3 - 1) + 5 * (5 - 1)];
  c[(5 - 1) + 5 * (1 - 1)] = c[(5 - 1) + 5 * (1 - 1)] - coeff * c[(3 - 1) + 5 * (1 - 1)];
  c[(5 - 1) + 5 * (2 - 1)] = c[(5 - 1) + 5 * (2 - 1)] - coeff * c[(3 - 1) + 5 * (2 - 1)];
  c[(5 - 1) + 5 * (3 - 1)] = c[(5 - 1) + 5 * (3 - 1)] - coeff * c[(3 - 1) + 5 * (3 - 1)];
  c[(5 - 1) + 5 * (4 - 1)] = c[(5 - 1) + 5 * (4 - 1)] - coeff * c[(3 - 1) + 5 * (4 - 1)];
  c[(5 - 1) + 5 * (5 - 1)] = c[(5 - 1) + 5 * (5 - 1)] - coeff * c[(3 - 1) + 5 * (5 - 1)];
  r[5 - 1] = r[5 - 1] - coeff * r[3 - 1];
  pivot = 1.00e0 / lhs[(4 - 1) + 5 * (4 - 1)];
  lhs[(4 - 1) + 5 * (5 - 1)] = lhs[(4 - 1) + 5 * (5 - 1)] * pivot;
  c[(4 - 1) + 5 * (1 - 1)] = c[(4 - 1) + 5 * (1 - 1)] * pivot;
  c[(4 - 1) + 5 * (2 - 1)] = c[(4 - 1) + 5 * (2 - 1)] * pivot;
  c[(4 - 1) + 5 * (3 - 1)] = c[(4 - 1) + 5 * (3 - 1)] * pivot;
  c[(4 - 1) + 5 * (4 - 1)] = c[(4 - 1) + 5 * (4 - 1)] * pivot;
  c[(4 - 1) + 5 * (5 - 1)] = c[(4 - 1) + 5 * (5 - 1)] * pivot;
  r[4 - 1] = r[4 - 1] * pivot;
  coeff = lhs[(1 - 1) + 5 * (4 - 1)];
  lhs[(1 - 1) + 5 * (5 - 1)] = lhs[(1 - 1) + 5 * (5 - 1)] - coeff * lhs[(4 - 1) + 5 * (5 - 1)];
  c[(1 - 1) + 5 * (1 - 1)] = c[(1 - 1) + 5 * (1 - 1)] - coeff * c[(4 - 1) + 5 * (1 - 1)];
  c[(1 - 1) + 5 * (2 - 1)] = c[(1 - 1) + 5 * (2 - 1)] - coeff * c[(4 - 1) + 5 * (2 - 1)];
  c[(1 - 1) + 5 * (3 - 1)] = c[(1 - 1) + 5 * (3 - 1)] - coeff * c[(4 - 1) + 5 * (3 - 1)];
  c[(1 - 1) + 5 * (4 - 1)] = c[(1 - 1) + 5 * (4 - 1)] - coeff * c[(4 - 1) + 5 * (4 - 1)];
  c[(1 - 1) + 5 * (5 - 1)] = c[(1 - 1) + 5 * (5 - 1)] - coeff * c[(4 - 1) + 5 * (5 - 1)];
  r[1 - 1] = r[1 - 1] - coeff * r[4 - 1];
  coeff = lhs[(2 - 1) + 5 * (4 - 1)];
  lhs[(2 - 1) + 5 * (5 - 1)] = lhs[(2 - 1) + 5 * (5 - 1)] - coeff * lhs[(4 - 1) + 5 * (5 - 1)];
  c[(2 - 1) + 5 * (1 - 1)] = c[(2 - 1) + 5 * (1 - 1)] - coeff * c[(4 - 1) + 5 * (1 - 1)];
  c[(2 - 1) + 5 * (2 - 1)] = c[(2 - 1) + 5 * (2 - 1)] - coeff * c[(4 - 1) + 5 * (2 - 1)];
  c[(2 - 1) + 5 * (3 - 1)] = c[(2 - 1) + 5 * (3 - 1)] - coeff * c[(4 - 1) + 5 * (3 - 1)];
  c[(2 - 1) + 5 * (4 - 1)] = c[(2 - 1) + 5 * (4 - 1)] - coeff * c[(4 - 1) + 5 * (4 - 1)];
  c[(2 - 1) + 5 * (5 - 1)] = c[(2 - 1) + 5 * (5 - 1)] - coeff * c[(4 - 1) + 5 * (5 - 1)];
  r[2 - 1] = r[2 - 1] - coeff * r[4 - 1];
  coeff = lhs[(3 - 1) + 5 * (4 - 1)];
  lhs[(3 - 1) + 5 * (5 - 1)] = lhs[(3 - 1) + 5 * (5 - 1)] - coeff * lhs[(4 - 1) + 5 * (5 - 1)];
  c[(3 - 1) + 5 * (1 - 1)] = c[(3 - 1) + 5 * (1 - 1)] - coeff * c[(4 - 1) + 5 * (1 - 1)];
  c[(3 - 1) + 5 * (2 - 1)] = c[(3 - 1) + 5 * (2 - 1)] - coeff * c[(4 - 1) + 5 * (2 - 1)];
  c[(3 - 1) + 5 * (3 - 1)] = c[(3 - 1) + 5 * (3 - 1)] - coeff * c[(4 - 1) + 5 * (3 - 1)];
  c[(3 - 1) + 5 * (4 - 1)] = c[(3 - 1) + 5 * (4 - 1)] - coeff * c[(4 - 1) + 5 * (4 - 1)];
  c[(3 - 1) + 5 * (5 - 1)] = c[(3 - 1) + 5 * (5 - 1)] - coeff * c[(4 - 1) + 5 * (5 - 1)];
  r[3 - 1] = r[3 - 1] - coeff * r[4 - 1];
  coeff = lhs[(5 - 1) + 5 * (4 - 1)];
  lhs[(5 - 1) + 5 * (5 - 1)] = lhs[(5 - 1) + 5 * (5 - 1)] - coeff * lhs[(4 - 1) + 5 * (5 - 1)];
  c[(5 - 1) + 5 * (1 - 1)] = c[(5 - 1) + 5 * (1 - 1)] - coeff * c[(4 - 1) + 5 * (1 - 1)];
  c[(5 - 1) + 5 * (2 - 1)] = c[(5 - 1) + 5 * (2 - 1)] - coeff * c[(4 - 1) + 5 * (2 - 1)];
  c[(5 - 1) + 5 * (3 - 1)] = c[(5 - 1) + 5 * (3 - 1)] - coeff * c[(4 - 1) + 5 * (3 - 1)];
  c[(5 - 1) + 5 * (4 - 1)] = c[(5 - 1) + 5 * (4 - 1)] - coeff * c[(4 - 1) + 5 * (4 - 1)];
  c[(5 - 1) + 5 * (5 - 1)] = c[(5 - 1) + 5 * (5 - 1)] - coeff * c[(4 - 1) + 5 * (5 - 1)];
  r[5 - 1] = r[5 - 1] - coeff * r[4 - 1];
  pivot = 1.00e0 / lhs[(5 - 1) + 5 * (5 - 1)];
  c[(5 - 1) + 5 * (1 - 1)] = c[(5 - 1) + 5 * (1 - 1)] * pivot;
  c[(5 - 1) + 5 * (2 - 1)] = c[(5 - 1) + 5 * (2 - 1)] * pivot;
  c[(5 - 1) + 5 * (3 - 1)] = c[(5 - 1) + 5 * (3 - 1)] * pivot;
  c[(5 - 1) + 5 * (4 - 1)] = c[(5 - 1) + 5 * (4 - 1)] * pivot;
  c[(5 - 1) + 5 * (5 - 1)] = c[(5 - 1) + 5 * (5 - 1)] * pivot;
  r[5 - 1] = r[5 - 1] * pivot;
  coeff = lhs[(1 - 1) + 5 * (5 - 1)];
  c[(1 - 1) + 5 * (1 - 1)] = c[(1 - 1) + 5 * (1 - 1)] - coeff * c[(5 - 1) + 5 * (1 - 1)];
  c[(1 - 1) + 5 * (2 - 1)] = c[(1 - 1) + 5 * (2 - 1)] - coeff * c[(5 - 1) + 5 * (2 - 1)];
  c[(1 - 1) + 5 * (3 - 1)] = c[(1 - 1) + 5 * (3 - 1)] - coeff * c[(5 - 1) + 5 * (3 - 1)];
  c[(1 - 1) + 5 * (4 - 1)] = c[(1 - 1) + 5 * (4 - 1)] - coeff * c[(5 - 1) + 5 * (4 - 1)];
  c[(1 - 1) + 5 * (5 - 1)] = c[(1 - 1) + 5 * (5 - 1)] - coeff * c[(5 - 1) + 5 * (5 - 1)];
  r[1 - 1] = r[1 - 1] - coeff * r[5 - 1];
  coeff = lhs[(2 - 1) + 5 * (5 - 1)];
  c[(2 - 1) + 5 * (1 - 1)] = c[(2 - 1) + 5 * (1 - 1)] - coeff * c[(5 - 1) + 5 * (1 - 1)];
  c[(2 - 1) + 5 * (2 - 1)] = c[(2 - 1) + 5 * (2 - 1)] - coeff * c[(5 - 1) + 5 * (2 - 1)];
  c[(2 - 1) + 5 * (3 - 1)] = c[(2 - 1) + 5 * (3 - 1)] - coeff * c[(5 - 1) + 5 * (3 - 1)];
  c[(2 - 1) + 5 * (4 - 1)] = c[(2 - 1) + 5 * (4 - 1)] - coeff * c[(5 - 1) + 5 * (4 - 1)];
  c[(2 - 1) + 5 * (5 - 1)] = c[(2 - 1) + 5 * (5 - 1)] - coeff * c[(5 - 1) + 5 * (5 - 1)];
  r[2 - 1] = r[2 - 1] - coeff * r[5 - 1];
  coeff = lhs[(3 - 1) + 5 * (5 - 1)];
  c[(3 - 1) + 5 * (1 - 1)] = c[(3 - 1) + 5 * (1 - 1)] - coeff * c[(5 - 1) + 5 * (1 - 1)];
  c[(3 - 1) + 5 * (2 - 1)] = c[(3 - 1) + 5 * (2 - 1)] - coeff * c[(5 - 1) + 5 * (2 - 1)];
  c[(3 - 1) + 5 * (3 - 1)] = c[(3 - 1) + 5 * (3 - 1)] - coeff * c[(5 - 1) + 5 * (3 - 1)];
  c[(3 - 1) + 5 * (4 - 1)] = c[(3 - 1) + 5 * (4 - 1)] - coeff * c[(5 - 1) + 5 * (4 - 1)];
  c[(3 - 1) + 5 * (5 - 1)] = c[(3 - 1) + 5 * (5 - 1)] - coeff * c[(5 - 1) + 5 * (5 - 1)];
  r[3 - 1] = r[3 - 1] - coeff * r[5 - 1];
  coeff = lhs[(4 - 1) + 5 * (5 - 1)];
  c[(4 - 1) + 5 * (1 - 1)] = c[(4 - 1) + 5 * (1 - 1)] - coeff * c[(5 - 1) + 5 * (1 - 1)];
  c[(4 - 1) + 5 * (2 - 1)] = c[(4 - 1) + 5 * (2 - 1)] - coeff * c[(5 - 1) + 5 * (2 - 1)];
  c[(4 - 1) + 5 * (3 - 1)] = c[(4 - 1) + 5 * (3 - 1)] - coeff * c[(5 - 1) + 5 * (3 - 1)];
  c[(4 - 1) + 5 * (4 - 1)] = c[(4 - 1) + 5 * (4 - 1)] - coeff * c[(5 - 1) + 5 * (4 - 1)];
  c[(4 - 1) + 5 * (5 - 1)] = c[(4 - 1) + 5 * (5 - 1)] - coeff * c[(5 - 1) + 5 * (5 - 1)];
  r[4 - 1] = r[4 - 1] - coeff * r[5 - 1];
  return;
}
```