triangular_russian.h File Reference

TRSM and TRTRI via Gray code tables. More...

#include <m4ri/mzd.h>

Go to the source code of this file.

Functions
void	_mzd_trsm_lower_left_russian (mzd_t const *L, mzd_t *B, int k)
	Solves L X = B with X and B matrices and L lower triangular using Gray code tables. More...
void	_mzd_trsm_upper_left_russian (mzd_t const *U, mzd_t *B, int k)
	Solves U X = B with X and B matrices and U upper triangular using Gray code tables. More...
mzd_t *	mzd_trtri_upper_russian (mzd_t *A, int k)
	Invert the upper triangular matrix A using Kronrod's method. More...

Detailed Description

TRSM and TRTRI via Gray code tables.

Author

Martin Albrecht marti.nosp@m.nral.nosp@m.brech.nosp@m.t@go.nosp@m.oglem.nosp@m.ail..nosp@m.com

Function Documentation

void _mzd_trsm_lower_left_russian	(	mzd_t const *	L,
		mzd_t *	B,
		int	k
	)

Solves L X = B with X and B matrices and L lower triangular using Gray code tables.

X is stored inplace on B.

Parameters

L	Input lower triangular matrix.
B	Input matrix, being overwritten by the solution matrix X
k	Size of Gray code tables or zero for automatic choice (recommended).

void _mzd_trsm_upper_left_russian	(	mzd_t const *	U,
		mzd_t *	B,
		int	k
	)

Solves U X = B with X and B matrices and U upper triangular using Gray code tables.

X is stored inplace on B.

Parameters

U	Input upper triangular matrix.
B	Input matrix, being overwritten by the solution matrix X
k	Size of Gray code tables or zero for automatic choice (recommended).

mzd_t* mzd_trtri_upper_russian	(	mzd_t *	A,
		int	k
	)

Invert the upper triangular matrix A using Kronrod's method.

Parameters

A	Matrix to be inverted (overwritten).
k	Table size parameter, may be 0 for automatic choice.

Returns

Inverse of A or throws an error

dummy offsets table for _mzd_ple_to_e

deal with the rest