#include "Matrix/MatrixSquareSymmetric.hpp"#include "Matrix/MatrixRectangular.hpp"#include "Matrix/AMatrixSquare.hpp"#include "Matrix/MatrixSquareGeneral.hpp"#include "Basic/VectorHelper.hpp"#include "Basic/AException.hpp"Macros | |
| #define | TRI(i) (((i) * ((i) + 1)) / 2) |
| #define | SQ(i, j, neq) ((j) * neq + (i)) |
| #define | AT(i, j) at[TRI(j)+(i)] /* for j >= i */ |
| #define | AL(i, j) al[SQ(i,j,neq)-TRI(j)] /* for i >= j */ |
| #define | BS(i, j) b[SQ(i,j,neq)] |
| #define | XS(i, j) x[SQ(i,j,neq)] |
| #define | AS(i, j) a[SQ(i,j,neq)] |
| #define | TL(i, j) tl[SQ(i,j,neq)-TRI(j)] /* for i >= j */ |
| #define | _TL(i, j) _tl[SQ(i,j,neq)-TRI(j)] /* for i >= j */ |
| #define | _XL(i, j) _xl[SQ(i,j,neq)-TRI(j)] /* for i >= j */ |
| #define | HA(i, j) ha[SQ(i,j,neq)] |
| #define AS | ( | i, | |
| j | |||
| ) | a[SQ(i,j,neq)] |
| #define AT | ( | i, | |
| j | |||
| ) | at[TRI(j)+(i)] /* for j >= i */ |
| #define BS | ( | i, | |
| j | |||
| ) | b[SQ(i,j,neq)] |
| #define HA | ( | i, | |
| j | |||
| ) | ha[SQ(i,j,neq)] |
| #define SQ | ( | i, | |
| j, | |||
| neq | |||
| ) | ((j) * neq + (i)) |
| #define TRI | ( | i | ) | (((i) * ((i) + 1)) / 2) |