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

```asm
<kernel_strsv_lt_inv_3_lib4>:
  endbr64
  push   %r12
  mov    %rsi,%rax
  mov    %edi,%r10d
  mov    %edx,%esi
  push   %rbp
  mov    %rcx,%rdx
  push   %rbx
  mov    0x20(%rsp),%rcx
  cmp    $0x4,%edi
  jle    L27
  movss  0xc(%r8),%xmm3
  movss  0xc(%rax),%xmm1
  pxor   %xmm0,%xmm0
  lea    0x0(,%rsi,4),%ebx
  movss  0x1c(%rax),%xmm4
  movaps %xmm0,%xmm2
  movslq %ebx,%rbx
  lea    0x10(%r8),%r11
  mulss  %xmm3,%xmm1
  shl    $0x2,%rbx
  mulss  %xmm3,%xmm4
  lea    (%rax,%rbx,1),%rdi
  mulss  0x2c(%rax),%xmm3
  subss  %xmm1,%xmm2
  movaps %xmm0,%xmm1
  subss  %xmm4,%xmm1
  subss  %xmm3,%xmm0
  cmp    $0x7,%r10d
  jle    L55
  lea    -0x8(%r10),%r12d
  shr    $0x2,%r12d
  mov    %r12d,%ebp
  lea    0x2(%rbp),%rsi
  shl    $0x4,%rsi
  add    %rsi,%r8
  mov    %rdi,%rsi
  nopl   (%rax)
L39:
  movss  (%r11),%xmm6
  movss  (%rsi),%xmm7
  add    $0x10,%r11
  movss  -0xc(%r11),%xmm5
  movss  -0x8(%r11),%xmm4
  mulss  %xmm6,%xmm7
  movss  -0x4(%r11),%xmm3
  subss  %xmm7,%xmm2
  movss  0x10(%rsi),%xmm7
  mulss  %xmm6,%xmm7
  mulss  0x20(%rsi),%xmm6
  subss  %xmm7,%xmm1
  subss  %xmm6,%xmm0
  movss  0x4(%rsi),%xmm6
  mulss  %xmm5,%xmm6
  subss  %xmm6,%xmm2
  movss  0x14(%rsi),%xmm6
  mulss  %xmm5,%xmm6
  mulss  0x24(%rsi),%xmm5
  subss  %xmm6,%xmm1
  subss  %xmm5,%xmm0
  movss  0x8(%rsi),%xmm5
  mulss  %xmm4,%xmm5
  subss  %xmm5,%xmm2
  movss  0x18(%rsi),%xmm5
  mulss  %xmm4,%xmm5
  mulss  0x28(%rsi),%xmm4
  subss  %xmm5,%xmm1
  subss  %xmm4,%xmm0
  movss  0xc(%rsi),%xmm4
  mulss  %xmm3,%xmm4
  subss  %xmm4,%xmm2
  movss  0x1c(%rsi),%xmm4
  mulss  %xmm3,%xmm4
  mulss  0x2c(%rsi),%xmm3
  add    %rbx,%rsi
  subss  %xmm4,%xmm1
  subss  %xmm3,%xmm0
  cmp    %r8,%r11
  jne    L39
  add    $0x1,%rbp
  lea    0x8(,%r12,4),%esi
  imul   %rbx,%rbp
  add    %rbp,%rdi
  cmp    %esi,%r10d
  jg     L33
L38:
  addss  0x8(%r9),%xmm0
  mulss  0x8(%rdx),%xmm0
  addss  (%r9),%xmm2
  addss  0x4(%r9),%xmm1
  movss  %xmm0,0x8(%rcx)
  movss  0x8(%rax),%xmm3
  mulss  %xmm0,%xmm3
  mulss  0x18(%rax),%xmm0
  subss  %xmm3,%xmm2
  subss  %xmm0,%xmm1
  mulss  0x4(%rdx),%xmm1
  movss  %xmm1,0x4(%rcx)
  mulss  0x4(%rax),%xmm1
  subss  %xmm1,%xmm2
  mulss  (%rdx),%xmm2
  movss  %xmm2,(%rcx)
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
  nopl   0x0(%rax,%rax,1)
L27:
  pxor   %xmm0,%xmm0
  lea    0xc(%rax),%rdi
  add    $0x4,%r8
  mov    $0x3,%esi
  movaps %xmm0,%xmm1
  movaps %xmm0,%xmm2
  jmp    L54
  nopl   0x0(%rax)
L33:
  movss  (%r8),%xmm3
  movss  (%rdi),%xmm4
  add    $0x4,%r8
  add    $0x1,%esi
  add    $0x4,%rdi
  mulss  %xmm3,%xmm4
  subss  %xmm4,%xmm2
  movss  0xc(%rdi),%xmm4
  mulss  %xmm3,%xmm4
  mulss  0x1c(%rdi),%xmm3
  subss  %xmm4,%xmm1
  subss  %xmm3,%xmm0
L54:
  cmp    %esi,%r10d
  jg     L33
  jmp    L38
L55:
  mov    %r11,%r8
  mov    $0x4,%esi
  jmp    L33
```

assistant: ```c
void kernel_strsv_lt_inv_3_lib4(int kmax, float *A, int sda, float *inv_diag_A, float *x, float *y, float *z) {
  const int bs = 4;
  int k;
  float *tA, *tx;
  tA = A;
  tx = x;
  float x_0, x_1, x_2, x_3, y_0 = 0, y_1 = 0, y_2 = 0;
  k = 3;
  if (kmax > 4) {
    x_3 = x[3];
    y_0 -= A[3 + bs * 0] * x_3;
    y_1 -= A[3 + bs * 1] * x_3;
    y_2 -= A[3 + bs * 2] * x_3;
    k = 4;
    A += 4 + (sda - 1) * bs;
    x += 4;
    for (; k < kmax - 3; k += 4) {
      x_0 = x[0];
      x_1 = x[1];
      x_2 = x[2];
      x_3 = x[3];
      y_0 -= A[0 + bs * 0] * x_0;
      y_1 -= A[0 + bs * 1] * x_0;
      y_2 -= A[0 + bs * 2] * x_0;
      y_0 -= A[1 + bs * 0] * x_1;
      y_1 -= A[1 + bs * 1] * x_1;
      y_2 -= A[1 + bs * 2] * x_1;
      y_0 -= A[2 + bs * 0] * x_2;
      y_1 -= A[2 + bs * 1] * x_2;
      y_2 -= A[2 + bs * 2] * x_2;
      y_0 -= A[3 + bs * 0] * x_3;
      y_1 -= A[3 + bs * 1] * x_3;
      y_2 -= A[3 + bs * 2] * x_3;
      A += sda * bs;
      x += 4;
    }
  } else {
    A += 3;
    x += 1;
  }
  for (; k < kmax; k++) {
    x_0 = x[0];
    y_0 -= A[0 + bs * 0] * x_0;
    y_1 -= A[0 + bs * 1] * x_0;
    y_2 -= A[0 + bs * 2] * x_0;
    A += 1;
    x += 1;
  }
  y_0 = y[0] + y_0;
  y_1 = y[1] + y_1;
  y_2 = y[2] + y_2;
  A = tA;
  x = tx;
  y_2 *= inv_diag_A[2];
  z[2] = y_2;
  y_0 -= A[2 + bs * 0] * y_2;
  y_1 -= A[2 + bs * 1] * y_2;
  y_1 *= inv_diag_A[1];
  z[1] = y_1;
  y_0 -= A[1 + bs * 0] * y_1;
  y_0 *= inv_diag_A[0];
  z[0] = y_0;
}
```