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

```asm
<mvnapp>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xe50,%rsp
  mov    %rdi,-0xe08(%rbp)
  mov    %rsi,-0xe10(%rbp)
  mov    %rdx,-0xe18(%rbp)
  mov    %rcx,-0xe20(%rbp)
  mov    %r8,-0xe28(%rbp)
  mov    %r9,-0xe30(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0xe38(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xe40(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0xe48(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0xe50(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0xe10(%rbp),%rax
  mov    (%rax),%eax
  mov    %eax,-0xde4(%rbp)
  movl   $0x1,-0xdf0(%rbp)
  jmp    L71
L13:
  mov    -0xdf0(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0xe18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xdf0(%rbp),%eax
  cltq
  movsd  %xmm0,-0xdd0(%rbp,%rax,8)
  mov    -0xdf0(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0xe20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xdf0(%rbp),%eax
  cltq
  movsd  %xmm0,-0xd30(%rbp,%rax,8)
  addl   $0x1,-0xdf0(%rbp)
L71:
  mov    -0xdf0(%rbp),%eax
  cmp    -0xde4(%rbp),%eax
  jle    L13
  movl   $0x1,-0xdf0(%rbp)
  movl   $0x0,-0xde8(%rbp)
  jmp    L4
L51:
  movl   $0x1,-0xdec(%rbp)
  jmp    L35
L65:
  mov    -0xde8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0xe28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xdec(%rbp),%eax
  movslq %eax,%rcx
  mov    -0xdf0(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  shl    $0x2,%rax
  add    %rdx,%rax
  shl    $0x2,%rax
  add    %rcx,%rax
  movsd  %xmm0,-0xc90(%rbp,%rax,8)
  addl   $0x1,-0xde8(%rbp)
  addl   $0x1,-0xdec(%rbp)
L35:
  mov    -0xdec(%rbp),%eax
  cmp    -0xde4(%rbp),%eax
  jle    L65
  addl   $0x1,-0xdf0(%rbp)
L4:
  mov    -0xdf0(%rbp),%eax
  cmp    -0xde4(%rbp),%eax
  jle    L51
  mov    -0xe08(%rbp),%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  jne    L68
  mov    -0xe38(%rbp),%rax
  mov    (%rax),%rdi
  mov    -0xe30(%rbp),%rax
  mov    (%rax),%r10d
  lea    -0xde0(%rbp),%r9
  lea    -0xc90(%rbp),%rcx
  lea    -0xd30(%rbp),%rdx
  lea    -0xdd0(%rbp),%rsi
  mov    -0xde4(%rbp),%eax
  lea    -0xdf4(%rbp),%r8
  push   %r8
  lea    -0xdd8(%rbp),%r8
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10d,%r8d
  mov    %eax,%edi
  call   <mvn1@plt>
  add    $0x10,%rsp
  jmp    L61
L68:
  mov    -0xe38(%rbp),%rax
  mov    (%rax),%rdi
  mov    -0xe30(%rbp),%rax
  mov    (%rax),%r10d
  lea    -0xde0(%rbp),%r9
  lea    -0xc90(%rbp),%rcx
  lea    -0xd30(%rbp),%rdx
  lea    -0xdd0(%rbp),%rsi
  mov    -0xde4(%rbp),%eax
  lea    -0xdf4(%rbp),%r8
  push   %r8
  lea    -0xdd8(%rbp),%r8
  push   %r8
  movq   %rdi,%xmm0
  mov    %r10d,%r8d
  mov    %eax,%edi
  call   <mvn2@plt>
  add    $0x10,%rsp
L61:
  movsd  -0xde0(%rbp),%xmm0
  mov    -0xe40(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0xdd8(%rbp),%xmm0
  mov    -0xe48(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0xdf4(%rbp),%edx
  mov    -0xe50(%rbp),%rax
  mov    %edx,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L22
  call   <__stack_chk_fail@plt>
L22:
  leave
  ret
```

assistant: ```c
void mvnapp(int *type, int *m0, double w0[], double x0[], double corr[], int *nsim, double *eps, double *pr, double *sd, int *ifail) {
  double r[20][20];
  double w[20], x[20];
  double p, perr;
  int i, m, j, k, ierr;
  void mvn1(int, double[], double[], double[][20], int, double, double *, double *, int *);
  void mvn2(int, double[], double[], double[][20], int, double, double *, double *, int *);
  m = *m0;
  for (i = 1; i <= m; i++) {
    w[i] = w0[i - 1];
    x[i] = x0[i - 1];
  }
  for (i = 1, k = 0; i <= m; i++) {
    for (j = 1; j <= m; j++) {
      r[i][j] = corr[k];
      k++;
    }
  }
  if (*type == 1)
    mvn1(m, w, x, r, *nsim, *eps, &p, &perr, &ierr);
  else
    mvn2(m, w, x, r, *nsim, *eps, &p, &perr, &ierr);
  *pr = p;
  *sd = perr;
  *ifail = ierr;
}
```