[[Imported from SVN: r6819]]

[[Imported from SVN: r6818]]

[[Imported from SVN: r6817]]

[[Imported from SVN: r6816]]

[[Imported from SVN: r6815]]

[[Imported from SVN: r6814]]

[[Imported from SVN: r6813]]

[[Imported from SVN: r6810]]

[[Imported from SVN: r6809]]

[[Imported from SVN: r6808]]

[[Imported from SVN: r6807]]

[[Imported from SVN: r6806]]

[[Imported from SVN: r6803]]

bugfix: use operator[] instead of find() if you want to force creation of a non-existant entry in a std::map

[[Imported from SVN: r6802]]

[[Imported from SVN: r6801]]

[[Imported from SVN: r6800]]

[[Imported from SVN: r6799]]

[[Imported from SVN: r6798]]

don't store the continuum interface boundary in the Steklov-Poincare step -- use the one from the complex instead

[[Imported from SVN: r6797]]

[[Imported from SVN: r6796]]

[[Imported from SVN: r6795]]

[[Imported from SVN: r6794]]

ake the interfaceBoundary a LeafBoundaryPatch, instead of a vector of LevelBoundaryPatches, and remove the hacks associated with the latter

[[Imported from SVN: r6793]]

[[Imported from SVN: r6792]]

[[Imported from SVN: r6779]]

After all, we are thinking 'Steklov-Poincare' here, which means
we are solving an equation on the interface.  The old code
taking into account the subdomain solutions wasn't really satisfying:
in particular, it was not symmetric wrt the subdomains.

[[Imported from SVN: r6762]]

[[Imported from SVN: r6761]]

[[Imported from SVN: r6760]]

[[Imported from SVN: r6759]]

[[Imported from SVN: r6758]]

[[Imported from SVN: r6757]]

[[Imported from SVN: r6756]]

[[Imported from SVN: r6755]]

[[Imported from SVN: r6754]]

[[Imported from SVN: r6753]]

[[Imported from SVN: r6752]]

[[Imported from SVN: r6751]]

[[Imported from SVN: r6750]]

[[Imported from SVN: r6749]]

[[Imported from SVN: r6748]]