PolyFEM
Loading...
Searching...
No Matches
polyfem::assembler::MooneyRivlin3ParamSymbolic Class Reference

#include <MooneyRivlin3ParamSymbolic.hpp>

Inheritance diagram for polyfem::assembler::MooneyRivlin3ParamSymbolic:
[legend]
Collaboration diagram for polyfem::assembler::MooneyRivlin3ParamSymbolic:
[legend]

Public Member Functions

 MooneyRivlin3ParamSymbolic ()
 
const GenericMatParamc1 () const
 
const GenericMatParamc2 () const
 
const GenericMatParamc3 () const
 
const GenericMatParamd1 () const
 
double compute_energy (const NonLinearAssemblerData &data) const override
 
Eigen::VectorXd assemble_gradient (const NonLinearAssemblerData &data) const override
 
Eigen::MatrixXd assemble_hessian (const NonLinearAssemblerData &data) const override
 
void add_multimaterial (const int index, const json &params, const Units &units) override
 
std::string name () const override
 
std::map< std::string, ParamFuncparameters () const override
 
void assign_stress_tensor (const OutputData &data, const int all_size, const ElasticityTensorType &type, Eigen::MatrixXd &all, const std::function< Eigen::MatrixXd(const Eigen::MatrixXd &)> &fun) const override
 
void compute_stress_grad_multiply_mat (const OptAssemblerData &data, const Eigen::MatrixXd &mat, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const override
 
void compute_stress_grad_multiply_stress (const OptAssemblerData &data, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const override
 
void compute_stress_grad_multiply_vect (const OptAssemblerData &data, const Eigen::MatrixXd &vect, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const override
 
template<int n_basis, int dim>
void compute_energy_aux_gradient_fast (const NonLinearAssemblerData &data, Eigen::Matrix< double, Eigen::Dynamic, 1 > &G_flattened) const
 
double assemble_energy (const bool is_volume, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev) const override
 
void assemble_gradient (const bool is_volume, const int n_basis, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, Eigen::MatrixXd &rhs) const override
 
virtual Eigen::VectorXd assemble_gradient (const NonLinearAssemblerData &data) const=0
 
void assemble_hessian (const bool is_volume, const int n_basis, const bool project_to_psd, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, utils::MatrixCache &mat_cache, StiffnessMatrix &grad) const override
 
virtual Eigen::MatrixXd assemble_hessian (const NonLinearAssemblerData &data) const=0
 
- Public Member Functions inherited from polyfem::assembler::NLAssembler
virtual ~NLAssembler ()=default
 
double assemble_energy (const bool is_volume, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev) const override
 
Eigen::VectorXd assemble_energy_per_element (const bool is_volume, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev) const override
 
void assemble_gradient (const bool is_volume, const int n_basis, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, Eigen::MatrixXd &rhs) const override
 
void assemble_hessian (const bool is_volume, const int n_basis, const bool project_to_psd, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, utils::MatrixCache &mat_cache, StiffnessMatrix &grad) const override
 
virtual bool is_linear () const override
 
- Public Member Functions inherited from polyfem::assembler::Assembler
virtual ~Assembler ()=default
 
int size () const
 
virtual void set_size (const int size)
 
virtual void assemble (const bool is_volume, const int n_basis, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, StiffnessMatrix &stiffness, const bool is_mass=false) const
 
virtual void compute_stiffness_value (const double t, const assembler::ElementAssemblyValues &vals, const Eigen::MatrixXd &local_pts, const Eigen::MatrixXd &displacement, Eigen::MatrixXd &tensor) const
 
virtual void compute_dstress_dmu_dlambda (const OptAssemblerData &data, Eigen::MatrixXd &dstress_dmu, Eigen::MatrixXd &dstress_dlambda) const
 
virtual void compute_stress_grad (const OptAssemblerData &data, const Eigen::MatrixXd &prev_grad_u_i, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const
 
virtual void compute_stress_prev_grad (const OptAssemblerData &data, const Eigen::MatrixXd &prev_grad_u_i, Eigen::MatrixXd &result) const
 
virtual VectorNd compute_rhs (const AutodiffHessianPt &pt) const
 
virtual Eigen::Matrix< AutodiffScalarGrad, Eigen::Dynamic, 1, 0, 3, 1 > kernel (const int dim, const AutodiffGradPt &rvect, const AutodiffScalarGrad &r) const
 
void set_materials (const std::vector< int > &body_ids, const json &body_params, const Units &units)
 
virtual void update_lame_params (const Eigen::MatrixXd &lambdas, const Eigen::MatrixXd &mus)
 
virtual bool is_fluid () const
 
- Public Member Functions inherited from polyfem::assembler::ElasticityAssembler
 ElasticityAssembler ()
 
virtual ~ElasticityAssembler ()=default
 
void compute_scalar_value (const OutputData &data, std::vector< NamedMatrix > &result) const override
 
void compute_tensor_value (const OutputData &data, std::vector< NamedMatrix > &result) const override
 
void compute_stress_tensor (const OutputData &data, const ElasticityTensorType &type, Eigen::MatrixXd &stresses) const
 
void compute_von_mises_stresses (const OutputData &data, Eigen::MatrixXd &stresses) const
 
bool is_solution_displacement () const override
 
bool is_tensor () const override
 

Private Member Functions

template<typename T >
compute_energy_aux (const NonLinearAssemblerData &data) const
 
template<int n_basis, int dim>
void compute_energy_hessian_aux_fast (const NonLinearAssemblerData &data, Eigen::MatrixXd &H) const
 
template<int n_basis, int dim>
void compute_energy_aux_gradient_fast (const NonLinearAssemblerData &data, Eigen::VectorXd &G_flattened) const
 
template<typename T >
elastic_energy (const RowVectorNd &p, const int el_id, const AutoDiffGradMat< T > &def_grad) const
 
void get_grad_hess_symbolic (const double c1, const double c2, const double c3, const double d1, const Eigen::MatrixXd &def_grad, Eigen::MatrixXd &grad, Eigen::MatrixXd &hess) const
 
void get_grad_hess_autodiff (const double c1, const double c2, const double c3, const double d1, const Eigen::MatrixXd &global_pts, const int el_id, const Eigen::MatrixXd &F, Eigen::MatrixXd &grad, Eigen::MatrixXd &hess) const
 

Private Attributes

GenericMatParam c1_
 
GenericMatParam c2_
 
GenericMatParam c3_
 
GenericMatParam d1_
 

Additional Inherited Members

- Public Types inherited from polyfem::assembler::Assembler
typedef std::pair< std::string, Eigen::MatrixXd > NamedMatrix
 
typedef std::function< double(const RowVectorNd &, const RowVectorNd &, double, int)> ParamFunc
 
- Protected Attributes inherited from polyfem::assembler::Assembler
int size_ = -1
 

Detailed Description

Definition at line 16 of file MooneyRivlin3ParamSymbolic.hpp.

Constructor & Destructor Documentation

◆ MooneyRivlin3ParamSymbolic()

polyfem::assembler::MooneyRivlin3ParamSymbolic::MooneyRivlin3ParamSymbolic ( )

Definition at line 8 of file MooneyRivlin3ParamSymbolic.cpp.

Member Function Documentation

◆ add_multimaterial()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::add_multimaterial ( const int  index,
const json params,
const Units units 
)
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 559 of file MooneyRivlin3ParamSymbolic.cpp.

References polyfem::assembler::GenericMatParam::add_multimaterial(), c1_, c2_, c3_, d1_, and polyfem::Units::stress().

Here is the call graph for this function:

◆ assemble_energy()

double polyfem::assembler::NLAssembler::assemble_energy ( const bool  is_volume,
const std::vector< basis::ElementBases > &  bases,
const std::vector< basis::ElementBases > &  gbases,
const AssemblyValsCache cache,
const double  t,
const double  dt,
const Eigen::MatrixXd &  displacement,
const Eigen::MatrixXd &  displacement_prev 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 244 of file Assembler.cpp.

◆ assemble_gradient() [1/3]

void polyfem::assembler::NLAssembler::assemble_gradient ( const bool  is_volume,
const int  n_basis,
const std::vector< basis::ElementBases > &  bases,
const std::vector< basis::ElementBases > &  gbases,
const AssemblyValsCache cache,
const double  t,
const double  dt,
const Eigen::MatrixXd &  displacement,
const Eigen::MatrixXd &  displacement_prev,
Eigen::MatrixXd &  rhs 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

◆ assemble_gradient() [2/3]

Eigen::VectorXd polyfem::assembler::MooneyRivlin3ParamSymbolic::assemble_gradient ( const NonLinearAssemblerData data) const
overridevirtual

Implements polyfem::assembler::NLAssembler.

Definition at line 14 of file MooneyRivlin3ParamSymbolic.cpp.

References polyfem::assembler::ElementAssemblyValues::basis_values, polyfem::assembler::Assembler::size(), and polyfem::assembler::NonLinearAssemblerData::vals.

Here is the call graph for this function:

◆ assemble_gradient() [3/3]

virtual Eigen::VectorXd polyfem::assembler::NLAssembler::assemble_gradient ( const NonLinearAssemblerData data) const
virtual

◆ assemble_hessian() [1/3]

void polyfem::assembler::NLAssembler::assemble_hessian ( const bool  is_volume,
const int  n_basis,
const bool  project_to_psd,
const std::vector< basis::ElementBases > &  bases,
const std::vector< basis::ElementBases > &  gbases,
const AssemblyValsCache cache,
const double  t,
const double  dt,
const Eigen::MatrixXd &  displacement,
const Eigen::MatrixXd &  displacement_prev,
utils::MatrixCache mat_cache,
StiffnessMatrix grad 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

◆ assemble_hessian() [2/3]

Eigen::MatrixXd polyfem::assembler::MooneyRivlin3ParamSymbolic::assemble_hessian ( const NonLinearAssemblerData data) const
overridevirtual

Implements polyfem::assembler::NLAssembler.

Definition at line 84 of file MooneyRivlin3ParamSymbolic.cpp.

References polyfem::assembler::ElementAssemblyValues::basis_values, polyfem::assembler::Assembler::size(), and polyfem::assembler::NonLinearAssemblerData::vals.

Here is the call graph for this function:

◆ assemble_hessian() [3/3]

virtual Eigen::MatrixXd polyfem::assembler::NLAssembler::assemble_hessian ( const NonLinearAssemblerData data) const
virtual

◆ assign_stress_tensor()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::assign_stress_tensor ( const OutputData data,
const int  all_size,
const ElasticityTensorType type,
Eigen::MatrixXd &  all,
const std::function< Eigen::MatrixXd(const Eigen::MatrixXd &)> &  fun 
) const
overridevirtual

◆ c1()

const GenericMatParam & polyfem::assembler::MooneyRivlin3ParamSymbolic::c1 ( ) const
inline

Definition at line 25 of file MooneyRivlin3ParamSymbolic.hpp.

References c1_.

Referenced by assign_stress_tensor(), compute_energy_aux_gradient_fast(), compute_energy_hessian_aux_fast(), compute_stress_grad_multiply_mat(), compute_stress_grad_multiply_stress(), compute_stress_grad_multiply_vect(), elastic_energy(), get_grad_hess_symbolic(), and parameters().

Here is the caller graph for this function:

◆ c2()

const GenericMatParam & polyfem::assembler::MooneyRivlin3ParamSymbolic::c2 ( ) const
inline

Definition at line 26 of file MooneyRivlin3ParamSymbolic.hpp.

References c2_.

Referenced by assign_stress_tensor(), compute_energy_aux_gradient_fast(), compute_energy_hessian_aux_fast(), compute_stress_grad_multiply_mat(), compute_stress_grad_multiply_stress(), compute_stress_grad_multiply_vect(), elastic_energy(), get_grad_hess_symbolic(), and parameters().

Here is the caller graph for this function:

◆ c3()

const GenericMatParam & polyfem::assembler::MooneyRivlin3ParamSymbolic::c3 ( ) const
inline

Definition at line 27 of file MooneyRivlin3ParamSymbolic.hpp.

References c3_.

Referenced by assign_stress_tensor(), compute_energy_aux_gradient_fast(), compute_energy_hessian_aux_fast(), compute_stress_grad_multiply_mat(), compute_stress_grad_multiply_stress(), compute_stress_grad_multiply_vect(), elastic_energy(), get_grad_hess_symbolic(), and parameters().

Here is the caller graph for this function:

◆ compute_energy()

double polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_energy ( const NonLinearAssemblerData data) const
overridevirtual

Implements polyfem::assembler::NLAssembler.

Definition at line 372 of file MooneyRivlin3ParamSymbolic.cpp.

◆ compute_energy_aux()

template<typename T >
T polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_energy_aux ( const NonLinearAssemblerData data) const
private

◆ compute_energy_aux_gradient_fast() [1/2]

template<int n_basis, int dim>
void polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_energy_aux_gradient_fast ( const NonLinearAssemblerData data,
Eigen::Matrix< double, Eigen::Dynamic, 1 > &  G_flattened 
) const

◆ compute_energy_aux_gradient_fast() [2/2]

template<int n_basis, int dim>
void polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_energy_aux_gradient_fast ( const NonLinearAssemblerData data,
Eigen::VectorXd &  G_flattened 
) const
private

◆ compute_energy_hessian_aux_fast()

template<int n_basis, int dim>
void polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_energy_hessian_aux_fast ( const NonLinearAssemblerData data,
Eigen::MatrixXd &  H 
) const
private

◆ compute_stress_grad_multiply_mat()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_stress_grad_multiply_mat ( const OptAssemblerData data,
const Eigen::MatrixXd &  mat,
Eigen::MatrixXd &  stress,
Eigen::MatrixXd &  result 
) const
overridevirtual

◆ compute_stress_grad_multiply_stress()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_stress_grad_multiply_stress ( const OptAssemblerData data,
Eigen::MatrixXd &  stress,
Eigen::MatrixXd &  result 
) const
overridevirtual

◆ compute_stress_grad_multiply_vect()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::compute_stress_grad_multiply_vect ( const OptAssemblerData data,
const Eigen::MatrixXd &  vect,
Eigen::MatrixXd &  stress,
Eigen::MatrixXd &  result 
) const
overridevirtual

◆ d1()

const GenericMatParam & polyfem::assembler::MooneyRivlin3ParamSymbolic::d1 ( ) const
inline

Definition at line 28 of file MooneyRivlin3ParamSymbolic.hpp.

References d1_.

Referenced by assign_stress_tensor(), compute_energy_aux_gradient_fast(), compute_energy_hessian_aux_fast(), compute_stress_grad_multiply_mat(), compute_stress_grad_multiply_stress(), compute_stress_grad_multiply_vect(), elastic_energy(), get_grad_hess_symbolic(), and parameters().

Here is the caller graph for this function:

◆ elastic_energy()

template<typename T >
T polyfem::assembler::MooneyRivlin3ParamSymbolic::elastic_energy ( const RowVectorNd p,
const int  el_id,
const AutoDiffGradMat< T > &  def_grad 
) const
private

Definition at line 317 of file MooneyRivlin3ParamSymbolic.cpp.

References c1(), c1_, c2(), c2_, c3(), c3_, d1(), d1_, polyfem::utils::determinant(), polyfem::first_invariant(), polyfem::second_invariant(), polyfem::assembler::Assembler::size(), and val.

Referenced by compute_energy_aux(), and get_grad_hess_autodiff().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ get_grad_hess_autodiff()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::get_grad_hess_autodiff ( const double  c1,
const double  c2,
const double  c3,
const double  d1,
const Eigen::MatrixXd &  global_pts,
const int  el_id,
const Eigen::MatrixXd &  F,
Eigen::MatrixXd &  grad,
Eigen::MatrixXd &  hess 
) const
private

Definition at line 352 of file MooneyRivlin3ParamSymbolic.cpp.

References elastic_energy(), polyfem::F, DiffScalarBase::setVariableCount(), and polyfem::assembler::Assembler::size().

Here is the call graph for this function:

◆ get_grad_hess_symbolic()

void polyfem::assembler::MooneyRivlin3ParamSymbolic::get_grad_hess_symbolic ( const double  c1,
const double  c2,
const double  c3,
const double  d1,
const Eigen::MatrixXd &  def_grad,
Eigen::MatrixXd &  grad,
Eigen::MatrixXd &  hess 
) const
private

Definition at line 342 of file MooneyRivlin3ParamSymbolic.cpp.

References c1(), c2(), c3(), d1(), polyfem::F, polyfem::autogen::generate_gradient(), and polyfem::autogen::generate_hessian().

Referenced by compute_stress_grad_multiply_mat(), compute_stress_grad_multiply_stress(), and compute_stress_grad_multiply_vect().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ name()

std::string polyfem::assembler::MooneyRivlin3ParamSymbolic::name ( ) const
inlineoverridevirtual

Implements polyfem::assembler::Assembler.

Definition at line 38 of file MooneyRivlin3ParamSymbolic.hpp.

◆ parameters()

std::map< std::string, Assembler::ParamFunc > polyfem::assembler::MooneyRivlin3ParamSymbolic::parameters ( ) const
overridevirtual

Implements polyfem::assembler::Assembler.

Definition at line 567 of file MooneyRivlin3ParamSymbolic.cpp.

References c1(), c2(), c3(), and d1().

Here is the call graph for this function:

Member Data Documentation

◆ c1_

◆ c2_

◆ c3_

◆ d1_


The documentation for this class was generated from the following files: