printing residual vector too

Author: shakedregev

examples/ScaleMicrogrid/ScaleMicrogrid.cpp

@@ -306,7 +306,8 @@ int test(index_type Nsize, real_type error_tol, bool debug_output)
 
   sys_model->initialize();
   sys_model->evaluateResidual();
-
+  // print the residual in matrix market format
+  sys_model->printResidualMatrixMarket("ScaleMicrogrid_Residual_N" + std::to_string(Nsize) + ".mtx", "ScaleMicrogrid Residual N" + std::to_string(Nsize));
   std::vector<real_type>& fres = sys_model->getResidual();
   if (debug_output)
   {

src/Model/PowerElectronics/SystemModelPowerElectronics.hpp

@@ -12,6 +12,46 @@
 
 namespace GridKit
 {
+  /**
+   * Writes a vector to a file in Matrix Market format
+   * 
+   * @param vec The vector to write
+   * @param filename The name of the output file
+   * @param header Additional header information/comments
+   * @return true if the write was successful, false otherwise
+   */
+  template <typename T>
+  bool writeVectorToMatrixMarket(const std::vector<T>& vec, const std::string& filename, const std::string& header) {
+      std::ofstream outFile(filename);
+      
+      if (!outFile.is_open()) {
+          std::cerr << "Error: Could not open file " << filename << " for writing." << std::endl;
+          return false;
+      }
+      
+      // Write Matrix Market header
+      outFile << "%%MatrixMarket vector array real general" << std::endl;
+      
+      // Write additional header information as comments
+      if (!header.empty()) {
+          outFile << "% " << header << std::endl;
+      }
+      
+      // Write the vector size
+      outFile << vec.size() << std::endl;
+      
+      // Write the vector elements
+      outFile << std::scientific << std::setprecision(16);
+      for (const auto& val : vec) {
+          outFile << val << std::endl;
+      }
+      
+      outFile.close();
+      
+      std::cout << "Vector successfully written to " << filename << std::endl;
+      return true;
+  }
+
 
   template <class ScalarT, typename IdxT>
   class PowerElectronicsModel : public CircuitComponent<ScalarT, IdxT>
@@ -324,6 +364,17 @@ namespace GridKit
       alpha_ = a;
     }
 
+    /**
+     * @brief print the system residual in COO format
+     *
+     * @param[in] filename
+     * @param[in] title
+     */
+    void printResidualMatrixMarket(std::string filename, std::string title)
+    {
+      writeVectorToMatrixMarket(f_, filename, title);
+    }
+
     /**
      * @brief print the system Jacobian in COO format
      *