user: 下面的汇编指令对应的C源代码是什么:

```asm
<gammabootstrap>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbx
  mov    %r8d,%ebx
  sub    $0x78,%rsp
  mov    %rdx,-0x88(%rbp)
  mov    %ecx,-0x48(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  test   %r8d,%r8d
  jle    L44
  lea    -0x1(%r8),%eax
  mov    %eax,-0x68(%rbp)
  cmp    $0x2,%eax
  jbe    L194
  mov    %r8d,%edx
  pxor   %xmm7,%xmm7
  lea    0x10(%rdi),%rcx
  mov    %rdi,%rax
  shr    $0x2,%edx
  movss  %xmm7,-0x64(%rbp)
  movaps %xmm7,%xmm1
  sub    $0x1,%edx
  shl    $0x4,%rdx
  add    %rcx,%rdx
L559:
  movups (%rax),%xmm0
  add    $0x10,%rax
  addss  %xmm0,%xmm1
  movaps %xmm0,%xmm2
  shufps $0x55,%xmm0,%xmm2
  addss  %xmm2,%xmm1
  movaps %xmm0,%xmm2
  unpckhps %xmm0,%xmm2
  shufps $0xff,%xmm0,%xmm0
  addss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  cmp    %rdx,%rax
  jne    L559
  mov    %ebx,%eax
  and    $0xfffffffc,%eax
  cmp    %eax,%ebx
  je     L641
L166:
  movslq %eax,%rdx
  addss  0x0(%r13,%rdx,4),%xmm1
  lea    0x0(,%rdx,4),%rcx
  lea    0x1(%rax),%edx
  cmp    %edx,%ebx
  jle    L279
  add    $0x2,%eax
  addss  0x4(%r13,%rcx,1),%xmm1
  cmp    %eax,%ebx
  jle    L279
  addss  0x8(%r13,%rcx,1),%xmm1
L279:
  pxor   %xmm7,%xmm7
  movaps %xmm1,%xmm3
  cmpl   $0x2,-0x68(%rbp)
  cvtsi2ss %ebx,%xmm7
  divss  %xmm7,%xmm3
  movss  %xmm7,-0x44(%rbp)
  jbe    L526
  lea    0x10(%r13),%rcx
L646:
  mov    %ebx,%edx
  movaps %xmm3,%xmm4
  pxor   %xmm1,%xmm1
  mov    %r13,%rax
  shr    $0x2,%edx
  shufps $0x0,%xmm4,%xmm4
  sub    $0x1,%edx
  shl    $0x4,%rdx
  add    %rcx,%rdx
  nopw   0x0(%rax,%rax,1)
L607:
  movups (%rax),%xmm0
  add    $0x10,%rax
  subps  %xmm4,%xmm0
  mulps  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  movaps %xmm0,%xmm2
  shufps $0x55,%xmm0,%xmm2
  addss  %xmm2,%xmm1
  movaps %xmm0,%xmm2
  unpckhps %xmm0,%xmm2
  shufps $0xff,%xmm0,%xmm0
  addss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  cmp    %rdx,%rax
  jne    L607
  mov    %ebx,%eax
  and    $0xfffffffc,%eax
  test   $0x3,%bl
  je     L444
L374:
  movslq %eax,%rdx
  movss  0x0(%r13,%rdx,4),%xmm0
  lea    0x0(,%rdx,4),%rcx
  lea    0x1(%rax),%edx
  subss  %xmm3,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %edx,%ebx
  jle    L444
  movss  0x4(%r13,%rcx,1),%xmm0
  add    $0x2,%eax
  subss  %xmm3,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %eax,%ebx
  jle    L444
  movss  0x8(%r13,%rcx,1),%xmm0
  subss  %xmm3,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
L444:
  divss  -0x44(%rbp),%xmm1
L329:
  movaps %xmm3,%xmm0
  movslq -0x48(%rbp),%rdx
  movslq %ebx,%rax
  mov    %rsp,%rsi
  mulss  %xmm3,%xmm0
  lea    0x0(,%rax,4),%rdi
  shl    $0x2,%rdx
  mov    %rdi,-0x50(%rbp)
  imul   %rdx,%rax
  divss  %xmm1,%xmm0
  add    $0xf,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rax
  sub    %rax,%rsi
  and    $0xfffffffffffffff0,%rcx
  divss  %xmm3,%xmm1
  movaps %xmm0,%xmm7
  unpcklps %xmm1,%xmm7
  movlps %xmm7,-0x80(%rbp)
  cmp    %rsi,%rsp
  je     L265
L351:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rsi,%rsp
  jne    L351
L265:
  and    $0xfff,%ecx
  sub    %rcx,%rsp
  test   %rcx,%rcx
  jne    L440
L242:
  lea    0xf(%rdx),%rdi
  mov    %rsp,%rdx
  mov    %rsp,%r14
  mov    %rdi,%rax
  mov    %rdi,-0xa0(%rbp)
  and    $0xfffffffffffff000,%rdi
  sub    %rdi,%rdx
  and    $0xfffffffffffffff0,%rax
  cmp    %rdx,%rsp
  je     L472
L404:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rdx,%rsp
  jne    L404
L472:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L616
L451:
  mov    -0xa0(%rbp),%rdi
  mov    %rsp,%rdx
  mov    %rsp,-0x70(%rbp)
  mov    %rdi,%rax
  and    $0xfffffffffffff000,%rdi
  sub    %rdi,%rdx
  and    $0xfffffffffffffff0,%rax
  cmp    %rdx,%rsp
  je     L465
L459:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rdx,%rsp
  jne    L459
L465:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L120
  orq    $0x0,-0x8(%rsp,%rax,1)
L120:
  lea    -0x3c(%rbp),%rdi
  mov    %rsp,-0x78(%rbp)
  call   <time@plt>
  mov    %eax,%edi
  call   <srand@plt>
  mov    -0x48(%rbp),%eax
  test   %eax,%eax
  jle    L563
  mov    -0x50(%rbp),%rax
  shr    $0x2,%rax
  mov    %rax,-0x60(%rbp)
  test   %ebx,%ebx
  jle    L394
  lea    -0x1(%rbx),%eax
  mov    $0xfffffffffffffffc,%rdx
  xor    %edi,%edi
  mov    %r14,-0x90(%rbp)
  mov    %eax,-0x68(%rbp)
  shl    $0x2,%rax
  sub    %rax,%rdx
  lea    0x4(%r14,%rax,1),%r12
  mov    %edi,%r14d
  mov    %rdx,-0x58(%rbp)
L61:
  mov    -0x58(%rbp),%rax
  lea    (%rax,%r12,1),%r15
  nopl   0x0(%rax,%rax,1)
L142:
  xor    %eax,%eax
  add    $0x4,%r15
  call   <rand@plt>
  pxor   %xmm0,%xmm0
  cltd
  idiv   %ebx
  movslq %edx,%rdx
  cvttss2si 0x0(%r13,%rdx,4),%eax
  cvtsi2ss %eax,%xmm0
  movss  %xmm0,-0x4(%r15)
  cmp    %r12,%r15
  jne    L142
  add    $0x1,%r14d
  add    -0x50(%rbp),%r12
  cmp    %r14d,-0x48(%rbp)
  jne    L61
  mov    -0x90(%rbp),%r14
L326:
  mov    %ebx,%eax
  mov    -0x48(%rbp),%r9d
  mov    %r14,-0x58(%rbp)
  mov    %ebx,%r8d
  shr    $0x2,%eax
  mov    -0x50(%rbp),%rdi
  mov    -0x60(%rbp),%r12
  and    $0xfffffffc,%r8d
  sub    $0x1,%eax
  lea    -0x1(%r9),%r13d
  mov    -0x70(%rbp),%r11
  movss  -0x64(%rbp),%xmm3
  shl    $0x4,%rax
  mov    %r13d,-0x94(%rbp)
  mov    %r14,%rsi
  xor    %ecx,%ecx
  mov    %rax,-0x90(%rbp)
  mov    -0x68(%rbp),%r13d
  lea    0x10(%rax),%r10
  nopl   0x0(%rax,%rax,1)
L551:
  movslq %ecx,%rdx
  movaps %xmm3,%xmm1
  test   %ebx,%ebx
  jle    L301
  cmp    $0x2,%r13d
  jbe    L633
  lea    (%r10,%rsi,1),%r15
  mov    %rsi,%rax
  movaps %xmm3,%xmm1
  nopw   0x0(%rax,%rax,1)
L87:
  movups (%rax),%xmm0
  add    $0x10,%rax
  addss  %xmm0,%xmm1
  movaps %xmm0,%xmm2
  shufps $0x55,%xmm0,%xmm2
  addss  %xmm2,%xmm1
  movaps %xmm0,%xmm2
  unpckhps %xmm0,%xmm2
  shufps $0xff,%xmm0,%xmm0
  addss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  cmp    %r15,%rax
  jne    L87
  mov    %r8d,%eax
  cmp    %r8d,%ebx
  je     L301
L565:
  imul   %r12,%rdx
  movslq %eax,%r15
  add    %rdx,%r15
  addss  (%r14,%r15,4),%xmm1
  lea    0x1(%rax),%r15d
  cmp    %r15d,%ebx
  jle    L301
  movslq %r15d,%r15
  add    $0x2,%eax
  add    %rdx,%r15
  addss  (%r14,%r15,4),%xmm1
  cmp    %eax,%ebx
  jle    L301
  cltq
  add    %rdx,%rax
  addss  (%r14,%rax,4),%xmm1
L301:
  divss  -0x44(%rbp),%xmm1
  movss  %xmm1,(%r11,%rcx,4)
  add    $0x1,%rcx
  add    %rdi,%rsi
  cmp    %r9,%rcx
  jne    L551
  mov    -0x90(%rbp),%r12
  mov    -0x94(%rbp),%r13d
  xor    %edx,%edx
  mov    -0x50(%rbp),%r10
  mov    -0x58(%rbp),%rdi
  mov    %r13d,-0x50(%rbp)
  lea    0x10(%r12),%rcx
  mov    -0x78(%rbp),%r11
  mov    -0x70(%rbp),%r12
  movss  -0x64(%rbp),%xmm5
  mov    -0x68(%rbp),%r13d
  nopl   0x0(%rax)
L474:
  movslq %edx,%rsi
  movaps %xmm5,%xmm1
  test   %ebx,%ebx
  jle    L548
  movss  (%r12,%rdx,4),%xmm4
  cmp    $0x2,%r13d
  jbe    L227
  movaps %xmm4,%xmm3
  lea    (%rcx,%rdi,1),%r15
  movaps %xmm5,%xmm1
  mov    %rdi,%rax
  shufps $0x0,%xmm3,%xmm3
  nopw   0x0(%rax,%rax,1)
L337:
  movups (%rax),%xmm0
  add    $0x10,%rax
  subps  %xmm3,%xmm0
  mulps  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  movaps %xmm0,%xmm2
  shufps $0x55,%xmm0,%xmm2
  addss  %xmm2,%xmm1
  movaps %xmm0,%xmm2
  unpckhps %xmm0,%xmm2
  shufps $0xff,%xmm0,%xmm0
  addss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  cmp    %r15,%rax
  jne    L337
  mov    %r8d,%eax
  cmp    %r8d,%ebx
  je     L548
L307:
  imul   -0x60(%rbp),%rsi
  movslq %eax,%r15
  add    %rsi,%r15
  movss  (%r14,%r15,4),%xmm0
  lea    0x1(%rax),%r15d
  subss  %xmm4,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %r15d,%ebx
  jle    L548
  movslq %r15d,%r15
  add    $0x2,%eax
  add    %rsi,%r15
  movss  (%r14,%r15,4),%xmm0
  subss  %xmm4,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %eax,%ebx
  jle    L548
  cltq
  add    %rsi,%rax
  movss  (%r14,%rax,4),%xmm0
  subss  %xmm4,%xmm0
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
L548:
  divss  -0x44(%rbp),%xmm1
  movss  %xmm1,(%r11,%rdx,4)
  add    $0x1,%rdx
  add    %r10,%rdi
  cmp    %r9,%rdx
  jne    L474
  mov    -0xa0(%rbp),%rdi
  mov    %rsp,%rdx
  mov    -0x50(%rbp),%r13d
  mov    %rdi,%rax
  and    $0xfffffffffffff000,%rdi
  sub    %rdi,%rdx
  and    $0xfffffffffffffff0,%rax
  cmp    %rdx,%rsp
  je     L535
L322:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rdx,%rsp
  jne    L322
L535:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L64
  orq    $0x0,-0x8(%rsp,%rax,1)
L64:
  mov    -0xa0(%rbp),%rdi
  lea    0x3(%rsp),%rax
  mov    %rsp,%rsi
  mov    %rax,%r9
  and    $0xfffffffffffffffc,%rax
  mov    %rdi,%rcx
  and    $0xfffffffffffff000,%rdi
  shr    $0x2,%r9
  mov    %rax,%r8
  sub    %rdi,%rsi
  and    $0xfffffffffffffff0,%rcx
  cmp    %rsi,%rsp
  je     L594
L464:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rsi,%rsp
  jne    L464
L594:
  and    $0xfff,%ecx
  sub    %rcx,%rsp
  test   %rcx,%rcx
  je     L387
  orq    $0x0,-0x8(%rsp,%rcx,1)
L387:
  lea    0x3(%rsp),%rdx
  mov    %rdx,%r14
  and    $0xfffffffffffffffc,%rdx
  shr    $0x2,%r14
  mov    %rdx,%rsi
  cmp    $0x2,%r13d
  jbe    L216
  mov    -0x48(%rbp),%edi
  mov    -0x70(%rbp),%r10
  xor    %ecx,%ecx
  mov    -0x78(%rbp),%r11
  shr    $0x2,%edi
  shl    $0x4,%rdi
  nopl   0x0(%rax,%rax,1)
L334:
  movups (%r10,%rcx,1),%xmm1
  movups (%r11,%rcx,1),%xmm6
  movaps %xmm1,%xmm0
  mulps  %xmm1,%xmm0
  divps  %xmm6,%xmm0
  movups %xmm0,(%rax,%rcx,1)
  movaps %xmm6,%xmm0
  divps  %xmm1,%xmm0
  movups %xmm0,(%rdx,%rcx,1)
  add    $0x10,%rcx
  cmp    %rcx,%rdi
  jne    L334
  mov    -0x48(%rbp),%ecx
  mov    %ecx,%edi
  and    $0xfffffffc,%edi
  and    $0x3,%ecx
  je     L91
L623:
  mov    -0x70(%rbp),%r11
  movslq %edi,%rcx
  mov    -0x78(%rbp),%rbx
  lea    0x1(%rdi),%r12d
  shl    $0x2,%rcx
  mov    -0x48(%rbp),%r15d
  add    %rcx,%r11
  add    %rcx,%rbx
  lea    (%rax,%rcx,1),%r10
  add    %rdx,%rcx
  movss  (%r11),%xmm1
  movss  (%rbx),%xmm0
  movaps %xmm1,%xmm2
  mulss  %xmm1,%xmm2
  divss  %xmm0,%xmm2
  divss  %xmm1,%xmm0
  movss  %xmm2,(%r10)
  movss  %xmm0,(%rcx)
  cmp    %r15d,%r12d
  jge    L163
  movss  0x4(%r11),%xmm1
  movss  0x4(%rbx),%xmm0
  add    $0x2,%edi
  movaps %xmm1,%xmm2
  mulss  %xmm1,%xmm2
  divss  %xmm0,%xmm2
  divss  %xmm1,%xmm0
  movss  %xmm2,0x4(%r10)
  movss  %xmm0,0x4(%rcx)
  cmp    %edi,%r15d
  jle    L96
  movss  0x8(%r11),%xmm1
  movss  0x8(%rbx),%xmm0
  movaps %xmm1,%xmm2
  mulss  %xmm1,%xmm2
  divss  %xmm0,%xmm2
  divss  %xmm1,%xmm0
  movss  %xmm2,0x8(%r10)
  movss  %xmm0,0x8(%rcx)
L96:
  cmp    $0x2,%r13d
  jbe    L438
L91:
  mov    -0x48(%rbp),%esi
  pxor   %xmm2,%xmm2
  xor    %ecx,%ecx
  movaps %xmm2,%xmm0
  shr    $0x2,%esi
  shl    $0x4,%rsi
  nopl   0x0(%rax,%rax,1)
L149:
  movups (%rax,%rcx,1),%xmm1
  addss  %xmm1,%xmm0
  movaps %xmm1,%xmm3
  shufps $0x55,%xmm1,%xmm3
  addss  %xmm0,%xmm3
  movaps %xmm1,%xmm0
  unpckhps %xmm1,%xmm0
  shufps $0xff,%xmm1,%xmm1
  addss  %xmm3,%xmm0
  addss  %xmm1,%xmm0
  movups (%rdx,%rcx,1),%xmm1
  add    $0x10,%rcx
  addss  %xmm1,%xmm2
  movaps %xmm1,%xmm3
  shufps $0x55,%xmm1,%xmm3
  addss  %xmm2,%xmm3
  movaps %xmm1,%xmm2
  unpckhps %xmm1,%xmm2
  shufps $0xff,%xmm1,%xmm1
  addss  %xmm3,%xmm2
  addss  %xmm1,%xmm2
  cmp    %rcx,%rsi
  jne    L149
  mov    -0x48(%rbp),%edi
  mov    %edi,%ecx
  and    $0xfffffffc,%ecx
  and    $0x3,%edi
  je     L384
L42:
  movslq %ecx,%rsi
  add    $0x1,%ecx
  shl    $0x2,%rsi
  lea    (%rax,%rsi,1),%r8
  add    %rdx,%rsi
  addss  (%r8),%xmm0
  addss  (%rsi),%xmm2
  cmp    -0x48(%rbp),%ecx
  jge    L384
L222:
  addss  0x4(%r8),%xmm0
  addss  0x4(%rsi),%xmm2
  add    $0x1,%ecx
  cmp    -0x48(%rbp),%ecx
  jge    L384
  addss  0x8(%r8),%xmm0
  addss  0x8(%rsi),%xmm2
L384:
  pxor   %xmm3,%xmm3
  movq   -0x80(%rbp),%xmm1
  mov    -0x88(%rbp),%rax
  cvtsi2ssl -0x48(%rbp),%xmm3
  addps  %xmm1,%xmm1
  divss  %xmm3,%xmm0
  divss  %xmm3,%xmm2
  unpcklps %xmm2,%xmm0
  subps  %xmm0,%xmm1
  movlps %xmm1,(%rax)
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L258
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L641:
  pxor   %xmm6,%xmm6
  movaps %xmm1,%xmm3
  cvtsi2ss %ebx,%xmm6
  divss  %xmm6,%xmm3
  movss  %xmm6,-0x44(%rbp)
  jmp    L646
L440:
  orq    $0x0,-0x8(%rsp,%rcx,1)
  jmp    L242
L44:
  pxor   %xmm7,%xmm7
  pxor   %xmm3,%xmm3
  cvtsi2ss %r8d,%xmm7
  movss  %xmm3,-0x64(%rbp)
  divss  %xmm7,%xmm3
  movss  %xmm7,-0x44(%rbp)
  movaps %xmm3,%xmm1
  jmp    L329
L194:
  pxor   %xmm6,%xmm6
  xor    %eax,%eax
  movss  %xmm6,-0x64(%rbp)
  movaps %xmm6,%xmm1
  jmp    L166
L526:
  xor    %eax,%eax
  pxor   %xmm1,%xmm1
  jmp    L374
L616:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L451
  nopl   0x0(%rax,%rax,1)
L227:
  xor    %eax,%eax
  movaps %xmm5,%xmm1
  jmp    L307
L633:
  xor    %eax,%eax
  movaps %xmm3,%xmm1
  jmp    L565
L563:
  pxor   %xmm2,%xmm2
  movaps %xmm2,%xmm0
  jmp    L384
L394:
  lea    -0x1(%rbx),%eax
  mov    %eax,-0x68(%rbp)
  jmp    L326
L216:
  xor    %edi,%edi
  jmp    L623
L438:
  movss  -0x64(%rbp),%xmm7
  mov    $0x1,%ecx
  movss  0x0(,%r9,4),%xmm0
  addss  %xmm7,%xmm0
  addss  0x0(,%r14,4),%xmm7
  movaps %xmm7,%xmm2
  jmp    L222
L163:
  cmp    $0x2,%r13d
  ja     L91
  pxor   %xmm2,%xmm2
  xor    %ecx,%ecx
  movaps %xmm2,%xmm0
  jmp    L42
L258:
  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;
}
```