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

```asm
<ep21R>:
  endbr64
  push   %r15
  imul   %r8d,%r9d
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rsi,0x18(%rsp)
  mov    %rcx,0x8(%rsp)
  test   %r8d,%r8d
  jle    L461
  mov    %rdi,%r15
  mov    %rdx,%rbp
  mov    %r8d,%r11d
  mov    %r9d,%r13d
  test   %r9d,%r9d
  jle    L85
  cmp    %r9d,%r8d
  mov    %r9d,%eax
  mov    %r9d,%edx
  mov    %rdi,%r14
  cmovle %r8d,%eax
  lea    0x8(%rbp),%r9
  mov    %rbp,%r8
  xor    %r12d,%r12d
  movq   D13(%rip),%xmm4
  pxor   %xmm3,%xmm3
  mov    %rsi,%r15
  mov    %rdi,%rcx
  mov    %eax,0x14(%rsp)
  xchg   %ax,%ax
L34:
  mov    %r12d,%ebx
  mov    %r12d,%esi
  cmp    $0x1,%r11d
  jne    L550
  cmp    %r12d,%r13d
  mov    %r11d,%edi
  cmovg  %edx,%edi
  jle    L624
  cmp    $0x1,%edx
  je     L624
  mov    %edi,%eax
  movapd %xmm3,%xmm0
  shr    %eax
  sub    $0x1,%eax
  shl    $0x4,%rax
  lea    0x10(%r8,%rax,1),%r10
  mov    %r8,%rax
  nop
L38:
  movupd (%rax),%xmm2
  add    $0x10,%rax
  mulpd  %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %rax,%r10
  jne    L38
  mov    %edi,%r10d
  and    $0xfffffffe,%r10d
  lea    (%r10,%rsi,1),%eax
  cmp    %r10d,%edi
  je     L84
L396:
  cltq
  movsd  0x0(%rbp,%rax,8),%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
L84:
  ucomisd %xmm0,%xmm3
  ja     L581
  sqrtsd %xmm0,%xmm0
L630:
  mov    0x8(%rsp),%rax
  movsd  (%rax,%r12,8),%xmm1
  movapd %xmm1,%xmm2
  andpd  %xmm4,%xmm2
  comisd %xmm2,%xmm0
  ja     L398
  comisd %xmm0,%xmm1
  jb     L91
  movsd  %xmm1,(%r15)
  cmp    %r12d,%r13d
  jle    L125
  cmp    $0x1,%r11d
  jne    L495
  lea    -0x1(%rdx),%r10d
  cmp    $0x2,%r10d
  jbe    L94
  cmp    %r9,%r14
  je     L94
  mov    %edx,%edi
  xor    %eax,%eax
  shr    %edi
  shl    $0x4,%rdi
  nopl   0x0(%rax,%rax,1)
L543:
  movupd (%r8,%rax,1),%xmm7
  movups %xmm7,(%r14,%rax,1)
  add    $0x10,%rax
  cmp    %rdi,%rax
  jne    L543
  mov    %edx,%eax
  and    $0xfffffffe,%eax
  add    %eax,%esi
  cmp    %eax,%edx
  je     L585
  movslq %esi,%rsi
  lea    0x1(%rbx),%eax
  add    $0x1,%r12
  add    $0x8,%r14
  movsd  0x0(%rbp,%rsi,8),%xmm0
  movsd  %xmm0,(%rcx,%rsi,8)
  lea    0x8(%r15),%rsi
L302:
  mov    %rsi,%r15
  mov    %r10d,%edx
  add    $0x8,%r8
  add    $0x8,%r9
  cmp    %r12d,0x14(%rsp)
  jg     L34
  mov    %rcx,%r15
  cmp    %eax,%r11d
  jle    L461
  cmp    $0x1,%r11d
  jne    L580
L457:
  mov    0x8(%rsp),%rbx
  movslq %eax,%rdx
  mov    0x18(%rsp),%rcx
  lea    0x0(,%rdx,8),%rdi
  movsd  (%rbx,%rdx,8),%xmm0
  add    %rdi,%rcx
  comisd %xmm3,%xmm0
  jb     L42
  movsd  %xmm0,(%rcx)
  cmp    %eax,%r13d
  jle    L461
  mov    %r13d,%r8d
  lea    (%r15,%rdi,1),%rsi
  lea    0x8(%rbp,%rdi,1),%rcx
  sub    %eax,%r8d
  cmp    %rcx,%rsi
  je     L152
  lea    -0x1(%r8),%ecx
  cmp    $0x2,%ecx
  jbe    L152
  mov    %r8d,%ecx
  add    %rbp,%rdi
  xor    %edx,%edx
  shr    %ecx
  shl    $0x4,%rcx
  nopw   0x0(%rax,%rax,1)
L295:
  movupd (%rdi,%rdx,1),%xmm6
  movups %xmm6,(%rsi,%rdx,1)
  add    $0x10,%rdx
  cmp    %rcx,%rdx
  jne    L295
  mov    %r8d,%edx
  and    $0xfffffffe,%edx
  add    %edx,%eax
  cmp    %edx,%r8d
  je     L461
  cltq
  movsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,(%r15,%rax,8)
L461:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L152:
  movsd  0x0(%rbp,%rdx,8),%xmm0
  movsd  %xmm0,(%r15,%rdx,8)
  add    $0x1,%rdx
  cmp    %edx,%r13d
  jg     L152
  jmp    L461
  nopw   0x0(%rax,%rax,1)
L91:
  movq   $0x0,(%r15)
  cmp    %r12d,%r13d
  jle    L125
  cmp    $0x1,%r11d
  jne    L143
  lea    -0x1(%rdx),%eax
  xor    %esi,%esi
  mov    %r14,%rdi
  mov    %r11d,0x40(%rsp)
  mov    %rcx,0x38(%rsp)
  add    $0x1,%r12
  add    $0x8,%r14
  lea    0x8(,%rax,8),%rdx
  mov    %r9,0x30(%rsp)
  mov    %r8,0x28(%rsp)
  movsd  %xmm3,0x20(%rsp)
  mov    %eax,0x48(%rsp)
  call   <memset@plt>
  mov    0x48(%rsp),%r10d
  movsd  0x20(%rsp),%xmm3
  lea    0x1(%rbx),%eax
  mov    0x28(%rsp),%r8
  mov    0x30(%rsp),%r9
  lea    0x8(%r15),%rsi
  movq   D13(%rip),%xmm4
  mov    0x38(%rsp),%rcx
  mov    0x40(%rsp),%r11d
  jmp    L302
  nopl   0x0(%rax,%rax,1)
L42:
  movq   $0x0,(%rcx)
  cmp    %eax,%r13d
  jle    L461
  sub    $0x1,%r13d
  add    $0x58,%rsp
  add    %r15,%rdi
  xor    %esi,%esi
  sub    %eax,%r13d
  pop    %rbx
  pop    %rbp
  lea    0x8(,%r13,8),%rdx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1060 <memset@plt>
L143:
  movslq %r11d,%rsi
  mov    %r12,%rax
  nopl   (%rax)
L465:
  movq   $0x0,(%rcx,%rax,8)
  add    %rsi,%rax
  cmp    %eax,%r13d
  jg     L465
L125:
  lea    0x1(%rbx),%eax
L271:
  add    $0x1,%r12
  lea    0x8(%r15),%rsi
  lea    -0x1(%rdx),%r10d
  add    $0x8,%r14
  jmp    L302
  nopl   0x0(%rax,%rax,1)
L398:
  addsd  %xmm0,%xmm1
  mulsd  D30(%rip),%xmm1
  movsd  %xmm1,(%r15)
  cmp    %r12d,%r13d
  jle    L125
  cmp    $0x1,%r11d
  jne    L493
  divsd  %xmm0,%xmm1
  lea    0x1(%rbx),%eax
  mulsd  (%r8),%xmm1
  movsd  %xmm1,(%r14)
  cmp    %eax,%r13d
  jle    L271
  add    $0x1,%r12
  lea    -0x1(%rdx),%esi
  add    $0x8,%r14
  mov    %rsi,%r10
  add    %r12,%rsi
  mov    %r12,%rdi
  lea    (%rcx,%rsi,8),%rsi
  cmp    %rsi,%r15
  lea    0x8(%r15),%rsi
  setae  %bl
  cmp    %rsi,%r14
  setae  0x48(%rsp)
  or     0x48(%rsp),%bl
  je     L22
  lea    0x10(%r8),%rbx
  cmp    %rbx,%r14
  setne  %bl
  cmp    $0x2,%edx
  setne  %dl
  test   %dl,%bl
  je     L22
  movsd  (%r15),%xmm2
  mov    %r10d,%edi
  xor    %edx,%edx
  shr    %edi
  divsd  %xmm0,%xmm2
  shl    $0x4,%rdi
  unpcklpd %xmm2,%xmm2
  nopw   0x0(%rax,%rax,1)
L268:
  movupd (%r9,%rdx,1),%xmm1
  mulpd  %xmm2,%xmm1
  movups %xmm1,(%r14,%rdx,1)
  add    $0x10,%rdx
  cmp    %rdx,%rdi
  jne    L268
  mov    %r10d,%edi
  and    $0xfffffffe,%edi
  lea    (%rdi,%rax,1),%edx
  cmp    %edi,%r10d
  je     L302
  movsd  (%r15),%xmm1
  movslq %edx,%rdx
  divsd  %xmm0,%xmm1
  mulsd  0x0(%rbp,%rdx,8),%xmm1
  movsd  %xmm1,(%rcx,%rdx,8)
  jmp    L302
  xchg   %ax,%ax
L22:
  movsd  (%r15),%xmm1
  divsd  %xmm0,%xmm1
  mulsd  0x0(%rbp,%rdi,8),%xmm1
  movsd  %xmm1,(%rcx,%rdi,8)
  add    $0x1,%rdi
  cmp    %edi,%r13d
  jg     L22
  jmp    L302
  nopw   0x0(%rax,%rax,1)
L550:
  movslq %r11d,%rdi
  mov    %r12,%rax
  movapd %xmm3,%xmm0
  nopw   0x0(%rax,%rax,1)
L424:
  movsd  0x0(%rbp,%rax,8),%xmm1
  add    %rdi,%rax
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  cmp    %eax,%r13d
  jg     L424
  jmp    L84
L94:
  mov    %r12,%rax
  cs nopw 0x0(%rax,%rax,1)
L488:
  movsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r13d
  jg     L488
L585:
  lea    0x1(%rbx),%eax
  add    $0x1,%r12
  lea    0x8(%r15),%rsi
  add    $0x8,%r14
  jmp    L302
  nopl   0x0(%rax,%rax,1)
L85:
  xor    %eax,%eax
  pxor   %xmm3,%xmm3
  cmp    $0x1,%r11d
  je     L457
L580:
  mov    0x8(%rsp),%rsi
  cltq
  mov    0x18(%rsp),%rdi
  movslq %r11d,%rcx
  movsd  (%rsi,%rax,8),%xmm0
  comisd %xmm3,%xmm0
  jb     L573
  nopw   0x0(%rax,%rax,1)
L450:
  movsd  %xmm0,(%rdi,%rax,8)
  cmp    %eax,%r13d
  jle    L480
  mov    %rax,%rdx
  nopl   (%rax)
L67:
  movsd  0x0(%rbp,%rdx,8),%xmm0
  movsd  %xmm0,(%r15,%rdx,8)
  add    %rcx,%rdx
  cmp    %edx,%r13d
  jg     L67
L480:
  add    $0x1,%rax
  cmp    %eax,%r11d
  jle    L461
L406:
  movsd  (%rsi,%rax,8),%xmm0
  comisd %xmm3,%xmm0
  jae    L450
L573:
  movq   $0x0,(%rdi,%rax,8)
  cmp    %eax,%r13d
  jle    L480
  mov    %rax,%rdx
  nopl   0x0(%rax)
L76:
  movq   $0x0,(%r15,%rdx,8)
  add    %rcx,%rdx
  cmp    %edx,%r13d
  jg     L76
  add    $0x1,%rax
  cmp    %eax,%r11d
  jg     L406
  jmp    L461
  xchg   %ax,%ax
L624:
  mov    %esi,%eax
  movapd %xmm3,%xmm0
  jmp    L396
L495:
  movslq %r11d,%rsi
  mov    %r12,%rax
  nopl   0x0(%rax)
L501:
  movsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    %rsi,%rax
  cmp    %eax,%r13d
  jg     L501
  jmp    L125
L493:
  divsd  %xmm0,%xmm1
  lea    (%r11,%rbx,1),%eax
  mulsd  (%r8),%xmm1
  movsd  %xmm1,(%r14)
  cmp    %eax,%r13d
  jle    L125
  movslq %r11d,%rsi
  cltq
L55:
  movsd  (%r15),%xmm1
  divsd  %xmm0,%xmm1
  mulsd  0x0(%rbp,%rax,8),%xmm1
  movsd  %xmm1,(%rcx,%rax,8)
  add    %rsi,%rax
  cmp    %eax,%r13d
  jg     L55
  jmp    L125
L581:
  mov    %r11d,0x4c(%rsp)
  mov    %rcx,0x40(%rsp)
  mov    %r9,0x38(%rsp)
  mov    %r8,0x30(%rsp)
  mov    %edx,0x28(%rsp)
  mov    %esi,0x48(%rsp)
  movsd  %xmm3,0x20(%rsp)
  call   <sqrt@plt>
  mov    0x4c(%rsp),%r11d
  mov    0x40(%rsp),%rcx
  mov    0x38(%rsp),%r9
  mov    0x30(%rsp),%r8
  mov    0x28(%rsp),%edx
  movq   D13(%rip),%xmm4
  movsd  0x20(%rsp),%xmm3
  mov    0x48(%rsp),%esi
  jmp    L630
```

tool_calls: name: read_data
arguments: data_tag: D13
type_guessed: qword
name: read_data
arguments: data_tag: D30
type_guessed: byte[8]

tool: D13: 9223372036854775807

tool: {"D30": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

assistant: ```c
void ep21R(double *x, double *t, double *u, double *v, int n, int k) {
  int i, j, tn = n * k;
  double temp;
  for (j = 0; j < n; j++) {
    temp = 0;
    for (i = j; i < tn; i += n)
      temp += u[i] * u[i];
    temp = sqrt(temp);
    if (temp > fabs(v[j])) {
      t[j] = (temp + v[j]) / 2;
      for (i = j; i < tn; i += n)
        x[i] = t[j] / temp * u[i];
    } else if (temp <= v[j]) {
      t[j] = v[j];
      for (i = j; i < tn; i += n)
        x[i] = u[i];
    } else {
      t[j] = 0;
      for (i = j; i < tn; i += n)
        x[i] = 0;
    }
  }
}
```