[GRASS-dev] Re: v.surf.rst

[Abhishek Shukla]

Hi list,
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.
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.


Please someone confirm that it is okay to solve matrix using Cholesky
decomposition.


---
Abhishek Shukla
Center for Security, Theory And Algorithm Research (CSTAR)
International Institute of Information Technology
Hyderabad
India
http://www.iiit.ac.in/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-dev/attachments/20100503/28b6ac4a/attachment.html