[petsc-dev] Backend independent VecGetArray for GPUs

Matthew Knepley knepley at gmail.com
Sun Oct 26 20:17:11 CDT 2014


On Sun, Oct 26, 2014 at 8:13 PM, Ashwin Srinath <ashwinsrnth at gmail.com>
wrote:

> Matthew, I think this fixes the scalar issue:
>
>
> https://github.com/ashwinsrnth/petsc-pycuda/commit/145b5f138ecbdddf3c8a2ab4caa86e132ece347f
>

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 <knepley at gmail.com>
> wrote:
>
>> On Sun, Oct 26, 2014 at 7:20 PM, Ashwin Srinath <ashwinsrnth at gmail.com>
>> 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 <dalcinl at gmail.com>
>>> wrote:
>>>
>>>> On 24 October 2014 03:56, Karl Rupp <rupp at iue.tuwien.ac.at> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20141026/31c16de9/attachment.html>


More information about the petsc-dev mailing list