Recreating fields change variable ordering ?
Consider the following skeleton:
SubdomainField< std::complex< double >, Form::Form0 > u("u", omega | gammaOut, FunctionSpaceTypeForm0::HierarchicalH1, order);
InterfaceField< std::complex< double >, Form::Form0 > g("g", sigma, FunctionSpaceTypeForm0::HierarchicalH1, order);
while(true)
{
gmshddm::problem::Formulation< std::complex< double > > formulation("Helmholtz", topology);
// Tell to the formulation that g is field that have to be exchanged between subdomains
formulation.addInterfaceField(g);
// Add terms, assemble, export matrix
}
If the field objects are recreated inside the loop instead of being reused, the matrices are different at each iteration. Why ? Is the ordering non deterministic ? At each run, the results are the same, but for each iteration of the loop the result is different (but appears deterministic)