Context Navigation

CylindricalGravity

v31	v32
59	59	=== Poisson Init routine ===
60	60
61		The first chunk of code in poisson.f90, '!PoissonInit', goes through and checks that the boundary conditions are set up correctly for the elliptic solver, and then sets up the 'Hypre' arrays (those arrays which are sent off to Hypre for solving Poisson's equation). In what follows, we'll '''focus on the ndim=2 case for illustration''' of the code.
	61	The first chunk of code in poisson.f90, '!PoissonInit', goes through and checks that the boundary conditions are set up correctly for the elliptic solver, and then sets up the 'Hypre' arrays (those arrays which are sent off to Hypre for solving Poisson's equation). In what follows, we'll '''focus on the ndim=2 case for illustration''' of the code, since cylindrical fixes amount to modifying the 2D stencil coefficients.
62	62
63	63	There are 2 main arrays in this section of the code: 1. 'offsets', and 2. 'stencil values'. (ientries seems less important).