Thanks so much, Dominic. I'll let you know soon how this branch works for my code.
Ashwin On Fri, Jul 18, 2014 at 7:33 PM, Dominic Meiser <[email protected]> wrote: > Hi Ashwin, > > I've been working on this issue the last couple of days and wanted to give > you an update. First off, thanks for the test case. This has helped a lot > in tracking down the issues. I've created a couple of unit tests based on > your test case. > > There were several independent bugs that prevented the VecView to work > with -dm_vec_type cusp. I've fixed some of them in this branch: > > https://bitbucket.org/dmeiser/petsc/branch/fix-cusp-dmda > > Viewing a vector created from a DM should now be possible. For instance > the test case you provided should work with this branch. > > There are still some issues remaining. DMGlobalToLocalBegin/End and > DMLocalToGlobalBegin/End are not working. I haven't been able to track > these down yet. I'll return to this next week. > > Cheers, > Dominic > > > > On 06/24/2014 08:15 PM, Ashwin Srinath wrote: > > Hello, petsc-users > > I'm having trouble *viewing* an mpicusp vector. Here's the simplest case > that reproduces the problem: > > int main(int argc, char** argv) { > > PetscInitialize(&argc, &argv, NULL, NULL); > > DM da; > Vec V; > > DMDACreate2d( PETSC_COMM_WORLD, > DM_BOUNDARY_NONE, DM_BOUNDARY_NONE, > DMDA_STENCIL_BOX, > 5, 5, > PETSC_DECIDE, PETSC_DECIDE, > 1, > 1, > NULL, NULL, > &da); > > DMCreateGlobalVector(da, &V); > > VecSet(V, 1); > VecView(V, PETSC_VIEWER_STDOUT_WORLD); > > PetscFinalize(); > return 0; > } > > I get the error: > [1]PETSC ERROR: Null argument, when expecting valid pointer > [0]PETSC ERROR: Trying to copy from a null pointer > > I executed with the following command: > mpiexec -n 2 ./main -dm_vec_type cusp -vec_type cusp > Both GPUs are attached to two different processes. > > This program works fine for vecmpi vectors, i.e., -dm_vec_type mpi and > -vec_type mpi. Also, I don't get an error unless I try to *view* the > vector. Can someone please point out what I'm doing wrong? > > Thanks for your time, > Ashwin Srinath > > > > > > > > -- > Dominic Meiser > Tech-X Corporation > 5621 Arapahoe Avenue > Boulder, CO 80303 > USA > Telephone: 303-996-2036 > Fax: 303-448-7756www.txcorp.com > >
