doxygen/libmesh/biharmonic_8h_source.html

#ifndef BIHARMONIC_H

#define BIHARMONIC_H


#include "libmesh/equation_systems.h"

#include "libmesh/replicated_mesh.h"

#include "libmesh/exodusII_io.h"

#include "libmesh/mesh_refinement.h"


// Bring in bits from the libMesh namespace.

// Just the bits we're using, since this is a header.

using libMesh::EquationSystems;

using libMesh::ExodusII_IO;

using libMesh::Point;

using libMesh::Real;

using libMesh::ReplicatedMesh;


#ifdef LIBMESH_ENABLE_AMR

using libMesh::MeshRefinement;

#endif


class Biharmonic : public EquationSystems

{

public:

  // Initial state enumeration

  enum InitialStateEnum {STRIP = 0,

                         ROD   = 1,

                         BALL  = 2};


  // Free energy enumeration

  enum FreeEnergyEnum {DOUBLE_WELL         = 1,

                       DOUBLE_OBSTACLE     = 2,

                       LOG_DOUBLE_WELL     = 3,

                       LOG_DOUBLE_OBSTACLE = 4};


  Biharmonic(ReplicatedMesh & mesh);


  ~Biharmonic()

  {

    // delete _meshRefinement;

  }


  // Misc. getters

  bool verbose()         { return _verbose; }

  Real dt0()             { return _dt0; }

  Real dt()              { return _dt; }


  // Public interface functions

  void viewParameters();

  void init();

  void step(const Real & dt = -1.0);

  void output(int timestep, const Real & t, Real & o_t, bool force = false);

  void run();


private:

  unsigned int  _dim, _N;

  Real _kappa, _theta, _theta_c;

  Real _tol;

  bool _growth, _degenerate, _cahn_hillard, _netforce;

  FreeEnergyEnum  _energy;

  int _log_truncation;

  bool _verbose;

  InitialStateEnum  _initialState;

  Point _initialCenter;

  Real _initialWidth;

  Real _dt0, _dt, _t0, _T, _t1;

  Real _cnWeight;

  //

  std::string  _ofile_base, _ofile;

  std::unique_ptr<ExodusII_IO> _exio;

  Real    _o_dt;

  int     _o_count;

  //

  friend class JR;

  class JR;       // forward

  ReplicatedMesh & _mesh;

  JR * _jr;

};


#endif // BIHARMONIC_H