stops creating plots with matplot++ if plotOn = false

empennage_design/src/taw/t_tail/low/low_t_tail_empennage_design_plot.cpp

@@ -13,12 +13,12 @@
 #include <aircraftGeometry2/processing/transform.h>
 #include <matplot/matplot.h>
 
+#include <algorithm>
 #include <array>
 #include <cstdlib>
 #include <deque>
 #include <map>
 #include <utility>
-#include <algorithm>
 
 #include "low_t_tail_empennage_design.h"
 
@@ -27,7 +27,8 @@
  *
  * \param surface aerodynamic surface object
  * \param object spar or control device object
- * \return std::tuple<std::vector<double>,std::vector<double>,std::string> coordinates and object name
+ * \return std::tuple<std::vector<double>,std::vector<double>,std::string>
+ * coordinates and object name
  */
 std::tuple<std::vector<double>, std::vector<double>, std::string> gen_polygon(
     const geom2::MultisectionSurface<geom2::AirfoilSection>& surface,
@@ -42,75 +43,96 @@ std::tuple<std::vector<double>, std::vector<double>, std::string> gen_polygon(
  * \param a alpha value between 0 (invisible) - 1 (solid)
  * \return std::array<float,4UL> matplot++ conform color format
  */
-//std::array<float, 4UL> rgba(float r, float b, float g, float a) { return {1 - a, r / 255, b / 255, g / 255}; }
+// std::array<float, 4UL> rgba(float r, float b, float g, float a) { return {1 -
+// a, r / 255, b / 255, g / 255}; }
 //
 ///* */
-//std::array<float, 4UL> spar_color = rgba(99, 110, 114, 1.0);
+// std::array<float, 4UL> spar_color = rgba(99, 110, 114, 1.0);
 //
-//std::map<std::string, std::array<float, 4UL>> control_device_colormap = {{"elevator", rgba(108, 92, 231, 1.0)},
-//                                                                         {"rudder", rgba(232, 67, 147, 1.0)}};
+// std::map<std::string, std::array<float, 4UL>> control_device_colormap =
+// {{"elevator", rgba(108, 92, 231, 1.0)},
+//                                                                          {"rudder",
+//                                                                          rgba(232,
+//                                                                          67,
+//                                                                          147, 1.0)}};
 //
-//std::tuple<std::vector<double>, std::vector<double>, std::string> gen_polygon(
-//    const geom2::MultisectionSurface<geom2::AirfoilSection>& surface,
-//    const geom2::MultisectionSurface<geom2::PolygonSection>& object) {
-//  double half_span = surface.is_symmetric ? 0.5 * geom2::measure::span(surface) : geom2::measure::span(surface);
-//  /* Spar coordinates positive values normally*/
-//  double from = object.sections.at(0).origin.z() * half_span;
-//  double to = object.sections.at(1).origin.z() * half_span;
+// std::tuple<std::vector<double>, std::vector<double>, std::string>
+// gen_polygon(
+//     const geom2::MultisectionSurface<geom2::AirfoilSection>& surface,
+//     const geom2::MultisectionSurface<geom2::PolygonSection>& object) {
+//   double half_span = surface.is_symmetric ? 0.5 *
+//   geom2::measure::span(surface) : geom2::measure::span(surface);
+//   /* Spar coordinates positive values normally*/
+//   double from = object.sections.at(0).origin.z() * half_span;
+//   double to = object.sections.at(1).origin.z() * half_span;
 //
-//  /* Chordwise coordinates */
-//  double chord_from_fwd = object.sections.at(0).get_contour().vertex(0).x();
-//  double chord_from_aft = object.sections.at(0).get_contour().vertex(1).x();
-//  double chord_to_fwd = object.sections.at(1).get_contour().vertex(0).x();
-//  double chord_to_aft = object.sections.at(1).get_contour().vertex(1).x();
+//   /* Chordwise coordinates */
+//   double chord_from_fwd = object.sections.at(0).get_contour().vertex(0).x();
+//   double chord_from_aft = object.sections.at(0).get_contour().vertex(1).x();
+//   double chord_to_fwd = object.sections.at(1).get_contour().vertex(0).x();
+//   double chord_to_aft = object.sections.at(1).get_contour().vertex(1).x();
 //
-//  auto f_chord_fwd = [chord_from_fwd, chord_to_fwd, half_span](double span_position) {
-//    return std::lerp(chord_from_fwd, chord_to_fwd, fabs(span_position) / half_span);
-//  };
-//  auto f_chord_aft = [chord_from_aft, chord_to_aft, half_span](double span_position) {
-//    return std::lerp(chord_from_aft, chord_to_aft, fabs(span_position) / half_span);
-//  };
-//  /* Add object from */
-//  std::vector<double> spanwise_coordinates{-fabs(from)};
+//   auto f_chord_fwd = [chord_from_fwd, chord_to_fwd, half_span](double
+//   span_position) {
+//     return std::lerp(chord_from_fwd, chord_to_fwd, fabs(span_position) /
+//     half_span);
+//   };
+//   auto f_chord_aft = [chord_from_aft, chord_to_aft, half_span](double
+//   span_position) {
+//     return std::lerp(chord_from_aft, chord_to_aft, fabs(span_position) /
+//     half_span);
+//   };
+//   /* Add object from */
+//   std::vector<double> spanwise_coordinates{-fabs(from)};
 //
-//  /* Add intermediate sections */
-//  for (auto& section : surface.sections) {
-//    if (fabs(section.origin.z()) > fabs(from) && fabs(section.origin.z()) < fabs(to)) {
-//      spanwise_coordinates.push_back(section.origin.z());
-//    }
-//  }
-//  /* Add object to */
-//  spanwise_coordinates.push_back(-fabs(to));
+//   /* Add intermediate sections */
+//   for (auto& section : surface.sections) {
+//     if (fabs(section.origin.z()) > fabs(from) && fabs(section.origin.z()) <
+//     fabs(to)) {
+//       spanwise_coordinates.push_back(section.origin.z());
+//     }
+//   }
+//   /* Add object to */
+//   spanwise_coordinates.push_back(-fabs(to));
 //
-//  /**/
-//  std::vector<double> x;
-//  std::vector<double> y;
-//  /* From in to out*/
-//  for (auto it = spanwise_coordinates.begin(); it != spanwise_coordinates.end(); ++it) {
-//    y.push_back(fabs(*it));
-//    x.push_back(geom2::measure::offset_LE(surface, *it).x() + f_chord_fwd(*it) * geom2::measure::chord(surface, *it));
-//  }
-//  /* From out to in*/
-//  for (auto rit = spanwise_coordinates.rbegin(); rit != spanwise_coordinates.rend(); ++rit) {
-//    y.push_back(fabs(*rit));
-//    x.push_back(geom2::measure::offset_LE(surface, *rit).x() +
-//                f_chord_aft(*rit) * geom2::measure::chord(surface, *rit));
-//  }
+//   /**/
+//   std::vector<double> x;
+//   std::vector<double> y;
+//   /* From in to out*/
+//   for (auto it = spanwise_coordinates.begin(); it !=
+//   spanwise_coordinates.end(); ++it) {
+//     y.push_back(fabs(*it));
+//     x.push_back(geom2::measure::offset_LE(surface, *it).x() +
+//     f_chord_fwd(*it) * geom2::measure::chord(surface, *it));
+//   }
+//   /* From out to in*/
+//   for (auto rit = spanwise_coordinates.rbegin(); rit !=
+//   spanwise_coordinates.rend(); ++rit) {
+//     y.push_back(fabs(*rit));
+//     x.push_back(geom2::measure::offset_LE(surface, *rit).x() +
+//                 f_chord_aft(*rit) * geom2::measure::chord(surface, *rit));
+//   }
 //
-//  /* Close loop */
-//  y.push_back(y.front());
-//  x.push_back(x.front());
+//   /* Close loop */
+//   y.push_back(y.front());
+//   x.push_back(x.front());
 //
-//  return {x, y, object.name};
-//}
+//   return {x, y, object.name};
+// }
 namespace taw {
 namespace low {
 
 fs::path T_Tail::plot_stabilizer_planform(
     const geom2::MultisectionSurface<geom2::AirfoilSection>& surface,
     const std::vector<geom2::MultisectionSurface<geom2::PolygonSection>>& spars,
-    const std::vector<geom2::MultisectionSurface<geom2::PolygonSection>>& devices) {
-  auto coordinates_xy = geom2::transform::outline_2d(surface, geom2::Direction_3(0, 0, 1));
+    const std::vector<geom2::MultisectionSurface<geom2::PolygonSection>>&
+        devices) {
+  /* Plot active */
+  if (!rtIO->plotOn) {
+    return "";
+  }
+  auto coordinates_xy =
+      geom2::transform::outline_2d(surface, geom2::Direction_3(0, 0, 1));
   std::deque<double> spanwise;
   std::deque<double> chordwise;
   auto span = geom2::measure::span(surface);
@@ -130,9 +152,10 @@ fs::path T_Tail::plot_stabilizer_planform(
   }
 
   // /* Get mac coordinates */
-  // auto mac_position_y = geom2::measure::mean_aerodynamic_chord_position(data->wing);
-  // auto mac_position = geom2::measure::offset_LE(data->wing,mac_position_y);
-  // auto mac = geom2::measure::chord(data->wing,mac_position_y);
+  // auto mac_position_y =
+  // geom2::measure::mean_aerodynamic_chord_position(data->wing); auto
+  // mac_position = geom2::measure::offset_LE(data->wing,mac_position_y); auto
+  // mac = geom2::measure::chord(data->wing,mac_position_y);
 
   /* Generate plot */
   std::string name = surface.name;
@@ -199,7 +222,8 @@ fs::path T_Tail::plot_stabilizer_planform(
   ax->legend()->box(true);
   ax->legend()->font_size(8);
 
-  fs::path plot_path = rtIO->getPlotDir() + "/" + surface.name + "_planfrom.svg";
+  fs::path plot_path =
+      rtIO->getPlotDir() + "/" + surface.name + "_planfrom.svg";
   matplot::save(f, plot_path.string(), "svg");
   std::this_thread::sleep_for(std::chrono::milliseconds(500));
   system("pkill gnuplot");
@@ -208,27 +232,36 @@ fs::path T_Tail::plot_stabilizer_planform(
 
 fs::path T_Tail::plot_thickness_and_twist_distribution(
     const geom2::MultisectionSurface<geom2::AirfoilSection>& surface) {
+  /* Plot active */
+  if (!rtIO->plotOn) {
+    return "";
+  }
   std::vector<double> eta;
   std::vector<double> thickness_to_chord_ratio;
   std::vector<double> twist;
-  double half_span = surface.is_symmetric ? geom2::measure::span(surface) * 0.5 : geom2::measure::span(surface);
+  double half_span = surface.is_symmetric ? geom2::measure::span(surface) * 0.5
+                                          : geom2::measure::span(surface);
   for (auto section : surface.sections) {
     eta.push_back(-section.origin.z() / half_span);
-    thickness_to_chord_ratio.push_back(geom2::measure::thickness_max(section) / section.get_chord_length());
+    thickness_to_chord_ratio.push_back(geom2::measure::thickness_max(section) /
+                                       section.get_chord_length());
     twist.push_back(section.get_twist_angle());
   }
-  const auto [t_to_c_min,t_to_c_max] = std::ranges::minmax_element(thickness_to_chord_ratio);
-  const auto [twist_min,twist_max] = std::ranges::minmax_element(twist);
+  const auto [t_to_c_min, t_to_c_max] =
+      std::ranges::minmax_element(thickness_to_chord_ratio);
+  const auto [twist_min, twist_max] = std::ranges::minmax_element(twist);
 
   auto f = matplot::figure(true);
 
   auto ax1 = f->add_subplot(2, 1, 0);
-  ax1->plot(eta, thickness_to_chord_ratio, "g-")->color("black").line_width(1.5);
+  ax1->plot(eta, thickness_to_chord_ratio, "g-")
+      ->color("black")
+      .line_width(1.5);
   // ax1->xlabel("\u03B7");
   ax1->x_axis().visible(true);
   ax1->grid(true);
   ax1->ylabel("t / c [1]");
-  ax1->ylim({*t_to_c_min - 0.01,*t_to_c_max + 0.01});
+  ax1->ylim({*t_to_c_min - 0.01, *t_to_c_max + 0.01});
 
   auto ax2 = f->add_subplot(2, 1, 1);
   ax2->plot(eta, twist, "b-")->color("black").line_style("--").line_width(1.5);
@@ -237,7 +270,8 @@ fs::path T_Tail::plot_thickness_and_twist_distribution(
   ax2->grid(true);
   ax2->ylim({*twist_min - 0.01, *twist_max + 0.01});
 
-  fs::path plot_path = rtIO->getPlotDir() + "/" + surface.name + "_thickness_and_twist_distribution.svg";
+  fs::path plot_path = rtIO->getPlotDir() + "/" + surface.name +
+                       "_thickness_and_twist_distribution.svg";
   matplot::save(f, plot_path.string(), "svg");
   std::this_thread::sleep_for(std::chrono::milliseconds(500));
   system("pkill gnuplot");
@@ -249,7 +283,12 @@ fs::path T_Tail::plot_thickness_and_twist_distribution(
  *
  * \return fs::path
  */
-std::vector<fs::path> T_Tail::plot_airfoils(const geom2::MultisectionSurface<geom2::AirfoilSection>& surface) {
+std::vector<fs::path> T_Tail::plot_airfoils(
+    const geom2::MultisectionSurface<geom2::AirfoilSection>& surface) {
+  /* Plot active */
+  if (!rtIO->plotOn) {
+    return {""};
+  }
   std::map<std::string, std::vector<geom2::Point_2>> airfoils_data{};
   std::vector<fs::path> airfoil_plots_path;
 
@@ -284,9 +323,11 @@ std::vector<fs::path> T_Tail::plot_airfoils(const geom2::MultisectionSurface<geo
       ax->xlim({0, 1});
       ax->axis(matplot::equal);
 
-      fs::path plot_path = rtIO->getPlotDir() + "/" + surface.name + "_" + airfoil_data.first + "_airfoil.svg";
+      fs::path plot_path = rtIO->getPlotDir() + "/" + surface.name + "_" +
+                           airfoil_data.first + "_airfoil.svg";
       matplot::save(f, plot_path.string(), "svg");
-      airfoil_plots_path.push_back(fs::relative(plot_path, fs::path(rtIO->getHtmlDir())));
+      airfoil_plots_path.push_back(
+          fs::relative(plot_path, fs::path(rtIO->getHtmlDir())));
       std::this_thread::sleep_for(std::chrono::milliseconds(500));
       system("pkill gnuplot");
     }