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

```asm
<mc_sampling>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xf0,%rsp
  mov    %edi,-0xa4(%rbp)
  mov    %esi,-0xa8(%rbp)
  mov    %edx,-0xac(%rbp)
  mov    %ecx,-0xb0(%rbp)
  mov    %r8d,-0xb4(%rbp)
  mov    %r9d,-0xb8(%rbp)
  movsd  %xmm0,-0xc0(%rbp)
  movsd  %xmm1,-0xc8(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xd0(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xd8(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0xe0(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D75(%rip),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xac(%rbp),%xmm0
  divsd  -0x68(%rbp),%xmm0
  movq   %xmm0,%rax
  mov    $0x0,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x60(%rbp)
  movsd  -0x60(%rbp),%xmm0
  movsd  D59(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  D33(%rip),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  D30(%rip),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0xc0(%rbp),%xmm0
  subsd  -0xc8(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movq   $0xffffffffffffffff,-0x88(%rbp)
  mov    -0xa4(%rbp),%ecx
  mov    -0xac(%rbp),%eax
  mov    $0x8,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <matrix@plt>
  mov    %rax,-0x40(%rbp)
  mov    -0xa4(%rbp),%ecx
  mov    -0xac(%rbp),%eax
  mov    $0x8,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <matrix@plt>
  mov    %rax,-0x38(%rbp)
  movl   $0x0,-0x90(%rbp)
  jmp    L118
L276:
  movl   $0x0,-0x8c(%rbp)
  jmp    L0
L1:
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x90(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0x40(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x8c(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movsd  (%rax),%xmm0
  movsd  %xmm0,(%rdx)
  addl   $0x1,-0x8c(%rbp)
L0:
  mov    -0x8c(%rbp),%eax
  cmp    -0xa4(%rbp),%eax
  jl     L1
  addl   $0x1,-0x90(%rbp)
L118:
  mov    -0x90(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jl     L276
  movl   $0x0,-0x90(%rbp)
  jmp    L27
L79:
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x90(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0xd8(%rbp),%rdx
  add    %rcx,%rdx
  movsd  (%rax),%xmm0
  movsd  %xmm0,(%rdx)
  addl   $0x1,-0x90(%rbp)
L27:
  mov    -0xb0(%rbp),%eax
  cmp    -0x90(%rbp),%eax
  jge    L79
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movl   $0x1,-0x98(%rbp)
  jmp    L216
L277:
  movsd  -0x70(%rbp),%xmm0
  addsd  -0xc8(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0x70(%rbp),%rax
  mov    $0x5,%edi
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movl   $0x0,-0x94(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movl   $0x0,-0x90(%rbp)
  jmp    L113
L15:
  movl   $0x0,-0x8c(%rbp)
  jmp    L52
L36:
  lea    -0x88(%rbp),%rax
  mov    %rax,%rdi
  call   <ran1@plt>
  movq   %xmm0,%rax
  mov    -0x90(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0x40(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x8c(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  movq   %rax,%xmm0
  mulsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,(%rdx)
  addl   $0x1,-0x8c(%rbp)
L52:
  mov    -0x8c(%rbp),%eax
  cmp    -0xa4(%rbp),%eax
  jl     L36
  addl   $0x1,-0x90(%rbp)
L113:
  mov    -0x90(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jl     L15
  mov    -0xac(%rbp),%edx
  mov    -0xa4(%rbp),%esi
  movsd  -0x60(%rbp),%xmm0
  mov    -0x28(%rbp),%rcx
  mov    -0x40(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rcx,%xmm0
  mov    %rax,%rdi
  call   <wave_function@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x80(%rbp)
  movl   $0x0,-0xa0(%rbp)
  jmp    L181
L54:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movl   $0x1,-0x9c(%rbp)
  jmp    L134
L104:
  movl   $0x0,-0x90(%rbp)
  jmp    L286
L178:
  movl   $0x0,-0x8c(%rbp)
  jmp    L132
L28:
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movsd  %xmm2,-0xe8(%rbp)
  lea    -0x88(%rbp),%rax
  mov    %rax,%rdi
  call   <ran1@plt>
  movq   %xmm0,%rax
  movsd  D29(%rip),%xmm1
  movq   %rax,%xmm0
  subsd  %xmm1,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  addsd  -0xe8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  jbe    L74
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  addsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L74:
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L229
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  subsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L229:
  addl   $0x1,-0x8c(%rbp)
L132:
  mov    -0x8c(%rbp),%eax
  cmp    -0xa4(%rbp),%eax
  jl     L28
  addl   $0x1,-0x90(%rbp)
L286:
  mov    -0x90(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jl     L178
  mov    -0xac(%rbp),%edx
  mov    -0xa4(%rbp),%esi
  movsd  -0x60(%rbp),%xmm0
  mov    -0x28(%rbp),%rcx
  mov    -0x38(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rcx,%xmm0
  mov    %rax,%rdi
  call   <wave_function@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  -0x80(%rbp),%xmm0
  cvttsd2si %xmm0,%edx
  mov    -0x10(%rbp),%rax
  mov    %edx,%esi
  movq   %rax,%xmm0
  lea    D78(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  lea    -0x88(%rbp),%rax
  mov    %rax,%rdi
  call   <ran1@plt>
  movsd  -0x10(%rbp),%xmm1
  divsd  -0x80(%rbp),%xmm1
  comisd %xmm0,%xmm1
  jb     L117
  movl   $0x0,-0x90(%rbp)
  jmp    L108
L289:
  movl   $0x0,-0x8c(%rbp)
  jmp    L212
L159:
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x90(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x40(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x8c(%rbp)
L212:
  mov    -0x8c(%rbp),%eax
  cmp    -0xa4(%rbp),%eax
  jl     L159
  addl   $0x1,-0x90(%rbp)
L108:
  mov    -0x90(%rbp),%eax
  cmp    -0xac(%rbp),%eax
  jl     L289
  movsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  addl   $0x1,-0x94(%rbp)
L117:
  mov    -0x9c(%rbp),%eax
  cmp    -0xb4(%rbp),%eax
  jle    L140
  movsd  -0x78(%rbp),%xmm0
  addsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x78(%rbp)
L140:
  addl   $0x1,-0x9c(%rbp)
L134:
  mov    -0xb8(%rbp),%edx
  mov    -0xb4(%rbp),%eax
  add    %edx,%eax
  cmp    %eax,-0x9c(%rbp)
  jle    L104
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xb8(%rbp),%xmm1
  movsd  -0x78(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x78(%rbp),%xmm0
  cvttsd2si %xmm0,%edx
  pxor   %xmm3,%xmm3
  cvtsi2sdl -0xa0(%rbp),%xmm3
  movq   %xmm3,%rax
  mov    %edx,%esi
  movq   %rax,%xmm0
  lea    D0(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  cmpl   $0x0,-0xa0(%rbp)
  je     L24
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x78(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L24:
  addl   $0x1,-0xa0(%rbp)
L181:
  mov    -0xa0(%rbp),%eax
  cmp    -0xa8(%rbp),%eax
  jle    L54
  mov    -0x94(%rbp),%edx
  mov    -0x70(%rbp),%rax
  mov    %edx,%esi
  movq   %rax,%xmm0
  lea    D74(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xa8(%rbp),%eax
  imul   -0xac(%rbp),%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xa8(%rbp),%eax
  imul   -0xac(%rbp),%eax
  imul   -0xac(%rbp),%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mov    -0x98(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x98(%rbp)
L216:
  mov    -0x98(%rbp),%eax
  cmp    -0xb0(%rbp),%eax
  jle    L277
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L248
  call   <__stack_chk_fail@plt>
L248:
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_label: D75
type: byte[8]
name: parse_data
arguments: data_label: D59
type: f32[2]
name: parse_data
arguments: data_label: D33
type: f64
name: parse_data
arguments: data_label: D30
type: f64
name: parse_data
arguments: data_label: D29
type: f64
name: parse_data
arguments: data_label: D78
type: string
name: parse_data
arguments: data_label: D0
type: string
name: parse_data
arguments: data_label: D74
type: string

tool: {"D75": [0xa5,0xda,0xa7,0xe3,0x31,0x03,0xd5,0x3f]}

tool: {"D59": [0.0,2.25]}

tool: D33: 9.27647

tool: D30: 40.88

tool: D29: 0.5

tool: D78: wfnew/wfold %12.5E%12.5E


tool: D0: Block %d energy %12.5E


tool: D74: variational parameter %12.5E accepted steps %d


assistant: ```c
void mc_sampling(int dimension, int number_blocks, int number_particles, int max_variations, int thermalization, int number_cycles, double initial_correlation, double delta_correlation, double *lowest_correlation, double *cumulative_e, double *cumulative_e2) {
  int blocks, cycles, variate, accept, i, j;
  long idum;
  double wfnew, wfold, correlation5, energy, energy2, delta_e, lowest_energy;
  double correlation, **r_old, **r_new;
  double l, kinetic_factor, density, scale, four_epsilon;
  density = 0.3648 * 0.9;
  l = pow((number_particles / density), 0.3333333);
  scale = l / 4;
  kinetic_factor = 9.27647;
  four_epsilon = 40.88;
  correlation = initial_correlation - delta_correlation;
  idum = -1;
  r_old = (double **)matrix(number_particles, dimension, sizeof(**r_old));
  r_new = (double **)matrix(number_particles, dimension, sizeof(**r_new));
  for (i = 0; i < number_particles; i++) {
    for (j = 0; j < dimension; j++) {
      r_old[i][j] = r_new[i][j] = 0;
    }
  }
  for (i = 0; i < max_variations + 1; i++) {
    cumulative_e[i] = cumulative_e2[i] = 0;
  }
  lowest_energy = 0.;
  for (variate = 1; variate <= max_variations; variate++) {
    correlation += delta_correlation;
    correlation5 = pow(correlation, 5);
    energy = energy2 = 0;
    accept = 0;
    delta_e = 0;
    for (i = 0; i < number_particles; i++) {
      for (j = 0; j < dimension; j++) {
        r_old[i][j] = l * ran1(&idum);
      }
    }
    wfold = wave_function(r_old, correlation5, l, dimension, number_particles);
    for (blocks = 0; blocks <= number_blocks; blocks++) {
      energy = 0;
      for (cycles = 1; cycles <= number_cycles + thermalization; cycles++) {
        for (i = 0; i < number_particles; i++) {
          for (j = 0; j < dimension; j++) {
            r_new[i][j] = r_old[i][j] + scale * (ran1(&idum) - 0.5);
            if (r_new[i][j] < 0.)
              r_new[i][j] += l;
            if (r_new[i][j] > 0.)
              r_new[i][j] -= l;
          }
        }
        wfnew = wave_function(r_new, correlation5, l, dimension, number_particles);
        printf("wfnew/wfold %12.5E%12.5E\n", wfnew, wfold);
        if (ran1(&idum) <= wfnew / wfold) {
          for (i = 0; i < number_particles; i++) {
            for (j = 0; j < dimension; j++) {
              r_old[i][j] = r_new[i][j];
            }
          }
          wfold = wfnew;
          accept = accept + 1;
        }
        if (cycles > thermalization) {
          energy += delta_e;
        }
      }
      energy /= number_cycles;
      printf("Block %d energy %12.5E\n", blocks, energy);
      if (blocks != 0) {
        cumulative_e[variate] += energy;
        cumulative_e2[variate] += energy * energy;
      }
    }
    printf("variational parameter %12.5E accepted steps %d\n", correlation, accept);
    cumulative_e[variate] /= number_blocks * number_particles;
    cumulative_e2[variate] /= number_blocks * number_particles * number_particles;
  }
}
```