TangentOutputModule.h

Go to the documentation of this file.

 
#pragma once
 
#include <jem/base/Error.h>
#include <jem/numeric/Sparse.h>
#include <jem/numeric/sparse/select.h>
#include <jive/algebra/AbstractMatrix.h>
#include <jive/algebra/ConstrainedMatrix.h>
#include <jive/algebra/SparseMatrixObject.h>
#include <jive/app/Module.h>
#include <jive/app/ModuleFactory.h>
#include <jive/fem/FEMatrixBuilder.h>
#include <jive/implict/Names.h>
#include <jive/implict/NonlinModule.h>
#include <jive/implict/SolverModule.h>
#include <jive/implict/utilities.h>
#include <jive/model/ModelFactory.h>
#include <jive/solver/Solver.h>
#include <jive/solver/declare.h>
#include <jive/util/FuncUtils.h>
 
#include "models/DirichletModel.h"
#include "models/PeriodicBCModel.h"
#include "modules/PBCGroupOutputModule.h"
 
using jem::Ref;
using jem::numeric::Function;
using jive::Properties;
using jive::SparseMatrix;
using jive::String;
using jive::algebra::AbstractMatrix;
using jive::algebra::ConstrainedMatrix;
using jive::app::Module;
using jive::app::ModuleFactory;
using jive::fem::FEMatrixBuilder;
using jive::implict::SolverModule;
using jive::model::Model;
using jive::util::FuncUtils;
 
class TangentOutputModule : public Module
{
public:
  JEM_DECLARE_CLASS(TangentOutputModule, Module);
 
  static const char *TYPE_NAME; 
 
  explicit TangentOutputModule(const String &name = "tangentOutput");
 
  virtual Status init(const Properties &conf, const Properties &props,
                      const Properties &globdat) override;
 
  virtual Status run(const Properties &globdat) override;
 
  virtual void shutdown(const Properties &globdat) override;
 
  static Ref<Module> makeNew(const String &name, const Properties &conf,
                             const Properties &props, const Properties &globdat);
 
  static void declare();
 
private:
 
  void readStrainStress_(const Vector &strains, const Vector &stresses,
                         const Properties &globdat);
 
  void reportStrainStress_(const Vector &H, const Vector &N);
 
  void getStrainStress_(const Matrix &strains, const Matrix &stresses,
                        const Properties &globdat);
 
  void storeTangentProps_(const Matrix &strains, const Matrix &stresses,
                          const Properties &globdat);
 
  void condenseMatrix_(const Matrix &strains, const Matrix &stresses,
                       const Properties &globdat);
 
protected:
  virtual ~TangentOutputModule();
 
private:
  String mode_;      
  idx_t rank_;       
  double thickness_; 
  double perturb_;   
 
  Ref<Model> masterModel_;             
  Ref<Function> sampleCond_;           
  Ref<GroupOutputModule> groupUpdate_; 
  Ref<SolverModule> solver_;           
  Ref<Constraints> cons_;              
 
  StringVector strains_;  
  StringVector stresses_; 
  StringVector sizes_;    
  IdxMatrix strainDofs_;  
};