WIP: Deal with imag prom

Author: phdum-a

palace/drivers/drivensolver.cpp

@@ -240,7 +240,6 @@ ErrorIndicator DrivenSolver::SweepAdaptive(SpaceOperator &space_op, int n_step,
   space_op.GetWavePortOp().SetSuppressOutput(
       true);  // Suppress wave port output for offline
 
-
   // Add ports to PROM if we do synthesis.
   if (iodata.solver.driven.adaptive_circuit_synthesis)
   {
@@ -493,22 +492,22 @@ void DrivenSolver::PrintPROMMatrices(const RomOperator &prom_op) const
   using VT = Units::ValueType;
 
   // PROM matrices should be real.
-  // TODO: Implications of periodic boundary conditions.
+  // TODO: Implications of periodic boundary conditions. Just print re & im separate.
   auto unit_H = iodata.units.GetScaleFactor<VT::INDUCTANCE>();
   Eigen::MatrixXd m_Linv = (1.0 / unit_H * v_d) * Kr.real() * v_d;
-  print_table(m_Linv, "prom-Linv.csv");
+  print_table(m_Linv, "prom-Linv-re.csv");
 
   auto unit_F = iodata.units.GetScaleFactor<VT::CAPACITANCE>();
   Eigen::MatrixXd m_C = (unit_F * v_d) * Mr.real() * v_d;
-  print_table(m_C, "prom-C.csv");
+  print_table(m_C, "prom-C-re.csv");
 
   // Cr need not exist, if no dissipation. Currently, Cr always exists since we need
   // dissipative ports for a  driven response, but this may change in the future.
   if (Cr.size() > 0)
   {
     auto unit_ohm = iodata.units.GetScaleFactor<VT::IMPEDANCE>();
     Eigen::MatrixXd m_Rinv = (1.0 / unit_ohm * v_d) * Cr.real() * v_d;
-    print_table(m_Rinv, "prom-Rinv.csv");
+    print_table(m_Rinv, "prom-Rinv-re.csv");
   }
 }
 }  // namespace palace