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

```asm
<form_LP1>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %r8,%r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x50(%rsp),%rax
  mov    0x78(%rsp),%rdx
  mov    0x68(%rsp),%rbp
  mov    %r9,-0x10(%rsp)
  movslq 0x38(%rsp),%r9
  mov    %edi,-0x14(%rsp)
  mov    %rcx,%rdi
  mov    (%rax,%r9,4),%eax
  lea    0x0(,%r9,4),%r15
  lea    0x0(,%r9,8),%rbx
  mov    %eax,-0x18(%rsp)
  mov    0x48(%rsp),%rax
  mov    -0x4(%rax,%r15,1),%r11d
  mov    0x40(%rsp),%rax
  movslq (%rax,%r9,4),%r14
  movl   $0x0,(%rdx,%r9,4)
  movslq (%rcx),%rdx
  mov    0x70(%rsp),%rcx
  mov    %r14,%rax
  mov    (%rcx,%r9,8),%rsi
  cmp    %edx,%r14d
  jle    L49
  mov    %r9,-0x8(%rsp)
  mov    (%r8,%r14,8),%r10
  lea    0x0(%rbp,%rbx,1),%r12
  mov    $0xffffffff,%ecx
  nop
L67:
  mov    (%r10,%rdx,4),%r9d
  mov    $0xffffffff,%r8d
  test   %r9d,%r9d
  je     L11
  mov    (%r12),%r8
  add    $0x1,%ecx
  movslq %ecx,%r9
  mov    %edx,(%r8,%r9,4)
  mov    %ecx,%r8d
L11:
  mov    %r8d,(%rsi,%rdx,4)
  add    $0x1,%rdx
  cmp    %edx,%eax
  jg     L67
  mov    -0x8(%rsp),%r9
  lea    0x1(%rcx),%r8d
L29:
  movl   $0xffffffff,(%rsi,%r14,4)
  add    $0x1,%eax
  cmp    0x4(%rdi),%eax
  jge    L36
  mov    0x0(%r13,%r14,8),%rdx
  cltq
  nopl   0x0(%rax)
L37:
  mov    (%rdx,%rax,4),%r12d
  mov    $0xffffffff,%r10d
  test   %r12d,%r12d
  je     L44
  mov    0x0(%rbp,%rbx,1),%r10
  add    $0x1,%ecx
  movslq %ecx,%r12
  mov    %eax,(%r10,%r12,4)
  mov    %ecx,%r10d
L44:
  mov    %r10d,(%rsi,%rax,4)
  add    $0x1,%rax
  cmp    %eax,0x4(%rdi)
  jg     L37
L36:
  mov    -0x18(%rsp),%eax
  test   %eax,%eax
  jne    L109
  mov    0x58(%rsp),%rax
  mov    %ecx,%esi
  sub    %r8d,%esi
  mov    (%rax,%r9,4),%edx
  mov    $0x1,%eax
  sub    $0x1,%edx
  cmp    %edx,%esi
  jl     L132
L109:
  add    0x80(%rsp),%r15
  mov    %ecx,(%r15)
  mov    %r15,%r12
  cmp    $0xffffffff,%ecx
  je     L18
  mov    -0x10(%rsp),%rax
  movslq -0x14(%rsp),%r13
  movslq %ecx,%rcx
  xor    %r10d,%r10d
  mov    0x0(%rbp,%r9,8),%rbp
  movslq %r11d,%r8
  mov    -0x8(%rax,%rbx,1),%rbx
  shl    $0x3,%r13
  mov    (%rax,%r9,8),%r15
  mov    (%rbx,%r14,8),%rdi
  lea    (%rdi,%r13,1),%r14
  nopl   0x0(%rax,%rax,1)
L62:
  movslq 0x0(%rbp,%r10,4),%rax
  mov    (%r15,%r10,8),%rsi
  mov    (%rbx,%rax,8),%rdx
  xor    %eax,%eax
  test   %r11d,%r11d
  jle    L122
  nopl   0x0(%rax)
L119:
  movsd  (%rdx,%rax,8),%xmm0
  subsd  (%rdi,%rax,8),%xmm0
  movsd  %xmm0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %r8,%rax
  jne    L119
L122:
  movsd  (%rdx,%r13,1),%xmm0
  subsd  (%r14),%xmm0
  lea    0x1(%r10),%rax
  movsd  %xmm0,(%rsi,%r13,1)
  cmp    %r10,%rcx
  je     L47
  mov    %rax,%r10
  jmp    L62
  nop
L47:
  mov    0x60(%rsp),%rax
  mov    (%rax,%r9,8),%rdx
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L108:
  movl   $0x1,(%rdx,%rax,4)
  add    $0x1,%rax
  cmp    %eax,(%r12)
  jge    L108
L18:
  xor    %eax,%eax
L132:
  mov    0x88(%rsp),%rdx
  mov    %eax,(%rdx)
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L49:
  xor    %r8d,%r8d
  mov    $0xffffffff,%ecx
  jmp    L29
```

assistant: ```c
void form_LP1(int nVar, int nSpt, int *SptType, int *SptIdx, int **RelTab, double ***A, int Lvl, int *Cell, int *Lvl2LPdim, int *FixLstPt, int *MinNumPt, int **PtIn, int **ToOrig, int **Pre2Cur, int *FrstPt, int *LstPt, int *info) {
  int i, j, j1, FixPt, nPtIn, Strt1Pt;
  const int Lvl1 = Lvl - 1;
  const int ynFixLstPt = FixLstPt[Lvl];
  const int PreLPdim = Lvl2LPdim[Lvl - 1];
  FixPt = Cell[Lvl];
  nPtIn = -1;
  FrstPt[Lvl] = 0;
  for (i = SptIdx[0]; i < FixPt; i++)
    if (RelTab[FixPt][i]) {
      ToOrig[Lvl][++nPtIn] = i;
      Pre2Cur[Lvl][i] = nPtIn;
    } else
      Pre2Cur[Lvl][i] = -1;
  Pre2Cur[Lvl][FixPt] = -1;
  Strt1Pt = nPtIn + 1;
  for (i = FixPt + 1; i < SptIdx[1]; i++)
    if (RelTab[FixPt][i]) {
      ToOrig[Lvl][++nPtIn] = i;
      Pre2Cur[Lvl][i] = nPtIn;
    } else
      Pre2Cur[Lvl][i] = -1;
  if (!ynFixLstPt && nPtIn - Strt1Pt < MinNumPt[Lvl] - 1) {
    *info = 1;
    return;
  }
  LstPt[Lvl] = nPtIn;
  for (i = 0; i <= LstPt[Lvl]; i++) {
    j = ToOrig[Lvl][i];
    for (j1 = 0; j1 < PreLPdim; j1++)
      A[Lvl][i][j1] = A[Lvl1][j][j1] - A[Lvl1][FixPt][j1];
    A[Lvl][i][nVar] = A[Lvl1][j][nVar] - A[Lvl1][FixPt][nVar];
  }
  for (i = 0; i <= LstPt[Lvl]; i++)
    PtIn[Lvl][i] = 1;
  *info = 0;
}
```