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

```asm
<gammabootstrap>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x68,%rsp
  mov    %rdi,%r14
  mov    %rdx,-0x78(%rbp)
  mov    %ecx,-0x48(%rbp)
  mov    %r8d,%r12d
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  test   %r8d,%r8d
  jle    L104
  mov    %rdi,%rax
  lea    -0x1(%r8),%edx
  lea    0x4(%rdi,%rdx,4),%rcx
  mov    %rdi,%rdx
  pxor   %xmm1,%xmm1
L180:
  addss  (%rdx),%xmm1
  add    $0x4,%rdx
  cmp    %rcx,%rdx
  jne    L180
  pxor   %xmm5,%xmm5
  cvtsi2ss %r12d,%xmm5
  movss  %xmm5,-0x44(%rbp)
  divss  %xmm5,%xmm1
  pxor   %xmm2,%xmm2
L333:
  movss  (%rax),%xmm0
  subss  %xmm1,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm2
  add    $0x4,%rax
  cmp    %rcx,%rax
  jne    L333
L30:
  divss  -0x44(%rbp),%xmm2
  movaps %xmm1,%xmm0
  mulss  %xmm1,%xmm0
  divss  %xmm2,%xmm0
  movss  %xmm0,-0x6c(%rbp)
  divss  %xmm1,%xmm2
  movss  %xmm2,-0x70(%rbp)
  movslq %r12d,%rax
  lea    0x0(,%rax,4),%rsi
  mov    %rsi,-0x50(%rbp)
  movslq -0x48(%rbp),%rdx
  lea    0x0(,%rdx,4),%rsi
  mov    %rsi,-0x68(%rbp)
  imul   %rsi,%rax
  add    $0xf,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffffff0,%rcx
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rdx
  sub    %rax,%rdx
L76:
  cmp    %rdx,%rsp
  je     L239
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L76
L239:
  mov    %rcx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L82
  orq    $0x0,-0x8(%rsp,%rax,1)
L82:
  mov    %rsp,%rbx
  mov    -0x68(%rbp),%rax
  add    $0xf,%rax
  mov    %rax,%rdx
  and    $0xfffffffffffffff0,%rdx
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L188:
  cmp    %rcx,%rsp
  je     L231
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L188
L231:
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L22
  orq    $0x0,-0x8(%rsp,%rax,1)
L22:
  mov    %rsp,-0x80(%rbp)
  mov    -0x68(%rbp),%rax
  lea    0xf(%rax),%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffffff0,%rax
  and    $0xfffffffffffff000,%rdx
  mov    %rsp,%rcx
  sub    %rdx,%rcx
L292:
  cmp    %rcx,%rsp
  je     L236
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L292
L236:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  test   %rdx,%rdx
  je     L226
  orq    $0x0,-0x8(%rsp,%rdx,1)
L226:
  mov    %rsp,-0x90(%rbp)
  lea    -0x3c(%rbp),%rdi
  call   <time@plt>
  mov    %eax,%edi
  call   <srand@plt>
  cmpl   $0x0,-0x48(%rbp)
  jle    L27
  lea    -0x1(%r12),%eax
  lea    0x4(%rbx,%rax,4),%r15
  neg    %rax
  shl    $0x2,%rax
  mov    %rax,-0x60(%rbp)
  mov    %r15,%r13
  mov    $0x0,%esi
  sub    $0x4,%rax
  mov    %rax,-0x58(%rbp)
  mov    %r15,-0x88(%rbp)
  mov    %esi,%r15d
  jmp    L11
L267:
  mov    $0x0,%eax
  call   <rand@plt>
  cltd
  idiv   %r12d
  movslq %edx,%rdx
  cvttss2si (%r14,%rdx,4),%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  movss  %xmm0,(%rbx)
  add    $0x4,%rbx
  cmp    %r13,%rbx
  jne    L267
L345:
  lea    0x1(%r15),%esi
  add    -0x50(%rbp),%r13
  cmp    %esi,-0x48(%rbp)
  je     L124
  mov    %esi,%r15d
L11:
  mov    -0x58(%rbp),%rax
  lea    (%rax,%r13,1),%rbx
  test   %r12d,%r12d
  jg     L267
  jmp    L345
L21:
  addss  (%rax),%xmm0
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L21
L204:
  divss  -0x44(%rbp),%xmm0
  movss  %xmm0,(%rdi)
  add    $0x4,%rdi
  add    -0x50(%rbp),%rdx
  cmp    %r9,%rdi
  je     L158
L280:
  lea    (%r8,%rdx,1),%rax
  pxor   %xmm0,%xmm0
  test   %r12d,%r12d
  jg     L21
  jmp    L204
L158:
  mov    $0x0,%edx
  mov    -0x60(%rbp),%rdi
  sub    $0x4,%rdi
  mov    -0x50(%rbp),%r8
  mov    -0x80(%rbp),%r10
  mov    -0x90(%rbp),%r9
  jmp    L399
L124:
  mov    %r15d,%eax
  mov    %r15d,%r11d
  mov    -0x88(%rbp),%r15
  mov    -0x80(%rbp),%rbx
  mov    %rbx,%rdi
  mov    %eax,%ecx
  lea    0x4(%rbx,%rcx,4),%r9
  mov    %r15,%rdx
  mov    -0x60(%rbp),%rax
  lea    -0x4(%rax),%r8
  jmp    L280
L121:
  movss  (%r10,%rdx,4),%xmm2
  lea    (%rdi,%r15,1),%rax
  pxor   %xmm1,%xmm1
L289:
  movss  (%rax),%xmm0
  subss  %xmm2,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  add    $0x4,%rax
  cmp    %rax,%r15
  jne    L289
L170:
  divss  -0x44(%rbp),%xmm1
  movss  %xmm1,(%r9,%rdx,4)
  lea    0x1(%rdx),%rax
  add    %r8,%r15
  cmp    %rcx,%rdx
  je     L269
  mov    %rax,%rdx
L399:
  pxor   %xmm1,%xmm1
  test   %r12d,%r12d
  jg     L121
  jmp    L170
L269:
  mov    -0x68(%rbp),%rax
  add    $0xf,%rax
  mov    %rax,%rdi
  and    $0xfffffffffffffff0,%rdi
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rdx
  sub    %rax,%rdx
L237:
  cmp    %rdx,%rsp
  je     L154
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L237
L154:
  mov    %rdi,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L215
  orq    $0x0,-0x8(%rsp,%rax,1)
L215:
  mov    %rsp,%rdx
  mov    -0x68(%rbp),%rax
  add    $0xf,%rax
  mov    %rax,%rdi
  and    $0xfffffffffffffff0,%rdi
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%r8
  sub    %rax,%r8
L138:
  cmp    %r8,%rsp
  je     L193
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L138
L193:
  mov    %rdi,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L167
  orq    $0x0,-0x8(%rsp,%rax,1)
L167:
  mov    %rsp,%rdi
  mov    $0x0,%eax
  mov    -0x80(%rbp),%r9
  mov    -0x90(%rbp),%r10
L14:
  movss  (%r9,%rax,4),%xmm1
  movss  (%r10,%rax,4),%xmm0
  movaps %xmm1,%xmm2
  mulss  %xmm1,%xmm2
  divss  %xmm0,%xmm2
  movss  %xmm2,(%rdx,%rax,4)
  divss  %xmm1,%xmm0
  movss  %xmm0,(%rdi,%rax,4)
  mov    %rax,%r8
  add    $0x1,%rax
  cmp    %rcx,%r8
  jne    L14
  mov    $0x0,%eax
  pxor   %xmm0,%xmm0
  movaps %xmm0,%xmm1
L146:
  addss  (%rdx,%rax,4),%xmm1
  addss  (%rdi,%rax,4),%xmm0
  mov    %rax,%r8
  add    $0x1,%rax
  cmp    %rcx,%r8
  jne    L146
  pxor   %xmm2,%xmm2
  cvtsi2ss %esi,%xmm2
  divss  %xmm2,%xmm1
  divss  %xmm2,%xmm0
  mov    $0x0,%eax
  mov    %r11d,%ecx
L248:
  mov    %eax,%edx
  add    $0x1,%eax
  cmp    %edx,%ecx
  jne    L248
L220:
  movss  -0x6c(%rbp),%xmm2
  addss  %xmm2,%xmm2
  subss  %xmm1,%xmm2
  mov    -0x78(%rbp),%rax
  movss  %xmm2,(%rax)
  movss  -0x70(%rbp),%xmm1
  addss  %xmm1,%xmm1
  subss  %xmm0,%xmm1
  movss  %xmm1,0x4(%rax)
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L64
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L27:
  pxor   %xmm0,%xmm0
  cvtsi2ssl -0x48(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  divss  %xmm0,%xmm1
  movaps %xmm1,%xmm0
  jmp    L220
L104:
  pxor   %xmm7,%xmm7
  cvtsi2ss %r8d,%xmm7
  movss  %xmm7,-0x44(%rbp)
  pxor   %xmm2,%xmm2
  movaps %xmm2,%xmm1
  divss  %xmm7,%xmm1
  jmp    L30
L64:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void gammabootstrap(float *surrogatenum, int *surrogatefreq, float *gammaparam, int numboot, int surrogatesize) {
  int samplesize = surrogatesize;
  float mainmean = 0;
  float mainvar = 0;
  int bootrep = numboot;
  for (int i = 0; i < surrogatesize; i++)
    mainmean = mainmean + surrogatenum[i];
  mainmean = mainmean / samplesize;
  for (int i = 0; i < surrogatesize; i++)
    mainvar = mainvar + ((surrogatenum[i] - mainmean) * (surrogatenum[i] - mainmean));
  mainvar = mainvar / samplesize;
  float mainalpha = (mainmean * mainmean) / mainvar;
  float maintheta = mainvar / mainmean;
  float bootsamples[bootrep][samplesize];
  float bootsampmean[bootrep];
  float bootsampvar[bootrep];
  time_t t;
  srand((unsigned)time(&t));
  for (int i = 0; i < bootrep; i++) {
    for (int j = 0; j < surrogatesize; j++) {
      int randomsample = surrogatenum[rand() % samplesize];
      bootsamples[i][j] = randomsample;
    }
  }
  for (int i = 0; i < bootrep; i++) {
    float tempmean = 0;
    for (int j = 0; j < samplesize; j++) {
      tempmean = tempmean + bootsamples[i][j];
    }
    tempmean = tempmean / samplesize;
    bootsampmean[i] = tempmean;
  }
  for (int i = 0; i < bootrep; i++) {
    float tempvar = 0;
    for (int j = 0; j < samplesize; j++) {
      tempvar = tempvar + ((bootsamples[i][j] - bootsampmean[i]) * (bootsamples[i][j] - bootsampmean[i]));
    }
    tempvar = tempvar / samplesize;
    bootsampvar[i] = tempvar;
  }
  float bootalpha[bootrep];
  float boottheta[bootrep];
  for (int i = 0; i < bootrep; i++) {
    bootalpha[i] = (bootsampmean[i] * bootsampmean[i]) / (bootsampvar[i]);
    boottheta[i] = (bootsampvar[i]) / (bootsampmean[i]);
  }
  float normalalphamean = 0;
  float normalalphavar = 0;
  float normalthetamean = 0;
  float normalthetavar = 0;
  for (int i = 0; i < bootrep; i++) {
    normalalphamean = normalalphamean + bootalpha[i];
    normalthetamean = normalthetamean + boottheta[i];
  }
  normalalphamean = normalalphamean / bootrep;
  normalthetamean = normalthetamean / bootrep;
  for (int i = 0; i < bootrep; i++) {
    normalalphavar = normalalphavar + ((bootalpha[i] - normalalphamean) * (bootalpha[i] - normalalphamean));
    normalthetavar = normalthetavar + ((boottheta[i] - normalthetamean) * (boottheta[i] - normalthetamean));
  }
  normalalphavar = normalalphavar / bootrep;
  normalthetavar = normalthetavar / bootrep;
  float alpha = (2 * mainalpha) - normalalphamean;
  float theta = (2 * maintheta) - normalthetamean;
  gammaparam[0] = alpha;
  gammaparam[1] = theta;
}
```