user: 请从给定的汇编指令推导出C源代码：

```asm
<dual_simplex_method>:
  endbr64
  push   %r15
  mov    %edx,%r15d
  push   %r14
  lea    0x1(%rsi),%r14d
  push   %r13
  movslq %r14d,%r14
  push   %r12
  lea    0x0(,%r14,8),%rax
  push   %rbp
  mov    %rdi,%rbp
  mov    %rax,%rdi
  push   %rbx
  mov    %esi,%ebx
  sub    $0x58,%rsp
  mov    %esi,0x48(%rsp)
  mov    %rcx,0x20(%rsp)
  mov    %r8,0x8(%rsp)
  mov    %rax,0x28(%rsp)
  call   <malloc@plt>
  mov    %rax,0x10(%rsp)
  test   %ebx,%ebx
  js     L148
  mov    %rax,%rcx
  mov    %rax,%r14
  movslq 0x48(%rsp),%rax
  lea    0x1(%r15),%r12d
  movslq %r12d,%rbx
  mov    %rax,0x18(%rsp)
  shl    $0x3,%rbx
  lea    0x8(%rcx,%rax,8),%r13
  nopl   0x0(%rax)
L117:
  mov    %rbx,%rdi
  add    $0x8,%r14
  call   <malloc@plt>
  mov    %rax,-0x8(%r14)
  cmp    %r14,%r13
  jne    L117
  mov    0x48(%rsp),%eax
  test   %eax,%eax
  jle    L148
  movslq %r15d,%rax
  movsd  D120(%rip),%xmm1
  xor    %r13d,%r13d
  xor    %ecx,%ecx
  mov    0x18(%rsp),%r8
  mov    %rax,0x30(%rsp)
  pxor   %xmm3,%xmm3
  lea    0x0(,%rax,8),%rsi
  mov    $0x1,%edi
  xor    %eax,%eax
  cs nopw 0x0(%rax,%rax,1)
L201:
  mov    0x0(%rbp,%rax,8),%rdx
  movsd  (%rdx,%rsi,1),%xmm0
  comisd %xmm0,%xmm1
  jbe    L267
  movapd %xmm0,%xmm2
  ucomisd %xmm0,%xmm3
  cmpnltsd %xmm3,%xmm2
  cmova  %eax,%r13d
  cmova  %edi,%ecx
  andpd  %xmm2,%xmm1
  andnpd %xmm0,%xmm2
  orpd   %xmm2,%xmm1
L267:
  add    $0x1,%rax
  cmp    %rax,%r8
  jne    L201
  mov    $0x1,%eax
  test   %ecx,%ecx
  je     L97
  test   %r15d,%r15d
  jle    L176
  mov    0x28(%rsp),%rax
  movslq %r13d,%rsi
  mov    %r15d,%edx
  xor    %ecx,%ecx
  mov    0x0(%rbp,%rsi,8),%r8
  movsd  D97(%rip),%xmm6
  pxor   %xmm5,%xmm5
  mov    -0x8(%rbp,%rax,1),%rdi
  movq   D80(%rip),%xmm4
  lea    -0x1(%r15),%eax
  mov    %rax,0x28(%rsp)
  mov    %r8,%rbx
  xor    %eax,%eax
  movapd %xmm6,%xmm1
  nopl   0x0(%rax)
L287:
  movsd  (%r8,%rax,8),%xmm2
  movsd  (%rdi,%rax,8),%xmm0
  divsd  %xmm2,%xmm0
  andpd  %xmm4,%xmm0
  comisd %xmm0,%xmm1
  jbe    L235
  movapd %xmm2,%xmm3
  ucomisd %xmm2,%xmm5
  cmpnltsd %xmm5,%xmm3
  cmova  %eax,%ecx
  andpd  %xmm3,%xmm1
  andnpd %xmm0,%xmm3
  orpd   %xmm3,%xmm1
L235:
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L287
  ucomisd %xmm6,%xmm1
  jp     L242
  je     L176
L242:
  movslq %ecx,%rax
  mov    %rbx,%r14
  mov    %r15d,0x4c(%rsp)
  mov    0x10(%rsp),%r9
  movsd  (%rbx,%rax,8),%xmm2
  mov    0x8(%rsp),%rbx
  mov    %rbp,%r8
  xor    %r10d,%r10d
  lea    0x0(,%rax,8),%r11
  mov    0x48(%rsp),%r15d
  mov    %rbp,0x38(%rsp)
  lea    (%rbx,%rsi,4),%rdx
  mov    0x20(%rsp),%rbx
  mov    %rbp,0x40(%rsp)
  mov    %r11,%rbp
  mov    (%rdx),%esi
  movq   D122(%rip),%xmm4
  lea    (%rbx,%rax,4),%rax
  movsd  D8(%rip),%xmm1
  mov    %r14,%rbx
  mov    (%rax),%edi
  mov    %edi,(%rdx)
  mov    %esi,(%rax)
  nopl   0x0(%rax)
L149:
  cmp    %r13d,%r10d
  sete   %dil
  xor    %edx,%edx
  xor    %eax,%eax
  jmp    L151
  nopl   (%rax)
L56:
  mov    (%r8),%r11
  mov    (%r9),%rsi
  movapd %xmm1,%xmm0
  divsd  (%r11,%rdx,1),%xmm0
  movsd  %xmm0,(%rsi,%rdx,1)
L137:
  add    $0x1,%eax
  add    $0x8,%rdx
  cmp    %eax,%r12d
  je     L41
L151:
  cmp    %ecx,%eax
  sete   %sil
  jne    L99
  test   %dil,%dil
  jne    L56
L99:
  cmp    %ecx,%eax
  setne  %r11b
  je     L277
  test   %dil,%dil
  je     L277
  mov    (%r8),%r11
  mov    (%r9),%rsi
  add    $0x1,%eax
  movsd  (%r11,%rdx,1),%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rsi,%rdx,1)
  add    $0x8,%rdx
  cmp    %eax,%r12d
  jne    L151
L41:
  lea    0x1(%r10),%eax
  add    $0x8,%r8
  add    $0x8,%r9
  cmp    %r10d,%r15d
  je     L260
  mov    %eax,%r10d
  jmp    L149
  nopl   0x0(%rax)
L277:
  cmp    %r13d,%r10d
  setne  %r14b
  test   %sil,%sil
  je     L163
  test   %r14b,%r14b
  je     L163
  mov    (%r8),%r11
  mov    (%r9),%rsi
  movsd  (%r11,%rdx,1),%xmm0
  xorpd  %xmm4,%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rsi,%rdx,1)
  jmp    L137
  nop
L163:
  test   %r11b,%r11b
  je     L137
  test   %r14b,%r14b
  je     L137
  mov    (%r8),%r11
  mov    (%r9),%rsi
  movsd  (%r11,%rdx,1),%xmm0
  movsd  (%r11,%rbp,1),%xmm3
  mulsd  (%rbx,%rdx,1),%xmm3
  mulsd  %xmm2,%xmm0
  subsd  %xmm3,%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%rsi,%rdx,1)
  jmp    L137
  nop
L21:
  lea    D12(%rip),%rdi
  xor    %r12d,%r12d
  lea    D110(%rip),%rbp
  call   <puts@plt>
  lea    D40(%rip),%r14
  nopl   (%rax)
L200:
  xor    %r13d,%r13d
  nopl   0x0(%rax,%rax,1)
L71:
  mov    (%rbx),%rax
  mov    %rbp,%rdi
  movsd  (%rax,%r13,8),%xmm0
  mov    $0x1,%eax
  add    $0x1,%r13
  call   <printf@plt>
  cmp    %r13d,%r15d
  jge    L71
  cmp    %r12d,0x48(%rsp)
  je     L349
  mov    0x8(%rsp),%rax
  movslq (%rax,%r12,4),%rax
  mov    %rax,%rsi
  imul   $0x66666667,%rax,%rax
  mov    %esi,%edx
  sar    $0x1f,%edx
  sar    $0x22,%rax
  sub    %edx,%eax
  cmp    $0x1,%eax
  je     L393
  cmp    $0x2,%eax
  jne    L378
  sub    $0x13,%esi
  lea    D44(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
L378:
  add    $0x1,%r12
  add    $0x8,%rbx
  cmp    %r12d,0x48(%rsp)
  jge    L200
  xor    %eax,%eax
L97:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L260:
  mov    0x18(%rsp),%rax
  mov    0x38(%rsp),%rbx
  xor    %r14d,%r14d
  mov    0x40(%rsp),%rbp
  mov    0x4c(%rsp),%r15d
  lea    0x8(,%rax,8),%r13
  mov    0x30(%rsp),%rax
  mov    %rbx,0x18(%rsp)
  mov    %rbp,%rbx
  mov    0x10(%rsp),%rbp
  lea    0x8(,%rax,8),%r12
  nopl   0x0(%rax)
L366:
  mov    (%rbx,%r14,1),%rdi
  mov    0x0(%rbp,%r14,1),%rsi
  mov    %r12,%rdx
  add    $0x8,%r14
  call   <memcpy@plt>
  cmp    %r13,%r14
  jne    L366
  lea    D135(%rip),%rdi
  mov    0x18(%rsp),%rbx
  lea    D88(%rip),%r13
  call   <puts@plt>
  mov    0x20(%rsp),%r14
  mov    0x28(%rsp),%rax
  lea    D104(%rip),%r12
  lea    0x4(%r14,%rax,4),%rbp
  jmp    L309
  nopw   0x0(%rax,%rax,1)
L314:
  lea    -0x14(%rsi),%eax
  cmp    $0x9,%eax
  jbe    L35
L223:
  add    $0x4,%r14
  cmp    %r14,%rbp
  je     L21
L309:
  mov    (%r14),%esi
  lea    -0xa(%rsi),%eax
  cmp    $0x9,%eax
  ja     L314
  movslq %esi,%rax
  mov    %esi,%edx
  mov    %r12,%rdi
  imul   $0x66666667,%rax,%rax
  sar    $0x1f,%edx
  sar    $0x22,%rax
  sub    %edx,%eax
  lea    (%rax,%rax,4),%eax
  add    %eax,%eax
  sub    %eax,%esi
  xor    %eax,%eax
  add    $0x1,%esi
  call   <printf@plt>
  jmp    L223
  nopl   0x0(%rax,%rax,1)
L35:
  movslq %esi,%rax
  mov    %esi,%edx
  mov    %r13,%rdi
  imul   $0x66666667,%rax,%rax
  sar    $0x1f,%edx
  sar    $0x22,%rax
  sub    %edx,%eax
  lea    (%rax,%rax,4),%eax
  add    %eax,%eax
  sub    %eax,%esi
  xor    %eax,%eax
  add    $0x1,%esi
  call   <printf@plt>
  jmp    L223
  nopl   0x0(%rax,%rax,1)
L393:
  sub    $0x9,%esi
  mov    %r14,%rdi
  xor    %eax,%eax
  call   <printf@plt>
  jmp    L378
  nopw   0x0(%rax,%rax,1)
L349:
  lea    D29(%rip),%rdi
  call   <puts@plt>
  jmp    L378
L148:
  add    $0x58,%rsp
  mov    $0x1,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L176:
  lea    D134(%rip),%rdi
  call   <puts@plt>
  mov    $0x2,%eax
  jmp    L97
```

tool_calls: name: read_data
arguments: data_label: D120
type_infered: f64
name: read_data
arguments: data_label: D97
type_infered: f64
name: read_data
arguments: data_label: D80
type_infered: qword
name: read_data
arguments: data_label: D122
type_infered: byte[8]
name: read_data
arguments: data_label: D8
type_infered: byte[8]
name: read_data
arguments: data_label: D12
type_infered: string
name: read_data
arguments: data_label: D110
type_infered: string
name: read_data
arguments: data_label: D40
type_infered: string
name: read_data
arguments: data_label: D44
type_infered: string
name: read_data
arguments: data_label: D135
type_infered: string
name: read_data
arguments: data_label: D88
type_infered: string
name: read_data
arguments: data_label: D104
type_infered: string
name: read_data
arguments: data_label: D29
type_infered: string
name: read_data
arguments: data_label: D134
type_infered: string

tool: D120: 100

tool: D97: 10000

tool: D80: 9223372036854775807

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

tool: {"D8": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: D12: 1

tool: D110: %lf	

tool: D40: X%d


tool: D44: Z%d


tool: D135: 
Table: 

tool: D88: Z%d		

tool: D104: X%d		

tool: D29: Z

tool: D134: The solution is unbounded

assistant: ```c
int dual_simplex_method(double **T, int no_of_eqns, int no_of_variables, int *column, int *row) {
  int flag = 0, mp = 0, i;
  double **A;
  A = (double **)malloc((no_of_eqns + 1) * sizeof(double *));
  for (i = 0; i < no_of_eqns + 1; i++) {
    A[i] = (double *)malloc((no_of_variables + 1) * sizeof(double));
  }
  double min = 100;
  for (i = 0; i < no_of_eqns; i++) {
    if (T[i][no_of_variables] < min && T[i][no_of_variables] < 0) {
      flag = 1;
      min = T[i][no_of_variables];
      mp = i;
    }
  }
  int most_negative_row = mp, pivot_column;
  double most_negative_value = min, pivot_value;
  if (flag == 0) {
    return 1;
  } else {
    int min_ratio_column = 0;
    double min_ratio = 10000;
    flag = 0;
    for (i = 0; i < no_of_variables; i++) {
      if (fabs(T[no_of_eqns][i] / T[most_negative_row][i]) < min_ratio && T[most_negative_row][i] < 0) {
        min_ratio = fabs(T[no_of_eqns][i] / T[most_negative_row][i]);
        min_ratio_column = i;
      }
    }
    if (min_ratio == 10000) {
      printf("The solution is unbounded\n");
      return 2;
    } else {
      pivot_column = min_ratio_column;
      pivot_value = T[most_negative_row][min_ratio_column];
      int temp1 = row[most_negative_row];
      row[most_negative_row] = column[min_ratio_column];
      column[min_ratio_column] = temp1;
      for (i = 0; i < no_of_eqns + 1; i++) {
        int j;
        for (j = 0; j < no_of_variables + 1; j++) {
          if (i == most_negative_row && j == min_ratio_column) {
            A[i][j] = 1 / T[i][j];
          } else if (i == most_negative_row && j != min_ratio_column) {
            A[i][j] = T[i][j] / pivot_value;
          } else if (j == min_ratio_column && i != most_negative_row) {
            A[i][j] = -T[i][j] / pivot_value;
          } else if (j != min_ratio_column && i != most_negative_row) {
            double temp = T[i][j] * pivot_value - T[most_negative_row][j] * T[i][min_ratio_column];
            A[i][j] = temp / ((double)pivot_value);
          }
        }
      }
    }
  }
  for (i = 0; i < no_of_eqns + 1; i++) {
    int j;
    for (j = 0; j < no_of_variables + 1; j++) {
      T[i][j] = A[i][j];
    }
  }
  printf("\nTable: \n");
  int k, l;
  for (l = 0; l < no_of_variables; l++) {
    if (column[l] / 10 == 1) {
      printf("X%d\t\t", (column[l] % 10) + 1);
    } else if (column[l] / 10 == 2) {
      printf("Z%d\t\t", (column[l] % 10) + 1);
    }
  }
  printf("1\n");
  for (k = 0; k < no_of_eqns + 1; k++) {
    for (l = 0; l < no_of_variables + 1; l++) {
      printf("%lf\t", T[k][l]);
    }
    if (k == no_of_eqns) {
      printf("Z\n");
    } else {
      switch (row[k] / 10) {
      case 1:
        printf("X%d\n", (row[k] % 10) + 1);
        break;
      case 2:
        printf("Z%d\n", (row[k] % 10) + 1);
        break;
      }
    }
  }
  return 0;
}
```