Context Navigation

CylindricalGravity

v33	v34
101	101	This routine is responsible for making changes to the stencil coefficients. Originally it was written to fix coefficients for boundary/ghost zones. This makes it a natural place to modify coefficients for differing geometry as well; this is where the fixes for 2.5D gravity will go.
102	102
103		Since ~~stencil values is a module~~ wide variable, we should not modify it itself in any routine. Instead, we define a new array, tempStencilValues, set it = !StencilValues, and then modify elements of this new array when necessary.
	103	Since "stencil values" is a module-wide variable, we should not modify it itself in any routine. Instead, we define a new array, tempStencilValues, set it = !StencilValues, and then modify elements of this new array when necessary.
104	104
105	105	The routine begins by allocating an array to hold all of the coefficients for all of the stencil's cells for each of the cells in the patch. This array is aptly called 'matrixvalues' and is mxmymz*num_stencil_values long (where mx,y,z are the number of cells in each direction of the patch, and num_stencil_values is self explanatory).