Wrote equivalent to OMEGA/FULLWV 9-point stencil. Still needs implementation of boundary conditions.
Attempt at writing free-surface boundary conditions. Not convinced this works.
Correct visual key for right-handed coordinate system in the new code. Also, fix some naming conventions.
Set up proper coefficients for the fixed boundary condition in the case that the free surface is not set.
Modify initialization of the grid terms in two main ways:
1. Adopt the convention used by S. Roecker, in which the edge entries use only the central buoyancy.
2. Reorder the terms in the coefficient calculation to make sure that they're in the same order as the rest of the equivalent lines of terms.
Fixed model parameter (didn't include inverse of density term). Also cleared up some of the notation, but added comments to make it clear how things relate back to the original code.
Changed name to avoid confusion.
Do not return a solver object; instead leave that for the calling program in case a different solver is appropriate.
Whole series of updates to the testing environment.
Several changes towards implementing the PML. In particular, the buoyancies are set up in a different way and the divisors are separated out.
Additionally, the padded arrays are now named as such to avoid accidentally using the wrong version.
First attempt at PML; pretty sure this doesn't work yet.
Fixed typo.
Semi-working stencil + PML.