user: What would the C source look like for the assembly code provided here:

```asm
<solvps>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,0x18(%rsp)
  test   %edx,%edx
  jle    L57
  mov    %rdi,%rbp
  pxor   %xmm3,%xmm3
  movslq %edx,%rax
  xor    %r9d,%r9d
  movsd  0x0(%rbp),%xmm0
  mov    %rax,0x20(%rsp)
  mov    %edx,%r14d
  mov    %rdi,%rbx
  lea    0x0(,%rax,8),%r12
  mov    %rsi,0x28(%rsp)
  mov    $0xffffffffffffffff,%r15
  mov    %r9d,%r13d
  comisd %xmm0,%xmm3
  lea    0x8(%r12),%rax
  mov    %rax,0x10(%rsp)
  jae    L64
  nopw   0x0(%rax,%rax,1)
L100:
  ucomisd %xmm0,%xmm3
  ja     L77
  sqrtsd %xmm0,%xmm0
L105:
  lea    0x1(%r13),%r8d
  movsd  %xmm0,0x0(%rbp)
  lea    0x0(%rbp,%r12,1),%rdi
  mov    %r8d,0x4(%rsp)
  cmp    %r8d,%r14d
  je     L46
  mov    %r15,%r9
  lea    (%rbx,%r12,1),%rax
  shr    $0x3,%r9
  mov    %rax,0x8(%rsp)
  add    $0x1,%r9
  mov    %r9,%r11
  mov    %r9,%rsi
  and    $0xfffffffffffffffe,%r11
  shr    %rsi
  lea    0x0(,%r11,8),%r10
  shl    $0x4,%rsi
  lea    (%rbx,%r10,1),%r13
  nopl   0x0(%rax,%rax,1)
L228:
  mov    %rax,%rdx
  cmp    %rbp,%rbx
  jae    L12
  cmp    $0x7,%r15
  jbe    L73
  xor    %edx,%edx
  movapd %xmm3,%xmm0
  nopl   0x0(%rax)
L101:
  movupd (%rbx,%rdx,1),%xmm2
  movupd (%rax,%rdx,1),%xmm4
  add    $0x10,%rdx
  mulpd  %xmm4,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %rdx,%rsi
  jne    L101
  lea    (%rax,%r10,1),%rdx
  mov    %r13,%rcx
  cmp    %r11,%r9
  je     L187
L95:
  movsd  (%rcx),%xmm1
  mulsd  (%rdx),%xmm1
  addsd  %xmm1,%xmm0
L187:
  movsd  (%rdi),%xmm1
  add    $0x1,%r8d
  add    %r12,%rax
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdi)
  divsd  0x0(%rbp),%xmm1
  movsd  %xmm1,(%rdi)
  add    %r12,%rdi
  cmp    %r8d,%r14d
  jne    L228
  mov    0x8(%rsp),%rbx
  add    0x10(%rsp),%rbp
  movsd  0x0(%rbp),%xmm1
  cmp    %rbx,%rbp
  jbe    L226
  mov    %rbx,%rax
  cs nopw 0x0(%rax,%rax,1)
L35:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  mulsd  %xmm0,%xmm0
  subsd  %xmm0,%xmm1
  movsd  %xmm1,0x0(%rbp)
  cmp    %rbp,%rax
  jb     L35
L226:
  movsd  0x0(%rbp),%xmm0
  add    $0x8,%r15
  mov    0x4(%rsp),%r13d
  comisd %xmm0,%xmm3
  jb     L100
L64:
  add    $0x38,%rsp
  mov    $0xffffffff,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  cs nopw 0x0(%rax,%rax,1)
L12:
  movapd %xmm3,%xmm0
  jmp    L187
L73:
  mov    %rbx,%rcx
  movapd %xmm3,%xmm0
  jmp    L95
L46:
  mov    0x18(%rsp),%rax
  mov    0x10(%rsp),%r8
  mov    %r13d,%r9d
  mov    0x28(%rsp),%r13
  mov    %rax,%rdi
  movsd  0x0(%r13),%xmm0
  add    %r8,%rdi
  divsd  (%rax),%xmm0
  movsd  %xmm0,0x0(%r13)
  cmp    $0x1,%r14d
  je     L61
  mov    0x20(%rsp),%rcx
  lea    0x8(%r13),%rdx
  lea    (%rax,%rcx,8),%rsi
  mov    $0x1,%ecx
L193:
  movsd  (%rdx),%xmm1
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L117:
  movsd  (%rsi,%rax,8),%xmm0
  mulsd  0x0(%r13,%rax,8),%xmm0
  add    $0x1,%rax
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdx)
  cmp    %eax,%ecx
  jg     L117
  divsd  (%rdi),%xmm1
  lea    0x1(%rcx),%eax
  movsd  %xmm1,(%rdx)
  add    %r8,%rdi
  add    $0x8,%rdx
  add    %r12,%rsi
  cmp    %ecx,%r9d
  je     L61
  mov    %eax,%ecx
  jmp    L193
L57:
  add    $0x38,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L61:
  mov    0x4(%rsp),%ecx
  mov    0x18(%rsp),%rdi
  mov    $0xfffffffffffffff8,%r8
  sub    %r12,%r8
  mov    %ecx,%eax
  imul   %ecx,%eax
  cltq
  lea    -0x8(%rdi,%rax,8),%rdi
  movslq %r9d,%rax
  lea    0x0(%r13,%rax,8),%rsi
L129:
  movsd  (%rsi),%xmm1
  lea    (%rdi,%r12,1),%rdx
  cmp    %ecx,%r14d
  jle    L26
  movslq %ecx,%rax
  nopw   0x0(%rax,%rax,1)
L22:
  movsd  0x0(%r13,%rax,8),%xmm0
  mulsd  (%rdx),%xmm0
  add    $0x1,%rax
  add    %r12,%rdx
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rsi)
  cmp    %eax,%r14d
  jg     L22
L26:
  divsd  (%rdi),%xmm1
  sub    $0x8,%rsi
  movsd  %xmm1,0x8(%rsi)
  add    %r8,%rdi
  lea    -0x1(%rcx),%eax
  sub    $0x2,%ecx
  js     L57
  mov    %eax,%ecx
  jmp    L129
L77:
  call   <sqrt@plt>
  pxor   %xmm3,%xmm3
  jmp    L105
```

assistant: ```c
int solvps(double *a, double *b, int n) {
  double *p, *q, *r, *s, t;
  int j, k;
  for (j = 0, p = a; j < n; ++j, p += n + 1) {
    for (q = a + j * n; q < p; ++q)
      *p -= *q * *q;
    if (*p <= 0.)
      return -1;
    *p = sqrt(*p);
    for (k = j + 1, q = p + n; k < n; ++k, q += n) {
      for (r = a + j * n, s = a + k * n, t = 0.; r < p;)
        t += *r++ * *s++;
      *q -= t;
      *q /= *p;
    }
  }
  for (j = 0, p = a; j < n; ++j, p += n + 1) {
    for (k = 0, q = a + j * n; k < j;)
      b[j] -= b[k++] * *q++;
    b[j] /= *p;
  }
  for (j = n - 1, p = a + n * n - 1; j >= 0; --j, p -= n + 1) {
    for (k = j + 1, q = p + n; k < n; q += n)
      b[j] -= b[k++] * *q;
    b[j] /= *p;
  }
  return 0;
}
```