user: What is the c source code of the assembly code below?

```asm
<qrupdate>:
  endbr64
  push   %r15
  mov    %rdx,%r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  push   %r12
  mov    %esi,%r12d
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rcx,0x20(%rsp)
  cmp    $0x1,%esi
  jle    L45
  lea    0x1(%rsi),%ecx
  mov    $0x1,%eax
  pxor   %xmm2,%xmm2
  movslq %ecx,%rcx
  shl    $0x3,%rcx
  lea    (%rdi,%rcx,1),%rdx
L27:
  add    $0x1,%eax
  movq   $0x0,-0x8(%rdx)
  add    %rcx,%rdx
  cmp    %eax,%r12d
  jne    L27
  lea    -0x1(%r12),%ebp
  movslq %ebp,%rax
  ucomisd (%r15,%rax,8),%xmm2
  lea    0x0(,%rax,8),%rdi
  mov    %rdi,0x28(%rsp)
  jp     L18
  jne    L18
  lea    -0x2(%r12),%eax
  xor    %edi,%edi
  cltq
  nopl   0x0(%rax,%rax,1)
L109:
  ucomisd (%r15,%rax,8),%xmm2
  mov    %eax,%ebp
  setnp  %dl
  cmovne %edi,%edx
  test   %eax,%eax
  setg   %cl
  sub    $0x1,%rax
  test   %cl,%dl
  jne    L109
L51:
  mov    %ebp,%r13d
  pxor   %xmm2,%xmm2
  sub    $0x1,%r13d
  js     L100
L146:
  movq   D26(%rip),%xmm5
  lea    0x1(%r13),%eax
  cltq
  lea    (%r15,%rax,8),%rbx
  movaps %xmm5,(%rsp)
  jmp    L49
  nopl   0x0(%rax)
L143:
  movsd  (%rbx),%xmm0
  andpd  D9(%rip),%xmm0
L132:
  sub    $0x1,%r13d
  movsd  %xmm0,-0x8(%rbx)
  sub    $0x8,%rbx
  cmp    $0xffffffff,%r13d
  je     L100
L49:
  movsd  (%rbx),%xmm1
  movsd  -0x8(%rbx),%xmm0
  mov    %r13d,%edx
  mov    %r12d,%esi
  xorpd  (%rsp),%xmm1
  mov    %r14,%rdi
  movsd  %xmm2,0x18(%rsp)
  call   <jrotate@plt>
  movsd  -0x8(%rbx),%xmm0
  movsd  0x18(%rsp),%xmm2
  ucomisd %xmm2,%xmm0
  jp     L82
  je     L143
L82:
  movsd  (%rbx),%xmm1
  mulsd  %xmm0,%xmm0
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm2
  ja     L11
  sqrtsd %xmm0,%xmm0
  jmp    L132
  nopl   0x0(%rax)
L100:
  movslq %r12d,%rdx
  xor    %eax,%eax
  test   %r12d,%r12d
  jle    L80
  nopw   0x0(%rax,%rax,1)
L24:
  mov    0x20(%rsp),%rsi
  movsd  (%r15),%xmm0
  mulsd  (%rsi,%rax,8),%xmm0
  addsd  (%r14,%rax,8),%xmm0
  movsd  %xmm0,(%r14,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L24
L80:
  test   %ebp,%ebp
  jle    L113
  mov    0x28(%rsp),%rax
  movq   D26(%rip),%xmm6
  movslq %r12d,%r15
  xor    %ebx,%ebx
  neg    %r15
  lea    0x10(%rax),%rdi
  lea    0x8(%r14,%rax,1),%r13
  shl    $0x3,%r15
  movaps %xmm6,(%rsp)
  mov    %rdi,0x18(%rsp)
  nopl   0x0(%rax,%rax,1)
L117:
  movsd  0x0(%r13),%xmm1
  xorpd  (%rsp),%xmm1
  mov    %ebx,%edx
  mov    %r12d,%esi
  movsd  0x0(%r13,%r15,1),%xmm0
  mov    %r14,%rdi
  add    $0x1,%ebx
  call   <jrotate@plt>
  add    0x18(%rsp),%r13
  cmp    %ebx,%ebp
  jne    L117
L113:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L45:
  lea    -0x1(%rsi),%ebp
  movslq %ebp,%rax
  shl    $0x3,%rax
  mov    %rax,0x28(%rsp)
  jmp    L51
  nopl   0x0(%rax)
L18:
  lea    -0x2(%r12),%r13d
  jmp    L146
L11:
  movsd  %xmm2,0x18(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm2
  jmp    L132
```

tool_calls: name: read_data
arguments: tag: D26
data_type: byte[8]
name: read_data
arguments: tag: D9
data_type: qword[2]

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

tool: {"D9": [0x7fffffffffffffff,0x0000000000000000]}

assistant: ```c
void qrupdate(double *R, int N, double *u, double *v) {
  int i, j, k, l;
  for (i = 1; i < N; ++i) {
    R[i * N + i - 1] = 0;
  }
  k = N - 1;
  while (u[k] == 0 && k > 0) {
    k--;
  }
  for (i = k - 1; i >= 0; --i) {
    jrotate(R, N, u[i], -u[i + 1], i);
    if (u[i] == 0) {
      u[i] = fabs(u[i + 1]);
    } else {
      u[i] = sqrt(u[i] * u[i] + u[i + 1] * u[i + 1]);
    }
  }
  for (i = 0; i < N; ++i) {
    R[i] += u[0] * v[i];
  }
  for (i = 0; i < k; ++i) {
    j = i * N;
    l = j + N;
    jrotate(R, N, R[j + i], -R[l + i], i);
  }
}
```