[petsc-users] Possible bug PETSc+Complex+CUDA

[Stefano Zampini]

Oh, there is also an issue I have recently noticed and did not have yet the time to fix it

With complex numbers, we use the definitions for complexes from thrust and this does not seem to be always compatible to whatever the C compiler uses
Matt, take a look at petscsytypes.h and you will see the issue https://gitlab.com/petsc/petsc/-/blob/master/include/petscsystypes.h#L208 <https://gitlab.com/petsc/petsc/-/blob/master/include/petscsystypes.h#L208>

For sure, you need to configure petsc with --with-clanguage=cxx, but even that does not to seem make it work on a CUDA box I have recently tried out (CUDA 10.1)
I believe the issue arise even if you call VecSet(v,0) on a VECCUDA

> On May 21, 2020, at 6:21 PM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Thu, May 21, 2020 at 10:53 AM Rui Silva <rui.silva at uam.es <mailto:rui.silva at uam.es>> wrote:
> Hello everyone,
> 
> I am trying to run PETSc (with complex numbers in the GPU). When I call 
> the VecWAXPY routine using the complex version of PETSc and mpicuda 
> vectors, the program fails with a segmentation fault. This problem does 
> not appear, if I run the complex version with mpi vectors or with the 
> real version using mpicuda vectors. Is there any problem using 
> CUDA+complex PETSc?
> 
> Furthermore, I use the -log_view option to run the complex+gpu code, 
> otherwise the program fails at the beggining.
> 
> What version of CUDA do you have? There are bugs in the versions before 10.2.
> 
>   Thanks,
> 
>     Matt
>  
> Best regards,
> 
> Rui Silva
> 
> -- 
>   Dr. Rui Emanuel Ferreira da Silva
>   Departamento de Física Teórica de la Materia Condensada
>   Universidad Autónoma de Madrid, Spain
>   https://ruiefdasilva.wixsite.com/ruiefdasilva <https://ruiefdasilva.wixsite.com/ruiefdasilva>
>   https://mmuscles.eu/ <https://mmuscles.eu/>
> 
> 
> 
> -- 
> 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
> 
> https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200521/78c64a67/attachment.html>