Grid_solvePoisson.F90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	Grid_solvePoisson (iSoln, iSrc, bcTypes, bcValues, poisfact)

Function/Subroutine Documentation

Grid_solvePoisson()

subroutine Grid_solvePoisson	(	integer	iSoln,
		integer	iSrc,
		integer	bcTypes,
		real	bcValues,
		real, intent(inout)	poisfact
	)

Driver routine for poisson solvers in the Grid

iSoln - the index for the solution variable (potential when used for self-gravity) iSrc - the index of the source variable (density when used for self-gravity) bcTypes - the boundary condition type; only the first entry is used. Only the first 2*NDIM elements are significant. They are interpreted in the order (X left, X right, Y left, Y right, Z left, Z right). Valid values are: GRID_PDE_BND_PERIODIC (1) GRID_PDE_BND_DIRICHLET (2) (homogeneous or constant Dirichlet) GRID_PDE_BND_NEUMANN (3) (homogeneous or constant Neumann) GRID_PDE_BND_ISOLATED (0) bcValues - the values to boundary conditions, currently not used (treated as 0) poisfact - scaling factor to be used in calculation

NOTES

The symbols listed above for bcTypes are declared as FORTRAN PARAMETERS in the module Grid_interfaces. Code using this interface should refer to that module with a USE statement, like this:

use Grid_interface, ONLY : GRID_PDE_BND_PERIODIC, GRID_PDE_BND_NEUMANN, & GRID_PDE_BND_ISOLATED, GRID_PDE_BND_DIRICHLET, & Grid_solvePoisson

Most implementations only support a limited subset of boundary condition types. Some implementations require that all significant elements of bcTypes are the same. (That is always the case for GRID_PDE_BND_ISOLATED.) Even if an implementation supports combinations of different boundary conditions on different sides of the domain, the types at left and right sides for the same axis direction will usually have to be the samme.

Support in some implementations provided with FLASH4:

GridSolvers/Multipole: GRID_PDE_BND_ISOLATED GridSolvers/Multipole_new: GRID_PDE_BND_ISOLATED GridSolvers/Multigrid (simple): GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET(hom.), GRID_PDE_BND_NEUMANN(hom.)(?) (same type in all directions), or GRID_PDE_BND_ISOLATED (requires Paramesh as Grid with NBlockX==NBlockY==NBlockZ==1) GridSolvers/MultigridMC (simple): GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET(hom.), GRID_PDE_BND_NEUMANN(hom.) (same type in all directions); GRID_PDE_BND_ISOLATED does not work properly. (requires Paramesh as Grid with NBlockX==NBlockY==NBlockZ==1) GridSolvers/Pfft: GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET(hom.), GRID_PDE_BND_NEUMANN(hom.) in various combinations, depending on GridSolvers/Pfft subdirectory (i.e. implementation configured in) (requires UG in pencil shape or Paramesh as Grid) GridSolvers/Multigrid hybrid with Pfft: GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET, GRID_PDE_BND_NEUMANN in various combinations, or GRID_PDE_BND_ISOLATED (requires Paramesh as Grid) GridSolvers/MultigridMC hybrid with Pfft: GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET, GRID_PDE_BND_NEUMANN in various combinations; GRID_PDE_BND_ISOLATED does not work properly. (requires Paramesh as Grid) GridSolvers/BiPCGStab preconditioned with MC multigrid: GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET, GRID_PDE_BND_NEUMANN in various combinations (requires Paramesh as Grid) GridSolvers/BiPCGStab preconditioned with (HG) multigrid: GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET, GRID_PDE_BND_NEUMANN in various combinations (requires Paramesh as Grid) GridSolvers/BiPCGStab preconditioned with HYPRE: GRID_PDE_BND_PERIODIC, GRID_PDE_BND_DIRICHLET, GRID_PDE_BND_NEUMANN in various combinations (requires Paramesh as Grid ?)

Definition at line 94 of file Grid_solvePoisson.F90.