KLUAdapter.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
 
extern "C" {
#include <klu.h>
}
 
#include <bitset>
#include <fstream>
#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 KLUAdapter : public DirectLinearSolver {
  std::vector<std::pair<UInt, UInt>> mChangedEntries;

  std::vector<Int> mVaryingColumns;
  std::vector<Int> mVaryingRows;

  klu_common mCommon;
  klu_numeric *mNumeric = nullptr;
  klu_symbolic *mSymbolic = nullptr;

  int mPreordering = AMD_ORDERING;

  int mPivotFaults = 0;
 
  PARTIAL_REFACTORIZATION_METHOD mPartialRefactorizationMethod =
      PARTIAL_REFACTORIZATION_METHOD::FACTORIZATION_PATH;

  Int nnz;
 
public:
  ~KLUAdapter() override;

  KLUAdapter();

  KLUAdapter(CPS::Logger::Log log);

  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;
 
protected:
  void printMatrixMarket(SparseMatrix &systemMatrix, int counter) const;

  void applyConfiguration() override;
};
} // namespace DPsim