Linear.hpp

Go to the documentation of this file.

//
// Created by Andrew Nazareth on 2/16/24.
//
 
#ifndef OASIS_LINEAR_HPP
#define OASIS_LINEAR_HPP
 
#include "Eigen/Dense"
#include <iostream>
#include <map>
#include <string>
#include <vector>
 
#include "Expression.hpp"
 
namespace Oasis {
 
typedef Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> MatrixXXD;
typedef Eigen::Matrix<double, Eigen::Dynamic, 1> Matrix1D;
 
auto SolveLinearSystems(std::vector<std::unique_ptr<Expression>>& exprs) -> std::map<std::string, double>;
 
auto SolveLinearSystems(MatrixXXD& matrix) -> Matrix1D;
 
auto SolveLinearSystems(MatrixXXD& matrixA, Matrix1D& matrixb) -> Matrix1D;
 
auto ConstructMatrices(const std::vector<std::unique_ptr<Expression>>& exprs)
    -> std::pair<std::pair<MatrixXXD, Matrix1D>, std::map<std::string, Eigen::Index>>;
 
template <class u, class v>
auto GetMapValue(std::map<u, v>& map, u key, v nextValue) -> std::pair<v, bool>
{
    v value;
    bool usedNextValue = false;
    if (auto p = map.find(key); p != map.end()) { // key is found
        value = p->second;
    } else { // key is not found
        map.insert(std::make_pair(key, nextValue));
        value = nextValue;
        usedNextValue = true;
    }
    return std::make_pair(value, usedNextValue);
}
 
}
 
#endif // OASIS_LINEAR_HPP