[petsc-users] DMCreateGlobalVector in fortran.

[Manuel Valera]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170204/1a2470ff/attachment.html>