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

```asm
<main>:
  endbr64
  push   %r12
  xor    %eax,%eax
  mov    %rsi,%r12
  push   %rbp
  push   %rbx
  mov    %edi,%ebx
  sub    $0x10,%rsp
  call   <omp_get_max_threads@plt>
  mov    %eax,%ebp
  cmp    $0x1,%ebx
  jg     L28
  pxor   %xmm1,%xmm1
  movsd  D17(%rip),%xmm0
  mov    $0x3b9aca00,%ebx
  lea    D59(%rip),%rdi
  cvtsi2sd %eax,%xmm1
  call   <printf@plt>
  xor    %eax,%eax
  call   <omp_get_wtime@plt>
  movsd  D58(%rip),%xmm3
  movsd  D31(%rip),%xmm4
  movsd  %xmm0,0x8(%rsp)
L27:
  movsd  D41(%rip),%xmm6
  movsd  D18(%rip),%xmm5
  pxor   %xmm1,%xmm1
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L13:
  pxor   %xmm0,%xmm0
  movapd %xmm5,%xmm2
  cvtsi2sd %rax,%xmm0
  add    $0x1,%rax
  addsd  %xmm6,%xmm0
  mulsd  %xmm3,%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm4,%xmm0
  divsd  %xmm0,%xmm2
  addsd  %xmm2,%xmm1
  cmp    %rbx,%rax
  jl     L13
L29:
  mulsd  %xmm1,%xmm3
  xor    %eax,%eax
  movsd  %xmm3,(%rsp)
  call   <omp_get_wtime@plt>
  movsd  (%rsp),%xmm3
  lea    D15(%rip),%rdi
  movapd %xmm0,%xmm1
  subsd  0x8(%rsp),%xmm1
  movapd %xmm3,%xmm0
  call   <printf@plt>
  add    $0x10,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
L28:
  mov    0x8(%r12),%rdi
  call   <atol@plt>
  pxor   %xmm2,%xmm2
  pxor   %xmm1,%xmm1
  lea    D59(%rip),%rdi
  cvtsi2sd %rax,%xmm2
  mov    %rax,%rbx
  cvtsi2sd %ebp,%xmm1
  movapd %xmm2,%xmm0
  movsd  %xmm2,(%rsp)
  call   <printf@plt>
  xor    %eax,%eax
  call   <omp_get_wtime@plt>
  movsd  D31(%rip),%xmm4
  movsd  (%rsp),%xmm2
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm4,%xmm3
  divsd  %xmm2,%xmm3
  test   %rbx,%rbx
  jg     L27
  pxor   %xmm1,%xmm1
  jmp    L29
```

tool_calls: name: disassemble_data
arguments: data_label: D17
type_guessed: f32[2]
name: disassemble_data
arguments: data_label: D59
type_guessed: string
name: disassemble_data
arguments: data_label: D58
type_guessed: byte[8]
name: disassemble_data
arguments: data_label: D31
type_guessed: f64
name: disassemble_data
arguments: data_label: D41
type_guessed: f64
name: disassemble_data
arguments: data_label: D18
type_guessed: f64
name: disassemble_data
arguments: data_label: D15
type_guessed: string

tool: {"D17": [0.0,25.725290298461914]}

tool: D59: Calculating PI using:
  %ld slices
  %d thread(s)


tool: {"D58": [0x95,0xd6,0x26,0xe8,0x0b,0x2e,0x11,0x3e]}

tool: D31: 1

tool: D41: 0.5

tool: D18: 4

tool: D15: Obtained value for PI: %.16g
Time taken:            %.16g seconds


assistant: ```c
int main(int argc, char **argv) {
  long i, num_steps = 1000000000;
  double step, x, sum, pi, taken;
  double start, stop;
  int num_threads = omp_get_max_threads();
  if (argc > 1) {
    num_steps = atol(argv[1]);
  }
  printf("Calculating PI using:\n"
         "  %ld slices\n"
         "  %d thread(s)\n",
         num_steps, num_threads);
  start = omp_get_wtime();
  sum = 0.0;
  step = 1.0 / num_steps;
#pragma omp parallel for private(x) reduction(+ : sum)
  for (i = 0; i < num_steps; i++) {
    x = (i + 0.5) * step;
    sum += 4.0 / (1.0 + x * x);
  }
  pi = sum * step;
  stop = omp_get_wtime();
  taken = ((double)(stop - start));
  printf("Obtained value for PI: %.16g\n"
         "Time taken:            %.16g seconds\n",
         pi, taken);
  return 0;
}
```