General initialization of the nested mesh hierarchy. More...

Functions
static int	dom_decomp (const int Nx, const int Ny, const int Nz, const int Np, int pNGx, int pNGy, int *pNGz)
	calls auto domain decomposition functions
static int	dom_decomp_2d (const int Nx, const int Ny, const int Np, int pNGx, int pNGy)
	finds optimum domain decomposition in 2D
static int	dom_decomp_3d (const int Nx, const int Ny, const int Nz, const int Np, int pNGx, int pNGy, int *pNGz)
	finds optimum domain decomposition in 3D
void	init_mesh (MeshS *pM)
	General initialization of the nested mesh hierarchy.
void	get_myGridIndex (DomainS pD, const int myID, int pi, int pj, int pk)
	Searches GData[][][] array to find i,j,k components of block being updated on this processor.

Detailed Description

General initialization of the nested mesh hierarchy.

PURPOSE: General initialization of the nested mesh hierarchy. Works for both nested and uniform meshes, on single and multiple processors. Each Mesh contains one or more levels, each level contains one or more Domains (regions with the same grid resolution), and depending on the parallelization each Domain contains one or more Grids (however there can never be more than one Grid per Domain being updated on any given processor). In the Mesh, this hierarchy is stored as an "array" of Domains indexed as Domain[nlevel][ndomain]. Since ndomain is different for each level, this "array" is not square, (really it is nlevel pointers, each to ndomain[nlevel] Domains).

Note for a uniform mesh on a single processor:

# of Mesh levels = # of Domains = # of Grids = 1 For a uniform mesh on multiple processors:
# of Mesh levels = # of Domains = 1; # of Grids = # of processors For a nested mesh on a single processor:
# of Domains = # of Grids

For a nested mesh on multiple processors, there is no relationship between these quantaties in general.

This function:

(1) sets properties of each Domain read from <domain> blocks in input file
(2) allocates and initializes the array of Domains,
(3) divides each Domain into one or more Grids depending on the parallelization.

This function supercedes init_domain() from v3.2. The init_grid() function initializes the data in each Grid structure in each Domain, including finding all child and parent Grids with SMR.

CONTAINS PUBLIC FUNCTIONS:

PRIVATE FUNCTION PROTOTYPES:

dom_decomp() - calls auto domain decomposition functions
dom_decomp_2d() - finds optimum domain decomposition in 2D
dom_decomp_3d() - finds optimum domain decomposition in 3D

Definition in file init_mesh.c.

Function Documentation

static int dom_decomp	(	const int	Nx,
		const int	Ny,
		const int	Nz,
		const int	Np,
		int *	pNGx,
		int *	pNGy,
		int *	pNGz
	)			`[static]`

calls auto domain decomposition functions

Calls apropriate 2D or 3D auto decomposition routines Functions written by T.A.G., added May 2007.

Definition at line 886 of file init_mesh.c.

References dom_decomp_2d(), and dom_decomp_3d().

Referenced by init_mesh().

Here is the call graph for this function:

Here is the caller graph for this function:

static int dom_decomp_2d	(	const int	Nx,
		const int	Ny,
		const int	Np,
		int *	pNGx,
		int *	pNGy
	)			`[static]`

finds optimum domain decomposition in 2D

Pptimizes domain decomposition in 2D.

The TOTAL amount of data communicated (summed over all processes and all INTERNAL boundaries) divided by 2*nghost (where the 2 is for two messages per internal interface) is computed and stored in the variable I, the minimum of which is I0. This assumes that in the x-direction we communicate only computational cells, while in the y-direction we comunicate the computational cells and x-direction ghost cells. Then Total x-communication is (rx - 1)*Ny*(2*nghost) Total y-communication is (ry - 1)*(Nx + rx*(2*nghost))*(2*nghost)

Definition at line 923 of file init_mesh.c.

References rx, and ry.

Referenced by dom_decomp(), and dom_decomp_3d().

Here is the caller graph for this function: