user: Converte this assembly snippet into its equivalent C code:

```asm
<main>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  jmp    L13
L40:
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movq   D63(%rip),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  -0x50(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  addsd  %xmm0,%xmm0
  movsd  D78(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x48(%rbp),%xmm0
  movq   D63(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x50(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  addsd  %xmm0,%xmm0
  movsd  D78(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x40(%rbp),%xmm0
  movq   D63(%rip),%xmm2
  xorpd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  -0x50(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  addsd  %xmm0,%xmm0
  movsd  D78(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x38(%rbp),%xmm0
  addsd  -0x30(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  movsd  D11(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x38(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x30(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L8
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x10(%rbp),%rax
  movq   %rax,%xmm0
  lea    D64(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  jmp    L13
L8:
  lea    D20(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
L13:
  lea    -0x40(%rbp),%rcx
  lea    -0x48(%rbp),%rdx
  lea    -0x50(%rbp),%rax
  mov    %rax,%rsi
  lea    D67(%rip),%rax
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <scanf@plt>
  cmp    $0x3,%eax
  je     L40
  mov    $0x0,%eax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L25
  call   <__stack_chk_fail@plt>
L25:
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_label: D63
type_guessed: byte[8]
name: parse_data
arguments: data_label: D78
type_guessed: byte[8]
name: parse_data
arguments: data_label: D11
type_guessed: byte[8]
name: parse_data
arguments: data_label: D64
type_guessed: string
name: parse_data
arguments: data_label: D20
type_guessed: string
name: parse_data
arguments: data_label: D67
type_guessed: string

tool: {"D63": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

tool: {"D78": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40]}

tool: {"D11": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x40]}

tool: D64: %.3lf


tool: D20: -1.000

tool: D67: %lf %lf %lf

assistant: ```c
int main() {
  double a, b, c, s, area, pre_area;
  double m1, m2, m3;
  while (scanf("%lf %lf %lf", &m1, &m2, &m3) == 3) {
    a = (2 * sqrt((2 * (m2) * (m2)) + (2 * (m3) * (m3)) + (-(m1) * (m1)))) / 3;
    b = (2 * sqrt((-(m2) * (m2)) + (2 * (m3) * (m3)) + (2 * (m1) * (m1)))) / 3;
    c = (2 * sqrt((2 * (m2) * (m2)) + (-(m3) * (m3)) + (2 * (m1) * (m1)))) / 3;
    s = (a + b + c) / 2;
    pre_area = s * (s - a) * (s - b) * (s - c);
    if (pre_area > 0) {
      area = sqrt(pre_area);
      printf("%.3lf\n", area);
    } else {
      printf("-1.000\n");
    }
  }
  return 0;
}
```