MatrixSquareSymmetric.cpp File Reference

#include "geoslib_old_f.h"
#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)]

Macro Definition Documentation

#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 AL	(		i,
			j
	)		al[SQ(i,j,neq)-TRI(j)] /* for i >= j */

#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 TL	(		i,
			j
	)		tl[SQ(i,j,neq)-TRI(j)] /* for i >= j */

#define TRI

(

i

)

(((i) * ((i) + 1)) / 2)

#define XS	(		i,
			j
	)		x[SQ(i,j,neq)]