user: Please show me the C implementation that matches the assembly code listed below:

```asm
<kernel_ssymv_2_lib4>:
  endbr64
  mov    %edi,%r10d
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x30(%rsp),%rdi
  mov    0x38(%rsp),%rbp
  mov    0x48(%rsp),%r13d
  test   %r10d,%r10d
  jle    L72
  mov    %esi,%r11d
  movss  (%r8),%xmm3
  movss  0x4(%r8),%xmm4
  mov    %rdx,%rsi
  cmp    $0x1,%r13d
  je     L111
  movss  D5(%rip),%xmm0
  xorps  %xmm0,%xmm3
  xorps  %xmm0,%xmm4
L111:
  cmpl   $0x1,0x40(%rsp)
  je     L4
  pxor   %xmm0,%xmm0
  test   %r11d,%r11d
  jle    L67
  xor    %edx,%edx
  movaps %xmm0,%xmm6
L113:
  lea    -0x1(%r11),%eax
  sub    %edx,%eax
  lea    0x4(,%rax,4),%rdx
  xor    %eax,%eax
  nopl   0x0(%rax)
L149:
  movss  (%rdi,%rax,1),%xmm2
  movss  (%rsi,%rax,1),%xmm5
  movaps %xmm3,%xmm7
  movss  0x10(%rsi,%rax,1),%xmm1
  mulss  %xmm5,%xmm7
  mulss  %xmm2,%xmm5
  addss  (%r9,%rax,1),%xmm7
  mulss  %xmm1,%xmm2
  mulss  %xmm4,%xmm1
  addss  %xmm5,%xmm6
  addss  %xmm2,%xmm0
  addss  %xmm7,%xmm1
  movss  %xmm1,(%r9,%rax,1)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L149
  add    %rax,%rsi
  add    %rax,%r9
  add    %rax,%rdi
L61:
  lea    -0x4(,%rcx,4),%eax
  cltq
  lea    (%rsi,%rax,4),%rsi
L47:
  lea    -0x3(%r10),%eax
  cmp    %eax,%r11d
  jge    L132
  lea    -0x4(%r10),%ebx
  shl    $0x2,%ecx
  mov    %rdi,%rdx
  mov    %rsi,%rax
  sub    %r11d,%ebx
  movslq %ecx,%rcx
  shr    $0x2,%ebx
  shl    $0x2,%rcx
  lea    0x1(%rbx),%r12d
  mov    %r12,%r14
  shl    $0x4,%r14
  lea    (%r9,%r14,1),%r8
L16:
  movss  (%rdx),%xmm1
  movss  (%rax),%xmm2
  movaps %xmm3,%xmm7
  add    $0x10,%r9
  movss  0x10(%rax),%xmm5
  add    $0x10,%rdx
  mulss  %xmm2,%xmm7
  mulss  %xmm1,%xmm2
  addss  -0x10(%r9),%xmm7
  mulss  %xmm5,%xmm1
  mulss  %xmm4,%xmm5
  addss  %xmm2,%xmm6
  movaps %xmm1,%xmm2
  addss  %xmm7,%xmm5
  addss  %xmm0,%xmm2
  movaps %xmm3,%xmm7
  movss  %xmm5,-0x10(%r9)
  movss  0x4(%rax),%xmm5
  movss  0x14(%rax),%xmm1
  movss  -0xc(%rdx),%xmm0
  mulss  %xmm5,%xmm7
  mulss  %xmm0,%xmm5
  addss  -0xc(%r9),%xmm7
  mulss  %xmm1,%xmm0
  mulss  %xmm4,%xmm1
  addss  %xmm5,%xmm6
  addss  %xmm0,%xmm2
  addss  %xmm7,%xmm1
  movaps %xmm3,%xmm7
  movss  %xmm1,-0xc(%r9)
  movss  0x8(%rax),%xmm5
  movss  0x18(%rax),%xmm1
  movss  -0x8(%rdx),%xmm0
  mulss  %xmm5,%xmm7
  mulss  %xmm0,%xmm5
  addss  -0x8(%r9),%xmm7
  mulss  %xmm1,%xmm0
  mulss  %xmm4,%xmm1
  addss  %xmm6,%xmm5
  addss  %xmm0,%xmm2
  addss  %xmm7,%xmm1
  movaps %xmm3,%xmm7
  movss  %xmm1,-0x8(%r9)
  movss  0xc(%rax),%xmm6
  movss  0x1c(%rax),%xmm1
  add    %rcx,%rax
  movss  -0x4(%rdx),%xmm0
  mulss  %xmm6,%xmm7
  mulss  %xmm0,%xmm6
  addss  -0x4(%r9),%xmm7
  mulss  %xmm1,%xmm0
  mulss  %xmm4,%xmm1
  addss  %xmm5,%xmm6
  addss  %xmm2,%xmm0
  addss  %xmm7,%xmm1
  movss  %xmm1,-0x4(%r9)
  cmp    %r8,%r9
  jne    L16
  imul   %r12,%rcx
  add    %r14,%rdi
  lea    0x4(%r11,%rbx,4),%r11d
  add    %rcx,%rsi
L9:
  cmp    %r11d,%r10d
  jle    L12
  sub    $0x1,%r10d
  xor    %eax,%eax
  sub    %r11d,%r10d
L27:
  movss  (%rdi,%rax,4),%xmm2
  movss  (%rsi,%rax,4),%xmm5
  movaps %xmm3,%xmm7
  mov    %rax,%rdx
  movss  0x10(%rsi,%rax,4),%xmm1
  mulss  %xmm5,%xmm7
  mulss  %xmm2,%xmm5
  addss  (%r8,%rax,4),%xmm7
  mulss  %xmm1,%xmm2
  mulss  %xmm4,%xmm1
  addss  %xmm5,%xmm6
  addss  %xmm2,%xmm0
  addss  %xmm7,%xmm1
  movss  %xmm1,(%r8,%rax,4)
  add    $0x1,%rax
  cmp    %r10,%rdx
  jne    L27
L12:
  movss  0x0(%rbp),%xmm2
  movss  0x4(%rbp),%xmm1
  cmp    $0x1,%r13d
  je     L30
  subss  %xmm6,%xmm2
  subss  %xmm0,%xmm1
  movss  %xmm2,0x0(%rbp)
  movss  %xmm1,0x4(%rbp)
L72:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax)
L4:
  movss  (%rdi),%xmm2
  mulss  (%rsi),%xmm2
  pxor   %xmm5,%xmm5
  add    $0x8,%rsi
  movss  0x4(%rdi),%xmm0
  movss  -0x4(%rsi),%xmm1
  add    $0x8,%r9
  add    $0x8,%rdi
  movaps %xmm0,%xmm6
  mulss  %xmm1,%xmm6
  mulss  0xc(%rsi),%xmm0
  addss  %xmm5,%xmm2
  mulss  %xmm3,%xmm1
  addss  -0x4(%r9),%xmm1
  addss  %xmm2,%xmm6
  movss  %xmm1,-0x4(%r9)
  addss  %xmm5,%xmm0
  cmp    $0x2,%r11d
  jle    L48
  mov    $0x2,%edx
  jmp    L113
  nopw   0x0(%rax,%rax,1)
L30:
  addss  %xmm6,%xmm2
  addss  %xmm0,%xmm1
  movss  %xmm2,0x0(%rbp)
  movss  %xmm1,0x4(%rbp)
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax,%rax,1)
L67:
  movaps %xmm0,%xmm6
  xor    %r11d,%r11d
  jmp    L47
  nopl   0x0(%rax,%rax,1)
L132:
  mov    %r9,%r8
  jmp    L9
L48:
  mov    $0x2,%r11d
  jmp    L61
```

tool_calls: name: disassemble_data
arguments: data_label: D5
data_type: byte[4]

tool: {"D5": [0x00,0x00,0x00,0x80]}

assistant: ```c
void kernel_ssymv_2_lib4(int kmax, int kna, float *A, int sda, float *x_n, float *y_n, float *x_t, float *y_t, int tri, int alg) {
  if (kmax <= 0)
    return;
  const int bs = 4;
  int k;
  float a_00, a_01, x_n_0, x_n_1, y_n_0, x_t_0, y_t_0, y_t_1;
  if (alg == 1) {
    x_n_0 = x_n[0];
    x_n_1 = x_n[1];
  } else {
    x_n_0 = -x_n[0];
    x_n_1 = -x_n[1];
  }
  y_t_0 = 0;
  y_t_1 = 0;
  k = 0;
  if (tri == 1) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    y_t_0 += a_00 * x_t_0;
    y_n[0] = y_n_0;
    y_n_0 = y_n[1];
    x_t_0 = x_t[1];
    a_00 = A[1 + bs * 0];
    a_01 = A[1 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_t_1 += a_01 * x_t_0;
    y_n[1] = y_n_0;
    A += 2;
    y_n += 2;
    x_t += 2;
    k += 2;
  }
  for (; k < kna; k++) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    a_01 = A[0 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n[0] = y_n_0;
    A += 1;
    y_n += 1;
    x_t += 1;
  }
  if (kna > 0 || tri == 1) {
    A += (sda - 1) * bs;
  }
  for (; k < kmax - 3; k += bs) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    a_01 = A[0 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n[0] = y_n_0;
    y_n_0 = y_n[1];
    x_t_0 = x_t[1];
    a_00 = A[1 + bs * 0];
    a_01 = A[1 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n[1] = y_n_0;
    y_n_0 = y_n[2];
    x_t_0 = x_t[2];
    a_00 = A[2 + bs * 0];
    a_01 = A[2 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n[2] = y_n_0;
    y_n_0 = y_n[3];
    x_t_0 = x_t[3];
    a_00 = A[3 + bs * 0];
    a_01 = A[3 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n[3] = y_n_0;
    A += sda * bs;
    y_n += 4;
    x_t += 4;
  }
  for (; k < kmax; k++) {
    y_n_0 = y_n[0];
    x_t_0 = x_t[0];
    a_00 = A[0 + bs * 0];
    a_01 = A[0 + bs * 1];
    y_n_0 += a_00 * x_n_0;
    y_t_0 += a_00 * x_t_0;
    y_n_0 += a_01 * x_n_1;
    y_t_1 += a_01 * x_t_0;
    y_n[0] = y_n_0;
    A += 1;
    y_n += 1;
    x_t += 1;
  }
  if (alg == 1) {
    y_t[0] += y_t_0;
    y_t[1] += y_t_1;
  } else {
    y_t[0] -= y_t_0;
    y_t[1] -= y_t_1;
  }
}
```