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

```asm
<CkNNStochasticProximityEmbedding>:
  endbr64
  push   %r15
  mov    %rdi,%r15
  xor    %edi,%edi
  push   %r14
  mov    %r8d,%r14d
  push   %r13
  mov    %rsi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  movq   %xmm0,%rbx
  sub    $0x88,%rsp
  mov    0xc0(%rsp),%eax
  mov    %edx,0x38(%rsp)
  mov    %ecx,0x58(%rsp)
  mov    %eax,0x3c(%rsp)
  mov    0xc8(%rsp),%eax
  mov    %r9d,0x68(%rsp)
  movsd  %xmm1,0x8(%rsp)
  mov    %eax,0x6c(%rsp)
  call   <time@plt>
  mov    %rax,%r12
  xor    %eax,%eax
  call   <getpid@plt>
  mov    %eax,%ebp
  xor    %eax,%eax
  call   <getpid@plt>
  mov    %ebp,%edi
  imul   %eax,%edi
  movslq %edi,%rdi
  add    %r12,%rdi
  call   <srand@plt>
  mov    0x58(%rsp),%ebp
  mov    0x38(%rsp),%edx
  mov    %r15,%rdi
  mov    %ebp,%esi
  call   <nearest_neighbours@plt>
  mov    %rax,0x60(%rsp)
  mov    %ebp,%eax
  imul   %r14d,%eax
  test   %eax,%eax
  jle    L112
  sub    $0x1,%eax
  mov    %r13,%r12
  lea    0x8(%r13,%rax,8),%rbp
L159:
  xor    %eax,%eax
  add    $0x8,%r12
  call   <rand@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x8(%r12)
  cmp    %rbp,%r12
  jne    L159
L112:
  mov    0x68(%rsp),%edx
  test   %edx,%edx
  jle    L60
  movq   %rbx,%xmm7
  subsd  0x8(%rsp),%xmm7
  mov    0x68(%rsp),%eax
  mov    %rbx,0x30(%rsp)
  movslq 0x38(%rsp),%rcx
  mov    0x60(%rsp),%rsi
  movl   $0x0,0x5c(%rsp)
  sub    $0x1,%eax
  movsd  %xmm7,0x78(%rsp)
  pxor   %xmm7,%xmm7
  cvtsi2sd %eax,%xmm7
  mov    %rcx,0x48(%rsp)
  lea    (%rsi,%rcx,4),%rcx
  mov    %rcx,0x50(%rsp)
  movsd  %xmm7,0x70(%rsp)
L18:
  mov    0x60(%rsp),%rax
  movl   $0x0,0x1c(%rsp)
  mov    %r14d,%r12d
  add    $0x4,%rax
  mov    %rax,0x40(%rsp)
  mov    0x3c(%rsp),%eax
  test   %eax,%eax
  jle    L62
  nopl   (%rax)
L90:
  xor    %eax,%eax
  call   <rand@plt>
  mov    0x38(%rsp),%ecx
  cltd
  idivl  0x58(%rsp)
  lea    0x1(%rdx),%eax
  mov    %edx,%ebp
  imul   %ecx,%eax
  cmp    %eax,%ecx
  jge    L74
  mov    %edx,%esi
  sub    0x38(%rsp),%eax
  mov    0x50(%rsp),%r14
  imul   %r12d,%esi
  sub    $0x1,%eax
  add    0x48(%rsp),%rax
  mov    %esi,0x20(%rsp)
  mov    0x40(%rsp),%rsi
  lea    (%rsi,%rax,4),%rax
  mov    %rax,0x10(%rsp)
  lea    -0x1(%r12),%eax
  mov    %eax,0x24(%rsp)
  lea    0x8(%r13),%rax
  mov    %rax,0x28(%rsp)
  mov    %r12d,%eax
  mov    %r14,%r12
  mov    %r13,%r14
  mov    %eax,%r13d
  jmp    L89
  nopl   0x0(%rax)
L42:
  add    $0x4,%r12
  cmp    0x10(%rsp),%r12
  je     L72
L89:
  mov    (%r12),%ebx
  mov    %r13d,%ecx
  mov    %ebp,%esi
  mov    %r14,%rdi
  mov    %ebx,%edx
  call   <ed@plt>
  mov    %ebx,%esi
  mov    %ebp,%edi
  movsd  %xmm0,0x8(%rsp)
  call   <trmIndex@plt>
  movsd  0x8(%rsp),%xmm0
  movsd  (%r15,%rax,8),%xmm1
  comisd %xmm1,%xmm0
  jbe    L42
  movsd  D8(%rip),%xmm3
  imul   %r13d,%ebx
  subsd  %xmm0,%xmm1
  mulsd  0x30(%rsp),%xmm3
  addsd  D7(%rip),%xmm0
  mulsd  %xmm1,%xmm3
  divsd  %xmm0,%xmm3
  test   %r13d,%r13d
  jle    L42
  mov    0x24(%rsp),%edx
  movslq %ebx,%rsi
  mov    0x28(%rsp),%rcx
  lea    (%r14,%rsi,8),%rax
  add    %rsi,%rdx
  lea    (%rcx,%rdx,8),%rdi
  movslq 0x20(%rsp),%rdx
  sub    %rsi,%rdx
  nopl   (%rax)
L32:
  movsd  (%rax,%rdx,8),%xmm1
  movapd %xmm1,%xmm0
  subsd  (%rax),%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rdx,8)
  movsd  (%rax),%xmm2
  add    $0x8,%rax
  movapd %xmm2,%xmm1
  subsd  %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rdi
  jne    L32
  add    $0x4,%r12
  cmp    0x10(%rsp),%r12
  jne    L89
L72:
  mov    %r13d,%r12d
  mov    %r14,%r13
L74:
  addl   $0x1,0x1c(%rsp)
  mov    0x1c(%rsp),%eax
  cmp    %eax,0x3c(%rsp)
  jne    L90
  mov    %r12d,%r14d
L62:
  mov    0x5c(%rsp),%eax
  mov    0x6c(%rsp),%ecx
  cltd
  idiv   %ecx
  test   %edx,%edx
  jne    L76
  mov    0x5c(%rsp),%ebx
  test   %ebx,%ebx
  je     L76
  test   %ecx,%ecx
  jle    L76
  mov    0x58(%rsp),%ecx
  mov    %r13,%rsi
  mov    %r15,%rdi
  mov    %r14d,%edx
  call   <stress@plt>
  movsd  0x30(%rsp),%xmm1
  mov    %ebx,%esi
  lea    D22(%rip),%rdi
  call   <printf@plt>
L76:
  movsd  0x30(%rsp),%xmm5
  movsd  0x78(%rsp),%xmm0
  divsd  0x70(%rsp),%xmm0
  addl   $0x1,0x5c(%rsp)
  mov    0x5c(%rsp),%eax
  subsd  %xmm0,%xmm5
  movsd  %xmm5,0x30(%rsp)
  cmp    %eax,0x68(%rsp)
  jne    L18
L60:
  mov    0x60(%rsp),%rdi
  call   <free@plt>
  mov    0x58(%rsp),%ecx
  mov    %r14d,%edx
  mov    %r13,%rsi
  add    $0x88,%rsp
  mov    %r15,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1110 <stress@plt>
```

tool_calls: name: disassemble_data
arguments: data_label: D8
type_infered: f64
name: disassemble_data
arguments: data_label: D7
type_infered: byte[8]
name: disassemble_data
arguments: data_label: D22
type_infered: string

tool: D8: 0.5

tool: {"D7": [0x3a,0x8c,0x30,0xe2,0x8e,0x79,0x45,0x3e]}

tool: D22: Cycle %d - Residual stress: %.3f, lambda %.3f


assistant: ```c
double CkNNStochasticProximityEmbedding(double *s, double *d_coords, int kn, int nelem, int dim, double maxlam, double minlam, int ncycle, int nstep, int stressfreq) {
  int *tmp;
  int *neighbours;
  int a = 0, b = 0, idx = 0, idxa = 0, idxb = 0, idxak = 0, idxbk = 0, idxab = 0;
  double dab = 0.0, rab = 0.0;
  double lam = maxlam;
  double t = 0.0;
  double finalstress = 0.0;
  srand(time(NULL) + getpid() * getpid());
  neighbours = nearest_neighbours(s, nelem, kn);
  for (int i = 0; i < nelem * dim; i++) {
    d_coords[i] = (double)rand();
  }
  for (int i = 0; i < ncycle; i++) {
    for (int j = 0; j < nstep; j++) {
      a = rand() % nelem;
      for (int k = kn; k < (a + 1) * kn; k++) {
        b = neighbours[k];
        dab = ed(d_coords, a, b, dim);
        rab = s[trmIndex(a, b)];
        if (dab > rab) {
          idxa = a * dim;
          idxb = b * dim;
          t = lam * 0.5 * (rab - dab) / (dab + 1e-8);
          for (int k = 0; k < dim; k++) {
            idxak = idxa + k;
            idxbk = idxb + k;
            d_coords[idxak] = d_coords[idxak] + t * (d_coords[idxak] - d_coords[idxbk]);
            d_coords[idxbk] = d_coords[idxbk] + t * (d_coords[idxbk] - d_coords[idxak]);
          }
        }
      }
    }
    if (i % stressfreq == 0 && i != 0 && stressfreq > 0)
      printf("Cycle %d - Residual stress: %.3f, lambda %.3f\n", i, stress(s, d_coords, dim, nelem), lam);
    lam = lam - (maxlam - minlam) / (double)(ncycle - 1);
  }
  free(neighbours);
  return (stress(s, d_coords, dim, nelem));
}
```