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

```asm
<umtz_cell_metric>:
  endbr64
  push   %r13
  pxor   %xmm13,%xmm13
  pxor   %xmm0,%xmm0
  push   %r12
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0xa8,%rsp
  lea    0x98(%rsp),%r12
  cvtss2sd (%rdi),%xmm13
  cvtss2sd 0xc(%rdi),%xmm0
  lea    0x90(%rsp),%r13
  mulsd  D34(%rip),%xmm0
  mov    %r13,%rsi
  mov    %r12,%rdi
  movsd  %xmm13,0x38(%rsp)
  divsd  D40(%rip),%xmm0
  call   <sincos@plt>
  mov    %r13,%rsi
  movsd  0x90(%rsp),%xmm5
  movsd  0x98(%rsp),%xmm7
  pxor   %xmm0,%xmm0
  mov    %r12,%rdi
  pxor   %xmm4,%xmm4
  cvtss2sd 0x10(%rbx),%xmm0
  cvtss2sd 0x4(%rbx),%xmm4
  mulsd  D34(%rip),%xmm0
  movsd  %xmm5,0x30(%rsp)
  movsd  %xmm7,0x28(%rsp)
  divsd  D40(%rip),%xmm0
  movsd  %xmm4,0x20(%rsp)
  call   <sincos@plt>
  pxor   %xmm0,%xmm0
  mov    %r13,%rsi
  mov    %r12,%rdi
  movsd  0x90(%rsp),%xmm12
  cvtss2sd 0x14(%rbx),%xmm0
  pxor   %xmm14,%xmm14
  movsd  0x98(%rsp),%xmm8
  mulsd  D34(%rip),%xmm0
  cvtss2sd 0x8(%rbx),%xmm14
  movsd  %xmm14,0x8(%rsp)
  movsd  %xmm12,0x18(%rsp)
  divsd  D40(%rip),%xmm0
  movsd  %xmm8,0x10(%rsp)
  call   <sincos@plt>
  movsd  0x30(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm12
  movsd  D57(%rip),%xmm2
  movsd  0x38(%rsp),%xmm13
  movapd %xmm5,%xmm0
  movsd  0x20(%rsp),%xmm4
  movsd  0x90(%rsp),%xmm11
  movapd %xmm12,%xmm9
  mulsd  %xmm5,%xmm0
  movapd %xmm13,%xmm1
  movsd  0x8(%rsp),%xmm14
  movapd %xmm5,%xmm10
  mulsd  %xmm4,%xmm1
  movapd %xmm11,%xmm3
  movsd  0x98(%rsp),%xmm6
  movsd  0x10(%rsp),%xmm8
  movsd  0x28(%rsp),%xmm7
  unpcklpd %xmm5,%xmm9
  unpcklpd %xmm12,%xmm10
  unpcklpd %xmm3,%xmm3
  subsd  %xmm0,%xmm2
  movapd %xmm12,%xmm0
  mulsd  %xmm12,%xmm0
  movapd %xmm1,%xmm15
  mulsd  %xmm14,%xmm15
  subsd  %xmm0,%xmm2
  movapd %xmm11,%xmm0
  mulsd  %xmm11,%xmm0
  subsd  %xmm0,%xmm2
  movapd %xmm5,%xmm0
  addsd  %xmm5,%xmm0
  mulsd  %xmm12,%xmm0
  mulsd  %xmm11,%xmm0
  addsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  ucomisd %xmm0,%xmm2
  ja     L17
  sqrtsd %xmm0,%xmm0
L16:
  mulsd  %xmm14,%xmm4
  mulsd  %xmm14,%xmm13
  mulsd  %xmm6,%xmm1
  unpcklpd %xmm6,%xmm6
  mulsd  %xmm15,%xmm0
  mulsd  %xmm12,%xmm5
  mulsd  %xmm7,%xmm4
  mulsd  %xmm8,%xmm13
  mulpd  %xmm10,%xmm3
  unpcklpd %xmm0,%xmm0
  subsd  %xmm11,%xmm5
  unpcklpd %xmm1,%xmm4
  divpd  %xmm0,%xmm4
  unpcklpd %xmm13,%xmm1
  subpd  %xmm9,%xmm3
  divpd  %xmm0,%xmm1
  movapd %xmm4,%xmm2
  addsd  %xmm4,%xmm2
  movapd %xmm1,%xmm0
  movapd %xmm1,%xmm13
  mulsd  %xmm1,%xmm13
  unpckhpd %xmm0,%xmm0
  addpd  %xmm1,%xmm1
  movapd %xmm0,%xmm14
  movapd %xmm4,%xmm0
  mulsd  %xmm14,%xmm2
  movapd %xmm14,%xmm15
  mulpd  %xmm4,%xmm1
  mulsd  %xmm4,%xmm0
  mulsd  %xmm14,%xmm15
  mulsd  %xmm5,%xmm2
  movapd %xmm7,%xmm5
  unpcklpd %xmm8,%xmm7
  mulsd  %xmm8,%xmm5
  mulpd  %xmm1,%xmm3
  mulpd  %xmm7,%xmm6
  unpcklpd %xmm15,%xmm0
  cvtpd2ps %xmm0,%xmm0
  divsd  %xmm5,%xmm2
  divpd  %xmm6,%xmm3
  unpcklpd %xmm2,%xmm13
  cvtpd2ps %xmm13,%xmm13
  movlhps %xmm13,%xmm0
  movups %xmm0,0x0(%rbp)
  cvtpd2ps %xmm3,%xmm3
  movlps %xmm3,0x10(%rbp)
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L17:
  movaps %xmm10,0x80(%rsp)
  movaps %xmm3,0x70(%rsp)
  movaps %xmm9,0x60(%rsp)
  movsd  %xmm6,0x58(%rsp)
  movsd  %xmm11,0x50(%rsp)
  movsd  %xmm8,0x48(%rsp)
  movsd  %xmm12,0x40(%rsp)
  movsd  %xmm7,0x38(%rsp)
  movsd  %xmm5,0x30(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm13,0x18(%rsp)
  movsd  %xmm15,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  call   <sqrt@plt>
  movapd 0x70(%rsp),%xmm3
  movapd 0x80(%rsp),%xmm10
  movapd 0x60(%rsp),%xmm9
  movsd  0x58(%rsp),%xmm6
  movsd  0x50(%rsp),%xmm11
  movsd  0x48(%rsp),%xmm8
  movsd  0x40(%rsp),%xmm12
  movsd  0x38(%rsp),%xmm7
  movsd  0x30(%rsp),%xmm5
  movsd  0x28(%rsp),%xmm14
  movsd  0x20(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm13
  movsd  0x10(%rsp),%xmm15
  movsd  0x8(%rsp),%xmm1
  jmp    L16
```

tool_calls: name: retrieve_data
arguments: label: D34
type: f64
name: retrieve_data
arguments: label: D40
type: f64
name: retrieve_data
arguments: label: D57
type: f64

tool: D34: 3.141592653589793

tool: D40: 180

tool: D57: 1

assistant: ```c
void umtz_cell_metric(const float *cell, float *metric) {
  double a, b, c, alph, beta, gamm, as, bs, cs, vol;
  a = cell[0];
  alph = cell[3] * (3.1415926535897932384626433832795028841972) / 180.0;
  b = cell[1];
  beta = cell[4] * (3.1415926535897932384626433832795028841972) / 180.0;
  c = cell[2];
  gamm = cell[5] * (3.1415926535897932384626433832795028841972) / 180.0;
  vol = a * b * c * sqrt(1.0 - cos(alph) * cos(alph) - cos(beta) * cos(beta) - cos(gamm) * cos(gamm) + 2.0 * cos(alph) * cos(beta) * cos(gamm));
  as = b * c * sin(alph) / vol;
  bs = c * a * sin(beta) / vol;
  cs = a * b * sin(gamm) / vol;
  metric[0] = (float)(as * as);
  metric[1] = (float)(bs * bs);
  metric[2] = (float)(cs * cs);
  metric[3] = (float)(2.0 * as * bs * (cos(beta) * cos(alph) - cos(gamm)) / (sin(alph) * sin(beta)));
  metric[4] = (float)(2.0 * cs * as * (cos(alph) * cos(gamm) - cos(beta)) / (sin(gamm) * sin(alph)));
  metric[5] = (float)(2.0 * bs * cs * (cos(gamm) * cos(beta) - cos(alph)) / (sin(beta) * sin(gamm)));
}
```