Domain decomposition works with two domains. See routine GEgenerate() written
by G. Erlebacher.

Next: output each domain to a file in Ellpack format. For now, not concerned with linking
them. The objective is to solve n derivative problems simultaneously and find out the
conditions under which one gets improved performance: more calculations gainst better
cache properties.