Simple interface - field without calculation

[kulikova.adriana]

Hello,

I would like to have a field "pFieldPrev" which only copied values of "pField" and is not calculated otherwise. If I add it in the following way, the KSP does not converge (I add the order, etc., and everything compiles):

CHKERR simpleInterface->addDomainField(pFieldPrev, H1, AINSWORTH_LEGENDRE_BASE, 1);

I assume it is trying to calculate the pFieldPrev field since it is included in the following:

CHKERR DMoFEMLoopFiniteElements(dM, simpleInterface->getDomainFEName(), feDomainLhsPtr);

Is there a specific way to do this in simple interface or should I create it on the side?

I would still like to be able to print the field and use it for calculations in the following iteration.

Thanks,

Adriana

[Karol Lewandowski]

I guess you want to add data field, not domain field
use this function from the simple 

  MoFEMErrorCode addDataField(const std::string &name, const FieldSpace space,
                              const FieldApproximationBase base,
                              const FieldCoefficientsNumber nb_of_coefficients,
                              const TagType tag_type = MB_TAG_SPARSE,
                              const enum MoFEMTypes bh = MF_ZERO,
                              int verb = -1);

​

[Karol Lewandowski]

The difference is that the domain field will add dofs to the problem and without filling them, yes KSP can struggle to converge because you added a bunch of zeroes to your problem. 

[kulikova.adriana]

Yes, this works nicely.

Thank you