user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<multiplets>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x7a0,%rsp
  mov    %rdi,-0x708(%rbp)
  mov    %rsi,-0x710(%rbp)
  mov    %rdx,-0x718(%rbp)
  mov    %rcx,-0x720(%rbp)
  mov    %r8,-0x728(%rbp)
  mov    %r9,-0x730(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x738(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0x740(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0x748(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0x750(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0x758(%rbp)
  mov    0x38(%rbp),%rax
  mov    %rax,-0x760(%rbp)
  mov    0x40(%rbp),%rax
  mov    %rax,-0x768(%rbp)
  mov    0x48(%rbp),%rax
  mov    %rax,-0x770(%rbp)
  mov    0x50(%rbp),%rax
  mov    %rax,-0x778(%rbp)
  mov    0x58(%rbp),%rax
  mov    %rax,-0x780(%rbp)
  mov    0x60(%rbp),%rax
  mov    %rax,-0x788(%rbp)
  mov    0x68(%rbp),%rax
  mov    %rax,-0x790(%rbp)
  mov    0x70(%rbp),%rax
  mov    %rax,-0x798(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D40(%rip),%xmm0
  movsd  %xmm0,-0x690(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x688(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x680(%rbp)
  movsd  D53(%rip),%xmm0
  movsd  %xmm0,-0x678(%rbp)
  movq   $0x0,-0x6d8(%rbp)
  mov    -0x710(%rbp),%rdx
  mov    -0x708(%rbp),%rax
  mov    %rdx,%rcx
  mov    $0x0,%edx
  mov    %rax,%rsi
  mov    $0x1,%edi
  call   <lmatrix@plt>
  mov    %rax,-0x6d0(%rbp)
  mov    -0x710(%rbp),%rax
  mov    $0x3,%ecx
  mov    $0x1,%edx
  mov    %eax,%esi
  mov    $0x1,%edi
  call   <dmatrix@plt>
  mov    %rax,-0x6c8(%rbp)
  mov    -0x710(%rbp),%rax
  mov    $0x3,%ecx
  mov    $0x1,%edx
  mov    %eax,%esi
  mov    $0x1,%edi
  call   <dmatrix@plt>
  mov    %rax,-0x6c0(%rbp)
  mov    $0x3,%ecx
  mov    $0x1,%edx
  mov    $0x3,%esi
  mov    $0x1,%edi
  call   <dmatrix@plt>
  mov    %rax,-0x6b8(%rbp)
  mov    $0x3,%ecx
  mov    $0x1,%edx
  mov    $0x64,%esi
  mov    $0x1,%edi
  call   <dmatrix@plt>
  mov    %rax,-0x6b0(%rbp)
  movq   $0x1,-0x700(%rbp)
  jmp    L279
L208:
  mov    -0x700(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x720(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x328,%rax
  mov    (%rax),%rax
  cmp    $0x1,%rax
  jle    L204
  mov    -0x700(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x720(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x328,%rax
  mov    (%rax),%rax
  add    $0x1,%rax
  mov    %rax,-0x6a8(%rbp)
  mov    -0x728(%rbp),%rax
  lea    0x8(%rax),%rdx
  mov    -0x700(%rbp),%rax
  mov    %rax,(%rdx)
  movq   $0x1,-0x6f0(%rbp)
  jmp    L240
L95:
  mov    -0x700(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x720(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6f0(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    -0x6f0(%rbp),%rdx
  add    $0x1,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0x728(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rax),%rax
  mov    %rax,(%rdx)
  addq   $0x1,-0x6f0(%rbp)
L240:
  mov    -0x700(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x720(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x328,%rax
  mov    (%rax),%rax
  cmp    %rax,-0x6f0(%rbp)
  jle    L95
  mov    -0x730(%rbp),%rdx
  mov    -0x728(%rbp),%rcx
  mov    -0x6a8(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <lsort@plt>
  movq   $0x0,-0x6f8(%rbp)
  movq   $0x1,-0x6f0(%rbp)
  jmp    L56
L149:
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6d0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    (%rax),%rax
  cmp    %rax,-0x6a8(%rbp)
  jne    L265
  movq   $0x1,-0x6e8(%rbp)
  jmp    L46
L254:
  mov    -0x6e8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x728(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x6d0(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x6e8(%rbp),%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  mov    (%rax),%rax
  cmp    %rax,%rdx
  jne    L22
  addq   $0x1,-0x6e8(%rbp)
L46:
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6d0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    (%rax),%rax
  cmp    %rax,-0x6e8(%rbp)
  jle    L254
  jmp    L164
L22:
  nop
L164:
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6d0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    (%rax),%rax
  cmp    %rax,-0x6e8(%rbp)
  jle    L265
  movq   $0x1,-0x6f8(%rbp)
L265:
  addq   $0x1,-0x6f0(%rbp)
L56:
  mov    -0x6f0(%rbp),%rax
  cmp    -0x6d8(%rbp),%rax
  jg     L276
  cmpq   $0x0,-0x6f8(%rbp)
  je     L149
L276:
  cmpq   $0x0,-0x6f8(%rbp)
  jne    L204
  addq   $0x1,-0x6d8(%rbp)
  mov    -0x6d8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6d0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    -0x6a8(%rbp),%rdx
  mov    %rdx,(%rax)
  movb   $0x0,-0x610(%rbp)
  movb   $0x0,-0x410(%rbp)
  mov    -0x6d8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x790(%rbp),%rax
  add    %rax,%rdx
  mov    -0x6a8(%rbp),%rax
  mov    %rax,(%rdx)
  movq   $0x1,-0x6f0(%rbp)
  jmp    L38
L301:
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x728(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,-0x698(%rbp)
  mov    -0x6d8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x798(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6f0(%rbp),%rax
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x698(%rbp),%rax
  mov    %rax,(%rdx)
  mov    -0x6d8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6d0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6f0(%rbp),%rax
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x698(%rbp),%rax
  mov    %rax,(%rdx)
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x778(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    0x8(%rax),%rax
  mov    %rax,-0x6e8(%rbp)
  mov    -0x6f0(%rbp),%rax
  cmp    -0x6a8(%rbp),%rax
  jne    L115
  lea    D51(%rip),%rax
  jmp    L215
L115:
  lea    D31(%rip),%rax
L215:
  mov    -0x698(%rbp),%rcx
  mov    -0x780(%rbp),%rdx
  add    %rcx,%rdx
  movzbl (%rdx),%edx
  movsbl %dl,%esi
  mov    -0x6e8(%rbp),%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0x750(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x6e8(%rbp),%rdi
  mov    -0x748(%rbp),%rdx
  add    %rdi,%rdx
  movzbl (%rdx),%edx
  movsbq %dl,%rdx
  lea    -0x210(%rbp),%rdi
  mov    %rax,%r9
  mov    %esi,%r8d
  lea    D30(%rip),%rax
  mov    %rax,%rsi
  mov    $0x0,%eax
  call   <sprintf@plt>
  lea    -0x210(%rbp),%rdx
  lea    -0x610(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <strcat@plt>
  mov    -0x698(%rbp),%rdx
  lea    -0x210(%rbp),%rax
  lea    D26(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <sprintf@plt>
  lea    -0x210(%rbp),%rdx
  lea    -0x410(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <strcat@plt>
  movq   $0x1,-0x6e8(%rbp)
  jmp    L37
L35:
  cmpq   $0x1,-0x6f0(%rbp)
  jle    L69
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x768(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  lea    0x30(%rax),%rdx
  mov    -0x6c8(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <dot@plt>
  pxor   %xmm1,%xmm1
  comisd %xmm0,%xmm1
  jbe    L69
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x768(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6e8(%rbp),%rax
  add    $0x6,%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6c8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6e8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movq   D0(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L272
L69:
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x768(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6e8(%rbp),%rax
  add    $0x6,%rax
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x6c8(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x6e8(%rbp),%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
L272:
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x770(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6e8(%rbp),%rax
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x6c0(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x6e8(%rbp),%rax
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x6e8(%rbp)
L37:
  cmpq   $0x3,-0x6e8(%rbp)
  jle    L35
  addq   $0x1,-0x6f0(%rbp)
L38:
  mov    -0x6f0(%rbp),%rax
  cmp    -0x6a8(%rbp),%rax
  jle    L301
  lea    -0x630(%rbp),%rdx
  mov    -0x6a8(%rbp),%rsi
  mov    -0x6c0(%rbp),%rax
  mov    %rdx,%rcx
  mov    $0x3,%edx
  mov    %rax,%rdi
  call   <ave_dmatrix@plt>
  lea    -0x650(%rbp),%rdx
  mov    -0x6a8(%rbp),%rsi
  mov    -0x6c8(%rbp),%rax
  mov    %rdx,%rcx
  mov    $0x3,%edx
  mov    %rax,%rdi
  call   <ave_dmatrix@plt>
  lea    -0x670(%rbp),%rdx
  lea    -0x690(%rbp),%rcx
  lea    -0x650(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <cross@plt>
  lea    -0x690(%rbp),%rdx
  lea    -0x650(%rbp),%rax
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <magang@plt>
  mov    %eax,%ecx
  mov    -0x6b8(%rbp),%rdx
  lea    -0x670(%rbp),%rax
  mov    %ecx,%esi
  mov    %rax,%rdi
  call   <arb_rotation@plt>
  movq   $0x0,-0x6a0(%rbp)
  movq   $0x1,-0x6f0(%rbp)
  jmp    L90
L217:
  mov    -0x6f0(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x728(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,-0x698(%rbp)
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x778(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    0x8(%rax),%rax
  mov    %rax,-0x6e8(%rbp)
  jmp    L286
L142:
  movq   $0x1,-0x6e0(%rbp)
  jmp    L193
L178:
  mov    -0x6e8(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x760(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x6e0(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x6e0(%rbp),%rax
  movsd  -0x630(%rbp,%rax,8),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x6e0(%rbp),%rax
  movsd  %xmm0,-0x650(%rbp,%rax,8)
  addq   $0x1,-0x6e0(%rbp)
L193:
  cmpq   $0x3,-0x6e0(%rbp)
  jle    L178
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x778(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x8,%rax
  mov    (%rax),%rdx
  mov    -0x6e8(%rbp),%rax
  sub    %rdx,%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x6b0(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rcx
  mov    -0x6b8(%rbp),%rdx
  lea    -0x650(%rbp),%rax
  mov    %rcx,%r9
  mov    $0x3,%r8d
  mov    $0x3,%ecx
  mov    $0x3,%esi
  mov    %rax,%rdi
  call   <multi_vec_Tmatrix@plt>
  addq   $0x1,-0x6e8(%rbp)
L286:
  mov    -0x698(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x778(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x10,%rax
  mov    (%rax),%rax
  cmp    %rax,-0x6e8(%rbp)
  jle    L142
  addq   $0x1,-0x6f0(%rbp)
L90:
  mov    -0x6f0(%rbp),%rax
  cmp    -0x6a8(%rbp),%rax
  jle    L217
L204:
  addq   $0x1,-0x700(%rbp)
L279:
  mov    -0x700(%rbp),%rax
  cmp    -0x718(%rbp),%rax
  jle    L208
  mov    -0x788(%rbp),%rax
  mov    -0x6d8(%rbp),%rdx
  mov    %rdx,(%rax)
  mov    -0x710(%rbp),%rcx
  mov    -0x708(%rbp),%rdx
  mov    -0x6d0(%rbp),%rax
  mov    %rcx,%r8
  mov    $0x0,%ecx
  mov    $0x1,%esi
  mov    %rax,%rdi
  call   <free_lmatrix@plt>
  mov    -0x710(%rbp),%rax
  mov    %eax,%edx
  mov    -0x6c8(%rbp),%rax
  mov    $0x3,%r8d
  mov    $0x1,%ecx
  mov    $0x1,%esi
  mov    %rax,%rdi
  call   <free_dmatrix@plt>
  mov    -0x710(%rbp),%rax
  mov    %eax,%edx
  mov    -0x6c0(%rbp),%rax
  mov    $0x3,%r8d
  mov    $0x1,%ecx
  mov    $0x1,%esi
  mov    %rax,%rdi
  call   <free_dmatrix@plt>
  mov    -0x6b8(%rbp),%rax
  mov    $0x3,%r8d
  mov    $0x1,%ecx
  mov    $0x3,%edx
  mov    $0x1,%esi
  mov    %rax,%rdi
  call   <free_dmatrix@plt>
  mov    -0x6b0(%rbp),%rax
  mov    $0x3,%r8d
  mov    $0x1,%ecx
  mov    $0x64,%edx
  mov    $0x1,%esi
  mov    %rax,%rdi
  call   <free_dmatrix@plt>
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L36
  call   <__stack_chk_fail@plt>
L36:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: label: D40
type_infered: byte[8]
name: retrieve_data
arguments: label: D53
type_infered: f64
name: retrieve_data
arguments: label: D51
type_infered: string
name: retrieve_data
arguments: label: D31
type_infered: string
name: retrieve_data
arguments: label: D30
type_infered: string
name: retrieve_data
arguments: label: D26
type_infered: string
name: retrieve_data
arguments: label: D0
type_infered: f32[2]

tool: {"D40": [0x85,0xeb,0x51,0xb8,0xfe,0x87,0xc3,0xc0]}

tool: D53: 1

tool: 

tool: D31:   +  

tool: D30: %c: %d %c%s

tool: D26: %d_

tool: {"D0": [0.0,-0.0]}

assistant: ```c
void multiplets(long num_ple, long max_ple, long num_residue, long **pair_info, long *ivec, long *idx1, char **AtomName, char **ResName, char *ChainID, long *ResSeq, char **Miscs, double **xyz, double **orien, double **org, long **seidx, char *bseq, long *num_multi, long *multi_idx, long **multi) {
  char pairstr[512], pairnum[512], tmp[512];
  double z[4] = {-9999.99, 0.0, 0.0, 1.0};
  double hinge[4], zave[4], pave[4], **ztot, **ptot, **rotmat, **xyz_residue;
  long i, inum_base, inum, is_exist, j, jr, k, m, n_unique = 0, **mtxple;
  mtxple = lmatrix(1, num_ple, 0, max_ple);
  ztot = dmatrix(1, max_ple, 1, 3);
  ptot = dmatrix(1, max_ple, 1, 3);
  rotmat = dmatrix(1, 3, 1, 3);
  xyz_residue = dmatrix(1, 100, 1, 3);
  for (i = 1; i <= num_residue; i++) {
    if (pair_info[i][101L] > 1) {
      inum_base = pair_info[i][101L] + 1;
      ivec[1] = i;
      for (j = 1; j <= pair_info[i][101L]; j++)
        ivec[j + 1] = pair_info[i][j];
      lsort(inum_base, ivec, idx1);
      is_exist = 0;
      for (j = 1; j <= n_unique && !is_exist; j++)
        if (inum_base == mtxple[j][0]) {
          for (k = 1; k <= mtxple[j][0]; k++)
            if (ivec[k] != mtxple[j][k])
              break;
          if (k > mtxple[j][0])
            is_exist = 1;
        }
      if (!is_exist) {
        mtxple[++n_unique][0] = inum_base;
        pairstr[0] = '\0';
        pairnum[0] = '\0';
        multi_idx[n_unique] = inum_base;
        for (j = 1; j <= inum_base; j++) {
          jr = ivec[j];
          multi[n_unique][j] = jr;
          mtxple[n_unique][j] = jr;
          k = seidx[jr][1];
          sprintf(tmp, "%c: %d %c%s", ChainID[k], ResSeq[k], bseq[jr], (j == inum_base) ? "" : "  +  ");
          strcat(pairstr, tmp);
          sprintf(tmp, "%d_", jr);
          strcat(pairnum, tmp);
          for (k = 1; k <= 3; k++) {
            if (j > 1 && dot(ztot[1], &orien[jr][6]) < 0.0)
              ztot[j][k] = -orien[jr][6 + k];
            else
              ztot[j][k] = orien[jr][6 + k];
            ptot[j][k] = org[jr][k];
          }
        }
        ave_dmatrix(ptot, inum_base, 3, pave);
        ave_dmatrix(ztot, inum_base, 3, zave);
        cross(zave, z, hinge);
        arb_rotation(hinge, magang(zave, z), rotmat);
        inum = 0;
        for (j = 1; j <= inum_base; j++) {
          jr = ivec[j];
          for (k = seidx[jr][1]; k <= seidx[jr][2]; k++) {
            for (m = 1; m <= 3; m++)
              zave[m] = xyz[k][m] - pave[m];
            multi_vec_Tmatrix(zave, 3, rotmat, 3, 3, xyz_residue[k - seidx[jr][1] + 1]);
          }
        }
      }
    }
  }
  *num_multi = n_unique;
  free_lmatrix(mtxple, 1, num_ple, 0, max_ple);
  free_dmatrix(ztot, 1, max_ple, 1, 3);
  free_dmatrix(ptot, 1, max_ple, 1, 3);
  free_dmatrix(rotmat, 1, 3, 1, 3);
  free_dmatrix(xyz_residue, 1, 100, 1, 3);
}
```