18 template <
class Model>
20 : cut_fraction_(cut_fraction)
24 template <
class Model>
43 template <
class Model>
46 return model_.conductivity_(p_head);
49 template <
class Model>
52 return model_.conductivity_(p_head);
55 template <
class Model>
58 return model_.water_content_(p_head);
61 template <
class Model>
64 return model_.water_content_(p_head);
73 : Bpar(0.5), Ppar(2), K_lower_limit(1.0E-20)
90 FFQs = pow(1 - pow(Qs_relative, 1/
m),
m);
127 T Kr, Q_unscaled, Q_cut_unscaled, FFQ;
130 Q_unscaled =
Q_rel(h);
132 FFQ = pow(1 - pow(Q_unscaled,1/
m),
m);
168 T Q = this->
Q_rel(h);
#define ASSERT_LT(a, b)
Definition of comparative assert macro (Less Than) only for debug mode.
#define ASSERT_GT(a, b)
Definition of comparative assert macro (Greater Than) only for debug mode.
fadbad::B< double > DiffDouble
auto water_content_diff(const DiffDouble &p_head) -> DiffDouble const override
auto conductivity_diff(const DiffDouble &p_head) -> DiffDouble const override
double conductivity(const double &p_head) const override
double water_content(const double &p_head) const override
void reset(SoilData data) override
SoilModelImplBase(double cut_fraction=0.999)
T conductivity_(const T &h) const
T water_content_(const T &h) const
void reset_(SoilData soil)
T Q_rel_inv(const T &q) const
T Q_rel(const T &h) const
T conductivity_(const T &h) const
k(h)=t(h)**(0.5)* (1- ((h)**n/(1+(h)**n)) **m)**2