<div dir="ltr">I'm using CUDA 10.1.105 and this looks like a logic bug and not a CUDA problem.<div><br></div><div>I added two MatSetType calls that I thought were the right thing to do. Your error is a type problem and test with DMPlex so the matrix creation goes through a different code path.</div><div><br></div><div>You could try removing them and see if it works.</div><div><br></div><div>That is really all there is in this branch at this point other than some ViennaCL stuff.</div><div><br></div><div>I would verify that a fresh pull from master works. I would also reconfigure. I have just pulled from master and pushed so you can do a fresh pull.</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jul 6, 2019 at 4:07 PM Xiangdong <<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Mark,</div><div><br></div>Yes, I am using your branch (mark/gamg-fix-viennacl-rebased) with Petsc Development GIT revision: v3.11.3-1064-g8b74309384. I compiled PETSc with CUDA 10.1.<div><br></div><div>Which version of CUDA are you using?</div><div><br></div><div>Thanks.</div><div><br></div><div>Xiangdong</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jul 6, 2019 at 1:43 PM Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I am not able to reproduce this error. THis is what I get. Are you pulling from mark/gamg-fix-viennacl-rebased ?<div><br></div><div>13:38  /gpfs/alpine/scratch/adams/geo127$ jsrun -n 1 -c 4 -a 4 -g 1 ./ex19 -da_refine 5 -snes_view -snes_monitor -ksp_monitor -mat_type aijcusparse -vec_type cuda -log_view<br>lid velocity = 0.000106281, prandtl # = 1., grashof # = 1.<br>  0 SNES Function norm 1.036007954337e-02 <br>    0 KSP Residual norm 9.144944502871e-02 <br>    1 KSP Residual norm 2.593759906204e-02 <br>    2 KSP Residual norm 1.669815200495e-02 <br>    3 KSP Residual norm 1.510777951698e-02 <br>    4 KSP Residual norm 1.458401237884e-02 <br>    5 KSP Residual norm 1.418635322926e-02 <br>    6 KSP Residual norm 1.377436725003e-02 <br>    7 KSP Residual norm 1.332236907186e-02 <br>    8 KSP Residual norm 1.288602527920e-02 <br>    9 KSP Residual norm 1.240018288138e-02 <br>   10 KSP Residual norm 1.186798872492e-02 <br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 5, 2019 at 11:07 PM Xiangdong <<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Mark,<div><br></div><div>After pulling your new commit, I got an error when running ex19 (even without nvprof). Same error with version 3.11.3. But ex19 run fine with old commit 64edb9194a9.</div><div><br></div><div>Below is a run time log error file. By the way, which version of cuda does PETSc support?</div><div><br></div><div>Thank you.</div><div><br></div><div>Xiangdong</div><div><br></div><div>lid velocity = 0.000106281, prandtl # = 1., grashof # = 1.<br>  0 SNES Function norm 1.036007954337e-02 <br>    0 KSP Residual norm 9.158537890035e-02 <br>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: Invalid argument<br>[0]PETSC ERROR: Object (seq) is not seqcuda or mpicuda<br>[0]PETSC ERROR: See <a href="https://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">https://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Development GIT revision: v3.11.3-1064-g8b74309384  GIT Date: 2019-07-03 10:24:59 -0400<br>[0]PETSC ERROR: ./extest on a arch-debug2 named new by eps Fri Jul  5 22:57:02 2019<br>[0]PETSC ERROR: Configure options PETSC_ARCH=arch-debug2 --with-debugging=1 --with-mpi-dir=/home/eps/MyLocal/mpi/mpich --with-blaslapack-dir=/home/eps/MyLocal/intel/mkl --with-cuda=1<br>[0]PETSC ERROR: #1 VecCUDAGetArrayRead() line 1304 in /home/eps/MyLocal/petsc-mark/src/vec/vec/impls/seq/seqcuda/<a href="http://veccuda2.cu" target="_blank">veccuda2.cu</a><br>[0]PETSC ERROR: #2 MatMultAdd_SeqAIJCUSPARSE() line 1415 in /home/eps/MyLocal/petsc-mark/src/mat/impls/aij/seq/seqcusparse/<a href="http://aijcusparse.cu" target="_blank">aijcusparse.cu</a><br>[0]PETSC ERROR: #3 MatMult_SeqAIJCUSPARSE() line 1347 in /home/eps/MyLocal/petsc-mark/src/mat/impls/aij/seq/seqcusparse/<a href="http://aijcusparse.cu" target="_blank">aijcusparse.cu</a><br>[0]PETSC ERROR: #4 MatMult() line 2368 in /home/eps/MyLocal/petsc-mark/src/mat/interface/matrix.c<br>[0]PETSC ERROR: #5 PCApplyBAorAB() line 662 in /home/eps/MyLocal/petsc-mark/src/ksp/pc/interface/precon.c<br>[0]PETSC ERROR: #6 KSP_PCApplyBAorAB() line 309 in /home/eps/MyLocal/petsc-mark/include/petsc/private/kspimpl.h<br>[0]PETSC ERROR: #7 KSPGMRESCycle() line 152 in /home/eps/MyLocal/petsc-mark/src/ksp/ksp/impls/gmres/gmres.c<br>[0]PETSC ERROR: #8 KSPSolve_GMRES() line 237 in /home/eps/MyLocal/petsc-mark/src/ksp/ksp/impls/gmres/gmres.c<br>[0]PETSC ERROR: #9 KSPSolve() line 764 in /home/eps/MyLocal/petsc-mark/src/ksp/ksp/interface/itfunc.c<br>[0]PETSC ERROR: #10 SNESSolve_NEWTONLS() line 225 in /home/eps/MyLocal/petsc-mark/src/snes/impls/ls/ls.c<br>[0]PETSC ERROR: #11 SNESSolve() line 4430 in /home/eps/MyLocal/petsc-mark/src/snes/interface/snes.c<br>[0]PETSC ERROR: #12 main() line 161 in /home/eps/MyCodes/petsctest/extest.c<br>[0]PETSC ERROR: PETSc Option Table entries:<br>[0]PETSC ERROR: -da_refine 5<br>[0]PETSC ERROR: -ksp_monitor<br>[0]PETSC ERROR: -mat_type aijcusparse<br>[0]PETSC ERROR: -snes_monitor<br>[0]PETSC ERROR: -snes_view<br>[0]PETSC ERROR: -vec_type cuda<br>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------<br>application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jul 3, 2019 at 10:27 AM Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">My branch mark/gamg-fix-viennacl-rebased fixes this problem for me. Please give it a try.<div>Thanks,</div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 2, 2019 at 10:41 PM Xiangdong via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello everyone,<div><br></div><div>When I run the ex19 with cuda like this:</div><div> mpiexec -np 4 ./ex19 -da_refine 5 -snes_view -snes_monitor -ksp_monitor -mat_type aijcusparse -vec_type cuda  -log_view</div><div><br></div><div>it worked fine and produced correct results.</div><div><br></div><div>However, when I tried to run this example with nvprof:</div><div>nvprof --profile-child-processes mpiexec -np 4 ./extest -da_refine 5 -snes_view -snes_monitor -ksp_monitor -mat_type aijcusparse -vec_type cuda  -log_view<br></div><div><br></div><div>I got errors like:</div><div>[3]PETSC ERROR: Error in external library<br>[3]PETSC ERROR: CUDA error 700<br>[3]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>[3]PETSC ERROR: Petsc Release Version 3.11.2, unknown <br>[3]PETSC ERROR: ./ex19 on a arch-opt named new by eps Tue Jul  2 22:26:01 2019<br>[3]PETSC ERROR: Configure options PETSC_ARCH=arch-opt --with-debugging=0 --with-mpi-dir=/home/eps/MyLocal/mpi/mpich --with-blaslapack-dir=/home/eps/MyLocal/intel/mkl --with-cuda-dir=/home/eps/MyLocal/cuda --with-hypre-dir=/home/eps/MyLocal/hypre-2.15.1/hypre-install --download-hdf5=1<br>[3]PETSC ERROR: #1 VecSet_SeqCUDA() line 785 in /home/eps/MyLocal/petsc/src/vec/vec/impls/seq/seqcuda/<a href="http://veccuda2.cu" target="_blank">veccuda2.cu</a><br>[3]PETSC ERROR: #2 VecSet() line 547 in /home/eps/MyLocal/petsc/src/vec/vec/interface/rvector.c<br>[3]PETSC ERROR: #3 VecCreate_MPICUDA() line 178 in /home/eps/MyLocal/petsc/src/vec/vec/impls/mpi/mpicuda/<a href="http://mpicuda.cu" target="_blank">mpicuda.cu</a><br>[3]PETSC ERROR: #4 VecSetType() line 51 in /home/eps/MyLocal/petsc/src/vec/vec/interface/vecreg.c<br>[3]PETSC ERROR: #5 VecCreate_CUDA() line 192 in /home/eps/MyLocal/petsc/src/vec/vec/impls/mpi/mpicuda/<a href="http://mpicuda.cu" target="_blank">mpicuda.cu</a><br>[3]PETSC ERROR: #6 VecSetType() line 51 in /home/eps/MyLocal/petsc/src/vec/vec/interface/vecreg.c<br>[3]PETSC ERROR: #7 MatCreateVecs() line 8996 in /home/eps/MyLocal/petsc/src/mat/interface/matrix.c<br>[3]PETSC ERROR: #8 MatFDColoringCreate() line 482 in /home/eps/MyLocal/petsc/src/mat/matfd/fdmatrix.c<br>[3]PETSC ERROR: #9 SNESComputeJacobian_DMDA() line 175 in /home/eps/MyLocal/petsc/src/snes/utils/dmdasnes.c<br>[3]PETSC ERROR: #10 SNESComputeJacobian() line 2718 in /home/eps/MyLocal/petsc/src/snes/interface/snes.c<br>[3]PETSC ERROR: #11 SNESSolve_NEWTONLS() line 222 in /home/eps/MyLocal/petsc/src/snes/impls/ls/ls.c<br>[3]PETSC ERROR: #12 SNESSolve() line 4560 in /home/eps/MyLocal/petsc/src/snes/interface/snes.c<br>[3]PETSC ERROR: #13 main() line 161 in /home/eps/MyCodes/petsctest/extest.c<br></div><div><br></div><div>The full run log is attached. </div><div><br></div><div>I am using  NVIDIA-SMI 418.56       Driver Version: 418.56       CUDA Version: 10.1.</div><div><br></div><div>I do not know why it is okay without nvprof but crashed with nvprof. Any suggestion to fix this?</div><div><br></div><div>Thank you.</div><div><br></div><div>Best,</div><div>Xiangdong</div><br></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>