Commit 68a73990 authored by Santiago Ospina's avatar Santiago Ospina

Modify the template definition for the CFL return type

Signed-off-by: Santiago Ospina De Los Ríos's avatarSantiago Ospina <santiago.ospina@iup.uni-heidelberg.de>
parent 3a9fee7f
......@@ -35,12 +35,12 @@ namespace Dorie{
* \f$||\beta||_{[L^\infty(T)]^d}h_T\f$.
*
* @tparam GF The grid function type.
* @tparam TF The time field. Type to represent time values.
* @tparam RF The range field. Type to represent CFL values.
*
* @return The CFL-condition (\f$\mathcal{CFL}\f$)
*/
template<class GF, class TF = double>
TF cfl_condition(const GF& gf, unsigned int intorder = 1)
template<class GF, class RF = double>
RF cfl_condition(const GF& gf, unsigned int intorder = 1)
{
using Traits = typename GF::Traits;
using DomainField = typename Traits::DomainFieldType;
......@@ -48,7 +48,7 @@ namespace Dorie{
const auto& gv = gf.getGridView();
TF suggested_dt = std::numeric_limits<TF>::max();
RF suggested_dt = std::numeric_limits<RF>::max();
// loop over the grid
for (const auto& cell : Dune::elements(gv,Dune::Partitions::interior))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment