From: Jens Lohne Eftang <jleftang@gm...>  20121120 22:32:08

All, According to this thread http://sourceforge.net/mailarchive/message.php?msg_id=24029850 in the devel list there were plans on changing the API for cholesky_solve and lu_solve such that the interface with two arguments did not change the arguments, and another interface with only one argument would. There seems to be only one implementation if cholesky_solve and lu_solve in the current libmesh version, which takes two arguments. Does it change the rhs input vector? (It seems to me that it doesn't but I just wanted to make sure). Also it seems that cholesky_solve() will change my matrix (due to the factorization I assume). Is it correct that subsequent solves will give the correct answer, but that if I need the matrix for matrix multiplication I will have to keep a copy? Thanks! Jens 