spatial_solver.h

Go to the documentation of this file.

#ifndef WMSPATIALSOLVER_H
#define WMSPATIALSOLVER_H
 
// Wm includes
#include "dfem/dfem.h"
#include "hostactions/wmtemporalsolver.h"
#include "patchprocesses/wmpatchprocess.h"
 
#include <memory>
 
// forward declarations
class WmMemoryObject;
 
namespace wxm
{
namespace dfem
{
namespace solver
{
 
class spatial_solver_t : public WmPatchProcess
{
protected:
    std::vector<size_t> rhs_idcs_;
    std::vector<size_t> all_idcs_;
 
public:
    bool test_dt = true;
 
    std::vector<std::unique_ptr<WmApplication>> _apps;
 
    typedef typename wxm::temporal_solver::variables_type variables_type;
 
    spatial_solver_t() = default;
 
    ~spatial_solver_t() override = default;
 
    void setup(const WxCryptSet& wxc) override;
 
    virtual const std::shared_ptr<TimestepConstraint>
    solve(const real time, const variables_type& input, variables_type& output) = 0;
 
    virtual void Barrier(const WxMsgBase& msg, const real time, variables_type& output)
    {
    }
 
    std::vector<wxm::array::patch_array_t*>
    get_patch_arrays(const variables_type& in) const;
 
    virtual const WxRange& element_periphery_range() const = 0;
 
    virtual const WxRange& element_interior_range() const = 0;
 
    virtual const WxRange& element_full_range() const = 0;
 
    // TODO: this is only here because zero_fluxes isn't designed to
    // take advantage of scopes yet. Need to fix in ndg_t
    virtual void zero_fluxes() = 0;
 
    const std::vector<size_t>& get_output_indices() const
    {
        return rhs_idcs_;
    }
 
    std::vector<size_t> convert_to_rhs_indexes(const std::vector<int>& indexes) const;
 
    virtual size_t get_num_dof() const = 0;
 
    // unused, but override just to satisfy WmPatchProcess interface
    // spatial_solver probably shouldn't be a patch process, but this is the easiest
    // way to get something which needs one object per local patch
    void process() override
    {
    }
};
 
} // namespace solver
} // namespace dfem
} // namespace wxm
 
#endif // WMSPATIALSOLVER_H