integrator.h

Go to the documentation of this file.

#ifndef WXM_INTEGRATOR_H
#define WXM_INTEGRATOR_H
 
#include "apps/integrate/integrate_app.h"
#include "basis/wmbasisarrayset.h"
#include "lib/geometry/unstructured_geometry/wmudggeometry.h"
#include "patchprocesses/wmpatchprocess.h"
 
#include "variable/distributed_variable.h"
#include "lib/wmunstructuredpatch.h"
 
#include "tools.h"
 
#include <memory>
 
namespace wxm
{
namespace dfem
{
class dg_integrator : public WmPatchProcess
{
protected:
    std::vector<std::unique_ptr<integrate_app>> apps;
    std::unique_ptr<WmUDGGeometry> dg_geometry;
    WmBasisArraySet basis_set;
    wxm::dfem::tools::scope_t element_scope;
 
    std::vector<size_t> in_idcs;
    std::vector<size_t> stage_idcs;
 
    std::vector<std::vector<real>> in_buf;
 
public:
    dg_integrator();
 
    void setup(const WxCryptSet& wxc) override;
 
    // unused, but override just to satisfy WmPatchProcess interface
    // this is the easiest way to get something which needs one object per local patch
    void process() override
    {
    }
 
    typedef std::vector<
        std::vector<std::pair<wxm::variable::distributed_variable_t*, size_t>>>
        variables_type;
 
    std::vector<std::vector<wxm::array::patch_array_t*>>
    get_patch_arrays(const variables_type& vars) const;
 
    virtual real integrate(const variables_type& vars);
 
    MPI_Op reduce_op;
};
} // namespace dfem
} // namespace wxm
#endif