user: Converte this assembly snippet into its equivalent C code:

```asm
<Mvm>:
  endbr64
  push   %r15
  movapd %xmm0,%xmm5
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edi,%ebx
  sub    $0x38,%rsp
  ucomisd D6(%rip),%xmm1
  mov    %esi,0x14(%rsp)
  mov    %rdx,0x18(%rsp)
  mov    %ecx,0x28(%rsp)
  mov    %r8,0x20(%rsp)
  mov    %r9,0x8(%rsp)
  jp     L17
  jne    L17
L125:
  test   %ebx,%ebx
  lea    0x3(%rbx),%edi
  mov    %ebx,%eax
  mov    0x14(%rsp),%esi
  cmovns %ebx,%edi
  sar    $0x1f,%eax
  shr    $0x1e,%eax
  mov    %esi,%edx
  lea    (%rbx,%rax,1),%r11d
  shr    $0x1f,%edx
  sar    $0x2,%edi
  and    $0x3,%r11d
  sub    %eax,%r11d
  mov    %esi,%eax
  shr    $0x1f,%eax
  add    %esi,%eax
  sar    %eax
  mov    %eax,(%rsp)
  lea    (%rsi,%rdx,1),%eax
  and    $0x1,%eax
  sub    %edx,%eax
  mov    %eax,0x2c(%rsp)
  cmp    $0x1,%esi
  jle    L53
  mov    0x20(%rsp),%rax
  movslq 0x28(%rsp),%r15
  lea    0x0(,%rdi,4),%r13d
  xor    %r14d,%r14d
  movslq %r13d,%r13
  mov    0x8(%rsp),%r8
  lea    0x8(%rax),%r12
  mov    0x18(%rsp),%rax
  shl    $0x3,%r15
  shl    $0x3,%r13
  lea    0x8(%rax),%rbp
  nopl   (%rax)
L88:
  movsd  -0x8(%r12),%xmm1
  movsd  (%r12),%xmm2
  mov    -0x8(%rbp),%r10
  mov    0x0(%rbp),%r9
  mulsd  %xmm5,%xmm1
  mulsd  %xmm5,%xmm2
  add    %r15,%r10
  add    %r15,%r9
  cmp    $0x3,%ebx
  jle    L118
  mov    %r8,%rax
  mov    %r10,%rcx
  mov    %r9,%rdx
  xor    %esi,%esi
  nopw   0x0(%rax,%rax,1)
L77:
  movsd  (%rdx),%xmm7
  movsd  (%rcx),%xmm0
  add    $0x1,%esi
  add    $0x20,%rax
  movsd  0x8(%rcx),%xmm6
  movsd  0x10(%rcx),%xmm4
  add    $0x20,%rdx
  add    $0x20,%rcx
  mulsd  %xmm2,%xmm7
  movsd  -0x8(%rcx),%xmm3
  mulsd  %xmm1,%xmm0
  addsd  -0x20(%rax),%xmm0
  mulsd  %xmm1,%xmm6
  addsd  -0x18(%rax),%xmm6
  mulsd  %xmm1,%xmm4
  addsd  -0x10(%rax),%xmm4
  mulsd  %xmm1,%xmm3
  addsd  %xmm7,%xmm0
  addsd  -0x8(%rax),%xmm3
  movsd  %xmm0,-0x20(%rax)
  movsd  -0x18(%rdx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm6,%xmm0
  movsd  %xmm0,-0x18(%rax)
  movsd  -0x10(%rdx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm4,%xmm0
  movsd  %xmm0,-0x10(%rax)
  movsd  -0x8(%rdx),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm3,%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %esi,%edi
  jg     L77
L118:
  add    %r13,%r8
  test   %r11d,%r11d
  jle    L127
  add    %r13,%r10
  add    %r13,%r9
  xor    %eax,%eax
L67:
  movsd  (%r10,%rax,8),%xmm0
  movsd  (%r9,%rax,8),%xmm3
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  addsd  (%r8,%rax,8),%xmm0
  movsd  %xmm0,(%r8,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r11d
  jg     L67
L127:
  add    $0x1,%r14d
  add    $0x10,%r12
  add    $0x10,%rbp
  cmp    %r14d,(%rsp)
  jg     L88
  mov    (%rsp),%eax
  mov    $0x1,%edx
  sub    $0x1,%eax
  add    $0x1,%rax
  cmpl   $0x1,0x14(%rsp)
  cmovle %rdx,%rax
  imul   %r13,%rax
  add    %rax,0x8(%rsp)
L53:
  cmpl   $0x1,0x2c(%rsp)
  jne    L104
  mov    (%rsp),%eax
  mov    0x20(%rsp),%rsi
  movslq 0x28(%rsp),%rdx
  add    %eax,%eax
  cltq
  mulsd  (%rsi,%rax,8),%xmm5
  mov    0x18(%rsp),%rsi
  mov    (%rsi,%rax,8),%rax
  xor    %esi,%esi
  lea    (%rax,%rdx,8),%rcx
  mov    0x8(%rsp),%rax
  mov    %rcx,%rdx
  cmp    $0x3,%ebx
  jle    L12
  cs nopw 0x0(%rax,%rax,1)
L8:
  movsd  (%rdx),%xmm0
  add    $0x1,%esi
  add    $0x20,%rax
  add    $0x20,%rdx
  mulsd  %xmm5,%xmm0
  addsd  -0x20(%rax),%xmm0
  movsd  %xmm0,-0x20(%rax)
  movsd  -0x18(%rdx),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  -0x18(%rax),%xmm0
  movsd  %xmm0,-0x18(%rax)
  movsd  -0x10(%rdx),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  -0x10(%rax),%xmm0
  movsd  %xmm0,-0x10(%rax)
  movsd  -0x8(%rdx),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %esi,%edi
  jg     L8
L12:
  test   %r11d,%r11d
  jle    L104
  lea    0x0(,%rdi,4),%eax
  lea    -0x1(%r11),%esi
  cltq
  lea    0x0(,%rax,8),%rdx
  lea    0x1(%rax,%rsi,1),%rax
  mov    0x8(%rsp),%rsi
  shl    $0x3,%rax
L87:
  movsd  (%rcx,%rdx,1),%xmm0
  mulsd  %xmm5,%xmm0
  addsd  (%rsi,%rdx,1),%xmm0
  movsd  %xmm0,(%rsi,%rdx,1)
  add    $0x8,%rdx
  cmp    %rdx,%rax
  jne    L87
L104:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L17:
  mov    0x8(%rsp),%rsi
  movapd %xmm1,%xmm0
  mov    %ebx,%edi
  mov    $0x1,%eax
  movsd  %xmm5,(%rsp)
  call   <Mscale@plt>
  movsd  (%rsp),%xmm5
  jmp    L125
```

tool_calls: name: read_data
arguments: data_tag: D6
data_type: f64

tool: D6: 1

assistant: ```c
void Mvm(m, n, alpha, A, j0, x, beta, y) int m, n, j0;
double alpha, beta;
double **A, *x, *y;
{
  register int i, j, m4, n2;
  register double *Aref;
  register double tmp;
  register double *Aref0, *Aref1;
  register double tmp0, tmp1;
  register double yval0, yval1, yval2, yval3;
  if (beta != 1.0)
    Mscale(m, beta, y);
  m4 = n2 = 0;
  m4 = m / 4;
  m = m % 4;
  n2 = n / 2;
  n = n % 2;
  for (j = 0; j < n2; j++) {
    tmp0 = alpha * x[2 * j];
    tmp1 = alpha * x[2 * j + 1];
    Aref0 = &(A[2 * j][j0]);
    Aref1 = &(A[2 * j + 1][j0]);
    for (i = 0; i < m4; i++) {
      yval0 = y[4 * i] + tmp0 * Aref0[4 * i];
      yval1 = y[4 * i + 1] + tmp0 * Aref0[4 * i + 1];
      yval2 = y[4 * i + 2] + tmp0 * Aref0[4 * i + 2];
      yval3 = y[4 * i + 3] + tmp0 * Aref0[4 * i + 3];
      y[4 * i] = yval0 + tmp1 * Aref1[4 * i];
      y[4 * i + 1] = yval1 + tmp1 * Aref1[4 * i + 1];
      y[4 * i + 2] = yval2 + tmp1 * Aref1[4 * i + 2];
      y[4 * i + 3] = yval3 + tmp1 * Aref1[4 * i + 3];
    }
    y += 4 * m4;
    Aref0 += 4 * m4;
    Aref1 += 4 * m4;
    for (i = 0; i < m; i++)
      y[i] += tmp0 * Aref0[i] + tmp1 * Aref1[i];
  }
  for (j = 0; j < n; j++) {
    tmp = alpha * x[2 * n2 + j];
    Aref = &(A[2 * n2 + j][j0]);
    for (i = 0; i < m4; i++) {
      y[4 * i] += tmp * Aref[4 * i];
      y[4 * i + 1] += tmp * Aref[4 * i + 1];
      y[4 * i + 2] += tmp * Aref[4 * i + 2];
      y[4 * i + 3] += tmp * Aref[4 * i + 3];
    }
    y += 4 * m4;
    Aref += 4 * m4;
    for (i = 0; i < m; i++)
      y[i] += tmp * Aref[i];
  }
}
```