matrix.hpp

Go to the documentation of this file.

#ifndef matrix_hpp
#define matrix_hpp
 
#include <functional>
#include <sstream>
 
namespace cpp_physics {
 
const std::out_of_range OOB = std::out_of_range("Row or column index out of bounds.");
const std::out_of_range SUBMAT_OOB = std::out_of_range("Submatrix exceeds bounds of original matrix.");
const std::invalid_argument SHAPE_MISMATCH = std::invalid_argument("Matrix-shapes don't match.");
const std::invalid_argument INNERDIM_MISMATCH = std::invalid_argument("Inner matrix-dimensions don't align.");
const std::invalid_argument DIM_MISMATCH = std::invalid_argument("Matrix-columns or -rows don't match.");
const std::invalid_argument TOO_FEW = std::invalid_argument("Too few matrices, minimum 1.");
const std::invalid_argument DIV_ZERO = std::invalid_argument("Division by 0.");
 
class Matrix {
private:
  std::vector<float> data;
  std::size_t rows;
  std::size_t cols;
 
public:
  Matrix(std::size_t rows, std::size_t cols = 1);
 
  static Matrix Identity(std::size_t size);
 
  std::pair<std::size_t, std::size_t> Shape() const;
 
  float Get(std::size_t row, std::size_t col) const;
 
  void applyFunction(std::function<float(float)> func);
 
  Matrix Submatrix(std::size_t rows, std::size_t cols, std::size_t row_start, std::size_t col_start) const;
 
  static Matrix Combine(std::vector<Matrix> matrices, bool vertical = false);
 
  float Sum() const;
 
  static Matrix memwise_Mul(const Matrix& lhs, const Matrix& rhs);
 
  Matrix memwise_iMul(const Matrix& other);
 
  static Matrix memwise_Div(const Matrix& lhs, const Matrix& rhs);
 
  Matrix memwise_iDiv(const Matrix& other);
 
  Matrix operator+=(const Matrix& other);
  Matrix operator+(const Matrix& other) const;
  Matrix operator-=(const Matrix& other);
  Matrix operator-(const Matrix& other) const;
  Matrix operator*=(const Matrix& other);
  Matrix operator*(const Matrix& other) const;
  Matrix operator*=(float scalar);
  Matrix operator*(float scalar) const;
  Matrix operator/=(float scalar);
  Matrix operator/(float scalar) const;
 
  Matrix operator=(const Matrix& other);
 
  Matrix operator-() const;
  Matrix operator!() const;
 
  float& operator[](std::size_t row, std::size_t col);
 
  bool operator==(const Matrix& other);
  bool operator!=(const Matrix& other);
};
 
Matrix operator*(float scalar, const Matrix& mat);
 
std::ostream& operator<<(std::ostream& flux, const Matrix& mat);
 
} // namespace cpp_physics
 
#endif