<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Oh, there is also an issue I have recently noticed and did not have yet the time to fix it<div class=""><br class=""><div class="">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</div><div class="">Matt, take a look at petscsytypes.h and you will see the issue <a href="https://gitlab.com/petsc/petsc/-/blob/master/include/petscsystypes.h#L208" class="">https://gitlab.com/petsc/petsc/-/blob/master/include/petscsystypes.h#L208</a></div><div class=""><br class=""></div><div class="">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)</div><div class="">I believe the issue arise even if you call VecSet(v,0) on a VECCUDA<br class=""><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">On May 21, 2020, at 6:21 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div dir="ltr" class="">On Thu, May 21, 2020 at 10:53 AM Rui Silva <<a href="mailto:rui.silva@uam.es" class="">rui.silva@uam.es</a>> wrote:<br class=""></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello everyone,<br class="">
<br class="">
I am trying to run PETSc (with complex numbers in the GPU). When I call <br class="">
the VecWAXPY routine using the complex version of PETSc and mpicuda <br class="">
vectors, the program fails with a segmentation fault. This problem does <br class="">
not appear, if I run the complex version with mpi vectors or with the <br class="">
real version using mpicuda vectors. Is there any problem using <br class="">
CUDA+complex PETSc?<br class="">
<br class="">
Furthermore, I use the -log_view option to run the complex+gpu code, <br class="">
otherwise the program fails at the beggining.<br class=""></blockquote><div class=""><br class=""></div><div class="">What version of CUDA do you have? There are bugs in the versions before 10.2.</div><div class=""><br class=""></div><div class=""> Thanks,</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Best regards,<br class="">
<br class="">
Rui Silva<br class="">
<br class="">
-- <br class="">
Dr. Rui Emanuel Ferreira da Silva<br class="">
Departamento de Física Teórica de la Materia Condensada<br class="">
Universidad Autónoma de Madrid, Spain<br class="">
<a href="https://ruiefdasilva.wixsite.com/ruiefdasilva" rel="noreferrer" target="_blank" class="">https://ruiefdasilva.wixsite.com/ruiefdasilva</a><br class="">
<a href="https://mmuscles.eu/" rel="noreferrer" target="_blank" class="">https://mmuscles.eu/</a><br class="">
<br class="">
</blockquote></div><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""><div dir="ltr" class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">-- Norbert Wiener</div><div class=""><br class=""></div><div class=""><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a><br class=""></div></div></div></div></div></div></div></div>
</div></blockquote></div><br class=""></div></div></div></body></html>