On Sun, Oct 26, 2014 at 8:13 PM, Ashwin Srinath <[email protected]> wrote:

Matthew, I think this fixes the scalar issue:

https://github.com/ashwinsrnth/petsc-pycuda/commit/145b5f138ecbdddf3c8a2ab4c...

Nice, you can do err checking using petsc4py as well: def getType(self): cdef const_char *cval = NULL CHKERR( KSPGetType(self.ksp, &cval) ) return bytes2str(cval) Thanks, Matt

On Sun, Oct 26, 2014 at 9:05 PM, Matthew Knepley <[email protected]> wrote:

...

On Sun, Oct 26, 2014 at 7:20 PM, Ashwin Srinath <[email protected]> wrote:

...

Hi,

I re-wrote the code to include the code from this branch - it's made everything a lot simpler: https://github.com/ashwinsrnth/petsc-pycuda/blob/master/GPUArray.pyx

Some notes:

1) You can check the return type of the PETSc calls for errors

2) Lisandro, do you have an automatic converter from PetscScalar to a numpy dtype? Using np.dtype.float64 is dicey.

Thanks,

Matt

...

If I can help with adding the functionality to `petsc4py`, I'd be glad to :)

Thanks, Ashwin

On Fri, Oct 24, 2014 at 2:43 PM, Lisandro Dalcin <[email protected]> wrote:

...

On 24 October 2014 03:56, Karl Rupp <[email protected]> wrote:

...

Both 1) and 2) are now pushed to the feature branch. Any feedback appreciated :-)

Thanks! Looks good. Ashwin, could you please re-write your code to test this branch? Once you confirm everything's working, we can follow discussing how to better support this in petsc4py.

-- Lisandro Dalcin ============ Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Numerical Porous Media Center (NumPor) King Abdullah University of Science and Technology (KAUST) http://numpor.kaust.edu.sa/

4700 King Abdullah University of Science and Technology al-Khawarizmi Bldg (Bldg 1), Office # 4332 Thuwal 23955-6900, Kingdom of Saudi Arabia http://www.kaust.edu.sa

Office Phone: +966 12 808-0459

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener