PeridigmForce.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#ifndef PERIDIGMFORCE_HPP
#define PERIDIGMFORCE_HPP

#include "Phalanx_ConfigDefs.hpp"
#include "Phalanx_Evaluator_WithBaseImpl.hpp"
#include "Phalanx_Evaluator_Derived.hpp"
#include "Phalanx_MDField.hpp"

#ifdef ALBANY_PERIDIGM
#include <Peridigm.hpp>
#include <Peridigm_AlbanyDiscretization.hpp>
#endif

namespace LCM {
template<typename EvalT, typename Traits>
class PeridigmForceBase : public PHX::EvaluatorWithBaseImpl<Traits>,
        public PHX::EvaluatorDerived<EvalT, Traits>  {

public:

  PeridigmForceBase(Teuchos::ParameterList& p,
        const Teuchos::RCP<Albany::Layouts>& dataLayout);

  void postRegistrationSetup(typename Traits::SetupData d,
                             PHX::FieldManager<Traits>& vm);

  void evaluateFields(typename Traits::EvalData d);

protected:

  Teuchos::RCP<Teuchos::ParameterList> peridigmParams;

  typedef typename EvalT::ScalarT ScalarT;
  typedef typename EvalT::MeshScalarT MeshScalarT;

  // Input:
  RealType density;
  PHX::MDField<ScalarT,Cell,QuadPoint,Dim> volume;
  PHX::MDField<ScalarT,Cell,QuadPoint,Dim> referenceCoordinates;
  PHX::MDField<ScalarT,Cell,QuadPoint,Dim> currentCoordinates;

  // Output:
  PHX::MDField<ScalarT,Cell,QuadPoint,Dim> force;
  PHX::MDField<ScalarT,Cell,QuadPoint,Dim> residual;

  unsigned int numQPs;
  unsigned int numDims;

#ifdef ALBANY_PERIDIGM
  // Peridigm objects
  Teuchos::RCP<PeridigmNS::Discretization> peridynamicDiscretization;
  Teuchos::RCP<PeridigmNS::Peridigm> peridigm;
#endif

};

// Inherted classes 
template<typename EvalT, typename Traits> class PeridigmForce;

// For all cases except those specialized below, just fall through to base class.
// The base class throws "Not Implemented" for evaluate fields.
template<typename EvalT, typename Traits>
class PeridigmForce : public PeridigmForceBase<EvalT, Traits> {
public:
  PeridigmForce(Teuchos::ParameterList& p, const Teuchos::RCP<Albany::Layouts>& dataLayout)
    : PeridigmForceBase<EvalT, Traits>(p, dataLayout) {};
};

// Template Specialization: Residual Eval calls Peridigm with doubles.
template<typename Traits>
class PeridigmForce<PHAL::AlbanyTraits::Residual, Traits> : public PeridigmForceBase<PHAL::AlbanyTraits::Residual, Traits> {
public:
  PeridigmForce(Teuchos::ParameterList& p, const Teuchos::RCP<Albany::Layouts>& dataLayout)
    : PeridigmForceBase<PHAL::AlbanyTraits::Residual,Traits>(p, dataLayout) {};
  void evaluateFields(typename Traits::EvalData d);
};

}

#endif