VolumeIntegralOperator.h

Go to the documentation of this file.

 
#ifndef included_AMP_VolumeIntegralOperator
#define included_AMP_VolumeIntegralOperator
 
#include "AMP/matrices/Matrix.h"
#include "AMP/operators/libmesh/NonlinearFEOperator.h"
#include "AMP/operators/libmesh/SourceNonlinearElement.h"
#include "AMP/operators/libmesh/VolumeIntegralOperatorParameters.h"
#include "AMP/vectors/MultiVariable.h"
 
#include <vector>
 
namespace AMP::Operator {
 
 
class VolumeIntegralOperator : public NonlinearFEOperator
{
public:
    explicit VolumeIntegralOperator( std::shared_ptr<const OperatorParameters> params );
 
    virtual ~VolumeIntegralOperator() {}
 
    std::string type() const override { return "VolumeIntegralOperator"; }
 
    void reset( std::shared_ptr<const OperatorParameters> ) override;
 
    void setAuxVariable( std::shared_ptr<AMP::LinearAlgebra::MultiVariable> var )
    {
        d_auxVariables = var;
    }
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> getInputVariable() const override
    {
        return d_inpVariables;
    }
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> getOutputVariable() const override
    {
        return d_outVariable;
    }
 
    std::shared_ptr<SourcePhysicsModel> getSourcePhysicsModel() { return d_sourcePhysicsModel; }
 
    std::shared_ptr<SourceNonlinearElement> getSourceElement() { return d_srcNonlinElem; }
 
    /*
       std::shared_ptr<AMP::LinearAlgebra::Matrix> getLinearizedVolumeIntegralOperator(
       std::shared_ptr<OperatorParameters> params);
     */
 
protected:
    explicit VolumeIntegralOperator( std::shared_ptr<const VolumeIntegralOperatorParameters> params,
                                     bool );
 
    std::shared_ptr<OperatorParameters>
    getJacobianParameters( AMP::LinearAlgebra::Vector::const_shared_ptr u ) override;
 
 
    void preAssembly( AMP::LinearAlgebra::Vector::const_shared_ptr u,
                      AMP::LinearAlgebra::Vector::shared_ptr r ) override;
 
    void postAssembly() override;
 
    void preElementOperation( const AMP::Mesh::MeshElement & ) override;
 
    void postElementOperation() override;
 
    void getNodeDofIndicesForCurrentElement();
 
    void init( std::shared_ptr<const VolumeIntegralOperatorParameters> params );
 
    std::string d_isInputType;
 
    std::vector<double> d_elementOutputVector; 
    std::shared_ptr<SourceNonlinearElement> d_srcNonlinElem; 
    std::shared_ptr<SourcePhysicsModel> d_sourcePhysicsModel; 
    std::vector<AMP::LinearAlgebra::Vector::const_shared_ptr>
        d_inVec; 
    std::vector<AMP::LinearAlgebra::Vector::const_shared_ptr>
        d_auxVec; 
    AMP::LinearAlgebra::Vector::shared_ptr d_multiAuxPtr;
 
    AMP::LinearAlgebra::Vector::shared_ptr d_outVec; 
    std::vector<size_t> d_dofIndices; 
    std::shared_ptr<AMP::Discretization::DOFManager> d_elementDofMap;
    std::shared_ptr<AMP::Discretization::DOFManager> d_nodeDofMap;
 
    AMP::Mesh::MeshElementVectorPtr d_currNodes;
 
private:
    std::shared_ptr<AMP::LinearAlgebra::MultiVariable> d_inpVariables; // Input Active variable
    std::shared_ptr<AMP::LinearAlgebra::MultiVariable> d_auxVariables; // Input Auxillary variable
    std::shared_ptr<AMP::LinearAlgebra::Variable> d_outVariable;       // Output variable
};
 
 
} // namespace AMP::Operator
 
#endif