SubchannelToCladGPMap.h

Go to the documentation of this file.

#ifndef included_AMP_SubchannelToCladGPMap
#define included_AMP_SubchannelToCladGPMap
 
#include "AMP/discretization/createLibmeshElements.h"
#include "AMP/operators/map/SubchannelToCladMap.h"
 
namespace AMP::Operator {
 
 
typedef SubchannelToCladMapParameters SubchannelToCladGPMapParameters;
 
 
class SubchannelToCladGPMap : public SubchannelToCladMap
{
public:
    std::string type() const override { return "SubchannelToCladGPMap"; }
 
    static bool validMapType( const std::string &s );
 
    explicit SubchannelToCladGPMap(
        std::shared_ptr<const AMP::Operator::OperatorParameters> params );
 
    virtual ~SubchannelToCladGPMap();
 
 
protected:
    // For a given subchannel, fill the elements of interest using the coordinates
    virtual void fillReturnVector( AMP::LinearAlgebra::Vector::shared_ptr vec,
                                   double range[4],
                                   std::shared_ptr<AMP::Mesh::Mesh> mesh,
                                   const std::vector<AMP::Mesh::MeshElementID> &ids,
                                   const std::vector<double> &z,
                                   const std::vector<double> &f ) override;
 
private:
    Discretization::createLibmeshElements libmeshElements;
 
    struct gaussPointZCoord {
        double z[4];
    };
    std::vector<gaussPointZCoord>
    getGaussPoints( std::shared_ptr<AMP::Mesh::Mesh> mesh,
                    const std::vector<AMP::Mesh::MeshElementID> &ids );
};
 
 
} // namespace AMP::Operator
 
#endif