user: 请将以下的汇编代码转换成C语言代码:

```asm
<c_sgbco>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %rdi,0x18(%rsp)
  mov    %esi,0x8(%rsp)
  mov    %edx,%r13d
  mov    %ecx,0x58(%rsp)
  mov    %r8d,0x30(%rsp)
  mov    %r9,0x40(%rsp)
  mov    0xe0(%rsp),%rax
  mov    %rax,0x68(%rsp)
  mov    0xe8(%rsp),%r15
  mov    %fs:0x28,%rax
  mov    %rax,0x98(%rsp)
  xor    %eax,%eax
  lea    0x1(%rcx),%ebp
  mov    %ebp,0x10(%rsp)
  mov    %ebp,%r12d
  add    %r8d,%r12d
  test   %edx,%edx
  jle    L353
  mov    $0x0,%r14d
  pxor   %xmm6,%xmm6
  movsd  %xmm6,0x28(%rsp)
  mov    $0x1,%ebx
  sub    %ecx,%edx
  mov    %edx,0x20(%rsp)
  mov    %r15,0x38(%rsp)
L101:
  lea    -0x1(%r12),%r15d
  lea    (%r14,%r15,1),%eax
  cltq
  mov    0x18(%rsp),%rcx
  lea    (%rcx,%rax,8),%rsi
  mov    %ebp,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  0x28(%rsp),%xmm7
  maxsd  %xmm0,%xmm7
  movsd  %xmm7,0x28(%rsp)
  cmp    %r12d,0x10(%rsp)
  cmovl  %r15d,%r12d
  cmp    %ebx,0x30(%rsp)
  setge  %al
  movzbl %al,%eax
  add    %eax,%ebp
  cmp    %ebx,0x20(%rsp)
  setle  %al
  movzbl %al,%eax
  sub    %eax,%ebp
  add    $0x1,%ebx
  add    0x8(%rsp),%r14d
  cmp    %ebx,%r13d
  jge    L101
  mov    0x38(%rsp),%r15
  sub    $0x8,%rsp
  lea    0x9c(%rsp),%rax
  push   %rax
  mov    0x50(%rsp),%r9
  mov    0x40(%rsp),%r12d
  mov    %r12d,%r8d
  mov    0x68(%rsp),%ebx
  mov    %ebx,%ecx
  mov    %r13d,%edx
  mov    0x18(%rsp),%ebp
  mov    %ebp,%esi
  mov    0x28(%rsp),%r14
  mov    %r14,%rdi
  mov    $0x0,%eax
  call   <c_sgbfa@plt>
  movslq %r13d,%rax
  mov    %rax,0x80(%rsp)
  shl    $0x3,%rax
  mov    %rax,0x88(%rsp)
  mov    %rax,%rdx
  mov    $0x0,%esi
  mov    %r15,%rdi
  call   <memset@plt>
  add    %r12d,%ebx
  mov    %ebx,0x94(%rsp)
  lea    0x1(%rbx),%eax
  mov    %eax,0x90(%rsp)
  movslq %ebp,%rax
  lea    0x0(,%rax,8),%rsi
  mov    %rsi,0x60(%rsp)
  movslq %ebx,%rdx
  mov    %rdx,0x98(%rsp)
  lea    (%r14,%rdx,8),%r12
  lea    -0x8(,%rax,8),%rax
  mov    %rax,0x30(%rsp)
  add    $0x10,%rsp
  mov    $0x1,%ebx
  movsd  D75(%rip),%xmm7
  movsd  %xmm7,0x10(%rsp)
  mov    $0x1,%eax
  mov    $0x0,%r14d
  mov    %ebp,0x48(%rsp)
  mov    %eax,%ebp
  jmp    L32
L306:
  pxor   %xmm1,%xmm1
  comisd %xmm2,%xmm1
  jae    L281
  movsd  0x10(%rsp),%xmm0
  andpd  D102(%rip),%xmm0
  xorpd  D236(%rip),%xmm0
  movsd  %xmm0,0x10(%rsp)
  jmp    L225
L281:
  movsd  0x10(%rsp),%xmm4
  andpd  D102(%rip),%xmm4
  movsd  %xmm4,0x10(%rsp)
  jmp    L225
L184:
  divsd  %xmm0,%xmm1
  mov    %r15,%rsi
  movsd  %xmm1,0x60(%rsp)
  movapd %xmm1,%xmm0
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  movsd  0x10(%rsp),%xmm3
  mulsd  0x60(%rsp),%xmm3
  movsd  %xmm3,0x10(%rsp)
  jmp    L398
L134:
  movapd %xmm3,%xmm6
  divsd  %xmm0,%xmm6
  movapd %xmm2,%xmm7
  divsd  %xmm0,%xmm7
L274:
  add    $0x1,%ebp
  mov    0x40(%rsp),%rcx
  mov    0x30(%rsp),%eax
  add    -0x4(%rcx,%rbx,4),%eax
  cmp    %r14d,%eax
  cmovl  %r14d,%eax
  cmp    %r13d,%eax
  cmovg  %r13d,%eax
  mov    %eax,%r14d
  cmp    %eax,%ebp
  jle    L267
L116:
  movsd  %xmm6,-0x8(%r15,%rbx,8)
  add    $0x1,%rbx
  mov    0x8(%rsp),%ecx
  add    %ecx,0x48(%rsp)
  add    0x50(%rsp),%r12
  cmp    %ebp,%r13d
  jl     L387
L32:
  movsd  -0x8(%r15,%rbx,8),%xmm2
  ucomisd D204(%rip),%xmm2
  jp     L306
  jne    L306
L225:
  movsd  0x10(%rsp),%xmm0
  subsd  %xmm2,%xmm0
  andpd  D102(%rip),%xmm0
  mov    %r12,0x38(%rsp)
  movsd  (%r12),%xmm1
  andpd  D102(%rip),%xmm1
  comisd %xmm1,%xmm0
  ja     L184
L398:
  movsd  -0x8(%r15,%rbx,8),%xmm0
  movsd  0x10(%rsp),%xmm5
  movapd %xmm5,%xmm3
  subsd  %xmm0,%xmm3
  xorpd  D236(%rip),%xmm5
  movapd %xmm5,%xmm2
  subsd  %xmm0,%xmm2
  mov    0x38(%rsp),%rax
  movsd  (%rax),%xmm0
  ucomisd D204(%rip),%xmm0
  jp     L134
  jne    L134
  movsd  D75(%rip),%xmm7
  movapd %xmm7,%xmm6
  jmp    L274
L267:
  movq   D102(%rip),%xmm0
  andpd  %xmm0,%xmm3
  andpd  %xmm0,%xmm2
  movslq %ebp,%rax
  lea    -0x8(%r15,%rax,8),%rcx
  movslq 0x80(%rsp),%rax
  movslq 0x48(%rsp),%rdx
  lea    -0x2(%rax,%rdx,1),%rax
  mov    0x18(%rsp),%rdx
  lea    (%rdx,%rax,8),%rdi
  mov    %r14d,%eax
  sub    %ebx,%eax
  lea    -0x1(%rax),%eax
  add    %rbx,%rax
  lea    0x8(%r15,%rax,8),%rsi
  mov    %rdi,%rdx
  mov    %rcx,%rax
  movapd %xmm0,%xmm5
L182:
  movsd  (%rax),%xmm4
  movsd  (%rdx),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm7,%xmm1
  addsd  %xmm4,%xmm1
  andpd  %xmm5,%xmm1
  addsd  %xmm1,%xmm2
  mulsd  %xmm6,%xmm0
  addsd  %xmm4,%xmm0
  movsd  %xmm0,(%rax)
  andpd  %xmm5,%xmm0
  addsd  %xmm0,%xmm3
  add    $0x8,%rax
  add    0x20(%rsp),%rdx
  cmp    %rsi,%rax
  jne    L182
  comisd %xmm3,%xmm2
  jbe    L116
  movapd %xmm7,%xmm1
  subsd  %xmm6,%xmm1
L298:
  movapd %xmm1,%xmm0
  mulsd  (%rdi),%xmm0
  addsd  (%rcx),%xmm0
  movsd  %xmm0,(%rcx)
  add    $0x8,%rcx
  add    0x20(%rsp),%rdi
  cmp    %rsi,%rcx
  jne    L298
  movapd %xmm7,%xmm6
  jmp    L116
L219:
  mov    0x20(%rsp),%eax
  sub    %ebx,%eax
  cltq
  mov    0x10(%rsp),%rdx
  lea    -0x8(%rdx,%rax,8),%rdx
  movslq %r15d,%rax
  mov    0x18(%rsp),%rcx
  lea    (%rcx,%rax,8),%rsi
  mov    %ebx,%edi
  mov    0x58(%rsp),%eax
  cmp    %ebx,%eax
  cmovle %eax,%edi
  mov    $0x0,%eax
  call   <c_sdot@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  addsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  jmp    L197
L162:
  movslq -0x4(%r12),%rax
  mov    0x10(%rsp),%rcx
  lea    -0x8(%rcx,%rax,8),%rax
  movsd  (%rax),%xmm0
  movsd  -0x8(%r14),%xmm1
  movsd  %xmm1,(%rax)
  movsd  %xmm0,-0x8(%r14)
  add    $0x1,%rbx
  sub    $0x4,%r12
  sub    $0x8,%rbp
  sub    0x8(%rsp),%r15d
  lea    0x1(%rbx),%eax
  cmp    %eax,%r13d
  jl     L58
L155:
  mov    %r13d,%eax
  sub    %ebx,%eax
  cmp    %eax,%r13d
  jg     L219
L197:
  mov    %rbp,%r14
  movsd  -0x8(%rbp),%xmm0
  andpd  D102(%rip),%xmm0
  comisd D75(%rip),%xmm0
  jbe    L162
  movsd  D75(%rip),%xmm1
  divsd  %xmm0,%xmm1
  mov    0x10(%rsp),%rsi
  movapd %xmm1,%xmm0
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  jmp    L162
L58:
  mov    0x30(%rsp),%r14
  mov    %rcx,%r15
  mov    %rcx,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D75(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  mov    0x80(%rsp),%ebp
  mov    $0x1,%ebx
  mov    D75(%rip),%rax
  mov    %rax,0x10(%rsp)
  mov    %r14,0x30(%rsp)
  mov    0x58(%rsp),%r14d
  mov    0x40(%rsp),%r12
  jmp    L53
L340:
  movsd  -0x8(%r15,%rbx,8),%xmm0
  andpd  D102(%rip),%xmm0
  comisd D75(%rip),%xmm0
  ja     L215
L6:
  add    $0x1,%rbx
  add    0x8(%rsp),%ebp
  cmp    %ebx,%r13d
  jl     L146
L53:
  movslq -0x4(%r12,%rbx,4),%rax
  lea    -0x8(%r15,%rax,8),%rax
  movsd  (%rax),%xmm0
  movsd  -0x8(%r15,%rbx,8),%xmm1
  movsd  %xmm1,(%rax)
  movsd  %xmm0,-0x8(%r15,%rbx,8)
  cmp    %ebx,%r13d
  jle    L340
  lea    0x1(%rbx),%eax
  cltq
  lea    -0x8(%r15,%rax,8),%rdx
  movslq %ebp,%rax
  mov    0x18(%rsp),%rcx
  lea    (%rcx,%rax,8),%rsi
  mov    %r13d,%edi
  sub    %ebx,%edi
  cmp    %r14d,%edi
  cmovg  %r14d,%edi
  mov    $0x1,%eax
  call   <c_saxpy@plt>
  jmp    L340
L215:
  movsd  D75(%rip),%xmm4
  divsd  %xmm0,%xmm4
  mov    %r15,%rsi
  movsd  %xmm4,0x20(%rsp)
  movapd %xmm4,%xmm0
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  movsd  0x10(%rsp),%xmm3
  mulsd  0x20(%rsp),%xmm3
  movsd  %xmm3,0x10(%rsp)
  jmp    L6
L146:
  mov    0x30(%rsp),%r14
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  D75(%rip),%xmm1
  divsd  %xmm0,%xmm1
  mov    %r15,%rsi
  movsd  %xmm1,0x20(%rsp)
  movapd %xmm1,%xmm0
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  movsd  0x20(%rsp),%xmm7
  mulsd  0x10(%rsp),%xmm7
  movsd  %xmm7,0x50(%rsp)
  mov    0x8(%rsp),%edx
  mov    %edx,%eax
  neg    %eax
  cltq
  shl    $0x3,%rax
  mov    %rax,0x48(%rsp)
  lea    -0x1(%r13),%r12d
  imul   %edx,%r12d
  movslq %r12d,%rax
  add    0x88(%rsp),%rax
  mov    0x18(%rsp),%rsi
  lea    (%rsi,%rax,8),%rbp
  movl   $0x0,0x10(%rsp)
  mov    %r13d,0x20(%rsp)
  mov    %r15,0x40(%rsp)
  mov    0x80(%rsp),%r15d
  jmp    L183
L20:
  divsd  %xmm0,%xmm1
  mov    0x40(%rsp),%rsi
  movsd  %xmm1,0x58(%rsp)
  movapd %xmm1,%xmm0
  mov    0x20(%rsp),%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  movsd  0x50(%rsp),%xmm6
  mulsd  0x58(%rsp),%xmm6
  movsd  %xmm6,0x50(%rsp)
  jmp    L287
L229:
  movsd  -0x8(%r13),%xmm0
  divsd  %xmm1,%xmm0
L212:
  movsd  %xmm0,-0x8(%r13)
  cmp    %r15d,%ebx
  mov    %r15d,%edi
  cmovle %ebx,%edi
  sub    $0x1,%edi
  sub    %edi,%ebx
  movslq %ebx,%rbx
  mov    0x40(%rsp),%rax
  lea    -0x8(%rax,%rbx,8),%rdx
  mov    %r15d,%eax
  sub    %edi,%eax
  mov    0x30(%rsp),%ecx
  lea    -0x1(%rcx,%rax,1),%eax
  cltq
  mov    0x18(%rsp),%rcx
  lea    (%rcx,%rax,8),%rsi
  xorpd  D236(%rip),%xmm0
  mov    $0x1,%eax
  call   <c_saxpy@plt>
  mov    0x10(%rsp),%eax
  lea    0x1(%rax),%edx
  sub    $0x8,%r14
  add    0x48(%rsp),%rbp
  sub    0x8(%rsp),%r12d
  add    $0x2,%eax
  cmp    %eax,0x20(%rsp)
  jl     L115
  mov    %edx,0x10(%rsp)
L183:
  mov    0x20(%rsp),%ebx
  sub    0x10(%rsp),%ebx
  mov    %r14,%r13
  movsd  -0x8(%r14),%xmm0
  andpd  D102(%rip),%xmm0
  mov    %r12d,0x30(%rsp)
  mov    %rbp,0x38(%rsp)
  movsd  0x0(%rbp),%xmm1
  andpd  D102(%rip),%xmm1
  comisd %xmm1,%xmm0
  ja     L20
L287:
  mov    0x38(%rsp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm6,%xmm6
  ucomisd %xmm6,%xmm1
  jp     L229
  movsd  D75(%rip),%xmm0
  je     L212
  jmp    L229
L115:
  mov    0x20(%rsp),%r13d
  mov    0x40(%rsp),%r15
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  D75(%rip),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rsp)
  mov    %r15,%rsi
  movapd %xmm1,%xmm0
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  movsd  0x28(%rsp),%xmm3
  ucomisd D204(%rip),%xmm3
  jp     L152
  pxor   %xmm0,%xmm0
  je     L74
L152:
  movsd  0x50(%rsp),%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  0x28(%rsp),%xmm0
L74:
  mov    0x68(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x98(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L258
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L387:
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D75(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  mov    0x40(%rsp),%rax
  mov    0x70(%rsp),%rdx
  lea    (%rax,%rdx,4),%r12
  mov    0x78(%rsp),%r14
  add    %r15,%r14
  lea    -0x1(%r13),%eax
  imul   0x8(%rsp),%eax
  mov    0x84(%rsp),%edx
  lea    0x1(%rdx,%rax,1),%eax
  mov    %r14,%rbp
  mov    $0x0,%ebx
  lea    0x1(%r13),%esi
  mov    %esi,0x20(%rsp)
  mov    %r14,0x30(%rsp)
  mov    %r15,0x10(%rsp)
  mov    %eax,%r15d
  jmp    L155
L353:
  sub    $0x8,%rsp
  lea    0x9c(%rsp),%rax
  push   %rax
  mov    0x50(%rsp),%r9
  mov    0x40(%rsp),%r8d
  mov    0x68(%rsp),%ecx
  mov    0x18(%rsp),%esi
  mov    0x28(%rsp),%rdi
  mov    $0x0,%eax
  call   <c_sgbfa@plt>
  movslq %r13d,%rdx
  shl    $0x3,%rdx
  mov    $0x0,%esi
  mov    %r15,%rdi
  call   <memset@plt>
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D75(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D75(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D75(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x0,%eax
  call   <c_sasum@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D75(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    %r15,%rsi
  mov    %r13d,%edi
  mov    $0x1,%eax
  call   <c_sscal@plt>
  add    $0x10,%rsp
  pxor   %xmm0,%xmm0
  jmp    L74
L258:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: label: D75
datatype: f64[2]
name: unpack_data
arguments: label: D102
datatype: qword[2]
name: unpack_data
arguments: label: D236
datatype: f32[4]
name: unpack_data
arguments: label: D204
datatype: f32[2]

tool: {"D75": [1.0,0.0]}

tool: {"D102": [0x7fffffffffffffff,0x0000000000000000]}

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

tool: {"D204": [0.0,0.0]}

assistant: ```c
void c_sgbco(double *abd, int lda, int n, int ml, int mu, int *ipvt, double *rcond, double *z) {
  int info;
  register int is, j, ju, k, kb, kp1, l, la, lm, lz, m, mm;
  double anorm, ek, s, sm, t, wk, wkm, ynorm;
  anorm = 0.;
  l = ml + 1;
  is = l + mu;
  for (j = 1; j <= n; j++) {
    anorm = ({
      const double _x = (double)(anorm);
      const double _y = (double)(c_sasum(l, &abd[is - 1 + (j - 1) * lda]));
      _x > _y ? _x : _y;
    });
    if (is > ml + 1) {
      is--;
    }
    if (j <= mu) {
      l++;
    }
    if (j >= n - ml) {
      l--;
    }
  }
  c_sgbfa(abd, lda, n, ml, mu, ipvt, &info);
  ek = 1.;
  memset(z, 0, n * sizeof(double));
  m = ml + mu + 1;
  ju = 0;
  for (k = 1; k <= n; k++) {
    if (z[k - 1] != 0.) {
      ek = ((-z[k - 1]) >= 0. ? fabs(ek) : -fabs(ek));
    }
    if (fabs(ek - z[k - 1]) > fabs(abd[m - 1 + (k - 1) * lda])) {
      s = fabs(abd[m - 1 + (k - 1) * lda]) / fabs(ek - z[k - 1]);
      c_sscal(n, s, z);
      ek *= s;
    }
    wk = ek - z[k - 1];
    wkm = -ek - z[k - 1];
    s = fabs(wk);
    sm = fabs(wkm);
    if (abd[m - 1 + (k - 1) * lda] != 0.) {
      wk /= abd[m - 1 + (k - 1) * lda];
      wkm /= abd[m - 1 + (k - 1) * lda];
    } else {
      wk = 1.;
      wkm = 1.;
    }
    kp1 = k + 1;
    ju = ({
      const int _i = (int)(({
        const int _i = (int)(ju);
        const int _j = (int)(mu + ipvt[k - 1]);
        _i > _j ? _i : _j;
      }));
      const int _j = (int)(n);
      _i < _j ? _i : _j;
    });
    mm = m;
    if (kp1 <= ju) {
      for (j = kp1; j <= ju; j++) {
        mm--;
        sm += fabs(z[j - 1] + wkm * abd[mm - 1 + (j - 1) * lda]);
        z[j - 1] += wk * abd[mm - 1 + (j - 1) * lda];
        s += fabs(z[j - 1]);
      }
      if (s < sm) {
        t = wkm - wk;
        wk = wkm;
        mm = m;
        for (j = kp1; j <= ju; j++) {
          mm--;
          z[j - 1] += t * abd[mm - 1 + (j - 1) * lda];
        }
      }
    }
    z[k - 1] = wk;
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  for (kb = 1; kb <= n; kb++) {
    k = n + 1 - kb;
    lm = ({
      const int _i = (int)(ml);
      const int _j = (int)(n - k);
      _i < _j ? _i : _j;
    });
    if (k < n) {
      z[k - 1] += c_sdot(lm, &abd[m + 1 - 1 + (k - 1) * lda], &z[k + 1 - 1]);
    }
    if (fabs(z[k - 1]) > 1.) {
      s = 1. / fabs(z[k - 1]);
      c_sscal(n, s, z);
    }
    l = ipvt[k - 1];
    t = z[l - 1];
    z[l - 1] = z[k - 1];
    z[k - 1] = t;
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  ynorm = 1.;
  for (k = 1; k <= n; k++) {
    l = ipvt[k - 1];
    t = z[l - 1];
    z[l - 1] = z[k - 1];
    z[k - 1] = t;
    lm = ({
      const int _i = (int)(ml);
      const int _j = (int)(n - k);
      _i < _j ? _i : _j;
    });
    if (k < n) {
      c_saxpy(lm, t, &abd[m + 1 - 1 + (k - 1) * lda], &z[k + 1 - 1]);
    }
    if (fabs(z[k - 1]) > 1.) {
      s = 1. / fabs(z[k - 1]);
      c_sscal(n, s, z);
      ynorm *= s;
    }
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  ynorm *= s;
  for (kb = 1; kb <= n; kb++) {
    k = n + 1 - kb;
    if (fabs(z[k - 1]) > fabs(abd[m - 1 + (k - 1) * lda])) {
      s = fabs(abd[m - 1 + (k - 1) * lda]) / fabs(z[k - 1]);
      c_sscal(n, s, z);
      ynorm *= s;
    }
    if (abd[m - 1 + (k - 1) * lda] != 0.) {
      z[k - 1] /= abd[m - 1 + (k - 1) * lda];
    } else {
      z[k - 1] = 1.;
    }
    lm = ({
           const int _i = (int)(k);
           const int _j = (int)(m);
           _i < _j ? _i : _j;
         }) -
         1;
    la = m - lm;
    lz = k - lm;
    t = -z[k - 1];
    c_saxpy(lm, t, &abd[la - 1 + (k - 1) * lda], &z[lz - 1]);
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  ynorm *= s;
  if (anorm != 0.) {
    *rcond = ynorm / anorm;
  } else {
    *rcond = 0.;
  }
  return;
}
```