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

```asm
<xorgqr>:
  endbr64
  push   %r15
  mov    %rdx,%r15
  push   %r14
  push   %r13
  push   %r12
  mov    %r8,%r12
  push   %rbp
  mov    %ecx,%ebp
  push   %rbx
  sub    $0x88,%rsp
  mov    %edi,0x14(%rsp)
  mov    %esi,0x20(%rsp)
  mov    %r9,0x50(%rsp)
  cmp    %esi,%edi
  jle    L249
  movslq 0x20(%rsp),%rdx
  mov    0x14(%rsp),%eax
  movslq %ecx,%rbx
  mov    %ebp,0x24(%rsp)
  lea    0x0(,%rbx,8),%rcx
  movsd  D42(%rip),%xmm7
  mov    %r8,0x28(%rsp)
  lea    0x0(,%rax,8),%edi
  mov    %edx,%eax
  lea    0x8(%rcx),%r10
  mov    %r15,0x18(%rsp)
  imul   %ebp,%eax
  movsd  %xmm7,0x8(%rsp)
  mov    %edx,%ebp
  mov    %rcx,%r13
  cltq
  lea    (%r15,%rax,8),%rbx
  add    %rdx,%rax
  lea    (%r15,%rax,8),%r14
  mov    %rbx,%r12
  mov    %edi,%ebx
  mov    %r14,%r15
  mov    %r10,%r14
  nop
L156:
  mov    %r12,%rdi
  mov    %ebx,%edx
  xor    %esi,%esi
  add    $0x1,%ebp
  call   <memset@plt>
  movsd  0x8(%rsp),%xmm4
  add    %r13,%r12
  movsd  %xmm4,(%r15)
  add    %r14,%r15
  cmp    %ebp,0x14(%rsp)
  jne    L156
  mov    0x18(%rsp),%r15
  mov    0x24(%rsp),%ebp
  mov    0x28(%rsp),%r12
L249:
  mov    0x20(%rsp),%esi
  lea    -0x1(%rsi),%eax
  test   %esi,%esi
  jle    L114
  imul   %ebp,%eax
  mov    %ebp,%ecx
  movslq %esi,%rbx
  mov    %ebp,%r13d
  not    %ecx
  movsd  D42(%rip),%xmm6
  neg    %r13d
  mov    %r15,%r8
  mov    %ecx,0x20(%rsp)
  movslq %ecx,%rcx
  movslq %r13d,%r13
  movq   D6(%rip),%xmm2
  lea    0x0(,%rcx,8),%rdi
  movslq %eax,%rdx
  mov    %rcx,0x40(%rsp)
  shl    $0x3,%r13
  mov    %rdi,0x48(%rsp)
  lea    (%rdx,%rbx,1),%rdi
  mov    0x14(%rsp),%edx
  lea    0x0(,%rdi,8),%r11
  mov    %rdi,0x28(%rsp)
  lea    (%rsi,%rax,1),%edi
  sub    %esi,%edx
  mov    %edi,%esi
  mov    %edi,0x30(%rsp)
  mov    %r11,%r15
  lea    -0x1(%rdx),%ecx
  lea    (%rdx,%rdi,1),%edi
  add    %ebp,%esi
  movsd  %xmm6,0x8(%rsp)
  mov    %ecx,0x34(%rsp)
  imul   %ebp,%ecx
  mov    %edi,0x24(%rsp)
  mov    %esi,0x38(%rsp)
  add    %ecx,%eax
  mov    %ecx,0x3c(%rsp)
  add    %ebp,%eax
  mov    %eax,0x70(%rsp)
  movslq %ebp,%rax
  mov    %rax,0x68(%rsp)
  shl    $0x3,%rax
  mov    %rax,0x60(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L80:
  cmp    %ebx,0x14(%rsp)
  jle    L165
  movsd  0x8(%rsp),%xmm5
  pxor   %xmm1,%xmm1
  movsd  %xmm5,-0x8(%r8,%r15,1)
  movsd  -0x8(%r12,%rbx,8),%xmm0
  ucomisd %xmm1,%xmm0
  jp     L104
  jne    L104
L115:
  mov    0x30(%rsp),%eax
  add    $0x1,%eax
  cmp    %eax,0x24(%rsp)
  jl     L165
  mov    0x34(%rsp),%edx
  cltq
  add    0x28(%rsp),%rdx
  lea    -0x8(%r8,%rax,8),%rax
  lea    0x8(%r8,%rdx,8),%rdx
L175:
  movsd  -0x8(%r12,%rbx,8),%xmm0
  add    $0x8,%rax
  xorpd  %xmm2,%xmm0
  mulsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L175
L165:
  movsd  0x8(%rsp),%xmm0
  subsd  -0x8(%r12,%rbx,8),%xmm0
  sub    $0x1,%rbx
  mov    0x20(%rsp),%ecx
  mov    0x40(%rsp),%rsi
  sub    %ebp,0x24(%rsp)
  add    %ecx,0x38(%rsp)
  add    %rsi,0x28(%rsp)
  add    %ebp,0x3c(%rsp)
  addl   $0x1,0x34(%rsp)
  movsd  %xmm0,-0x8(%r8,%r15,1)
  test   %ebx,%ebx
  jle    L114
  lea    -0x1(%rbx),%eax
  xor    %esi,%esi
  mov    %r8,0x18(%rsp)
  lea    0x8(,%rax,8),%rdx
  neg    %rax
  lea    -0x10(%r15,%rax,8),%rdi
  add    %r8,%rdi
  call   <memset@plt>
  mov    0x20(%rsp),%ecx
  mov    0x18(%rsp),%r8
  add    %ecx,0x30(%rsp)
  movq   D6(%rip),%xmm2
  add    0x48(%rsp),%r15
  jmp    L80
  xchg   %ax,%ax
L104:
  mov    0x14(%rsp),%eax
  mov    %ebx,%ecx
  mov    %eax,%edi
  lea    0x1(%rax),%r14d
  movslq 0x24(%rsp),%rax
  sub    %ebx,%edi
  sub    %ebx,%r14d
  lea    (%r8,%rax,8),%rax
  nopw   0x0(%rax,%rax,1)
L14:
  ucomisd -0x8(%rax),%xmm1
  jp     L91
  jne    L91
  sub    $0x8,%rax
  sub    $0x1,%r14d
  jne    L14
L91:
  movslq 0x38(%rsp),%r11
  movslq 0x3c(%rsp),%rdx
  lea    -0x1(%r14),%eax
  add    0x70(%rsp),%ecx
  add    %r11,%rdx
  lea    (%rcx,%r14,1),%esi
  add    %rax,%rdx
  neg    %rax
  shl    $0x3,%rax
  lea    (%r8,%rdx,8),%rdx
  mov    %rax,0x18(%rsp)
  lea    -0x8(%rax),%r9
  nopl   0x0(%rax,%rax,1)
L281:
  mov    %edi,%r10d
  sub    $0x1,%edi
  cmp    %esi,%ecx
  jge    L207
  lea    (%r9,%rdx,1),%rax
  cs nopw 0x0(%rax,%rax,1)
L278:
  ucomisd (%rax),%xmm1
  jp     L60
  jne    L60
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L278
L207:
  sub    %ebp,%esi
  sub    %ebp,%ecx
  add    %r13,%rdx
  test   %edi,%edi
  jne    L281
  xor    %r10d,%r10d
  test   %r14d,%r14d
  jne    L198
  jmp    L115
  nopl   (%rax)
L60:
  mov    %r11,0x58(%rsp)
  test   %r14d,%r14d
  je     L115
  mov    0x50(%rsp),%rdi
  lea    0x0(,%r10,8),%edx
  xor    %esi,%esi
  mov    %r8,0x78(%rsp)
  mov    %r10d,0x74(%rsp)
  call   <memset@plt>
  mov    0x74(%rsp),%r10d
  mov    0x38(%rsp),%r9d
  mov    %r14d,0x74(%rsp)
  mov    0x58(%rsp),%r11
  mov    0x78(%rsp),%r8
  mov    %r15,0x58(%rsp)
  pxor   %xmm1,%xmm1
  lea    -0x1(%r10),%eax
  mov    0x28(%rsp),%rdx
  lea    (%r14,%r9,1),%ecx
  mov    %r10d,0x78(%rsp)
  imul   %ebp,%eax
  lea    -0x8(%r8,%r11,8),%rdi
  mov    0x50(%rsp),%rsi
  mov    %ecx,%r10d
  sub    %r11,%rdx
  mov    0x60(%rsp),%r14
  mov    $0x8,%r11d
  mov    0x68(%rsp),%r15
  sub    0x18(%rsp),%r11
  mov    %r13,0x18(%rsp)
  add    %r9d,%eax
  mov    %eax,%r13d
  test   %ebp,%ebp
  js     L225
L32:
  cmp    %r13d,%r9d
  jg     L126
L271:
  cmp    %r10d,%r9d
  jge    L116
  lea    (%r11,%rdi,1),%rcx
  mov    %rdi,%rax
  movapd %xmm1,%xmm3
L162:
  movsd  (%rax),%xmm0
  mulsd  (%rax,%rdx,8),%xmm0
  add    $0x8,%rax
  addsd  %xmm0,%xmm3
  cmp    %rax,%rcx
  jne    L162
L59:
  addsd  (%rsi),%xmm3
  add    %ebp,%r9d
  add    $0x8,%rsi
  add    %ebp,%r10d
  add    %r14,%rdi
  sub    %r15,%rdx
  movsd  %xmm3,-0x8(%rsi)
  test   %ebp,%ebp
  jns    L32
L225:
  cmp    %r13d,%r9d
  jge    L271
L126:
  mov    0x18(%rsp),%r13
  mov    0x58(%rsp),%r15
  mov    0x74(%rsp),%r14d
  mov    0x78(%rsp),%r10d
  movsd  -0x8(%r12,%rbx,8),%xmm0
L198:
  sub    $0x8,%rsp
  xorpd  D6(%rip),%xmm0
  mov    %r10d,%esi
  mov    %r14d,%edi
  push   %rbp
  mov    0x48(%rsp),%r9d
  mov    0x60(%rsp),%rcx
  mov    0x40(%rsp),%edx
  mov    %r8,0x28(%rsp)
  call   <xgerc@plt>
  pop    %rax
  movq   D6(%rip),%xmm2
  pop    %rdx
  mov    0x18(%rsp),%r8
  jmp    L115
  nopw   0x0(%rax,%rax,1)
L114:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L116:
  movapd %xmm1,%xmm3
  jmp    L59
```

tool_calls: name: retrieve_data
arguments: data_label: D42
data_type: f64
name: retrieve_data
arguments: data_label: D6
data_type: f32[4]

tool: D42: 1

tool: {"D6": [0.0,-0.0,0.0,0.0]}

assistant: ```c
void xorgqr(int n, int k, double A_data[], int lda, const double tau_data[], double *work_data) {
  double c;
  int b_i;
  int c_i;
  int exitg1;
  int i;
  int i1;
  int ia;
  int iac;
  int iaii;
  int ic0;
  int itau;
  int ix;
  int lastc;
  int lastv;
  char exitg2;
  i = n - 1;
  for (b_i = k; b_i <= i; b_i++) {
    ia = b_i * lda;
    i1 = n - 1;
    memset(&A_data[ia], 0, (i1 + 1) * sizeof(double));
    A_data[ia + b_i] = 1.0;
  }
  itau = k - 1;
  for (c_i = k; c_i >= 1; c_i--) {
    iaii = c_i + (c_i - 1) * lda;
    if (c_i < n) {
      A_data[iaii - 1] = 1.0;
      ic0 = iaii + lda;
      if (tau_data[itau] != 0.0) {
        lastv = (n - c_i) + 1;
        b_i = (iaii + n) - c_i;
        while ((lastv > 0) && (A_data[b_i - 1] == 0.0)) {
          lastv--;
          b_i--;
        }
        lastc = n - c_i;
        exitg2 = 0;
        while ((!exitg2) && (lastc > 0)) {
          b_i = ic0 + (lastc - 1) * lda;
          ia = b_i;
          do {
            exitg1 = 0;
            if (ia <= (b_i + lastv) - 1) {
              if (A_data[ia - 1] != 0.0)
                exitg1 = 1;
              else
                ia++;
            } else {
              lastc--;
              exitg1 = 2;
            }
          } while (exitg1 == 0);
          if (exitg1 == 1)
            exitg2 = 1;
        }
      } else {
        lastv = 0;
        lastc = 0;
      }
      if (lastv > 0) {
        if (lastc != 0) {
          memset(&work_data[0], 0, lastc * sizeof(double));
          b_i = 0;
          i = ic0 + lda * (lastc - 1);
          for (iac = ic0; lda < 0 ? iac >= i : iac <= i; iac += lda) {
            ix = iaii;
            c = 0.0;
            i1 = (iac + lastv) - 1;
            for (ia = iac; ia <= i1; ia++) {
              c += A_data[ia - 1] * A_data[ix - 1];
              ix++;
            }
            work_data[b_i] += c;
            b_i++;
          }
        }
        xgerc(lastv, lastc, -tau_data[itau], iaii, work_data, A_data, ic0, lda);
      }
    }
    if (c_i < n) {
      b_i = iaii + 1;
      i = (iaii + n) - c_i;
      for (ix = b_i; ix <= i; ix++)
        A_data[ix - 1] *= -tau_data[itau];
    }
    A_data[iaii - 1] = 1.0 - tau_data[itau];
    for (b_i = 0; b_i <= c_i - 2; b_i++)
      A_data[(iaii - b_i) - 2] = 0.0;
    itau--;
  }
}
```