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

```asm
<sweep>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %r9,%r13
  push   %r12
  mov    %esi,%r12d
  push   %rbp
  mov    %edi,%ebp
  push   %rbx
  sub    $0x48,%rsp
  mov    0x88(%rsp),%ebx
  mov    %rdx,0x38(%rsp)
  mov    %r8d,0x2c(%rsp)
  mov    0x80(%rsp),%r14
  movsd  %xmm0,0x18(%rsp)
  movsd  %xmm1,0x30(%rsp)
  test   %ebx,%ebx
  je     L94
  mov    %edi,%eax
  cltd
  idiv   %ebx
  mov    %eax,0x24(%rsp)
L27:
  mov    %r12d,%eax
  cltd
  idiv   %ebx
  mov    %eax,0xc(%rsp)
  lea    0x1(%rcx),%eax
  mov    %eax,0x28(%rsp)
  cmp    0x2c(%rsp),%eax
  jg     L77
  mov    %ebx,0x88(%rsp)
  nopw   0x0(%rax,%rax,1)
L29:
  mov    0x24(%rsp),%eax
  xor    %edx,%edx
  mov    %r14,%rbx
  mov    %edx,%r15d
  test   %eax,%eax
  jle    L20
  nopl   0x0(%rax)
L91:
  mov    0xc(%rsp),%edi
  xor    %r14d,%r14d
  test   %edi,%edi
  jle    L85
  mov    %rbx,%rax
  mov    %r14d,%ebx
  mov    %rax,%r14
  nopl   0x0(%rax)
L104:
  sub    $0x8,%rsp
  mov    %r14,%r9
  mov    %r13,%r8
  mov    %ebx,%ecx
  mov    0x90(%rsp),%eax
  mov    %r15d,%edx
  mov    %r12d,%esi
  mov    %ebp,%edi
  add    $0x1,%ebx
  push   %rax
  call   <copy_block@plt>
  pop    %r8
  pop    %r9
  cmp    %ebx,0xc(%rsp)
  jne    L104
  mov    %r14,%rbx
L85:
  add    $0x1,%r15d
  cmp    0x24(%rsp),%r15d
  jne    L91
  mov    %r12d,0x20(%rsp)
  mov    %rbx,%r14
  movsd  0x30(%rsp),%xmm1
  xor    %r15d,%r15d
  mov    0x38(%rsp),%rbx
  nopw   0x0(%rax,%rax,1)
L63:
  mov    0xc(%rsp),%edx
  xor    %r12d,%r12d
  test   %edx,%edx
  jle    L65
  mov    %r12d,%eax
  mov    %r15d,%r12d
  mov    %eax,%r15d
  nopl   0x0(%rax)
L87:
  mov    0x88(%rsp),%eax
  mov    %r15d,%esi
  mov    %r14,%rcx
  mov    %ebp,%r9d
  mov    %rbx,%r8
  mov    %r13,%rdx
  mov    %r12d,%edi
  add    $0x1,%r15d
  push   %rax
  mov    0x28(%rsp),%eax
  push   %rax
  movsd  0x28(%rsp),%xmm0
  movsd  %xmm1,0x20(%rsp)
  call   <compute_estimate@plt>
  pop    %rcx
  pop    %rsi
  cmp    %r15d,0xc(%rsp)
  movsd  0x10(%rsp),%xmm1
  jne    L87
  mov    %r12d,%r15d
L65:
  add    $0x1,%r15d
  cmp    0x24(%rsp),%r15d
  jne    L63
  mov    0x20(%rsp),%r12d
L20:
  addl   $0x1,0x28(%rsp)
  mov    0x28(%rsp),%eax
  cmp    %eax,0x2c(%rsp)
  jge    L29
L77:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L94:
  movl   $0x1,0x24(%rsp)
  mov    %edi,%ebx
  jmp    L27
```

assistant: ```c
void sweep(int nx, int ny, double dx, double dy, double *f_, int itold, int itnew, double *u_, double *unew_, int block_size) {
  int it;
  int block_x, block_y;
  if (block_size == 0)
    block_size = nx;
  int max_blocks_x = (nx / block_size);
  int max_blocks_y = (ny / block_size);
#pragma omp parallel shared(u_, unew_, f_, max_blocks_x, max_blocks_y, nx, ny, dx, dy, itold, itnew, block_size) private(it, block_x, block_y)
  {
#pragma omp single
    {
      for (it = itold + 1; it <= itnew; it++) {
        for (block_x = 0; block_x < max_blocks_x; block_x++) {
          for (block_y = 0; block_y < max_blocks_y; block_y++) {
#pragma omp task shared(u_, unew_, nx, ny, block_size) firstprivate(block_x, block_y)
            copy_block(nx, ny, block_x, block_y, u_, unew_, block_size);
          }
        }
#pragma omp taskwait
        for (block_x = 0; block_x < max_blocks_x; block_x++) {
          for (block_y = 0; block_y < max_blocks_y; block_y++) {
#pragma omp task default(none) shared(u_, unew_, f_, dx, dy, nx, ny, block_size) firstprivate(block_x, block_y)
            compute_estimate(block_x, block_y, u_, unew_, f_, dx, dy, nx, ny, block_size);
          }
        }
#pragma omp taskwait
      }
    }
  }
}
```