setSolutionIntoInterfaceFields

Hi again,

I would like to look at the eigenvectors of the \mathcal{F}=(\mathcal{I}-\mathcal{A}) operator of the Schwarz procedure. To this end, I proceed with the following:

1. get the discrete operator (associated with \mathcal{F}) via gmshddm::problem::Formulation::computeMatrix(),
2. compute its eigenvalues and eigenvectors with gmshfem::algebra::eigenvalues<>(),
3. insert each eigenvector (i.e. a g_{ij} vector) in the gmshddm::problem::Formulation and solve the volume PDEs.

Would it be a better/cleaner way?

If not, then a method to implement the last step is required, if I'm not mistaken. To this end, I suggest simply a clean wrapper around gmshddm::problem::Formulation::_fillG(). The naming gmshddm::problem::Formulation::setSolutionIntoInterfaceFields() is an echo of gmshfem::problem::Formulation::setSolutionIntoFields() in gmshfem.

Cheers, Nicolas.