Add Power Electronics sparsity pattern computation into allocate

CHANGELOG.md

@@ -37,6 +37,7 @@
 - Added support for DependencyTracking::Variable in PowerElectronics models.
 - Updated Jacobian value storage from `ScalarT` to `RealT`.
 - Added a header file defining constants to be used throughout the code.
+- Added Jacobian sparsity pattern computation into `PowerElectronics`.
 
 ## v0.1
 

GridKit/LinearAlgebra/SparseMatrix/COO_Matrix.hpp

@@ -292,11 +292,6 @@ namespace GridKit
     template <typename RealT, typename IdxT>
     inline void COO_Matrix<RealT, IdxT>::axpy(RealT alpha, COO_Matrix<RealT, IdxT>& a)
     {
-      if (alpha == 0)
-      {
-        return;
-      }
-
       if (!this->sorted_)
       {
         this->sortSparse();
@@ -372,9 +367,6 @@ namespace GridKit
     template <typename RealT, typename IdxT>
     inline void COO_Matrix<RealT, IdxT>::axpy(RealT alpha, std::vector<IdxT> r, std::vector<IdxT> c, std::vector<RealT> v)
     {
-      if (alpha == 0)
-        return;
-
       if (!this->sorted_)
       {
         this->sortSparse();

GridKit/LinearAlgebra/SparseMatrix/CsrMatrix.cpp

@@ -818,6 +818,26 @@ namespace GridKit
       }
     }
 
+    /**
+     * @brief Changes the size of the matrix. De-allocates the matrix, if needed, since the
+     *        old indices corresponded to a differently-size matrix and no longer make sense.
+     *        Does not re-allocate any memory.
+     * @param n The new number of rows
+     * @param m The new number of columns
+     * @return 0 if succesful
+     */
+    template <typename RealT, typename IdxT>
+    int CsrMatrix<RealT, IdxT>::resize(IdxT n, IdxT m)
+    {
+      destroyMatrixData(memory::HOST);
+      destroyMatrixData(memory::DEVICE);
+
+      n_ = n;
+      m_ = m;
+
+      return 0;
+    }
+
     template class CsrMatrix<double, long int>;
     template class CsrMatrix<double, size_t>;
   } // namespace LinearAlgebra

GridKit/LinearAlgebra/SparseMatrix/CsrMatrix.hpp

@@ -76,6 +76,8 @@ namespace GridKit
       virtual int setValuesPointer(RealT*              new_vals,
                                    memory::MemorySpace memspace);
 
+      int resize(IdxT n, IdxT m);
+
     private:
       IdxT n_{0};   ///< number of rows
       IdxT m_{0};   ///< number of columns

GridKit/Model/PowerElectronics/DistributedGenerator/DistributedGenerator.cpp

@@ -39,9 +39,9 @@ namespace GridKit
   {
     // internals [\delta_i, Pi, Qi, phi_di, phi_qi, gamma_di, gamma_qi, il_di, il_qi, vo_di, vo_qi, io_di, io_qi]
     // externals [\omega_ref, vba_out, vbb_out]
-    size_           = 16;
-    n_intern_       = 13;
-    n_extern_       = 3;
+    size_           = NUM_INTERNALS + NUM_EXTERNALS;
+    n_intern_       = NUM_INTERNALS;
+    n_extern_       = NUM_EXTERNALS;
     extern_indices_ = {0, 1, 2};
     idc_            = id;
   }

GridKit/Model/PowerElectronics/DistributedGenerator/DistributedGenerator.hpp

@@ -64,6 +64,9 @@ namespace GridKit
     using CircuitComponent<ScalarT, IdxT>::n_intern_;
 
   public:
+    static constexpr size_t NUM_INTERNALS = 13;
+    static constexpr size_t NUM_EXTERNALS = 3;
+
     DistributedGenerator(IdxT                                        id,
                          DistributedGeneratorParameters<RealT, IdxT> parm,
                          bool                                        reference_frame);

GridKit/Model/PowerElectronics/MicrogridBusDQ/MicrogridBusDQ.cpp

@@ -24,9 +24,9 @@ namespace GridKit
     : RN_(RN)
   {
     // externals [vbus_d, vbus_q]
-    size_           = 2;
-    n_intern_       = 0;
-    n_extern_       = 2;
+    size_           = NUM_INTERNALS + NUM_EXTERNALS;
+    n_intern_       = NUM_INTERNALS;
+    n_extern_       = NUM_EXTERNALS;
     extern_indices_ = {0, 1};
     idc_            = id;
   }

GridKit/Model/PowerElectronics/MicrogridBusDQ/MicrogridBusDQ.hpp

@@ -43,6 +43,9 @@ namespace GridKit
     using CircuitComponent<ScalarT, IdxT>::n_intern_;
 
   public:
+    static constexpr size_t NUM_INTERNALS = 0;
+    static constexpr size_t NUM_EXTERNALS = 2;
+
     MicrogridBusDQ(IdxT id, RealT RN);
     virtual ~MicrogridBusDQ();
 

GridKit/Model/PowerElectronics/MicrogridLine/MicrogridLine.cpp

@@ -28,9 +28,9 @@ namespace GridKit
   {
     // internals [id, iq]
     // externals [\omegaref, vbd_in, vbq_in, vbd_out, vbq_out]
-    size_           = 7;
-    n_intern_       = 2;
-    n_extern_       = 5;
+    size_           = NUM_INTERNALS + NUM_EXTERNALS;
+    n_intern_       = NUM_INTERNALS;
+    n_extern_       = NUM_EXTERNALS;
     extern_indices_ = {0, 1, 2, 3, 4};
     idc_            = id;
   }

GridKit/Model/PowerElectronics/MicrogridLine/MicrogridLine.hpp

@@ -44,6 +44,9 @@ namespace GridKit
     using CircuitComponent<ScalarT, IdxT>::n_intern_;
 
   public:
+    static constexpr size_t NUM_INTERNALS = 2;
+    static constexpr size_t NUM_EXTERNALS = 5;
+
     MicrogridLine(IdxT id, RealT R, RealT L);
     virtual ~MicrogridLine();
 

GridKit/Model/PowerElectronics/MicrogridLoad/MicrogridLoad.cpp

@@ -26,9 +26,9 @@ namespace GridKit
   {
     // internals [id, iq]
     // externals [\omegaref, vbd_out, vbq_out]
-    size_           = 5;
-    n_intern_       = 2;
-    n_extern_       = 3;
+    size_           = NUM_INTERNALS + NUM_EXTERNALS;
+    n_intern_       = NUM_INTERNALS;
+    n_extern_       = NUM_EXTERNALS;
     extern_indices_ = {0, 1, 2};
     idc_            = id;
   }

GridKit/Model/PowerElectronics/MicrogridLoad/MicrogridLoad.hpp

@@ -44,6 +44,9 @@ namespace GridKit
     using CircuitComponent<ScalarT, IdxT>::n_intern_;
 
   public:
+    static constexpr size_t NUM_INTERNALS = 2;
+    static constexpr size_t NUM_EXTERNALS = 3;
+
     MicrogridLoad(IdxT id, RealT R, RealT L);
     virtual ~MicrogridLoad();