[petsc-users] DMCreateGlobalVector in fortran.

[Barry Smith]

   See:  http://www.mcs.anl.gov/petsc/documentation/changes/dev.html  search for DMDACreate  and update your code.

  Barry


> On Feb 4, 2017, at 6:02 PM, Manuel Valera <mvalera at mail.sdsu.edu> wrote:
> 
> 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 petsc
> use petscdm
> use petscdmda
> 
> IMPLICIT NONE
> 
> Vec                                        :: xpdmda,yp,zp,xcp,ycp,zcp
> DM                                         :: da,dac,dau,dav,daw
> 
> CONTAINS
> SUBROUTINE InitGridP
> 
> IMPLICIT 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