def_arnoldi Subroutine

public subroutine def_arnoldi(V, H, k, f2c, CurrentGrid)

Arnoldi prosess of DEF_fullgmres Attention: h_(i,j)=(w,v_i)=v^H*w, for complex value, so the code should be dot_product(v_i,w)

Arguments

Type IntentOptional Attributes Name
complex(kind=realdp), intent(inout), dimension(1-LAP:f2c%nxc+LAP,1-LAP:f2c%nyc+LAP,def_mg_miter+1) :: V
complex(kind=realdp), intent(inout), dimension(def_mg_miter+1,def_mg_miter) :: H
integer, intent(in) :: k
type(TwoGrids), intent(inout) :: f2c
type(Gridpara) :: CurrentGrid

Calls

proc~~def_arnoldi~~CallsGraph proc~def_arnoldi def_arnoldi proc~ex Ex proc~def_arnoldi->proc~ex proc~mg_dot_prod mg_dot_prod proc~def_arnoldi->proc~mg_dot_prod proc~mg_norm mg_norm proc~def_arnoldi->proc~mg_norm proc~mgcslp_invmhx MGCSLP_invMHx proc~def_arnoldi->proc~mgcslp_invmhx proc~ex->proc~mgcslp_invmhx proc~coarsegridpara CoarseGridpara proc~ex->proc~coarsegridpara proc~helm_a2hx Helm_A2hx proc~ex->proc~helm_a2hx proc~helm_ahx Helm_Ahx proc~ex->proc~helm_ahx proc~ztx ZTx proc~ex->proc~ztx proc~ztzx ZTZx proc~ex->proc~ztzx proc~zx Zx proc~ex->proc~zx mpi_allreduce mpi_allreduce proc~mg_dot_prod->mpi_allreduce proc~mg_norm->mpi_allreduce proc~f_cycle F_cycle proc~mgcslp_invmhx->proc~f_cycle proc~grid_destroy grid_destroy proc~mgcslp_invmhx->proc~grid_destroy proc~gridsys_define gridsys_define proc~mgcslp_invmhx->proc~gridsys_define proc~v_cycle V_cycle proc~mgcslp_invmhx->proc~v_cycle proc~f_cycle->proc~f_cycle proc~f_cycle->proc~grid_destroy proc~f_cycle->proc~v_cycle proc~coarsegrid_create coarsegrid_create proc~f_cycle->proc~coarsegrid_create proc~cslp_op_bc CSLP_OP_BC proc~f_cycle->proc~cslp_op_bc proc~damp_jacobi_smoother Damp_Jacobi_smoother proc~f_cycle->proc~damp_jacobi_smoother proc~mg_bicgstab mg_bicgstab proc~f_cycle->proc~mg_bicgstab proc~mg_fullgmres mg_fullgmres proc~f_cycle->proc~mg_fullgmres proc~prolongation_en_correct prolongation_en_correct proc~f_cycle->proc~prolongation_en_correct proc~restriction restriction proc~f_cycle->proc~restriction proc~helmholtz2d_bc_mg Helmholtz2d_BC_mg proc~helm_a2hx->proc~helmholtz2d_bc_mg proc~helmholtz2d_o4cmpct Helmholtz2d_O4cmpct proc~helm_a2hx->proc~helmholtz2d_o4cmpct proc~helmholtz2d_red_glk Helmholtz2d_ReD_Glk proc~helm_a2hx->proc~helmholtz2d_red_glk proc~helm_ahx->proc~helmholtz2d_bc_mg proc~v_cycle->proc~grid_destroy proc~v_cycle->proc~v_cycle proc~v_cycle->proc~coarsegrid_create proc~v_cycle->proc~cslp_op_bc proc~v_cycle->proc~damp_jacobi_smoother proc~v_cycle->proc~mg_bicgstab proc~v_cycle->proc~mg_fullgmres proc~v_cycle->proc~prolongation_en_correct proc~v_cycle->proc~restriction proc~mg_check_xy2d mg_check_xy2d proc~ztx->proc~mg_check_xy2d proc~ztzx->proc~ztx proc~ztzx->proc~zx proc~zx->proc~mg_check_xy2d proc~cslp_op_bc->proc~mg_check_xy2d proc~cslp2d_stencils cslp2d_stencils proc~cslp_op_bc->proc~cslp2d_stencils proc~if_bcs if_BCs proc~cslp_op_bc->proc~if_bcs proc~damp_jacobi_smoother->proc~mg_check_xy2d proc~damp_jacobi_smoother->proc~cslp2d_stencils proc~damp_jacobi_smoother->proc~if_bcs proc~helmholtz2d_bc_mg->proc~mg_check_xy2d proc~helmholtz2d_stencils Helmholtz2d_stencils proc~helmholtz2d_bc_mg->proc~helmholtz2d_stencils proc~helmholtz2d_bc_mg->proc~if_bcs proc~helmholtz2d_o4cmpct->proc~mg_check_xy2d proc~helmholtz2d_red_glk->proc~mg_check_xy2d proc~extrpltghostbcs ExtrpltGhostBCs proc~helmholtz2d_red_glk->proc~extrpltghostbcs proc~helmholtz2d_red_glk->proc~helmholtz2d_stencils proc~red_glk_stencils ReD_Glk_stencils proc~helmholtz2d_red_glk->proc~red_glk_stencils proc~mg_bicgstab->proc~mg_dot_prod proc~mg_bicgstab->proc~mg_norm proc~mg_bicgstab->proc~cslp_op_bc proc~mg_check_x2d mg_check_x2d proc~mg_check_xy2d->proc~mg_check_x2d proc~mg_check_y2d mg_check_y2d proc~mg_check_xy2d->proc~mg_check_y2d proc~mg_fullgmres->proc~mg_norm proc~mg_fullgmres->proc~cslp_op_bc proc~mg_apply_givens_rotation mg_apply_givens_rotation proc~mg_fullgmres->proc~mg_apply_givens_rotation proc~mg_arnoldi mg_arnoldi proc~mg_fullgmres->proc~mg_arnoldi proc~mg_back_substitute mg_back_substitute proc~mg_fullgmres->proc~mg_back_substitute proc~prolongation_en_correct->proc~mg_check_xy2d proc~restriction->proc~mg_check_xy2d proc~mg_checkreal_xy2d mg_checkreal_xy2d proc~restriction->proc~mg_checkreal_xy2d proc~dirichlet Dirichlet proc~if_bcs->proc~dirichlet proc~first_order_neum first_order_Neum proc~if_bcs->proc~first_order_neum cdsqrt cdsqrt proc~mg_apply_givens_rotation->cdsqrt proc~mg_arnoldi->proc~mg_dot_prod proc~mg_arnoldi->proc~mg_norm proc~mg_arnoldi->proc~cslp_op_bc mpi_sendrecv mpi_sendrecv proc~mg_check_x2d->mpi_sendrecv proc~mg_check_y2d->mpi_sendrecv proc~mg_checkreal_x2d mg_checkreal_x2d proc~mg_checkreal_xy2d->proc~mg_checkreal_x2d proc~mg_checkreal_y2d mg_checkreal_y2d proc~mg_checkreal_xy2d->proc~mg_checkreal_y2d proc~mg_checkreal_x2d->mpi_sendrecv proc~mg_checkreal_y2d->mpi_sendrecv

Called by

proc~~def_arnoldi~~CalledByGraph proc~def_arnoldi def_arnoldi proc~def_fullgmres DEF_fullgmres proc~def_fullgmres->proc~def_arnoldi proc~invey invEy proc~invey->proc~def_fullgmres proc~qx Qx proc~qx->proc~invey proc~p_adef1x P_ADEF1x proc~p_adef1x->proc~qx proc~p_defx P_DEFx proc~p_defx->proc~qx proc~p_tlkmx P_TLKMx proc~p_tlkmx->proc~qx proc~px Px proc~px->proc~qx proc~def_px DEF_Px proc~def_px->proc~p_adef1x proc~def_px->proc~p_defx proc~def_px->proc~p_tlkmx proc~precond_x Precond_x proc~precond_x->proc~def_px proc~full_pgcr full_pgcr proc~full_pgcr->proc~precond_x proc~full_pgmres full_pgmres proc~full_pgmres->proc~precond_x proc~pfgmres pfgmres proc~pfgmres->proc~precond_x proc~pre_fullgmres Pre_fullgmres proc~pre_fullgmres->proc~precond_x proc~prearnoldi Prearnoldi proc~pre_fullgmres->proc~prearnoldi proc~pre_restartgmres Pre_restartgmres proc~pre_restartgmres->proc~precond_x proc~pre_restartgmres->proc~prearnoldi proc~prearnoldi->proc~precond_x program~helmholtz_2d helmholtz_2d program~helmholtz_2d->proc~full_pgcr program~helmholtz_2d->proc~full_pgmres program~helmholtz_2d->proc~pfgmres program~helmholtz_2d->proc~pre_fullgmres program~helmholtz_2d->proc~pre_restartgmres