From f6b8e0eae4374f339487a33e3e4fe5462d5816e1 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?J=C3=B6rg=20Frings-F=C3=BCrst?= Date: Sat, 25 Nov 2017 10:16:00 +0100 Subject: New upstream version 2.0.0 --- numlib/ludecomp.h | 14 ++++++++++++++ 1 file changed, 14 insertions(+) mode change 100644 => 100755 numlib/ludecomp.h (limited to 'numlib/ludecomp.h') diff --git a/numlib/ludecomp.h b/numlib/ludecomp.h old mode 100644 new mode 100755 index cdb4382..1075a5b --- a/numlib/ludecomp.h +++ b/numlib/ludecomp.h @@ -69,6 +69,8 @@ int *pivx /* Pivoting row permutations record */ ); /* Invert a matrix A using lu decomposition */ +/* NOTE that it returns transposed inverse by normal convention. */ +/* Use sym_matrix_trans() to fix this, or use matrix_trans_mult() */ /* Return 1 if the matrix is singular, 0 if OK */ int lu_invert( @@ -76,7 +78,19 @@ double **a, /* A[][] input matrix, returns inversion of A */ int n /* Dimensionality */ ); +/* Invert a matrix A using lu decomposition, and polish it. */ +/* NOTE that it returns transposed inverse by normal convention. */ +/* Use sym_matrix_trans() to fix this, or use matrix_trans_mult() */ +/* Return 1 if the matrix is singular, 0 if OK */ +int +lu_polished_invert( +double **a, /* A[][] input matrix, returns inversion of A */ +int n /* Dimensionality */ +); + /* Pseudo-Invert matrix A using lu decomposition */ +/* NOTE that it returns transposed inverse by normal convention. */ +/* Use matrix_trans() to fix this, or use matrix_trans_mult(). */ /* Return 1 if the matrix is singular, 0 if OK */ int lu_psinvert( -- cgit v1.2.3