user: Please show me the C implementation that matches the assembly code listed below:

```asm
<data_from_contour>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %esi,-0x4c(%rbp)
  mov    %rdx,-0x58(%rbp)
  mov    %rcx,-0x60(%rbp)
  mov    %r8,-0x68(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movl   $0x0,-0x30(%rbp)
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x48(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm3
  movapd %xmm0,%xmm2
  subsd  %xmm3,%xmm2
  mov    -0x48(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    -0x60(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x60(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%rax
  movsd  %xmm0,(%rax)
  movl   $0x0,-0x34(%rbp)
  jmp    L7
L13:
  mov    -0x34(%rbp),%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x34(%rbp),%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x34(%rbp),%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm3
  movapd %xmm0,%xmm2
  subsd  %xmm3,%xmm2
  mov    -0x34(%rbp),%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm3
  subsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x58(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x28(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L33
  mov    -0x34(%rbp),%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x58(%rbp),%rax
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L33:
  mov    -0x60(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd -0x28(%rbp),%xmm0
  jbe    L3
  mov    -0x34(%rbp),%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x60(%rbp),%rax
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L3:
  addl   $0x1,-0x34(%rbp)
L7:
  mov    -0x34(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L13
  mov    -0x30(%rbp),%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x20(%rbp),%xmm1
  movapd %xmm0,%xmm2
  subsd  %xmm1,%xmm2
  mov    -0x30(%rbp),%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x18(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movapd %xmm2,%xmm1
  movq   %rax,%xmm0
  call   <atan2@plt>
  mov    $0x5a,%esi
  mov    %eax,%edi
  call   <fmod@plt>
  movq   %xmm0,%rax
  mov    -0x68(%rbp),%rdx
  mov    %rax,(%rdx)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L48
  call   <__stack_chk_fail@plt>
L48:
  leave
  ret
```

assistant: ```c
void data_from_contour(double *cont, int anz, double *a, double *b, double *tet) {
  int i, ai, bi;
  double d, center[2];
  ai = 0;
  *a = *b = sqrt(((cont[0] - center[0]) * (cont[0] - center[0])) + ((cont[1] - center[1]) * (cont[1] - center[1])));
  for (i = 0; i < anz; i++) {
    d = sqrt(((cont[i * 2] - center[0]) * (cont[i * 2] - center[0])) + ((cont[i * 2 + 1] - center[1]) * (cont[i * 2 + 1] - center[1])));
    if (*a < d) {
      ai = i;
      *a = d;
    }
    if (*b > d) {
      bi = i;
      *b = d;
    }
  }
  *tet = fmod(atan2(cont[ai * 2 + 1] - center[1], cont[ai * 2] - center[0]), 90);
}
```