user: Translate the below assembly code into C programming language:

```asm
<gauss_make_triangular>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rdi,0x28(%rsp)
  mov    %rsi,0x48(%rsp)
  mov    %edx,%r15d
  movslq %edx,%rax
  lea    -0x1(%rdx),%edx
  mov    %edx,-0x4(%rsi,%rax,4)
  test   %r15d,%r15d
  jle    L59
  mov    %rdi,%rbp
  lea    0x8(,%rax,8),%rbx
  mov    %rbx,0x50(%rsp)
  shl    $0x3,%rax
  mov    %rax,0x8(%rsp)
  mov    %rdi,0x30(%rsp)
  mov    %r15d,0x38(%rsp)
  movq   $0x1,0x20(%rsp)
  movq   $0x0,0x18(%rsp)
  lea    (%r15,%r15,1),%eax
  mov    %eax,0x5c(%rsp)
  jmp    L35
L50:
  mov    0x18(%rsp),%r13d
  jmp    L89
L6:
  cmp    0x3c(%rsp),%r15d
  jle    L15
L144:
  mov    0x8(%rsp),%r8
  lea    (%r8,%rbp,1),%rdi
  mov    %r8,%rcx
  add    0x30(%rsp),%rcx
  mov    %rdi,%rdx
  mov    0x20(%rsp),%r9
L117:
  mov    %r9,%rax
L67:
  movsd  (%rcx),%xmm0
  mulsd  0x0(%rbp,%rax,8),%xmm0
  addsd  (%rdx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r15d
  jg     L67
  add    $0x1,%r12d
  add    %r8,%rdx
  add    %r8,%rcx
  cmp    %r12d,%r15d
  jg     L117
  addq   $0x1,0x18(%rsp)
  addq   $0x1,0x20(%rsp)
  add    0x5c(%rsp),%esi
  mov    %esi,0x38(%rsp)
  mov    0x50(%rsp),%rbx
  add    %rbx,0x30(%rsp)
  mov    %rdi,%rbp
L35:
  mov    0x18(%rsp),%rdi
  mov    %edi,0x14(%rsp)
  mov    0x20(%rsp),%rbx
  mov    %ebx,%r12d
  mov    %ebx,0x3c(%rsp)
  cmp    %ebx,%r15d
  jle    L50
  movslq 0x38(%rsp),%rax
  add    %rdi,%rax
  mov    0x28(%rsp),%rsi
  lea    (%rsi,%rax,8),%r14
  mov    %edi,%r13d
  mov    %ebx,0x58(%rsp)
  mov    %rbp,0x40(%rsp)
  mov    %rsi,%r12
L94:
  movsd  (%r14),%xmm0
  call   <fabs@plt>
  mov    %rax,%rbp
  mov    %r15d,%eax
  imul   %r13d,%eax
  add    0x14(%rsp),%eax
  cltq
  movsd  (%r12,%rax,8),%xmm0
  call   <fabs@plt>
  cmp    %rax,%rbp
  cmovg  %ebx,%r13d
  add    $0x1,%ebx
  add    0x8(%rsp),%r14
  cmp    %ebx,%r15d
  jne    L94
  mov    0x58(%rsp),%r12d
  mov    0x40(%rsp),%rbp
L89:
  mov    0x48(%rsp),%rax
  mov    0x18(%rsp),%rdi
  mov    %r13d,(%rax,%rdi,4)
  mov    %r15d,%ecx
  imul   %r13d,%ecx
  mov    0x14(%rsp),%eax
  add    %ecx,%eax
  cltq
  mov    0x28(%rsp),%rbx
  lea    (%rbx,%rax,8),%rax
  movsd  (%rax),%xmm1
  mov    0x38(%rsp),%esi
  sub    %r15d,%esi
  mov    0x30(%rsp),%rbx
  movsd  (%rbx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  %xmm1,(%rbx)
  pxor   %xmm2,%xmm2
  ucomisd %xmm2,%xmm1
  jp     L70
  je     L33
L70:
  cmp    0x3c(%rsp),%r15d
  jle    L15
  movslq 0x38(%rsp),%rax
  add    0x18(%rsp),%rax
  mov    0x28(%rsp),%rdi
  lea    (%rdi,%rax,8),%rax
  mov    %r12d,%edx
L24:
  movsd  (%rax),%xmm0
  xorpd  D3(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%edx
  add    0x8(%rsp),%rax
  cmp    %edx,%r15d
  jne    L24
  cmp    0x18(%rsp),%r13d
  je     L6
  movslq %ecx,%rcx
  mov    0x28(%rsp),%rax
  lea    (%rax,%rcx,8),%rdx
  mov    0x20(%rsp),%rax
L103:
  movsd  (%rdx,%rax,8),%xmm0
  movsd  0x0(%rbp,%rax,8),%xmm1
  movsd  %xmm1,(%rdx,%rax,8)
  movsd  %xmm0,0x0(%rbp,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r15d
  jg     L103
  jmp    L144
L59:
  mov    $0x1,%eax
  jmp    L8
L33:
  mov    $0x0,%eax
L8:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L15:
  mov    $0x1,%eax
  jmp    L8
```

tool_calls: name: read_data
arguments: data_tag: D3
datatype: byte[16]

tool: {"D3": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
int gauss_make_triangular(double *A, int *p, int n) {
  p[n - 1] = n - 1;
  for (int k = 0; k < n; ++k) {
    int m = k;
    for (int i = k + 1; i < n; ++i)
      if (fabs(A[k + n * i]) > fabs(A[k + n * m]))
        m = i;
    p[k] = m;
    double t1 = A[k + n * m];
    A[k + n * m] = A[k + n * k];
    A[k + n * k] = t1;
    if (t1 != 0) {
      for (int i = k + 1; i < n; ++i)
        A[k + n * i] /= -t1;
      if (k != m)
        for (int i = k + 1; i < n; ++i) {
          double t2 = A[i + n * m];
          A[i + n * m] = A[i + n * k];
          A[i + n * k] = t2;
        }
      for (int j = k + 1; j < n; ++j)
        for (int i = k + 1; i < n; ++i)
          A[i + n * j] += A[k + j * n] * A[i + k * n];
    } else
      return 0;
  }
  return 1;
}
```