Matrix<3, 3> m;                               // Create a 3x3 matrix filled with zeros per default
Matrix<3, 3> m({0, 1, 2, 3, 4, 5, 6, 7, 8});  // Create 3x3 matrix from array

// Sepcial types
Matrix<3, 3> m(ZEROS);    // 3x3 matrix explicitly filled with zeros
Matrix<3, 3> m(ONES);     // 3x3 matrix explicitly filled with ones
Matrix<3, 3> m(IDENTITY); // 3x3 identity matrix (ones along diagonal)

typedef Matrix<2, 2> Matrix2;
typedef Matrix<3, 3> Matrix3;
typedef Matrix<4, 4> Matrix4;

const Matrix<2, 2> MatrixI2(IDENTITY);
const Matrix<3, 3> MatrixI3(IDENTITY);
const Matrix<4, 4> MatrixI4(IDENTITY);

typedef Matrix<2, 1> VectorCol2;
typedef Matrix<3, 1> VectorCol3;
typedef Matrix<4, 1> VectorCol4;

typedef Matrix<1, 2> VectorRow2;
typedef Matrix<1, 3> VectorRow3;
typedef Matrix<1, 4> VectorRow4;

typedef VectorCol2 Vector2;
typedef VectorCol4 Vector3;
typedef VectorCol4 Vector4;

m.get(1, 2);   // Get the matrix element at row 1 column 2 (count starts at 0)
m.get(8);      // Get the matrix element at absolute index 8 (count starts at 0 and goes rowwise)
m[8];          // Get the matrix element at absolute index 8 (returns value and not a reference, thus setting value is not allowed)    

m.set(1, 2, 10.0f);   // Set the matrix element at row 1 column 2 to 10.0
m.set(8, 10.0f);      // Set the matrix element at absolute index 8 to 10.0

m.get_row(1);       // Get all elements of row 1 as array
m.get_col(1);       // Get all elements of column 1 as array

m.set_row(1, {1, 2, 3});       // Set all elements of row 1 from array
m.set_col(1, {1, 2, 3});       // Set all elements of column 1 from array

m.rows();       // Get the number of matrix rows
m.cols();       // Get the number of matrix columns

Matrix<2, 2> sm = m.submatrix(1, 2);  // Delete row 1 and column 2 from matrix to get submatrix

Matrix<2, 2> m1, m2;

Matrix<2, 2> sum = m1 + m2;      // Matrix elementwise addition
Matrix<2, 2> dif = m1 - m2;      // Matrix elementwise subtraction
Matrix<2, 2> pro = m1 * m2;      // Matrix multiplication
Matrix<2, 2> sca = m1 * 2.0f;    // (Elemetwise) multiplication by scalar
Matrix<2, 2> sca = m1 / 2.0f;    // (Elemetwise) division by scalar

m.transposed(); // Get transposed matrix

m.adjugate();   // Get adjugate matrix

m.inverse();    // Get inverse matrix

m.trace(); // Get trace of the matrix 

m.det(); // Get determinant of the matrix 

m1 == m2   // Returns true if the matrices have same dimensions and all values are the same (difference < 1e-9)
m1 != m2   // Opposite of ==

v.norm();                       // Get vector norm
v.norm2();                      // Get vector norm squared
v.normalized();                 // Get normalized copy of vector
v.normalize();                  // Normalize vector

v.dot_product(v2);              // Calculate dot product v * v2
v.cross_product(v2);            // Calculate cross product v x v2 (3x1 vectors only)
v.cross_product_normalized(v2); // Get cross product as normalized vector (3x1 vectors only)

v.angle_between_vectors(v2);    // Get angle between v and v2 in rad (2x1 or 3x1 vectors only)
v.rotation_axis(v2);            // Get normalized rotation axis to rotate v to v2 (3x1 vectors only)

std::cout << m << std::endl;  // Print using << operator

m.print();                    // Same as above