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

```asm
<calcSC_ak1>:
  endbr64
  push   %r15
  push   %r14
  mov    %r9,%r14
  push   %r13
  mov    %rdx,%r13
  push   %r12
  push   %rbp
  mov    %r8,%rbp
  mov    %r9,%r8
  push   %rbx
  movslq %edi,%rbx
  mov    %ebx,%edi
  shl    $0x3,%rbx
  sub    $0x38,%rsp
  mov    0x90(%rsp),%r12
  mov    0x98(%rsp),%rdx
  mov    %rsi,0x20(%rsp)
  mov    %rcx,0x28(%rsp)
  mov    %rbp,%rcx
  mov    0x88(%rsp),%r15
  movsd  %xmm0,0x18(%rsp)
  mov    %r12,%rsi
  movapd %xmm2,%xmm0
  movsd  %xmm1,0x10(%rsp)
  movapd %xmm3,%xmm1
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <QuadSD_ak1@plt>
  movsd  0x0(%rbp),%xmm6
  movsd  (%rsp),%xmm2
  movq   D14(%rip),%xmm0
  movsd  -0x8(%r12,%rbx,1),%xmm3
  movsd  D7(%rip),%xmm7
  movsd  0x8(%rsp),%xmm1
  andpd  %xmm0,%xmm3
  andpd  %xmm0,%xmm6
  movsd  0x10(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm5
  mulsd  %xmm7,%xmm3
  comisd %xmm6,%xmm3
  jb     L11
  movsd  -0x10(%r12,%rbx,1),%xmm3
  movsd  (%r14),%xmm6
  mov    $0x3,%eax
  andpd  %xmm0,%xmm3
  andpd  %xmm0,%xmm6
  mulsd  %xmm7,%xmm3
  comisd %xmm6,%xmm3
  jae    L1
L11:
  mulsd  %xmm4,%xmm1
  movsd  %xmm1,(%r15)
  movsd  (%r14),%xmm3
  movsd  0x0(%rbp),%xmm1
  movapd %xmm3,%xmm6
  andpd  %xmm0,%xmm6
  andpd  %xmm1,%xmm0
  comisd %xmm0,%xmm6
  movapd %xmm5,%xmm0
  jae    L9
  divsd  %xmm1,%xmm0
  mov    0x70(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x78(%rsp),%rax
  movsd  (%r14),%xmm0
  divsd  0x0(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x70(%rsp),%rax
  mulsd  (%rax),%xmm2
  mov    0x80(%rsp),%rax
  movsd  %xmm2,(%rax)
  mov    0x70(%rsp),%rax
  movsd  (%r15),%xmm0
  divsd  0x0(%rbp),%xmm0
  addsd  %xmm0,%xmm2
  movsd  (%rax),%xmm0
  mov    0x20(%rsp),%rax
  mulsd  %xmm4,%xmm2
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x0(%r13)
  movsd  (%r14),%xmm0
  divsd  0x0(%rbp),%xmm0
  mulsd  %xmm5,%xmm0
  subsd  %xmm0,%xmm4
  movsd  %xmm4,(%rax)
  mov    0x80(%rsp),%rax
  movsd  (%r15),%xmm1
  movsd  (%rax),%xmm0
  mov    0x78(%rsp),%rax
  mulsd  (%r14),%xmm0
  mulsd  (%rax),%xmm1
  mov    0x28(%rsp),%rax
  addsd  %xmm1,%xmm0
  addsd  %xmm5,%xmm0
  movsd  %xmm0,(%rax)
  mov    $0x1,%eax
L1:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L9:
  divsd  %xmm3,%xmm0
  mov    0x70(%rsp),%rax
  movapd %xmm4,%xmm1
  movsd  %xmm0,(%rax)
  mov    0x78(%rsp),%rax
  movsd  0x0(%rbp),%xmm0
  divsd  (%r14),%xmm0
  movsd  %xmm0,(%rax)
  movapd %xmm2,%xmm0
  mov    0x80(%rsp),%rax
  mulsd  %xmm4,%xmm0
  movsd  %xmm0,(%rax)
  mov    0x70(%rsp),%rax
  divsd  (%r14),%xmm1
  addsd  %xmm5,%xmm0
  mulsd  (%r15),%xmm1
  mulsd  (%rax),%xmm0
  mov    0x78(%rsp),%rax
  addsd  %xmm0,%xmm1
  movsd  %xmm1,0x0(%r13)
  mulsd  (%rax),%xmm4
  mov    0x20(%rsp),%rax
  subsd  %xmm5,%xmm4
  movsd  %xmm4,(%rax)
  mov    0x78(%rsp),%rax
  addsd  (%rax),%xmm2
  mov    0x28(%rsp),%rax
  mulsd  %xmm5,%xmm2
  addsd  (%r15),%xmm2
  movsd  %xmm2,(%rax)
  add    $0x38,%rsp
  mov    $0x2,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D14
type: qword
name: disassemble_data
arguments: data_label: D7
type: byte[8]

tool: D14: 9223372036854775807

tool: {"D7": [0x00,0x00,0x00,0x00,0x00,0x00,0x19,0x3d]}

assistant: ```c
int calcSC_ak1(int N, double a, double b, double *a1, double *a3, double *a7, double *c, double *d, double *e, double *f, double *g, double *h, double K[100 + 1], double u, double v, double qk[100 + 1]) {
  int dumFlag = 3;
  QuadSD_ak1(N, u, v, K, qk, c, d);
  if (fabs((*c)) <= (100.0 * 2.2204460492503131e-016 * fabs(K[N - 1]))) {
    if (fabs((*d)) <= (100.0 * 2.2204460492503131e-016 * fabs(K[N - 2])))
      return dumFlag;
  }
  *h = v * b;
  if (fabs((*d)) >= fabs((*c))) {
    dumFlag = 2;
    *e = a / (*d);
    *f = (*c) / (*d);
    *g = u * b;
    *a3 = (*e) * ((*g) + a) + (*h) * (b / (*d));
    *a1 = -a + (*f) * b;
    *a7 = (*h) + ((*f) + u) * a;
  } else {
    dumFlag = 1;
    *e = a / (*c);
    *f = (*d) / (*c);
    *g = (*e) * u;
    *a3 = (*e) * a + ((*g) + (*h) / (*c)) * b;
    *a1 = -(a * ((*d) / (*c))) + b;
    *a7 = (*g) * (*d) + (*h) * (*f) + a;
  }
  return dumFlag;
}
```