matrix inversion BUG FIX
Status: Alpha
Brought to you by:
tramm
The file /efis/src/matrix.c has a bug. It performs the LU
decomposition right. invL is also fine but invU is wrong.
I implemented and optimized the calculation of invL
and invU without any function calls and it works fine with
the test I have done (12x12 and 13x13 matrix
inversions)
I also added some code from the /sim/src/mat/old-
matlib.cpp for the n=1 and n=2 cases.
I hope it will be useful for you.
I sent this e-mail to Mr. Hudson, but I had no response
and i don't know if he received it.
Santiago Fuster
Universidad Politecnica de Valencia
/efis/src/matrix.c