Dear Andreas, Mohsen, Thank you for the pointers. I will try the scikits approach as you suggested.
Andreas, PyCUDA looks like a great project. I saw one of your talks online Andreas and it sparked my interest. At the moment the performance of my code is not really an issue, but I anticipate I will want to spend a bit of time looking into optimisation at some point in the near future. I was just curious if significant speed-up was possible with relatively little changes to the existing code. It seems this may be possible with scikits. Although with my specific problem being sparse and tri-diagonal probably better performance can be had using PyCUDA with CUSP or cuSPARSE wrappers as you suggested. So right now I don't have any specific questions regarding the SparseSolve code, when I revisit this and have invested a bit of time learning CUDA I probably will. Thank you for the help. Best wishes, Dan On 12 March 2013 00:56, Andreas Kloeckner <[email protected]> wrote: > Dear Dan, > > "[email protected]" <[email protected]> writes: > > Is it possible to solve a linear system of the form Ax=d using PyCUDA? > > The actual system I wish to solve is sparse and tridiagonal. > > As long as it's positive definite and you don't mind installing pymetis > and Cython, the answer is yes. > > That said, writing a CUSP wrapper is likely a saner approach than using > this code. It's undocumented for a reason. :) > > > I noticed the sparse module in the Docs, > > http://wiki.tiker.net/PyCuda/Examples/SparseSolve, but I could > > understand the code (I'm very familiar with Python, but not with > > CUDA). Is there an example with comments or a tutorial would be even > > better. At the moment my code is using scipy and numpy. > > What about it is not making sense? If you've got a specific question, > I'd be happy to try and answer. > > Andreas >
_______________________________________________ PyCUDA mailing list [email protected] http://lists.tiker.net/listinfo/pycuda
