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

```asm
<jacobian_fd>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  mov    %rdi,%rbp
  lea    0x0(,%rdx,8),%edi
  push   %rbx
  mov    %edx,%ebx
  sub    $0x58,%rsp
  mov    %rcx,0x20(%rsp)
  call   <malloc@plt>
  mov    %rax,%r14
  test   %ebx,%ebx
  jle    L56
  lea    0x8(%r12),%rax
  cmp    %rax,%r14
  je     L67
  lea    -0x1(%rbx),%eax
  cmp    $0x2,%eax
  jbe    L67
  mov    %ebx,%r13d
  mov    %ebx,%eax
  mov    %r12,0x18(%rsp)
  movsd  D11(%rip),%xmm4
  shr    %r13d
  shl    $0x3,%rax
  shl    $0x4,%r13
  mov    %rax,0x28(%rsp)
  lea    0x0(%r13,%r14,1),%rax
  mov    %r12,%r13
  movsd  %xmm4,0x38(%rsp)
  mov    %rax,0x48(%rsp)
  mov    %ebx,%eax
  sub    %r14,%r13
  and    $0xfffffffe,%eax
  mov    %eax,%r15d
  mov    %eax,0x44(%rsp)
  shl    $0x3,%r15
  lea    (%r14,%r15,1),%rax
  lea    (%r12,%r15,1),%rcx
  mov    %ebx,%r12d
  xor    %r15d,%r15d
  mov    %rax,0x30(%rsp)
  mov    %rcx,%rbx
  xchg   %ax,%ax
L74:
  mov    0x18(%rsp),%rax
  movsd  0x38(%rsp),%xmm2
  mov    %r14,%rdx
  mulsd  (%rax,%r15,1),%xmm2
  mov    0x48(%rsp),%rax
  nopl   0x0(%rax)
L65:
  movupd (%rdx,%r13,1),%xmm3
  add    $0x10,%rdx
  movups %xmm3,-0x10(%rdx)
  cmp    %rax,%rdx
  jne    L65
  cmp    0x44(%rsp),%r12d
  je     L37
  movsd  (%rbx),%xmm0
  mov    0x30(%rsp),%rax
  movsd  %xmm0,(%rax)
L37:
  lea    (%r14,%r15,1),%rdx
  movapd %xmm2,%xmm0
  mov    %r12d,%esi
  mov    %r14,%rdi
  addsd  (%rdx),%xmm0
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm0,(%rdx)
  call   *%rbp
  mov    0x18(%rsp),%rdi
  mov    %r12d,%esi
  movsd  %xmm0,0x8(%rsp)
  call   *%rbp
  movsd  0x8(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm2
  mov    0x20(%rsp),%rax
  subsd  %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax,%r15,1)
  add    $0x8,%r15
  cmp    0x28(%rsp),%r15
  jne    L74
L56:
  add    $0x58,%rsp
  mov    %r14,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1060 <free@plt>
L67:
  mov    %ebx,%r13d
  mov    %r12,%r15
  movsd  D11(%rip),%xmm5
  lea    0x0(,%r13,8),%rax
  sub    %r14,%r15
  xor    %r13d,%r13d
  lea    (%r14,%rax,1),%rdx
  mov    %rax,0x28(%rsp)
  mov    %r15,%rax
  mov    %r13,%r15
  movsd  %xmm5,0x38(%rsp)
  mov    %r14,%r13
  mov    %rdx,%r14
  mov    %rax,%rdx
  nopl   0x0(%rax,%rax,1)
L24:
  movsd  0x38(%rsp),%xmm2
  mulsd  (%r12,%r15,1),%xmm2
  mov    %r13,%rax
  nopw   0x0(%rax,%rax,1)
L12:
  movsd  (%rdx,%rax,1),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,-0x8(%rax)
  cmp    %r14,%rax
  jne    L12
  lea    0x0(%r13,%r15,1),%rax
  movapd %xmm2,%xmm0
  mov    %rdx,0x18(%rsp)
  mov    %ebx,%esi
  addsd  (%rax),%xmm0
  movsd  %xmm2,0x10(%rsp)
  mov    %r13,%rdi
  movsd  %xmm0,(%rax)
  call   *%rbp
  mov    %ebx,%esi
  mov    %r12,%rdi
  movsd  %xmm0,0x8(%rsp)
  call   *%rbp
  movsd  0x8(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm2
  mov    0x20(%rsp),%rax
  mov    0x18(%rsp),%rdx
  subsd  %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax,%r15,1)
  add    $0x8,%r15
  cmp    0x28(%rsp),%r15
  jne    L24
  mov    %r13,%r14
  jmp    L56
```

tool_calls: name: read_data
arguments: data_tag: D11
type_infered: byte[8]

tool: {"D11": [0x40,0xeb,0xd6,0x89,0x33,0xfa,0x60,0x3e]}

assistant: ```c
void jacobian_fd(double (*funcpt)(double *, int), double *x, int N, double *f) {
  int i, j;
  double step, fd;
  double *xi;
  fd = sqrt((double)1e-15);
  xi = (double *)malloc(sizeof(double) * N);
  for (i = 0; i < N; ++i) {
    step = fd * x[i];
    for (j = 0; j < N; ++j) {
      xi[j] = x[j];
    }
    xi[i] += step;
    f[i] = (funcpt(xi, N) - funcpt(x, N)) / step;
  }
  free(xi);
}
```