sigma_physics.h

Go to the documentation of this file.

#include "libmesh/fem_system.h"
#include "libmesh/mesh.h"
#include "libmesh/mesh_function.h"
#include "libmesh/fem_physics.h"
#include "libmesh/system.h"
 
// Using namespace libMesh
using namespace libMesh;
 
// FEMSystem, TimeSolver and  NewtonSolver will handle most tasks,
// but we must specify element residuals
class SigmaPhysics : public FEMPhysics
{
public:
  // Constructor
  SigmaPhysics() : FEMPhysics() {}
 
  // Special functions
  SigmaPhysics (SigmaPhysics &&) = default;
  SigmaPhysics (const SigmaPhysics &) = default;
  SigmaPhysics & operator= (const SigmaPhysics &) = default;
  SigmaPhysics & operator= (SigmaPhysics &&) = default;
  virtual ~SigmaPhysics () = default;
 
  Real & k() { return _k; }
  bool & analytic_jacobians() { return _analytic_jacobians; }
 
  // System initialization
  virtual void init_data (System & sys);
 
  virtual std::unique_ptr<DifferentiablePhysics> clone_physics() override;
 
protected:
 
  // Context initialization
  virtual void init_context (DiffContext & context) override;
 
  // Element residual and jacobian calculations
  // Time dependent parts
  virtual bool element_time_derivative (bool request_jacobian,
                                        DiffContext & context) override;
 
  // The parameters to solve for
  Real _k;
 
  // Variables to hold the computed QoIs
  Number computed_QoI[1];
 
  // The FE type to use
  std::string _fe_family;
  unsigned int _fe_order;
 
  // Index for T variable
  unsigned int T_var;
 
  // Calculate Jacobians analytically or not?
  bool _analytic_jacobians;
};