user: 下面的汇编指令对应的C源代码是什么：

```asm
<main>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x128,%rsp
  mov    %edi,-0xc4(%rbp)
  mov    %rsi,-0xd0(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
  mov    %rax,%r13
  movl   $0x1000,-0x94(%rbp)
  movl   $0x1000,-0x90(%rbp)
  movl   $0x3e8,-0x8c(%rbp)
  movss  D93(%rip),%xmm0
  movss  %xmm0,-0x88(%rbp)
  movss  D25(%rip),%xmm0
  movss  %xmm0,-0x84(%rbp)
  movss  D77(%rip),%xmm0
  movss  %xmm0,-0xbc(%rbp)
  mov    -0x90(%rbp),%esi
  mov    -0x94(%rbp),%edi
  movslq %esi,%rax
  sub    $0x1,%rax
  mov    %rax,-0x80(%rbp)
  movslq %esi,%rax
  mov    %rax,%rcx
  mov    $0x0,%ebx
  movslq %esi,%rax
  lea    0x0(,%rax,4),%r12
  movslq %edi,%rax
  sub    $0x1,%rax
  mov    %rax,-0x78(%rbp)
  movslq %esi,%rax
  mov    %rax,%r10
  mov    $0x0,%r11d
  movslq %edi,%rax
  mov    %rax,%r8
  mov    $0x0,%r9d
  mov    %r11,%rdx
  imul   %r8,%rdx
  mov    %r9,%rax
  imul   %r10,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r10,%rax
  mul    %r8
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rax
  mov    %rax,-0xe0(%rbp)
  movq   $0x0,-0xd8(%rbp)
  movslq %edi,%rax
  mov    %rax,%r14
  mov    $0x0,%r15d
  mov    -0xe0(%rbp),%r9
  mov    -0xd8(%rbp),%r10
  mov    %r10,%rdx
  imul   %r14,%rdx
  mov    %r9,%rax
  imul   %r15,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r9,%rax
  mul    %r14
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rdx
  movslq %edi,%rax
  imul   %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%ebx
  mov    $0x0,%edx
  div    %rbx
  imul   $0x10,%rax,%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L227:
  cmp    %rcx,%rsp
  je     L161
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L227
L161:
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  mov    %rdx,%rax
  and    $0xfff,%eax
  test   %rax,%rax
  je     L244
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L244:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x70(%rbp)
  mov    -0x90(%rbp),%esi
  mov    -0x94(%rbp),%edi
  movslq %esi,%rax
  sub    $0x1,%rax
  mov    %rax,-0x68(%rbp)
  movslq %esi,%rax
  mov    %rax,-0x130(%rbp)
  movq   $0x0,-0x128(%rbp)
  movslq %esi,%rax
  lea    0x0(,%rax,4),%rbx
  movslq %edi,%rax
  sub    $0x1,%rax
  mov    %rax,-0x60(%rbp)
  movslq %esi,%rax
  mov    %rax,-0xf0(%rbp)
  movq   $0x0,-0xe8(%rbp)
  movslq %edi,%rax
  mov    %rax,-0x100(%rbp)
  movq   $0x0,-0xf8(%rbp)
  mov    -0xf0(%rbp),%r8
  mov    -0xe8(%rbp),%r9
  mov    %r9,%rdx
  mov    -0x100(%rbp),%r10
  mov    -0xf8(%rbp),%r11
  imul   %r10,%rdx
  mov    %r11,%rax
  imul   %r8,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r8,%rax
  mul    %r10
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rax
  mov    %rax,-0x110(%rbp)
  movq   $0x0,-0x108(%rbp)
  movslq %edi,%rax
  mov    %rax,-0x120(%rbp)
  movq   $0x0,-0x118(%rbp)
  mov    -0x110(%rbp),%r8
  mov    -0x108(%rbp),%r9
  mov    %r9,%rdx
  mov    -0x120(%rbp),%r14
  mov    -0x118(%rbp),%r15
  imul   %r14,%rdx
  mov    %r15,%rax
  imul   %r8,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r8,%rax
  mul    %r14
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rdx
  movslq %edi,%rax
  imul   %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L131:
  cmp    %rdx,%rsp
  je     L3
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L131
L3:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L50
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L50:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x58(%rbp)
  mov    -0x94(%rbp),%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x50(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x140(%rbp)
  movq   $0x0,-0x138(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0x150(%rbp)
  movq   $0x0,-0x148(%rbp)
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L260:
  cmp    %rdx,%rsp
  je     L112
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L260
L112:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L48
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L48:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x48(%rbp)
  mov    -0x94(%rbp),%eax
  imul   -0x90(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edx
  mov    -0x70(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  movl   $0x0,-0xb8(%rbp)
  jmp    L119
L211:
  mov    -0x70(%rbp),%rax
  mov    -0xb8(%rbp),%edx
  movslq %edx,%rdx
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax,%rdx,4)
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0x94(%rbp),%eax
  lea    -0x1(%rax),%esi
  mov    -0x70(%rbp),%rax
  mov    -0xb8(%rbp),%edx
  movslq %edx,%rdi
  movslq %esi,%rdx
  imul   %rcx,%rdx
  add    %rdi,%rdx
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax,%rdx,4)
  addl   $0x1,-0xb8(%rbp)
L119:
  mov    -0xb8(%rbp),%eax
  cmp    -0x90(%rbp),%eax
  jl     L211
  movl   $0x0,-0xb4(%rbp)
  jmp    L24
L280:
  pxor   %xmm0,%xmm0
  cvtsi2ssl -0xb4(%rbp),%xmm0
  mulss  -0x88(%rbp),%xmm0
  mov    -0x94(%rbp),%eax
  sub    $0x1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2ss %eax,%xmm1
  divss  %xmm1,%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <sinf@plt>
  movd   %xmm0,%eax
  mov    -0x48(%rbp),%rdx
  mov    -0xb4(%rbp),%ecx
  movslq %ecx,%rcx
  mov    %eax,(%rdx,%rcx,4)
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0x48(%rbp),%rax
  mov    -0xb4(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm0
  mov    -0x70(%rbp),%rax
  mov    -0xb4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  mov    -0x48(%rbp),%rax
  mov    -0xb4(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm2
  movss  %xmm2,-0xe0(%rbp)
  movss  -0x88(%rbp),%xmm0
  movss  D88(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <expf@plt>
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0x90(%rbp),%eax
  lea    -0x1(%rax),%edx
  mulss  -0xe0(%rbp),%xmm0
  mov    -0x70(%rbp),%rax
  movslq %edx,%rsi
  mov    -0xb4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  addl   $0x1,-0xb4(%rbp)
L24:
  mov    -0xb4(%rbp),%eax
  cmp    -0x94(%rbp),%eax
  jl     L280
  pxor   %xmm3,%xmm3
  cvtsi2sdl -0x94(%rbp),%xmm3
  movq   %xmm3,%rax
  mov    -0x90(%rbp),%edx
  mov    %edx,%esi
  movq   %rax,%xmm0
  lea    D89(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  mov    $0x0,%eax
  call   <StartTimer@plt>
  movl   $0x0,-0xb0(%rbp)
  movl   $0x1,-0xac(%rbp)
  jmp    L79
L169:
  mov    -0x58(%rbp),%rax
  mov    -0xac(%rbp),%edx
  movslq %edx,%rdx
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax,%rdx,4)
  mov    %rbx,%rcx
  shr    $0x2,%rcx
  mov    -0x94(%rbp),%eax
  lea    -0x1(%rax),%esi
  mov    -0x58(%rbp),%rax
  mov    -0xac(%rbp),%edx
  movslq %edx,%rdi
  movslq %esi,%rdx
  imul   %rcx,%rdx
  add    %rdi,%rdx
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax,%rdx,4)
  addl   $0x1,-0xac(%rbp)
L79:
  mov    -0xac(%rbp),%eax
  cmp    -0x90(%rbp),%eax
  jl     L169
  movl   $0x1,-0xa8(%rbp)
  jmp    L29
L257:
  mov    %rbx,%rcx
  shr    $0x2,%rcx
  mov    -0x48(%rbp),%rax
  mov    -0xa8(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm0
  mov    -0x58(%rbp),%rax
  mov    -0xa8(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  mov    -0x48(%rbp),%rax
  mov    -0xa8(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm4
  movss  %xmm4,-0xe0(%rbp)
  movss  -0x88(%rbp),%xmm0
  movss  D88(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <expf@plt>
  mov    %rbx,%rcx
  shr    $0x2,%rcx
  mov    -0x90(%rbp),%eax
  lea    -0x1(%rax),%edx
  mulss  -0xe0(%rbp),%xmm0
  mov    -0x58(%rbp),%rax
  movslq %edx,%rsi
  mov    -0xa8(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  addl   $0x1,-0xa8(%rbp)
L29:
  mov    -0xa8(%rbp),%eax
  cmp    -0x94(%rbp),%eax
  jl     L257
  jmp    L28
L30:
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0xbc(%rbp)
  movl   $0x1,-0xa4(%rbp)
  jmp    L168
L243:
  movl   $0x1,-0xa0(%rbp)
  jmp    L123
L262:
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0xa0(%rbp),%eax
  lea    0x1(%rax),%edx
  mov    -0x70(%rbp),%rax
  movslq %edx,%rsi
  mov    -0xa4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  (%rax,%rdx,4),%xmm1
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0xa0(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x70(%rbp),%rax
  movslq %edx,%rsi
  mov    -0xa4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  (%rax,%rdx,4),%xmm0
  addss  %xmm0,%xmm1
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0xa4(%rbp),%eax
  lea    -0x1(%rax),%esi
  mov    -0x70(%rbp),%rax
  mov    -0xa0(%rbp),%edx
  movslq %edx,%rdi
  movslq %esi,%rdx
  imul   %rcx,%rdx
  add    %rdi,%rdx
  movss  (%rax,%rdx,4),%xmm0
  addss  %xmm0,%xmm1
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0xa4(%rbp),%eax
  lea    0x1(%rax),%esi
  mov    -0x70(%rbp),%rax
  mov    -0xa0(%rbp),%edx
  movslq %edx,%rdi
  movslq %esi,%rdx
  imul   %rcx,%rdx
  add    %rdi,%rdx
  movss  (%rax,%rdx,4),%xmm0
  addss  %xmm0,%xmm1
  mov    %rbx,%rcx
  shr    $0x2,%rcx
  movss  D47(%rip),%xmm0
  mulss  %xmm1,%xmm0
  mov    -0x58(%rbp),%rax
  mov    -0xa0(%rbp),%edx
  movslq %edx,%rsi
  mov    -0xa4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  mov    %rbx,%rcx
  shr    $0x2,%rcx
  mov    -0x58(%rbp),%rax
  mov    -0xa0(%rbp),%edx
  movslq %edx,%rsi
  mov    -0xa4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  (%rax,%rdx,4),%xmm0
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0x70(%rbp),%rax
  mov    -0xa0(%rbp),%edx
  movslq %edx,%rsi
  mov    -0xa4(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  (%rax,%rdx,4),%xmm1
  subss  %xmm1,%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <fabsf@plt>
  mov    %eax,%edx
  mov    -0xbc(%rbp),%eax
  mov    %edx,%edi
  movd   %eax,%xmm0
  call   <fmaxf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0xbc(%rbp)
  addl   $0x1,-0xa0(%rbp)
L123:
  mov    -0x90(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0xa0(%rbp)
  jl     L262
  addl   $0x1,-0xa4(%rbp)
L168:
  mov    -0x94(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0xa4(%rbp)
  jl     L243
  movl   $0x1,-0x9c(%rbp)
  jmp    L271
L190:
  movl   $0x1,-0x98(%rbp)
  jmp    L115
L90:
  mov    %rbx,%rsi
  shr    $0x2,%rsi
  mov    %r12,%rcx
  shr    $0x2,%rcx
  mov    -0x58(%rbp),%rax
  mov    -0x98(%rbp),%edx
  movslq %edx,%rdi
  mov    -0x9c(%rbp),%edx
  movslq %edx,%rdx
  imul   %rsi,%rdx
  add    %rdi,%rdx
  movss  (%rax,%rdx,4),%xmm0
  mov    -0x70(%rbp),%rax
  mov    -0x98(%rbp),%edx
  movslq %edx,%rsi
  mov    -0x9c(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  addl   $0x1,-0x98(%rbp)
L115:
  mov    -0x90(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0x98(%rbp)
  jl     L90
  addl   $0x1,-0x9c(%rbp)
L271:
  mov    -0x94(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0x9c(%rbp)
  jl     L190
  mov    -0xb0(%rbp),%edx
  movslq %edx,%rax
  imul   $0x51eb851f,%rax,%rax
  shr    $0x20,%rax
  sar    $0x5,%eax
  mov    %edx,%ecx
  sar    $0x1f,%ecx
  sub    %ecx,%eax
  imul   $0x64,%eax,%ecx
  mov    %edx,%eax
  sub    %ecx,%eax
  test   %eax,%eax
  jne    L240
  pxor   %xmm0,%xmm0
  cvtss2sd -0xbc(%rbp),%xmm0
  pxor   %xmm5,%xmm5
  cvtsi2sdl -0xb0(%rbp),%xmm5
  movq   %xmm5,%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  lea    D87(%rip),%rax
  mov    %rax,%rdi
  mov    $0x2,%eax
  call   <printf@plt>
L240:
  addl   $0x1,-0xb0(%rbp)
L28:
  movss  -0xbc(%rbp),%xmm0
  comiss -0x84(%rbp),%xmm0
  jbe    L47
  mov    -0xb0(%rbp),%eax
  cmp    -0x8c(%rbp),%eax
  jl     L30
L47:
  mov    $0x0,%eax
  call   <GetTimer@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  movsd  D64(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  lea    D49(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  mov    %r13,%rsp
  mov    $0x0,%eax
  mov    -0x38(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L61
  call   <__stack_chk_fail@plt>
L61:
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

tool_calls: name: retrieve_data
arguments: data_tag: D93
datatype: byte[4]
name: retrieve_data
arguments: data_tag: D25
datatype: byte[4]
name: retrieve_data
arguments: data_tag: D77
datatype: byte[4]
name: retrieve_data
arguments: data_tag: D88
datatype: byte[4]
name: retrieve_data
arguments: data_tag: D89
datatype: string
name: retrieve_data
arguments: data_tag: D47
datatype: byte[4]
name: retrieve_data
arguments: data_tag: D87
datatype: string
name: retrieve_data
arguments: data_tag: D64
datatype: byte[8]
name: retrieve_data
arguments: data_tag: D49
datatype: string

tool: {"D93": [0xdb,0x0f,0x49,0x40]}

tool: {"D25": [0xac,0xc5,0x27,0x37]}

tool: {"D77": [0x00,0x00,0x80,0x3f]}

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

tool: D89: Jacobi relaxation Calculation: %d x %d mesh


tool: {"D47": [0x00,0x00,0x80,0x3e]}

tool: D87: %5d, %0.6f


tool: {"D64": [0x00,0x00,0x00,0x00,0x00,0x40,0x8f,0x40]}

tool: D49:  total: %f s


assistant: ```c
int main(int argc, char **argv) {
  int n = 4096;
  int m = 4096;
  int iter_max = 1000;
  const float pi = 2.0f * asinf(1.0f);
  const float tol = 1.0e-5f;
  float error = 1.0f;
  float A[n][m];
  float Anew[n][m];
  float y0[n];
  memset(A, 0, n * m * sizeof(float));
  for (int i = 0; i < m; i++) {
    A[0][i] = 0.f;
    A[n - 1][i] = 0.f;
  }
  for (int j = 0; j < n; j++) {
    y0[j] = sinf(pi * j / (n - 1));
    A[j][0] = y0[j];
    A[j][m - 1] = y0[j] * expf(-pi);
  }
  printf("Jacobi relaxation Calculation: %d x %d mesh\n", n, m);
  StartTimer();
  int iter = 0;
#pragma omp parallel for shared(Anew)
  for (int i = 1; i < m; i++) {
    Anew[0][i] = 0.f;
    Anew[n - 1][i] = 0.f;
  }
#pragma omp parallel for shared(Anew)
  for (int j = 1; j < n; j++) {
    Anew[j][0] = y0[j];
    Anew[j][m - 1] = y0[j] * expf(-pi);
  }
  while (error > tol && iter < iter_max) {
    error = 0.f;
#pragma omp parallel for shared(m, n, Anew, A)
#pragma acc kernels
    for (int j = 1; j < n - 1; j++) {
      for (int i = 1; i < m - 1; i++) {
        Anew[j][i] = 0.25f * (A[j][i + 1] + A[j][i - 1] + A[j - 1][i] + A[j + 1][i]);
        error = fmaxf(error, fabsf(Anew[j][i] - A[j][i]));
      }
    }
#pragma omp parallel for shared(m, n, Anew, A)
#pragma acc kernels
    for (int j = 1; j < n - 1; j++) {
      for (int i = 1; i < m - 1; i++) {
        A[j][i] = Anew[j][i];
      }
    }
    if (iter % 100 == 0)
      printf("%5d, %0.6f\n", iter, error);
    iter++;
  }
  double runtime = GetTimer();
  printf(" total: %f s\n", runtime / 1000.f);
}
```