**NAME**

Grid_solvePoisson

**SYNOPSIS**

call Grid_solvePoisson(integer(IN) :: iSoln, integer(IN) :: iSrc, integer(IN) :: bcTypes(6), real(IN) :: bcValues(2,6), real(INOUT) :: poisfact)

**DESCRIPTION**

Driver routine for poisson solvers in the Grid

**ARGUMENTS**

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 ?)

High-Energy Density Physics

Thermonuclear-Powered Supernovae

Fluid-Structure Interaction

Implicit Solvers

Code Requests

Download

Contributed code

Flash Test downloads

Visualization downloads

Documentation and User Support

Publications

FLASH internal

The Flash Center for Computational Science is based at the University of Chicago and is supported by U.S. DOE and NSF.