TSElasticPreAndPostStep.hpp

Go to the documentation of this file.

/** @file TSElasticPreAndPostStep.hpp
 * @brief Elastic TS pre/post-step registration and H1 transfer lifecycle.
 */
 
#ifndef ESHELBIAN_TS_ELASTIC_PRE_AND_POST_STEP_HPP
#define ESHELBIAN_TS_ELASTIC_PRE_AND_POST_STEP_HPP
 
namespace EshelbianPlasticity {
 
struct TSElasticPreAndPostStep {
 
  /** Register coefficient reset before solves, with optional H1 bookkeeping. */
  static MoFEMErrorCode setPreStep(EshelbianCore *ep_ptr, TS ts,
                                   bool h1_update = false);
 
  /** Own projection resources for one solve; also clean up on early exits. */
  struct ProjectionScope {
    explicit ProjectionScope(EshelbianCore *ep_ptr);
    ~ProjectionScope();
    ProjectionScope(const ProjectionScope &) = delete;
    ProjectionScope &operator=(const ProjectionScope &) = delete;
    MoFEMErrorCode close();
 
  private:
    MPI_Comm comm;
    bool active = false;
  };
 
  static MoFEMErrorCode preStepFun(TS ts);
  static MoFEMErrorCode postStepFun(TS ts);
 
private:
  static MoFEMErrorCode postStepInitialise(EshelbianCore *ep_ptr);
  static MoFEMErrorCode postStepDestroy();
 
  static inline SmartPetscObj<KSP> prjKsp;
  static inline SmartPetscObj<Vec> prjD;
  static inline SmartPetscObj<Vec> prjF;
  static inline SmartPetscObj<DM> dmPrjSpatial;
  static inline EshelbianCore *epPtr;
};
} // namespace EshelbianPlasticity
 
#endif // ESHELBIAN_TS_ELASTIC_PRE_AND_POST_STEP_HPP