Corrected a variety of errors relating to subdoains. Both ELL and ELL_MULTI formats
now will work with method_8a_multi (in the former case, there is only a single subdomain).
If I use method_8a_multi, I use the stencil values obtained from the rbffd coe. If using
method_8a() (only works with ELL and MM formats), I generate my own stencil values
(based on a variety of strategies not related to actual derivative stencils.)
Next step: Cuthill Mcgee.
Currently: I get same Gflop rate for one or 4 domains, as expected for large bandwidth
problems. Cuthill McGee should hopefully, alleviate this as the grids get smaller (that
is the hypothesis).