residual_box Subroutine

private subroutine residual_box(mg, id, nc)

Arguments

Type IntentOptional Attributes Name
type(mg_t), intent(inout) :: mg
integer, intent(in) :: id
integer, intent(in) :: nc

Called by

proc~~residual_box~2~~CalledByGraph proc~residual_box~2 residual_box proc~max_residual_lvl~2 max_residual_lvl proc~max_residual_lvl~2->proc~residual_box~2 proc~update_coarse~2 update_coarse proc~update_coarse~2->proc~residual_box~2 proc~mg_fas_fmg~2 mg_fas_fmg proc~mg_fas_fmg~2->proc~update_coarse~2 proc~mg_fas_vcycle~2 mg_fas_vcycle proc~mg_fas_fmg~2->proc~mg_fas_vcycle~2 proc~mg_fas_vcycle~2->proc~max_residual_lvl~2 proc~mg_fas_vcycle~2->proc~update_coarse~2 proc~diffusion_solve_acoeff~2 diffusion_solve_acoeff proc~diffusion_solve_acoeff~2->proc~mg_fas_fmg~2 proc~diffusion_solve_acoeff~2->proc~mg_fas_vcycle~2 proc~diffusion_solve_vcoeff~2 diffusion_solve_vcoeff proc~diffusion_solve_vcoeff~2->proc~mg_fas_fmg~2 proc~diffusion_solve_vcoeff~2->proc~mg_fas_vcycle~2 proc~diffusion_solve~2 diffusion_solve proc~diffusion_solve~2->proc~mg_fas_fmg~2 proc~diffusion_solve~2->proc~mg_fas_vcycle~2

Source Code

  subroutine residual_box(mg, id, nc)
    type(mg_t), intent(inout) :: mg
    integer, intent(in)       :: id
    integer, intent(in)       :: nc

    call mg%box_op(mg, id, nc, mg_ires)

    mg%boxes(id)%cc(1:nc, mg_ires) = mg%boxes(id)%cc(1:nc,&
        mg_irhs) - mg%boxes(id)%cc(1:nc, mg_ires)
  end subroutine residual_box