Commit fe349522 authored by Thibault Hallouin's avatar Thibault Hallouin
Browse files

Merge remote-tracking branch 'origin/dev' into dev

1 merge request!3release v0.1.0
Pipeline #42712 passed with stage
in 2 minutes and 15 seconds
Showing with 11 additions and 11 deletions
+11 -11
...@@ -216,7 +216,7 @@ namespace evalhyd ...@@ -216,7 +216,7 @@ namespace evalhyd
// else if m_cdt provided, use them to generate t_msk // else if m_cdt provided, use them to generate t_msk
else if (m_cdt.size() > 0) else if (m_cdt.size() > 0)
{ {
xt::xtensor<bool, 2> c_msk = xt::zeros<bool>({n_msk, n_tim}); B2 c_msk = xt::zeros<bool>({n_msk, n_tim});
for (int m = 0; m < n_msk; m++) for (int m = 0; m < n_msk; m++)
xt::view(c_msk, m) = xt::view(c_msk, m) =
...@@ -228,7 +228,7 @@ namespace evalhyd ...@@ -228,7 +228,7 @@ namespace evalhyd
} }
// if neither t_msk nor m_cdt provided, generate dummy mask // if neither t_msk nor m_cdt provided, generate dummy mask
else else
return xt::xtensor<bool, 2>{xt::ones<bool>({std::size_t{1}, n_tim})}; return B2({xt::ones<bool>({std::size_t{1}, n_tim})});
}; };
auto msk = gen_msk(); auto msk = gen_msk();
...@@ -242,7 +242,7 @@ namespace evalhyd ...@@ -242,7 +242,7 @@ namespace evalhyd
} }
else if ( transform == "sqrt" ) else if ( transform == "sqrt" )
{ {
return xt::eval(xt::sqrt(q)); return D2(xt::sqrt(q));
} }
else if ( transform == "inv" ) else if ( transform == "inv" )
{ {
...@@ -250,7 +250,7 @@ namespace evalhyd ...@@ -250,7 +250,7 @@ namespace evalhyd
// determine an epsilon value to avoid zero divide // determine an epsilon value to avoid zero divide
epsilon = xt::mean(q_obs_)() * 0.01; epsilon = xt::mean(q_obs_)() * 0.01;
return xt::eval(1. / (q + epsilon)); return D2(1. / (q + epsilon));
} }
else if ( transform == "log" ) else if ( transform == "log" )
{ {
...@@ -258,7 +258,7 @@ namespace evalhyd ...@@ -258,7 +258,7 @@ namespace evalhyd
// determine an epsilon value to avoid log zero // determine an epsilon value to avoid log zero
epsilon = xt::mean(q_obs_)() * 0.01; epsilon = xt::mean(q_obs_)() * 0.01;
return xt::eval(xt::log(q + epsilon)); return D2(xt::log(q + epsilon));
} }
else if ( transform == "pow" ) else if ( transform == "pow" )
{ {
...@@ -268,11 +268,11 @@ namespace evalhyd ...@@ -268,11 +268,11 @@ namespace evalhyd
// determine an epsilon value to avoid zero divide // determine an epsilon value to avoid zero divide
epsilon = xt::mean(q_obs_)() * 0.01; epsilon = xt::mean(q_obs_)() * 0.01;
return xt::eval(xt::pow(q + epsilon, exponent)); return D2(xt::pow(q + epsilon, exponent));
} }
else else
{ {
return xt::eval(xt::pow(q, exponent)); return D2(xt::pow(q, exponent));
} }
} }
else else
...@@ -283,8 +283,8 @@ namespace evalhyd ...@@ -283,8 +283,8 @@ namespace evalhyd
} }
}; };
auto obs = q_transform(q_obs_); const D2 obs = q_transform(q_obs_);
auto prd = q_transform(q_prd_); const D2 prd = q_transform(q_prd_);
// generate bootstrap experiment if requested // generate bootstrap experiment if requested
std::vector<xt::xkeep_slice<int>> exp; std::vector<xt::xkeep_slice<int>> exp;
......
...@@ -246,7 +246,7 @@ namespace evalhyd ...@@ -246,7 +246,7 @@ namespace evalhyd
// generate masks from conditions if provided // generate masks from conditions if provided
auto gen_msk = [&]() { auto gen_msk = [&]() {
xt::xtensor<bool, 4> c_msk = xt::zeros<bool>({n_sit, n_ltm, n_msk, n_tim}); B4 c_msk = xt::zeros<bool>({n_sit, n_ltm, n_msk, n_tim});
if (m_cdt.size() > 0) if (m_cdt.size() > 0)
for (int s = 0; s < n_sit; s++) for (int s = 0; s < n_sit; s++)
for (int l = 0; l < n_ltm; l++) for (int l = 0; l < n_ltm; l++)
...@@ -259,7 +259,7 @@ namespace evalhyd ...@@ -259,7 +259,7 @@ namespace evalhyd
); );
return c_msk; return c_msk;
}; };
const xt::xtensor<bool, 4> c_msk = gen_msk(); const B4 c_msk = gen_msk();
// generate bootstrap experiment if requested // generate bootstrap experiment if requested
std::vector<xt::xkeep_slice<int>> b_exp; std::vector<xt::xkeep_slice<int>> b_exp;
......
Supports Markdown
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