richards_fs_height.h

Go to the documentation of this file.

//
// Created by julian on 2/26/26.
//
 
#ifndef UG4_RICHARDS_FS_HEIGHT_H
#define UG4_RICHARDS_FS_HEIGHT_H
#include "../../level_set_pos.h"
#include "lib_disc/function_spaces/grid_function_global_user_data.h"
namespace ug {
namespace d3f {
 
template <typename TGridFunction>
class RichardsFreeSurfaceHeight
:   public StdDependentUserData<RichardsFreeSurfaceHeight<TGridFunction>, number, TGridFunction::dim> {
 
public:
    typedef TGridFunction gf_type;
    static const int dim = gf_type::dim;
    typedef RichardsFreeSurfaceHeight<gf_type> this_type;
    typedef StdDependentUserData<this_type, number, dim> base_type;
    typedef typename gf_type::domain_type domain_type;
    typedef typename domain_type::position_accessor_type pos_acc_type;
 
    typedef LSPositionZ<TGridFunction> measurer_type;
 
    typedef GlobalGridFunctionNumberData<TGridFunction> glob_userdata_type;
 
    typedef Attachment<size_t> num_a_type;
 
    typedef Attachment<number> val_a_type;
 
    typedef Grid::VertexAttachmentAccessor<num_a_type> num_aa_type;
 
    typedef Grid::VertexAttachmentAccessor<val_a_type> val_aa_type;
 
    RichardsFreeSurfaceHeight(
        SmartPtr<gf_type> spGf,
        const char * subset = "top",
        const char * cmp = "p"
        );
 
    void reinit();
 
    virtual bool continuous () const {return true;}
 
    virtual bool requires_grid_fct () const {return true;}
 
    template <int refDim>
    void eval_and_deriv
    (
        number vValue[], 
        const MathVector<dim> vGlobIP[],
        number time,
        int si,
        GridObject* elem,
        const MathVector<dim> vCornerCoords[],
        const MathVector<refDim> vLocIP[],
        const size_t nip,
        LocalVector* u,
        bool bDeriv,
        int s,
        std::vector<std::vector<number> > vvvDeriv[],
        const MathMatrix<refDim, dim>* vJT = NULL
    );
private:
    struct t_pnt_data
    {
        MathVector<dim-1> xy; 
        number z; 
        bool valid; 
 
        t_pnt_data (MathVector<dim-1> & the_xy)
        : xy (the_xy), valid(false) {};
    };
    std::vector<t_pnt_data> m_measurement_points;
    SubsetGroup m_SsGrp;
    SmartPtr<gf_type> m_spGf;
    const char * m_cmp;
    glob_userdata_type m_glob_userdata;
    num_a_type m_a_vrt_num;
    val_a_type m_a_vrt_val;
    SmartPtr<std::vector<number>> m_sp_vrt_val;
    measurer_type m_sp_ls_pos_z;
};
 
 
 
 
}
}
 
#include "richards_fs_height_impl.h"
#endif //UG4_USERDATA_VAL_SURF_PROJ_H