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

```asm
<bpflow>:
  endbr64
  push   %r15
  mov    %esi,%eax
  mov    %rcx,%r15
  push   %r14
  push   %r13
  mov    %rdx,%r13
  push   %r12
  lea    (%rdi,%rsi,1),%r12d
  push   %rbp
  push   %rbx
  mov    %r9,%rbx
  sub    $0x28,%rsp
  mov    0x60(%rsp),%r14
  mov    0x68(%rsp),%r9
  mov    %edi,(%rsp)
  test   %r12d,%r12d
  jle    L171
  lea    0x4(%r8),%rdi
  mov    %r9,%rcx
  mov    %rbx,%r10
  mov    %r8,%rsi
  sub    %rdi,%rcx
  cmp    $0x8,%rcx
  seta   %cl
  sub    %rdi,%r10
  cmp    $0x8,%r10
  lea    0x4(%rbx),%r10
  seta   %dil
  and    %edi,%ecx
  mov    %r9,%rdi
  sub    %r10,%rdi
  cmp    $0x8,%rdi
  seta   %dil
  test   %dil,%cl
  je     L82
  cmp    $0x1,%r12d
  je     L82
  lea    -0x1(%r12),%ecx
  mov    %r12d,%r11d
  cmp    $0x2,%ecx
  jbe    L0
  mov    %r12d,%edi
  xor    %ecx,%ecx
  pxor   %xmm0,%xmm0
  shr    $0x2,%edi
  shl    $0x4,%rdi
  nopl   0x0(%rax)
L239:
  movdqu (%rsi,%rcx,1),%xmm1
  movups %xmm1,(%rbx,%rcx,1)
  movups %xmm0,(%r9,%rcx,1)
  add    $0x10,%rcx
  cmp    %rdi,%rcx
  jne    L239
  mov    %r12d,%r10d
  movslq %r12d,%rdi
  and    $0xfffffffc,%r10d
  mov    %r10d,%ecx
  cmp    %r12d,%r10d
  je     L265
  mov    %r12d,%r11d
  sub    %r10d,%r11d
  cmp    $0x1,%r11d
  je     L2
L176:
  mov    %r10d,%edi
  mov    (%rsi,%rdi,4),%r10
  mov    %r10,(%rbx,%rdi,4)
  mov    %r11d,%r10d
  and    $0xfffffffe,%r10d
  movq   $0x0,(%r9,%rdi,4)
  movslq %r12d,%rdi
  add    %r10d,%ecx
  cmp    %r10d,%r11d
  je     L265
L2:
  movslq %ecx,%rcx
  movslq %r12d,%rdi
  mov    (%rsi,%rcx,4),%esi
  mov    %esi,(%rbx,%rcx,4)
  movl   $0x0,(%r9,%rcx,4)
L265:
  movslq (%rsp),%rdx
  lea    0x0(%r13,%rdi,4),%rsi
  lea    0x0(,%rdx,4),%r10
  mov    %rdx,0x10(%rsp)
  lea    0x0(%r13,%r10,1),%rbp
  movslq 0x0(%rbp),%rcx
  cmp    (%rsi),%ecx
  jge    L135
  nopl   0x0(%rax)
L128:
  movl   $0x0,(%r14,%rcx,4)
  add    $0x1,%rcx
  cmp    %ecx,(%rsi)
  jg     L128
L135:
  cmp    %r12d,(%rsp)
  jge    L42
  sub    $0x1,%eax
  add    0x10(%rsp),%rax
  add    %rbx,%r10
  lea    0x4(%rbx,%rax,4),%rdx
  cs nopw 0x0(%rax,%rax,1)
L27:
  movslq 0x0(%rbp),%rax
  mov    0x4(%rbp),%r8d
  cmp    %r8d,%eax
  jl     L9
  jmp    L132
  nop
L169:
  cmp    %ecx,%esi
  cmovle %esi,%ecx
  sub    %ecx,%esi
  mov    %esi,(%rdi)
  sub    %ecx,(%r10)
  mov    %ecx,(%r14,%rax,4)
  mov    0x4(%rbp),%r8d
  add    $0x1,%rax
  cmp    %eax,%r8d
  jle    L132
L9:
  mov    (%r10),%ecx
L216:
  test   %ecx,%ecx
  je     L132
  movslq (%r15,%rax,4),%rsi
  lea    (%rbx,%rsi,4),%rdi
  mov    (%rdi),%esi
  test   %esi,%esi
  jne    L169
  add    $0x1,%rax
  cmp    %eax,%r8d
  jg     L216
L132:
  add    $0x4,%r10
  add    $0x4,%rbp
  cmp    %r10,%rdx
  jne    L27
L42:
  lea    0x0(,%r12,4),%edi
  mov    %r9,0x18(%rsp)
  call   <smalloc@plt>
  mov    (%rsp),%edi
  mov    %rax,0x8(%rsp)
  test   %edi,%edi
  jle    L194
  xor    %ebp,%ebp
  mov    %rbp,(%rsp)
  mov    0x18(%rsp),%rbp
  nopl   0x0(%rax)
L10:
  mov    (%rsp),%rax
  mov    %eax,%r12d
  mov    (%rbx,%rax,4),%eax
  test   %eax,%eax
  jne    L34
  jmp    L173
  nopl   0x0(%rax,%rax,1)
L84:
  test   %eax,%eax
  je     L173
L34:
  sub    $0x8,%rsp
  mov    %rbx,%rcx
  mov    %r15,%rdx
  mov    %r13,%rsi
  push   0x10(%rsp)
  mov    %rbp,%r9
  mov    %r14,%r8
  mov    %r12d,%edi
  call   <augment@plt>
  pop    %rdx
  pop    %rcx
  mov    (%rsp),%rdx
  mov    (%rbx,%rdx,4),%esi
  test   %esi,%esi
  jne    L84
L173:
  addq   $0x1,(%rsp)
  mov    (%rsp),%rax
  cmp    0x10(%rsp),%rax
  jne    L10
L194:
  mov    0x8(%rsp),%rdi
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <sfree@plt>
L82:
  movslq %r12d,%rdi
  xor    %ecx,%ecx
  nopl   (%rax)
L228:
  mov    (%rsi,%rcx,4),%edx
  mov    %edx,(%rbx,%rcx,4)
  movl   $0x0,(%r9,%rcx,4)
  add    $0x1,%rcx
  cmp    %rdi,%rcx
  jne    L228
  jmp    L265
L171:
  movslq %r12d,%rdi
  jmp    L265
L0:
  xor    %ecx,%ecx
  xor    %r10d,%r10d
  jmp    L176
```

assistant: ```c
void bpflow(int n_left, int n_right, int *pointers, int *indices, int *vweight, int *resid, int *flow, int *touched) {
  int *seen;
  int neighbor;
  int flow1;
  int modified;
  int i, j;
  for (i = 0; i < n_left + n_right; i++) {
    resid[i] = vweight[i];
    touched[i] = 0;
  }
  for (i = pointers[n_left]; i < pointers[n_left + n_right]; i++) {
    flow[i] = 0;
  }
  for (i = n_left; i < n_left + n_right; i++) {
    for (j = pointers[i]; j < pointers[i + 1] && resid[i]; j++) {
      neighbor = indices[j];
      if (resid[neighbor]) {
        flow1 = ((resid[i]) < (resid[neighbor]) ? (resid[i]) : (resid[neighbor]));
        resid[neighbor] -= flow1;
        resid[i] -= flow1;
        flow[j] = flow1;
      }
    }
  }
  seen = smalloc((n_left + n_right) * sizeof(int));
  for (i = 0; i < n_left; i++) {
    modified = 1;
    while (resid[i] && modified) {
      modified = augment(i, pointers, indices, resid, flow, touched, seen);
    }
  }
  sfree(seen);
}
```