ConstraintsEliminationOperator.h

Go to the documentation of this file.

 
#ifndef included_AMP_ConstraintsEliminationOperator
#define included_AMP_ConstraintsEliminationOperator
 
#include "AMP/operators/Operator.h"
#include "AMP/vectors/Variable.h"
#include "AMP/vectors/Vector.h"
#include <vector>
 
namespace AMP::Operator {
 
class ConstraintsEliminationOperator : public Operator
{
 
public:
    explicit ConstraintsEliminationOperator( std::shared_ptr<const OperatorParameters> params )
        : Operator( params )
    {
        AMP_INSIST( params->d_db->keyExists( "InputVariable" ), "key not found" );
        std::string inpVarName = params->d_db->getString( "InputVariable" );
        d_InputVariable.reset( new AMP::LinearAlgebra::Variable( inpVarName ) );
 
        AMP_INSIST( params->d_db->keyExists( "OutputVariable" ), "key not found" );
        std::string outVarName = params->d_db->getString( "OutputVariable" );
        d_OutputVariable.reset( new AMP::LinearAlgebra::Variable( outVarName ) );
    }
 
    virtual ~ConstraintsEliminationOperator() {}
 
    std::string type() const override { return "ConstraintsEliminationOperator"; }
 
    void apply( AMP::LinearAlgebra::Vector::const_shared_ptr u,
                AMP::LinearAlgebra::Vector::shared_ptr f ) override;
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> getInputVariable() const override;
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> getOutputVariable() const override;
 
    virtual void addSlaveToMaster( AMP::LinearAlgebra::Vector::shared_ptr u ) = 0;
    virtual void copyMasterToSlave( AMP::LinearAlgebra::Vector::shared_ptr u ) = 0;
    virtual void setSlaveToZero( AMP::LinearAlgebra::Vector::shared_ptr u );
    virtual void addShiftToSlave( AMP::LinearAlgebra::Vector::shared_ptr u );
 
protected:
    std::vector<size_t> d_SlaveIndices;
    std::vector<double> d_SlaveShift;
 
private:
    std::shared_ptr<AMP::LinearAlgebra::Variable> d_InputVariable;  
    std::shared_ptr<AMP::LinearAlgebra::Variable> d_OutputVariable; 
};
} // namespace AMP::Operator
 
#endif