SparseLUAdapter.h

/* Copyright 2017-2021 Institute for Automation of Complex Power Systems,
 *                     EONERC, RWTH Aachen University
 *
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at https://mozilla.org/MPL/2.0/.
 *********************************************************************************/
 
#pragma once
 
#include <bitset>
#include <iostream>
#include <list>
#include <memory>
#include <unordered_map>
#include <vector>
 
#include <dpsim/Config.h>
#include <dpsim/Definitions.h>
#include <dpsim/DirectLinearSolver.h>
 
namespace DPsim {
class SparseLUAdapter : public DirectLinearSolver {
  Eigen::SparseLU<CPS::SparseMatrixRow, Eigen::COLAMDOrdering<int>>
      LUFactorizedSparse;
 
public:
  using DirectLinearSolver::DirectLinearSolver;

  ~SparseLUAdapter() override;

  void preprocessing(SparseMatrix &systemMatrix,
                     std::vector<std::pair<UInt, UInt>>
                         &listVariableSystemMatrixEntries) override;

  void factorize(SparseMatrix &systemMatrix) override;

  void refactorize(SparseMatrix &systemMatrix) override;

  void partialRefactorize(SparseMatrix &systemMatrix,
                          std::vector<std::pair<UInt, UInt>>
                              &listVariableSystemMatrixEntries) override;

  Matrix solve(Matrix &rightSideVector) override;
};
} // namespace DPsim