EMT_SSNComp.h

// SPDX-FileCopyrightText: 2026 Institute for Automation of Complex Power Systems, EONERC, RWTH Aachen University
// SPDX-License-Identifier: MPL-2.0
 
#pragma once
 
#include <array>
#include <vector>
 
#include <dpsim-models/MNASimPowerComp.h>
#include <dpsim-models/Solver/MNAInterface.h>
 
namespace CPS {
namespace EMT {

class SSNComp : public MNASimPowerComp<Real> {
private:
  const Int mInputSize;
  const Int mOutputSize;
 
protected:
  Real mTimeStep;
 
  Matrix mW;
  Matrix mYHist;
 
  Matrix mA;
  Matrix mB;
  Matrix mC;
  Matrix mD;
 
  Matrix mdA;
  Matrix mdB;
 
  const Attribute<Matrix>::Ptr mX;
 
  SSNComp(String uid, String name, Int inputSize, Int outputSize,
          Logger::Level logLevel = Logger::Level::off);
 
  virtual Matrix calculateHistoryVector() const;
 
  virtual MatrixComp calculateSteadyStateStateFromInput(const MatrixComp &u,
                                                        Real frequency) const;
  virtual MatrixComp
  calculateSteadyStateOutputFromInput(const MatrixComp &x,
                                      const MatrixComp &u) const;
 
  virtual void updateState(const Matrix &uOld, const Matrix &uNew);

  virtual void updateLogAttributes(const Matrix &u) const;
 
  virtual void recomputeDiscreteModel();

  virtual void updateStateSpaceModel();
 
  virtual Attribute<Matrix>::Ptr inputAttribute() const = 0;
  virtual Attribute<Matrix>::Ptr outputAttribute() const = 0;
 
public:
  struct LocalAbcStateBlock {
    std::array<UInt, 3> indices;

    String name;
  };

  UInt getStateCount() const;

  virtual std::vector<String> getLocalStateNames() const;

  virtual std::vector<LocalAbcStateBlock> getLocalAbcStateBlocks() const;

  const Matrix &getDiscreteA() const;

  const Matrix &getDiscreteB() const;

  const Matrix &getC() const;
 
  void setParameters(const Matrix &A, const Matrix &B, const Matrix &C,
                     const Matrix &D);
 
  void mnaCompInitialize(Real omega, Real timeStep,
                         Attribute<Matrix>::Ptr leftVector) override final;
 
  void mnaCompAddPreStepDependencies(
      AttributeBase::List &prevStepDependencies,
      AttributeBase::List &attributeDependencies,
      AttributeBase::List &modifiedAttributes) override final;
 
  void mnaCompPreStep(Real time, Int timeStepCount) override final;
 
  void mnaCompAddPostStepDependencies(
      AttributeBase::List &prevStepDependencies,
      AttributeBase::List &attributeDependencies,
      AttributeBase::List &modifiedAttributes,
      Attribute<Matrix>::Ptr &leftVector) override final;
};
 
} // namespace EMT
} // namespace CPS