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

```asm
<daxpy_vec>:
  endbr64
  push   %r15
  mov    %rdi,%r10
  push   %r14
  push   %r13
  mov    %r8d,%r13d
  push   %r12
  movslq %ecx,%r12
  push   %rbp
  movslq %edx,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x18,%rsp
  ucomisd D13(%rip),%xmm1
  jp     L99
  jne    L99
L6:
  lea    -0x3(%r13),%r15d
  xor    %edx,%edx
  test   %r15d,%r15d
  jle    L125
L46:
  movslq %ebp,%r11
  movslq %r12d,%rax
  mov    %rbx,%rcx
  mov    %r10,%rdx
  mov    %rax,%r9
  lea    (%r11,%r11,2),%r14
  xor    %edi,%edi
  lea    0x0(,%rax,8),%rsi
  shl    $0x4,%rax
  lea    0x0(,%r11,8),%r8
  shl    $0x5,%r9
  shl    $0x3,%r14
  add    %r10,%rax
  shl    $0x5,%r11
  nopl   0x0(%rax)
L119:
  movsd  (%rcx),%xmm1
  add    $0x4,%edi
  mulsd  %xmm0,%xmm1
  addsd  (%rdx),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  (%rcx,%r8,1),%xmm1
  mulsd  %xmm0,%xmm1
  addsd  (%rdx,%rsi,1),%xmm1
  movsd  %xmm1,(%rdx,%rsi,1)
  movsd  (%rcx,%r8,2),%xmm1
  add    %r9,%rdx
  mulsd  %xmm0,%xmm1
  addsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  movsd  (%rcx,%r14,1),%xmm1
  add    %r11,%rcx
  mulsd  %xmm0,%xmm1
  addsd  (%rax,%rsi,1),%xmm1
  movsd  %xmm1,(%rax,%rsi,1)
  add    %r9,%rax
  cmp    %edi,%r15d
  jg     L119
  lea    -0x4(%r13),%eax
  shr    $0x2,%eax
  add    $0x1,%eax
  imul   %rax,%r9
  lea    0x0(,%rax,4),%edx
  imul   %r11,%rax
  add    %r9,%r10
  add    %rax,%rbx
L125:
  cmp    %edx,%r13d
  je     L132
L145:
  lea    -0x1(%r13),%eax
  cmp    %edx,%eax
  jle    L98
  movsd  (%rbx),%xmm1
  shl    $0x3,%rbp
  shl    $0x3,%r12
  add    $0x2,%edx
  add    %rbp,%rbx
  mulsd  %xmm0,%xmm1
  addsd  (%r10),%xmm1
  movsd  %xmm1,(%r10)
  movsd  (%rbx),%xmm1
  add    %r12,%r10
  add    %rbp,%rbx
  mulsd  %xmm0,%xmm1
  addsd  (%r10),%xmm1
  movsd  %xmm1,(%r10)
  add    %r12,%r10
L98:
  cmp    %r13d,%edx
  jge    L132
  mulsd  (%rbx),%xmm0
  addsd  (%r10),%xmm0
  movsd  %xmm0,(%r10)
L132:
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L99:
  ucomisd D15(%rip),%xmm1
  jp     L30
  je     L156
L30:
  test   %r13d,%r13d
  jle    L6
  movslq %r13d,%rax
  lea    (%r10,%rax,8),%rdx
  mov    %r10,%rax
  nopw   0x0(%rax,%rax,1)
L32:
  movsd  (%rax),%xmm2
  add    $0x8,%rax
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L32
  lea    -0x3(%r13),%r15d
  xor    %edx,%edx
  test   %r15d,%r15d
  jg     L46
  jmp    L145
  nopw   0x0(%rax,%rax,1)
L156:
  test   %r13d,%r13d
  jle    L6
  movslq %r13d,%rdx
  mov    %r10,%rdi
  xor    %esi,%esi
  movsd  %xmm0,0x8(%rsp)
  shl    $0x3,%rdx
  lea    -0x3(%r13),%r15d
  call   <memset@plt>
  movsd  0x8(%rsp),%xmm0
  xor    %edx,%edx
  mov    %rax,%r10
  test   %r15d,%r15d
  jg     L46
  jmp    L145
```

tool_calls: name: unpack_data
arguments: data_tag: D13
datatype: f64
name: unpack_data
arguments: data_tag: D15
datatype: f64

tool: D13: 1

tool: D15: 0

assistant: ```c
void daxpy_vec(double *Y, const double *X, double alpha, double beta, int incX, int incY, int N) {
  register int i;
  register double *y0;
  register const double *x0;
  if (beta != 1.0) {
    y0 = Y;
    if (beta == 0.0) {
      for (i = 0; i < N; i++) {
        y0[0] = 0.0;
        y0++;
      }
    } else {
      for (i = 0; i < N; i++) {
        y0[0] *= beta;
        y0++;
      }
    }
  }
  y0 = Y;
  x0 = X;
  for (i = 0; i < N - 3; i += 4) {
    y0[0] += alpha * x0[0];
    x0 += incX;
    y0 += incY;
    y0[0] += alpha * x0[0];
    x0 += incX;
    y0 += incY;
    y0[0] += alpha * x0[0];
    x0 += incX;
    y0 += incY;
    y0[0] += alpha * x0[0];
    x0 += incX;
    y0 += incY;
  }
  if (i == N) {
    return;
  }
  if (i < N - 1) {
    y0[0] += alpha * x0[0];
    x0 += incX;
    y0 += incY;
    y0[0] += alpha * x0[0];
    x0 += incX;
    y0 += incY;
    i += 2;
  }
  if (i < N) {
    y0[0] += alpha * x0[0];
  }
  return;
}
```