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

```asm
<dd_exp>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  sub    $0x158,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x148(%rsp)
  xor    %eax,%eax
  mov    D466(%rip),%rax
  mov    %rax,0x50(%rsp)
  mov    D222(%rip),%rax
  mov    %rax,0x58(%rsp)
  mov    D213(%rip),%rax
  mov    %rax,0x60(%rsp)
  mov    D287(%rip),%rax
  mov    %rax,0x68(%rsp)
  mov    D400(%rip),%rax
  mov    %rax,0x70(%rsp)
  mov    D292(%rip),%rax
  mov    %rax,0x78(%rsp)
  mov    D111(%rip),%rax
  mov    %rax,0x80(%rsp)
  mov    D418(%rip),%rax
  mov    %rax,0x88(%rsp)
  mov    D207(%rip),%rax
  mov    %rax,0x90(%rsp)
  mov    D148(%rip),%rax
  mov    %rax,0x98(%rsp)
  mov    D403(%rip),%rax
  mov    %rax,0xa0(%rsp)
  mov    D333(%rip),%rax
  mov    %rax,0xa8(%rsp)
  mov    D77(%rip),%rax
  mov    %rax,0xb0(%rsp)
  mov    D110(%rip),%rax
  mov    %rax,0xb8(%rsp)
  mov    D443(%rip),%rax
  mov    %rax,0xc0(%rsp)
  mov    D118(%rip),%rax
  mov    %rax,0xc8(%rsp)
  mov    D261(%rip),%rax
  mov    %rax,0xd0(%rsp)
  mov    D444(%rip),%rax
  mov    %rax,0xd8(%rsp)
  mov    D48(%rip),%rax
  mov    %rax,0xe0(%rsp)
  mov    D147(%rip),%rax
  mov    %rax,0xe8(%rsp)
  mov    D315(%rip),%rax
  mov    %rax,0xf0(%rsp)
  mov    D262(%rip),%rax
  mov    %rax,0xf8(%rsp)
  mov    D413(%rip),%rax
  mov    %rax,0x100(%rsp)
  mov    D331(%rip),%rax
  mov    %rax,0x108(%rsp)
  mov    D339(%rip),%rax
  mov    %rax,0x110(%rsp)
  mov    D374(%rip),%rax
  mov    %rax,0x118(%rsp)
  mov    D417(%rip),%rax
  mov    %rax,0x120(%rsp)
  mov    D314(%rip),%rax
  mov    %rax,0x128(%rsp)
  mov    D53(%rip),%rax
  movsd  (%rdi),%xmm0
  movsd  D27(%rip),%xmm1
  mov    %rax,0x130(%rsp)
  mov    D154(%rip),%rax
  comisd %xmm0,%xmm1
  mov    %rax,0x138(%rsp)
  jae    L46
  comisd D104(%rip),%xmm0
  jb     L84
  movsd  D457(%rip),%xmm0
  pxor   %xmm1,%xmm1
L41:
  movsd  %xmm0,0x0(%rbp)
  movsd  %xmm1,0x8(%rbp)
L113:
  mov    0x148(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L67
  add    $0x158,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L84:
  mov    %rdi,%r12
  call   <dd_is_zero@plt>
  cmp    $0x1,%eax
  je     L40
  mov    %r12,%rdi
  call   <dd_is_one@plt>
  cmp    $0x1,%eax
  je     L93
  movsd  (%r12),%xmm0
  movsd  D13(%rip),%xmm2
  movsd  D432(%rip),%xmm1
  movsd  D112(%rip),%xmm3
  divsd  %xmm2,%xmm0
  addsd  D352(%rip),%xmm0
  movapd %xmm0,%xmm5
  movsd  %xmm0,0x8(%rsp)
  andpd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm3
  jbe    L4
  cvttsd2si %xmm5,%rax
  pxor   %xmm0,%xmm0
  movsd  D264(%rip),%xmm6
  andnpd %xmm5,%xmm1
  cvtsi2sd %rax,%xmm0
  movapd %xmm0,%xmm3
  cmpnlesd %xmm5,%xmm3
  andpd  %xmm6,%xmm3
  subsd  %xmm3,%xmm0
  orpd   %xmm1,%xmm0
  movsd  %xmm0,0x8(%rsp)
L4:
  lea    0x10(%rsp),%r14
  movsd  0x8(%rsp),%xmm0
  lea    0x40(%rsp),%rbx
  mov    D290(%rip),%rax
  mov    %r14,%rdi
  movsd  %xmm2,0x10(%rsp)
  lea    0x20(%rsp),%r15
  lea    0x60(%rsp),%r13
  mov    %rax,0x18(%rsp)
  call   <dd_mlt_d@plt>
  mov    %r14,%rdx
  mov    %r12,%rdi
  mov    %r14,%rsi
  call   <dd_sub@plt>
  mov    %r14,%rsi
  mov    %r14,%rdi
  lea    0x30(%rsp),%r12
  movsd  D146(%rip),%xmm0
  call   <dd_mul_pwr2@plt>
  mov    %rbx,%rsi
  mov    %r14,%rdi
  call   <dd_sqr@plt>
  mov    D352(%rip),%rax
  mov    %r15,%rsi
  mov    %rbx,%rdi
  movq   %rax,%xmm0
  call   <dd_mul_pwr2@plt>
  mov    %r14,%rsi
  mov    %r15,%rdi
  call   <dd_inc@plt>
  mov    %r14,%rsi
  mov    %rbx,%rdi
  call   <dd_mlt@plt>
  lea    0x50(%rsp),%rsi
  mov    %r12,%rdx
  mov    %rbx,%rdi
  call   <dd_mul@plt>
  movl   $0x0,0x4(%rsp)
L99:
  mov    %r12,%rsi
  mov    %r15,%rdi
  call   <dd_inc@plt>
  mov    %r14,%rsi
  mov    %rbx,%rdi
  call   <dd_mlt@plt>
  mov    %r13,%rsi
  mov    %r12,%rdx
  mov    %rbx,%rdi
  addl   $0x1,0x4(%rsp)
  add    $0x10,%r13
  call   <dd_mul@plt>
  movsd  0x30(%rsp),%xmm0
  andpd  D432(%rip),%xmm0
  comisd D304(%rip),%xmm0
  jbe    L13
  cmpl   $0x5,0x4(%rsp)
  jne    L99
L13:
  mov    %r12,%rsi
  mov    %r15,%rdi
  mov    $0x9,%r13d
  call   <dd_inc@plt>
  nopl   0x0(%rax)
L17:
  mov    D4(%rip),%rax
  mov    %rbx,%rsi
  mov    %r15,%rdi
  movq   %rax,%xmm0
  call   <dd_mul_pwr2@plt>
  mov    %r12,%rsi
  mov    %r15,%rdi
  call   <dd_sqr@plt>
  mov    %r15,%rdx
  mov    %r12,%rsi
  mov    %rbx,%rdi
  call   <dd_add@plt>
  sub    $0x1,%r13d
  jne    L17
  mov    D264(%rip),%rax
  mov    %rbp,%rsi
  mov    %r15,%rdi
  movq   %rax,%xmm0
  call   <dd_add_dd_d@plt>
  cvttsd2si 0x8(%rsp),%esi
  mov    %rbp,%rdi
  call   <dd_ldexp@plt>
  jmp    L113
L46:
  pxor   %xmm0,%xmm0
  movapd %xmm0,%xmm1
  jmp    L41
  nopl   (%rax)
L93:
  movsd  D197(%rip),%xmm0
  movsd  D62(%rip),%xmm1
  jmp    L41
  nopl   (%rax)
L40:
  movsd  D264(%rip),%xmm0
  pxor   %xmm1,%xmm1
  jmp    L41
L67:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D466
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D222
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D213
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D287
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D400
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D292
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D111
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D418
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D207
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D148
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D403
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D333
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D77
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D110
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D443
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D118
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D261
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D444
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D48
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D147
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D315
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D262
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D413
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D331
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D339
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D374
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D417
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D314
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D53
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D27
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D154
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D104
data_type: f64[2]
name: extract_symbolic_data
arguments: label: D457
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D13
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D432
data_type: qword[2]
name: extract_symbolic_data
arguments: label: D112
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D352
data_type: f64
name: extract_symbolic_data
arguments: label: D264
data_type: f64
name: extract_symbolic_data
arguments: label: D290
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D146
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D304
data_type: byte[16]
name: extract_symbolic_data
arguments: label: D4
data_type: f64
name: extract_symbolic_data
arguments: label: D197
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D62
data_type: byte[8]

tool: {"D466": [0x55,0x55,0x55,0x55,0x55,0x55,0xc5,0x3f]}

tool: {"D222": [0x55,0x55,0x55,0x55,0x55,0x55,0x65,0x3c]}

tool: {"D213": [0x55,0x55,0x55,0x55,0x55,0x55,0xa5,0x3f]}

tool: {"D287": [0x55,0x55,0x55,0x55,0x55,0x55,0x45,0x3c]}

tool: {"D400": [0x11,0x11,0x11,0x11,0x11,0x11,0x81,0x3f]}

tool: {"D292": [0x11,0x11,0x11,0x11,0x11,0x11,0x01,0x3c]}

tool: {"D111": [0x17,0x6c,0xc1,0x16,0x6c,0xc1,0x56,0x3f]}

tool: {"D418": [0x9f,0xf4,0x49,0x9f,0xf4,0x49,0xef,0xbb]}

tool: {"D207": [0x1a,0xa0,0x01,0x1a,0xa0,0x01,0x2a,0x3f]}

tool: {"D148": [0x1a,0xa0,0x01,0x1a,0xa0,0x01,0x6a,0x3b]}

tool: {"D403": [0x1a,0xa0,0x01,0x1a,0xa0,0x01,0xfa,0x3e]}

tool: {"D333": [0x1a,0xa0,0x01,0x1a,0xa0,0x01,0x3a,0x3b]}

tool: {"D77": [0x34,0xc7,0x56,0xa5,0xe3,0x1d,0xc7,0x3e]}

tool: {"D110": [0x00,0x6c,0xdc,0x8d,0x4f,0x15,0x6c,0xbb]}

tool: {"D443": [0x5c,0x9f,0x78,0xb7,0x4f,0x7e,0x92,0x3e]}

tool: {"D118": [0x9a,0xa9,0x4f,0x5b,0xc0,0xbb,0x3c,0x3b]}

tool: {"D261": [0xe4,0x44,0xf5,0x67,0x45,0xe6,0x5a,0x3e]}

tool: {"D444": [0x09,0xf2,0xd1,0x06,0x2e,0x06,0xfc,0xba]}

tool: {"D48": [0x98,0xd8,0xf8,0xef,0xd8,0xee,0x21,0x3e]}

tool: {"D147": [0x06,0x4c,0xe1,0x59,0xc9,0xae,0xc2,0xba]}

tool: {"D315": [0x09,0x6d,0xa8,0x13,0x46,0x12,0xe6,0x3d]}

tool: {"D262": [0xbe,0x8e,0x74,0xcc,0xe0,0x28,0x8f,0x3a]}

tool: {"D413": [0x9d,0x7c,0xc0,0xa8,0x74,0x39,0xa9,0x3d]}

tool: {"D331": [0x1f,0xfd,0x2e,0x8a,0x6f,0x5d,0x30,0x3a]}

tool: {"D339": [0x1f,0xb8,0x33,0xe7,0xf3,0xe7,0x6a,0x3d]}

tool: {"D374": [0xcb,0xe8,0x0e,0x6b,0x65,0xd8,0xe1,0x39]}

tool: {"D417": [0x1f,0xb8,0x33,0xe7,0xf3,0xe7,0x2a,0x3d]}

tool: {"D314": [0xcb,0xe8,0x0e,0x6b,0x65,0xd8,0xa1,0x39]}

tool: {"D53": [0x4a,0xad,0x30,0x70,0xc7,0x52,0xe9,0x3c]}

tool: {"D27": [0.0,-4.1923828125]}

tool: {"D154": [0x6c,0xdc,0x5d,0x46,0x81,0xc9,0x8a,0x39]}

tool: {"D104": [709.0,0.6931471805599453]}

tool: {"D457": [0.0,-1.875]}

tool: {"D13": [0xef,0x39,0xfa,0xfe,0x42,0x2e,0xe6,0x3f]}

tool: {"D432": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D112": [0.0,176.0]}

tool: D352: 0.5

tool: D264: 1

tool: {"D290": [0x3f,0x80,0x39,0x3b,0x9e,0xbc,0x7a,0x3c]}

tool: {"D146": [0.0,0.875]}

tool: {"D304": [0xf9,0xff,0xff,0xff,0xff,0xff,0xdf,0x38,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x40]}

tool: D4: 2

tool: {"D197": [0x69,0x57,0x14,0x8b,0x0a,0xbf,0x05,0x40]}

tool: {"D62": [0x3a,0x01,0xb1,0xe2,0x7e,0xd5,0xa4,0x3c]}

assistant: ```c
void dd_exp(const double *a, double *b) {
  const double k = 512.0;
  const double inv_k = 1.0 / k;
  const double exp1hi = 2.718281828459045091e+00;
  const double exp1lo = 1.445646891729250158e-16;
  const double log2hi = 6.931471805599452862e-01;
  const double log2lo = 2.319046813846299558e-17;
  const double dd_eps = 4.93038065763132e-32;
  double m, tol;
  double r[2];
  double s[2];
  double t[2];
  double p[2];
  int i;
  double *inv_fac;
  double i_fac[30];
  i_fac[0] = 1.66666666666666657e-01;
  i_fac[1] = 9.25185853854297066e-18;
  i_fac[2] = 4.16666666666666644e-02;
  i_fac[3] = 2.31296463463574266e-18;
  i_fac[4] = 8.33333333333333322e-03;
  i_fac[5] = 1.15648231731787138e-19;
  i_fac[6] = 1.38888888888888894e-03;
  i_fac[7] = -5.30054395437357706e-20;
  i_fac[8] = 1.98412698412698413e-04;
  i_fac[9] = 1.72095582934207053e-22;
  i_fac[10] = 2.48015873015873016e-05;
  i_fac[11] = 2.15119478667758816e-23;
  i_fac[12] = 2.75573192239858925e-06;
  i_fac[13] = -1.85839327404647208e-22;
  i_fac[14] = 2.75573192239858883e-07;
  i_fac[15] = 2.37677146222502973e-23;
  i_fac[16] = 2.50521083854417202e-08;
  i_fac[17] = -1.44881407093591197e-24;
  i_fac[18] = 2.08767569878681002e-09;
  i_fac[19] = -1.20734505911325997e-25;
  i_fac[20] = 1.60590438368216133e-10;
  i_fac[21] = 1.25852945887520981e-26;
  i_fac[22] = 1.14707455977297245e-11;
  i_fac[23] = 2.06555127528307454e-28;
  i_fac[24] = 7.64716373181981641e-13;
  i_fac[25] = 7.03872877733453001e-30;
  i_fac[26] = 4.77947733238738525e-14;
  i_fac[27] = 4.39920548583408126e-31;
  i_fac[28] = 2.81145725434552060e-15;
  i_fac[29] = 1.65088427308614326e-31;
  if (a[0] <= -709.0) {
    b[0] = 0.0;
    b[1] = 0.0;
    return;
  }
  if (a[0] >= 709.0) {
    b[0] = -1.0;
    b[1] = 0.0;
    return;
  }
  if (dd_is_zero(a) == 1) {
    b[0] = 1.0;
    b[1] = 0.0;
    return;
  }
  if (dd_is_one(a) == 1) {
    b[0] = exp1hi;
    b[1] = exp1lo;
    return;
  }
  m = floor(a[0] / log2hi + 0.5);
  r[0] = log2hi;
  r[1] = log2lo;
  dd_mlt_d(r, m);
  dd_sub(a, r, r);
  dd_mul_pwr2(r, inv_k, r);
  dd_sqr(r, p);
  dd_mul_pwr2(p, 0.5, s);
  dd_inc(s, r);
  dd_mlt(p, r);
  inv_fac = &i_fac[0];
  dd_mul(p, inv_fac, t);
  i = 0;
  tol = inv_k * dd_eps;
  do {
    dd_inc(s, t);
    dd_mlt(p, r);
    inv_fac = &i_fac[2 * (++i)];
    dd_mul(p, inv_fac, t);
  } while (fabs(t[0]) > tol && i < 5);
  dd_inc(s, t);
  for (i = 0; i < 9; i++) {
    dd_mul_pwr2(s, 2.0, p);
    dd_sqr(s, t);
    dd_add(p, t, s);
  }
  dd_add_dd_d(s, 1.0, b);
  i = (int)m;
  dd_ldexp(b, i);
}
```