user: 下面的汇编指令对应的C源代码是什么:

```asm
<getT>:
  endbr64
  push   %r15
  mov    %edi,%eax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xd8,%rsp
  mov    %rsi,0x90(%rsp)
  mov    0x118(%rsp),%r15
  mov    %rdx,0xb0(%rsp)
  mov    %rcx,0x88(%rsp)
  mov    %r8,0x78(%rsp)
  mov    %r9,0xb8(%rsp)
  test   %edi,%edi
  jle    L28
  mov    %rdx,%rdi
  movslq %eax,%rdx
  mov    %eax,%eax
  mov    %r9,0x60(%rsp)
  lea    0x0(,%rdx,8),%rsi
  mov    %rax,0xc8(%rsp)
  shl    $0x3,%rax
  xor    %r14d,%r14d
  mov    %rsi,0xa8(%rsp)
  mov    0x148(%rsp),%rsi
  mov    %rdi,0x70(%rsp)
  mov    %rsi,0xa0(%rsp)
  mov    0x110(%rsp),%rsi
  mov    %rax,0xc0(%rsp)
  mov    %rsi,0x68(%rsp)
L62:
  mov    0x90(%rsp),%rax
  mov    0xa0(%rsp),%r12
  mov    %r15,0x30(%rsp)
  mov    0xb0(%rsp),%r13
  mov    %rax,0x20(%rsp)
  mov    (%rax,%r14,4),%ecx
  mov    0x88(%rsp),%rax
  mov    %rax,0x28(%rsp)
  mov    0x78(%rsp),%rax
  mov    %rax,0x48(%rsp)
  mov    0x128(%rsp),%rax
  mov    %rax,0x40(%rsp)
  mov    0x120(%rsp),%rax
  mov    %rax,0x38(%rsp)
  mov    0xc0(%rsp),%rax
  add    %r12,%rax
  mov    %rax,0x98(%rsp)
  mov    0xb8(%rsp),%rax
  mov    %rax,0x50(%rsp)
  mov    0x110(%rsp),%rax
  mov    %rax,0x58(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L52:
  movq   $0x0,(%r12)
  test   %ecx,%ecx
  jle    L1
  movl   $0x0,0x84(%rsp)
  mov    0x20(%rsp),%rax
  mov    (%rax),%edx
  nopw   0x0(%rax,%rax,1)
L14:
  mov    0x88(%rsp),%rax
  mov    0x84(%rsp),%ebp
  add    (%rax,%r14,4),%ebp
  test   %edx,%edx
  jle    L1
  mov    0x130(%rsp),%rax
  movslq %ebp,%rbp
  mov    %r13,%rsi
  xor    %ebx,%ebx
  shl    $0x3,%rbp
  add    %rbp,%rax
  mov    %rax,0x10(%rsp)
  mov    0x138(%rsp),%rax
  add    %rbp,%rax
  add    0x140(%rsp),%rbp
  mov    %rax,0x18(%rsp)
  mov    %r15,%rax
  mov    %rbp,%r13
  mov    %r12,%r15
  mov    %rsi,%rbp
  mov    %rax,%r12
  cs nopw 0x0(%rax,%rax,1)
L8:
  mov    0x28(%rsp),%rax
  movsd  0x0(%rbp),%xmm0
  mov    (%r12,%r14,4),%ecx
  movsd  0x0(%r13),%xmm4
  mov    (%rax),%edx
  mov    0x70(%rsp),%rax
  mulsd  (%rax),%xmm0
  add    %ebx,%edx
  mov    0x130(%rsp),%rax
  add    $0x1,%ebx
  movslq %edx,%rdx
  mulsd  (%rax,%rdx,8),%xmm0
  mov    0x10(%rsp),%rax
  mulsd  (%rax),%xmm0
  mov    0x138(%rsp),%rax
  mulsd  (%rax,%rdx,8),%xmm0
  mov    0x18(%rsp),%rax
  mulsd  (%rax),%xmm0
  mov    0x140(%rsp),%rax
  movsd  %xmm0,0x8(%rsp)
  movsd  (%rax,%rdx,8),%xmm0
  mov    0x68(%rsp),%rax
  movsd  (%rax),%xmm7
  mov    0x60(%rsp),%rax
  movsd  (%rax),%xmm6
  mov    0x78(%rsp),%rax
  movsd  (%rax,%r14,8),%xmm5
  mov    0x128(%rsp),%rax
  mov    (%rax,%r14,4),%r9d
  mov    0x120(%rsp),%rax
  mov    (%rax,%r14,4),%r8d
  mov    0x58(%rsp),%rax
  movsd  (%rax),%xmm3
  mov    0x50(%rsp),%rax
  movsd  (%rax),%xmm2
  mov    0x48(%rsp),%rax
  movsd  (%rax),%xmm1
  mov    0x40(%rsp),%rax
  mov    (%rax),%edx
  mov    0x38(%rsp),%rax
  mov    (%rax),%esi
  mov    0x30(%rsp),%rax
  mov    (%rax),%edi
  call   <kinetic@plt>
  mov    0x20(%rsp),%rax
  mulsd  0x8(%rsp),%xmm0
  addsd  (%r15),%xmm0
  mov    (%rax),%edx
  movsd  %xmm0,(%r15)
  cmp    %ebx,%edx
  jg     L8
  mov    %r12,%rax
  mov    %r15,%r12
  addl   $0x1,0x84(%rsp)
  mov    %rbp,%r13
  mov    %rax,%r15
  mov    0x90(%rsp),%rax
  mov    (%rax,%r14,4),%ecx
  mov    0x84(%rsp),%eax
  cmp    %ecx,%eax
  jl     L14
L1:
  addq   $0x4,0x20(%rsp)
  add    $0x8,%r12
  add    $0x8,%r13
  addq   $0x4,0x28(%rsp)
  addq   $0x8,0x58(%rsp)
  addq   $0x8,0x50(%rsp)
  addq   $0x8,0x48(%rsp)
  addq   $0x4,0x40(%rsp)
  addq   $0x4,0x38(%rsp)
  addq   $0x4,0x30(%rsp)
  cmp    %r12,0x98(%rsp)
  jne    L52
  mov    0xa8(%rsp),%rdi
  addq   $0x8,0x70(%rsp)
  add    $0x1,%r14
  add    %rdi,0xa0(%rsp)
  addq   $0x8,0x68(%rsp)
  addq   $0x8,0x60(%rsp)
  cmp    %r14,0xc8(%rsp)
  jne    L62
L28:
  add    $0xd8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void getT(int nbf, int *nprim, double *normc, int *istart, double *xcenter, double *ycenter, double *zcenter, int *lpower, int *mpower, int *npower, double *normp, double *coef, double *alpha, double *T) {
  int i, j, iprim, jprim, jindex, iindex;
  for (j = 0; j < nbf; j++) {
    for (i = 0; i < nbf; i++) {
      T[i + j * nbf] = 0.;
      for (jprim = 0; jprim < nprim[j]; jprim++) {
        jindex = istart[j] + jprim;
        for (iprim = 0; iprim < nprim[i]; iprim++) {
          iindex = istart[i] + iprim;
          T[i + j * nbf] += normc[i] * normc[j] * normp[iindex] * normp[jindex] * coef[iindex] * coef[jindex] * kinetic(alpha[iindex], lpower[i], mpower[i], npower[i], xcenter[i], ycenter[i], zcenter[i], alpha[jindex], lpower[j], mpower[j], npower[j], xcenter[j], ycenter[j], zcenter[j]);
        }
      }
    }
  }
}
```