complete_support_map Subroutine

public subroutine complete_support_map(comm, equi, mesh, mesh_staggered, vol, qsupp_fwd, qsupp_bwd, psupp_fwd, psupp_bwd)

Completes the operators for the support operator method - Turns map matrices qsupp into parallel gradient matrices - Establishes parallel divergence matrices psupp

Arguments

Type	Intent	Optional	Attributes		Name
integer,	intent(in)			::	comm	MPI communicator
class(equilibrium_t),	intent(inout)			::	equi	Equilibrium
type(mesh_cart_t),	intent(in)			::	mesh	Mesh
type(mesh_cart_t),	intent(in)			::	mesh_staggered	Mesh at toroidally staggered (fwd) position
real(kind=FP),	intent(in),		dimension(mesh%get_n_points())	::	vol	Flux box volumes
type(csrmat_t),	intent(inout)			::	qsupp_fwd	Map matrix on input, parallel gradient matrix on output (fwd direction)
type(csrmat_t),	intent(inout)			::	qsupp_bwd	Map matrix on input, parallel gradient matrix on output (bwd direction)
type(csrmat_t),	intent(out),		allocatable	::	psupp_fwd	Parallel divergence matrix (fwd direction)
type(csrmat_t),	intent(out),		allocatable	::	psupp_bwd	Parallel divergence matrix (bwd direction)