El 13/02/2013, a las 15:17, Barry Smith escribió:

Thanks for noting this. Is the problem the (*right)->map operation in?

ierr = PetscLayoutReference(mat->cmap,&(*right)->map);CHKERRQ(ierr);

This can be corrected easily by replacing the PetscLayoutReference() here with VecSetLayout(). Does that resolve the problem?

General note: We should really move the PetscLayout operations out of Vec mpi impl location.

Barry

Yes it works. https://bitbucket.org/petsc/petsc-dev/commits/8264b0c102238be7a957e713db98c3... By the way, I am getting a warning which I think is new - when compiling with nvcc in optimized builds: /home/jroman/soft/petsc-dev/include/petsc-private/petscimpl.h(201): warning: "cc" clobber ignored Jose