user: Demonstrate the C code that generates the assembly output below:

```asm
<xorgqr>:
  endbr64
  push   %r15
  mov    %rdx,%r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %ecx,%ebp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %edi,0x20(%rsp)
  mov    %esi,0x8(%rsp)
  mov    %r8,0x10(%rsp)
  mov    %r9,0x70(%rsp)
  cmp    %esi,%edi
  jle    L88
  movslq %esi,%rdx
  movslq %ebp,%rbx
  movsd  D22(%rip),%xmm7
  lea    0x0(,%rdi,8),%r14d
  mov    %edx,%eax
  lea    0x0(,%rbx,8),%rdi
  mov    %r15,0x18(%rsp)
  imul   %ebp,%eax
  lea    0x8(%rdi),%rsi
  mov    %ebp,0x24(%rsp)
  mov    %edx,%ebp
  movsd  %xmm7,(%rsp)
  mov    %rsi,%r12
  cltq
  lea    (%r15,%rax,8),%rbx
  add    %rdx,%rax
  lea    (%r15,%rax,8),%r13
  mov    %rbx,%r15
  mov    %rdi,%rbx
  nopl   0x0(%rax)
L225:
  mov    %r15,%rdi
  mov    %r14d,%edx
  xor    %esi,%esi
  add    $0x1,%ebp
  call   <memset@plt>
  movsd  (%rsp),%xmm4
  add    %rbx,%r15
  movsd  %xmm4,0x0(%r13)
  add    %r12,%r13
  cmp    %ebp,0x20(%rsp)
  jne    L225
  mov    0x18(%rsp),%r15
  mov    0x24(%rsp),%ebp
L88:
  mov    0x8(%rsp),%edi
  lea    -0x1(%rdi),%edx
  test   %edi,%edi
  jle    L315
  mov    %ebp,%eax
  movslq %edi,%r14
  mov    0x10(%rsp),%rbx
  mov    %ebp,%r12d
  not    %eax
  neg    %r12d
  movsd  D22(%rip),%xmm7
  mov    %r15,%r8
  mov    %eax,0x68(%rsp)
  cltq
  movslq %r12d,%r12
  pxor   %xmm1,%xmm1
  mov    %rax,0x48(%rsp)
  shl    $0x3,%rax
  mov    %rax,0x50(%rsp)
  mov    %edx,%eax
  movslq %edx,%rdx
  imul   %ebp,%eax
  lea    (%rbx,%rdx,8),%rbx
  movsd  %xmm7,(%rsp)
  movslq %eax,%rcx
  lea    (%rax,%rdi,1),%edx
  lea    (%rcx,%r14,1),%r9
  mov    0x20(%rsp),%ecx
  mov    %edx,0x24(%rsp)
  lea    0x0(,%r9,8),%r11
  mov    %r9,0x40(%rsp)
  sub    %edi,%ecx
  mov    %r9,0x30(%rsp)
  mov    %r11,%r13
  mov    %ecx,%esi
  add    %edx,%esi
  mov    %esi,0x18(%rsp)
  lea    (%rdx,%rbp,1),%esi
  lea    -0x1(%rcx),%edx
  imul   %ebp,%edx
  mov    %esi,0x38(%rsp)
  lea    -0x1(%r12),%rsi
  shl    $0x3,%r12
  mov    %rsi,0x58(%rsp)
  shl    $0x3,%rsi
  mov    %rsi,0x60(%rsp)
  lea    (%r15,%r11,1),%rsi
  mov    %ecx,%r11d
  add    %edx,%eax
  mov    %rsi,0x28(%rsp)
  add    %ebp,%eax
  mov    %edx,0x3c(%rsp)
  mov    %eax,0x6c(%rsp)
  movslq %ebp,%rax
  shl    $0x3,%rax
  mov    %rax,0x78(%rsp)
  nopl   0x0(%rax)
L199:
  cmp    %r14d,0x20(%rsp)
  jle    L274
  movsd  (%rsp),%xmm5
  movsd  %xmm5,-0x8(%r8,%r13,1)
  movsd  (%rbx),%xmm0
  ucomisd %xmm1,%xmm0
  jp     L345
  jne    L345
L181:
  mov    0x24(%rsp),%eax
  lea    0x1(%rax),%edx
  cmp    %edx,0x18(%rsp)
  jl     L274
  mov    %r11d,%eax
  add    0x30(%rsp),%rax
  mov    0x28(%rsp),%rsi
  lea    0x8(%rbx),%rcx
  lea    (%r8,%rax,8),%rax
  cmp    %rax,%rbx
  setae  %al
  cmp    %rcx,%rsi
  setae  %cl
  or     %cl,%al
  je     L342
  cmp    $0x1,%r11d
  je     L342
  mov    %r11d,%eax
  movsd  (%rbx),%xmm2
  xorpd  D49(%rip),%xmm2
  shr    %eax
  sub    $0x1,%eax
  unpcklpd %xmm2,%xmm2
  shl    $0x4,%rax
  lea    0x10(%rsi,%rax,1),%rcx
  mov    %rsi,%rax
  nopl   0x0(%rax,%rax,1)
L208:
  movupd (%rax),%xmm0
  add    $0x10,%rax
  mulpd  %xmm2,%xmm0
  movups %xmm0,-0x10(%rax)
  cmp    %rcx,%rax
  jne    L208
  mov    %r11d,%ecx
  and    $0xfffffffe,%ecx
  lea    (%rcx,%rdx,1),%eax
  cmp    %r11d,%ecx
  je     L274
  cltq
  movsd  (%rbx),%xmm0
  xorpd  D49(%rip),%xmm0
  lea    -0x8(%r8,%rax,8),%rax
  mulsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
L274:
  movsd  (%rsp),%xmm0
  subsd  (%rbx),%xmm0
  sub    $0x1,%r14
  sub    $0x8,%rbx
  mov    0x50(%rsp),%rax
  mov    0x68(%rsp),%esi
  add    $0x1,%r11d
  mov    0x60(%rsp),%rcx
  mov    0x58(%rsp),%rdi
  add    %rcx,0x28(%rsp)
  lea    (%rax,%r13,1),%r15
  mov    0x48(%rsp),%rcx
  add    %esi,0x24(%rsp)
  sub    %ebp,0x18(%rsp)
  add    %esi,0x38(%rsp)
  add    %ebp,0x3c(%rsp)
  add    %rcx,0x40(%rsp)
  add    %rdi,0x30(%rsp)
  movsd  %xmm0,-0x8(%r8,%r13,1)
  test   %r14d,%r14d
  jle    L315
  lea    -0x1(%r14),%eax
  xor    %esi,%esi
  mov    %r11d,0x10(%rsp)
  lea    0x8(,%rax,8),%rdx
  neg    %rax
  mov    %r8,0x8(%rsp)
  lea    -0x10(%r13,%rax,8),%rdi
  mov    %r15,%r13
  add    %r8,%rdi
  call   <memset@plt>
  mov    0x8(%rsp),%r8
  mov    0x10(%rsp),%r11d
  pxor   %xmm1,%xmm1
  jmp    L199
  nopl   0x0(%rax)
L345:
  mov    0x20(%rsp),%eax
  mov    %r14d,%ecx
  mov    %eax,%edi
  lea    0x1(%rax),%r10d
  movslq 0x18(%rsp),%rax
  sub    %r14d,%edi
  sub    %r14d,%r10d
  lea    (%r8,%rax,8),%rax
  nopl   0x0(%rax)
L174:
  ucomisd -0x8(%rax),%xmm1
  jp     L317
  jne    L317
  sub    $0x8,%rax
  sub    $0x1,%r10d
  jne    L174
L317:
  movslq 0x38(%rsp),%rdx
  movslq 0x3c(%rsp),%rax
  lea    -0x1(%r10),%r9d
  add    0x6c(%rsp),%ecx
  add    %rdx,%rax
  mov    %rdx,0x8(%rsp)
  lea    (%rcx,%r10,1),%esi
  add    %r9,%rax
  not    %r9
  lea    (%r8,%rax,8),%rdx
  shl    $0x3,%r9
  nopw   0x0(%rax,%rax,1)
L57:
  mov    %edi,%r15d
  sub    $0x1,%edi
  cmp    %esi,%ecx
  jge    L349
  lea    (%r9,%rdx,1),%rax
  xchg   %ax,%ax
L172:
  ucomisd (%rax),%xmm1
  jp     L28
  jne    L28
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L172
L349:
  sub    %ebp,%esi
  sub    %ebp,%ecx
  add    %r12,%rdx
  test   %edi,%edi
  jne    L57
  xor    %r15d,%r15d
  test   %r10d,%r10d
  jne    L10
  jmp    L181
  nopl   (%rax)
L28:
  test   %r10d,%r10d
  mov    %r10d,0x10(%rsp)
  je     L181
  mov    0x70(%rsp),%rdi
  lea    0x0(,%r15,8),%edx
  xor    %esi,%esi
  mov    %r8,0x80(%rsp)
  mov    %r11d,0x9c(%rsp)
  call   <memset@plt>
  lea    -0x1(%r15),%eax
  mov    0x10(%rsp),%r10d
  mov    0x38(%rsp),%esi
  imul   %ebp,%eax
  mov    0x8(%rsp),%rcx
  mov    0x80(%rsp),%r8
  mov    %r12,0x8(%rsp)
  lea    (%r10,%rsi,1),%edx
  mov    %r13,0x10(%rsp)
  mov    0x70(%rsp),%rdi
  pxor   %xmm1,%xmm1
  mov    %r10d,0x80(%rsp)
  mov    0x24(%rsp),%r13d
  lea    -0x8(%r8,%rcx,8),%rcx
  mov    %edx,%r10d
  add    %esi,%eax
  mov    %r14,0x88(%rsp)
  mov    0x9c(%rsp),%r11d
  mov    %rbx,0x90(%rsp)
  mov    0x28(%rsp),%r14
  mov    %eax,%r12d
  mov    %r15d,0x98(%rsp)
  mov    0x78(%rsp),%rbx
  test   %ebp,%ebp
  js     L50
L177:
  cmp    %r12d,%esi
  jg     L78
L110:
  cmp    %esi,%r10d
  jle    L167
  mov    %r10d,%r9d
  sub    %esi,%r9d
  cmp    $0x1,%r9d
  je     L146
  mov    %r9d,%edx
  lea    -0x8(%r14),%r15
  movapd %xmm1,%xmm0
  xor    %eax,%eax
  shr    %edx
  shl    $0x4,%rdx
  nopw   0x0(%rax,%rax,1)
L320:
  movupd (%r15,%rax,1),%xmm3
  movupd (%rcx,%rax,1),%xmm6
  add    $0x10,%rax
  mulpd  %xmm6,%xmm3
  addsd  %xmm3,%xmm0
  unpckhpd %xmm3,%xmm3
  addsd  %xmm3,%xmm0
  cmp    %rax,%rdx
  jne    L320
  mov    %r9d,%r15d
  and    $0xfffffffe,%r15d
  lea    (%rsi,%r15,1),%eax
  lea    0x0(%r13,%r15,1),%edx
  cmp    %r9d,%r15d
  je     L214
L19:
  movslq %edx,%rdx
  cltq
  movsd  -0x8(%r8,%rdx,8),%xmm2
  mulsd  -0x8(%r8,%rax,8),%xmm2
  addsd  %xmm2,%xmm0
L214:
  addsd  (%rdi),%xmm0
  add    %ebp,%esi
  add    $0x8,%rdi
  add    %ebp,%r10d
  add    %rbx,%rcx
  movsd  %xmm0,-0x8(%rdi)
  test   %ebp,%ebp
  jns    L177
L50:
  cmp    %r12d,%esi
  jge    L110
L78:
  mov    0x90(%rsp),%rbx
  mov    0x8(%rsp),%r12
  mov    0x10(%rsp),%r13
  mov    0x80(%rsp),%r10d
  mov    0x98(%rsp),%r15d
  mov    0x88(%rsp),%r14
  movsd  (%rbx),%xmm0
L10:
  mov    %r11d,0x10(%rsp)
  sub    $0x8,%rsp
  mov    %r15d,%esi
  mov    %r10d,%edi
  push   %rbp
  xorpd  D49(%rip),%xmm0
  mov    0x48(%rsp),%r9d
  mov    0x80(%rsp),%rcx
  mov    0x34(%rsp),%edx
  mov    %r8,0x18(%rsp)
  call   <xgerc@plt>
  pop    %rax
  pop    %rdx
  mov    0x8(%rsp),%r8
  mov    0x10(%rsp),%r11d
  jmp    L181
  nopl   0x0(%rax,%rax,1)
L315:
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L342:
  movslq %edx,%rax
  lea    -0x1(%r11),%edx
  add    0x40(%rsp),%rdx
  lea    -0x8(%r8,%rax,8),%rax
  lea    0x8(%r8,%rdx,8),%rdx
  nopl   0x0(%rax,%rax,1)
L337:
  movsd  (%rbx),%xmm0
  xorpd  D49(%rip),%xmm0
  add    $0x8,%rax
  mulsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L337
  jmp    L274
L167:
  movapd %xmm1,%xmm0
  jmp    L214
L146:
  mov    %r13d,%edx
  mov    %esi,%eax
  movapd %xmm1,%xmm0
  jmp    L19
```

tool_calls: name: read_data
arguments: label: D22
type_infered: f64
name: read_data
arguments: label: D49
type_infered: f32[4]

tool: D22: 1

tool: {"D49": [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--;
  }
}
```