PFSolverPowerPolarSparse.h

// SPDX-FileCopyrightText: 2026 Institute for Automation of Complex Power Systems, EONERC, RWTH Aachen University
// SPDX-License-Identifier: MPL-2.0
 
#pragma once
 
#include <dpsim/DirectLinearSolver.h>
#include <dpsim/PFSolverPowerPolar.h>
 
namespace DPsim {
class PFSolverPowerPolarSparse : public PFSolverPowerPolar {
protected:
  CPS::SparseMatrixRow mJsparse;
  std::shared_ptr<DirectLinearSolver> mLinearSolver;
  std::vector<std::pair<CPS::UInt, CPS::UInt>> mVariableSystemMatrixEntries;

  void setUpJacobianStorage() override;
  void calculateJacobian() override;
  void solveJacobianSystem() override;

  void buildJacobianPattern();
  bool isConnected(CPS::UInt k, CPS::UInt j);
 
public:
  PFSolverPowerPolarSparse(CPS::String name, const CPS::SystemTopology &system,
                           CPS::Real timeStep, CPS::Logger::Level logLevel);
  virtual ~PFSolverPowerPolarSparse(){};
};
} // namespace DPsim