user: Translate the below assembly code into C programming language:

```asm
<rcomp>:
  endbr64
  push   %rbp
  mov    %rdi,%rbp
  sub    $0x10,%rsp
  movsd  (%rdi),%xmm2
  comisd D130(%rip),%xmm2
  movsd  (%rsi),%xmm1
  jb     L47
  divsd  %xmm2,%xmm1
  pxor   %xmm2,%xmm2
  ucomisd %xmm2,%xmm1
  movsd  %xmm1,D115(%rip)
  jp     L22
  movapd %xmm2,%xmm0
  jne    L22
  add    $0x10,%rsp
  pop    %rbp
  ret
  nopl   0x0(%rax)
L22:
  movsd  D100(%rip),%xmm4
  movsd  0x0(%rbp),%xmm3
  lea    D115(%rip),%rdi
  movsd  D109(%rip),%xmm0
  movapd %xmm4,%xmm1
  movsd  %xmm3,(%rsp)
  divsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm1
  mulsd  %xmm1,%xmm0
  subsd  %xmm4,%xmm0
  mulsd  %xmm1,%xmm0
  addsd  D42(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm3,%xmm1
  mulsd  D113(%rip),%xmm1
  subsd  D137(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,D117(%rip)
  call   <rlog@plt>
  movsd  (%rsp),%xmm3
  movsd  0x0(%rbp),%xmm1
  pxor   %xmm2,%xmm2
  movsd  D117(%rip),%xmm4
  mulsd  %xmm0,%xmm3
  ucomisd %xmm1,%xmm2
  subsd  %xmm3,%xmm4
  movsd  %xmm4,D65(%rip)
  ja     L54
  sqrtsd %xmm1,%xmm1
L55:
  movapd %xmm4,%xmm0
  movsd  %xmm1,(%rsp)
  call   <exp@plt>
  movsd  (%rsp),%xmm1
  mulsd  D65(%rip),%xmm1
  add    $0x10,%rsp
  movapd %xmm0,%xmm2
  pop    %rbp
  movapd %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  ret
  nopl   0x0(%rax,%rax,1)
L47:
  movapd %xmm1,%xmm0
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  call   <log@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  (%rsp),%xmm1
  comisd D100(%rip),%xmm2
  mulsd  %xmm2,%xmm0
  subsd  %xmm1,%xmm0
  jae    L53
  movsd  %xmm2,(%rsp)
  call   <exp@plt>
  movsd  (%rsp),%xmm2
  mov    %rbp,%rdi
  mulsd  %xmm0,%xmm2
  movsd  %xmm2,(%rsp)
  call   <gam1@plt>
  addsd  D100(%rip),%xmm0
  mulsd  (%rsp),%xmm0
  add    $0x10,%rsp
  pop    %rbp
  ret
  nopl   0x0(%rax)
L53:
  call   <exp@plt>
  mov    %rbp,%rdi
  movsd  %xmm0,(%rsp)
  call   <Xgamm@plt>
  movsd  (%rsp),%xmm1
  add    $0x10,%rsp
  pop    %rbp
  divsd  %xmm0,%xmm1
  movapd %xmm1,%xmm0
  ret
L54:
  movapd %xmm1,%xmm0
  movsd  %xmm4,(%rsp)
  call   <sqrt@plt>
  movsd  (%rsp),%xmm4
  movapd %xmm0,%xmm1
  jmp    L55
```

tool_calls: name: disassemble_data
arguments: tag: D130
type_infered: f64[2]
name: disassemble_data
arguments: tag: D100
type_infered: f64[2]
name: disassemble_data
arguments: tag: D109
type_infered: f64
name: disassemble_data
arguments: tag: D42
type_infered: f64
name: disassemble_data
arguments: tag: D113
type_infered: f64
name: disassemble_data
arguments: tag: D137
type_infered: f64
name: disassemble_data
arguments: tag: D65
type_infered: byte[8]

tool: {"D130": [20.0,1.0]}

tool: {"D100": [1.0,0.75]}

tool: D109: 0.75

tool: D42: 3.5

tool: D113: 1260

tool: D137: 105

tool: {"D65": [0x56,0x36,0xd4,0x33,0x45,0x88,0xd9,0x3f]}

assistant: ```c
double rcomp(double *a, double *x) {
  static double rt2pin = .398942280401433e0;
  static double rcomp, t, t1, u;
  rcomp = 0.0e0;
  if (*a >= 20.0e0)
    goto S20;
  t = *a * log(*x) - *x;
  if (*a >= 1.0e0)
    goto S10;
  rcomp = *a * exp(t) * (1.0e0 + gam1(a));
  return rcomp;
S10:
  rcomp = exp(t) / Xgamm(a);
  return rcomp;
S20:
  u = *x / *a;
  if (u == 0.0e0)
    return rcomp;
  t = pow(1.0e0 / *a, 2.0);
  t1 = (((0.75e0 * t - 1.0e0) * t + 3.5e0) * t - 105.0e0) / (*a * 1260.0e0);
  t1 -= (*a * rlog(&u));
  rcomp = rt2pin * sqrt(*a) * exp(t1);
  return rcomp;
}
```