Hi list,<br>I want to ask why LU decomposition has been used for matrix solving in RST library when matrix being solved is symmetric. I am planning to solve it using Cholesky decomposition. <br>I want to confirm with grass developers that there is no issue if I solve it using Cholesky decomposition. My goal is not to change current serial implementation, I am implementing  (started, at least) v.surf.rst on GPU and planning to use Cholesky instead of LU as it is expected to show better performance over LU.<br>

<br><div><br>Please someone confirm that it is okay to solve matrix using Cholesky decomposition.<br><br><br>---<br>Abhishek Shukla<br>Center for Security, Theory And Algorithm Research (CSTAR)<br>International Institute of Information Technology<br>

Hyderabad<br>India<br><a href="http://www.iiit.ac.in/">http://www.iiit.ac.in/</a></div>