[DLG] Added pi_div_2 to constants

juanmanzanero · Dec 24, 2023 · b7faf4b · b7faf4b
1 parent f6f6999
commit b7faf4b
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 1 deletion.
diff --git a/lion/foundation/constants.h b/lion/foundation/constants.h
@@ -8,6 +8,9 @@ enum Axis{ X, Y, Z };
 template<typename T = double>
 constexpr T pi_T = T{ 3.1415926535897932385L };
 
+template<typename T = double>
+constexpr T pi_div_2_T = T{ 1.570796326794896558L };
+
 template<typename T = double>
 constexpr T eps_T = std::numeric_limits<T>::epsilon();
 
@@ -20,6 +23,7 @@ constexpr T nan_T = std::numeric_limits<T>::quiet_NaN();
 
 // handy names when using double, our usual case
 constexpr auto pi{ pi_T<double> };
+constexpr auto pi_div_2{ pi_div_2_T<double> };
 constexpr auto eps{ eps_T<double> };
 constexpr auto inf{ inf_T<double> };
 constexpr auto Inf{ inf };

diff --git a/lion/io/Xml_element.h b/lion/io/Xml_element.h
@@ -399,15 +399,18 @@ void Xml_element::add_children_and_values_from_map(const MapType &m)
 
     for (const auto &mi : m) {
         std::ostringstream ss;
-        ss << std::setprecision(17);
 
         if constexpr (map_of_vectors) {
+            ss << std::setprecision(std::numeric_limits<typename mapped_type::value_type>::max_digits10);
+
             ss << mi.second.front();
             for (auto mis = std::next(mi.second.cbegin()); mis != mi.second.cend(); ++mis) {
                 ss << ", " << *mis;
             }
         }
         else {
+            ss << std::setprecision(std::numeric_limits<mapped_type>::max_digits10);
+
             ss << mi.second;
         }
         add_child(mi.first).set_value(ss.str());

diff --git a/lion/math/ipopt_optimize_nlp_finite_differences.h b/lion/math/ipopt_optimize_nlp_finite_differences.h
@@ -156,6 +156,7 @@ struct Ipopt_optimize_NLP_finite_differences
 
 
         // configure ipopt's options
+        app->Options()->SetStringValue("sb", "yes");
         app->Options()->SetStringValue("linear_solver", opts.linear_solver);
         if (opts.linear_solver == "ma97") {
             app->Options()->SetStringValue("ma97_order", "metis");