Hello devs, Ive been creating the framework to distribute the arrays from my model using the global/local vectors and the DMDAs, for this i created a small subroutine that is responsible for making the necessary calls:
SUBROUTINE DAs(da,globalv,localv) ! use PetscObjects, only :: ierr ! Umbrella program to update and communicate the arrays in a ! distributed fashion using the DMDA objects from PETSc. ! Manuel Valera 1/20/17 ! Arguments: ! da = DMDA array either 1d,2d or 3d, already created and setup ! globalv = global vector to be operated into ! localv = local chunk each processor works in. Vec,intent(inout) :: globalv,localv DM,intent(inout) :: da call DMCreateGlobalVector(da,globalv,ierr) call DMCreateLocalVector(da,localv,ierr) call DMGlobalToLocalBegin(da,globalv,INSERT_VALUES,localv,ierr) call DMGlobalToLocalEnd(da,globalv,INSERT_VALUES,localv,ierr) END SUBROUTINE DAs I call this from another program where i have declared a vector as: #include "petsc/finclude/petscdm.h"#include "petsc/finclude/petscdmda.h"use petscuse petscdmuse petscdmda IMPLICIT NONE Vec :: xpdmda,yp,zp,xcp,ycp,zcp DM :: da,dac,dau,dav,daw CONTAINSSUBROUTINE InitGridPIMPLICIT NONE!PetscInt :: ierr! CREATE DMDA:call CDA3(da,IMax-1,JMax-1,KMax-1)!xpdmda,yp,zp have the same dimensions than x,y,z!xpl,ypl,zpl are the local versions of xp,yp,zp!We create the global and local vectors and communicate their values:call DAs(da,xpdmda,xpl,x) And at this call it crashes with the following message: [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Null argument, when expecting valid pointer [0]PETSC ERROR: Null Object: Parameter # 2 [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. [0]PETSC ERROR: Petsc Development GIT revision: v3.7.5-2949-gc3c607c GIT Date: 2017-01-20 17:34:16 -0600 [0]PETSC ERROR: ./ucmsLEP on a arch-linux2-c-debug named ocean by valera Sat Feb 4 15:53:53 2017 [0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-fblaslapack --download-mpich --download-hdf5 --download-netcdf --download-hypre --download-metis --download-parmetis --download-trillinos --with-debugging=1 [0]PETSC ERROR: #3 VecSetLocalToGlobalMapping() line 79 in /usr/dataC/home/valera/petsc/src/vec/vec/interface/vector.c [0]PETSC ERROR: #4 DMCreateGlobalVector_DA() line 41 in /usr/dataC/home/valera/petsc/src/dm/impls/da/dadist.c [0]PETSC ERROR: #5 DMCreateGlobalVector() line 840 in /usr/dataC/home/valera/petsc/src/dm/interface/dm.c So my question is, do i have to setup the vectors in any additional way? i followed the fortran DMDA examples closely and i didnt notice anything extra, Thanks, Manuel
