[petsc-users] DMPlexMatSetClosure from Fortran

[Dharmendar Reddy]

Hello,
         I had my code working earlier, but after an update to
 Petsc Development GIT revision: b45093d79c132b908a729253318f8991
55602205  GIT Date: 2013-11-11 16:19:16 -0600

I am seeing errors related to some DMPlex routines. I fixed some based
on earlier queries. Looks like i have found one more.

I call DMPlexMatSetClosure by passing petsc_null_object for default
and global sections. I do the same for DMPlexVecSetClosure.  The
residual part of the code which uses vecsetclosure works fine. But the
Matsetclosure gives an error:

due to finction call:
    call DMPlexMatSetClosure(dm, PETSC_NULL_OBJECT, PETSC_NULL_OBJECT,Jac, &
         cellId,this%Kelm_Array, ADD_VALUES,ierr)

 [0]PETSC ERROR: --------------------- Error Message
----------------------------
--------
[0]PETSC ERROR: Null argument, when expecting valid pointer!
[0]PETSC ERROR: Null Object: Parameter # 2!
[0]PETSC ERROR: ----------------------------------------------------------------
--------
[0]PETSC ERROR: Petsc Development GIT revision: b45093d79c132b908a729253318f8991
55602205  GIT Date: 2013-11-11 16:19:16 -0600
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: ----------------------------------------------------------------
--------
[0]PETSC ERROR: PoisTest on a mpi_rScalar_Debug named login4.stampede.tacc.utexa
s.edu by Reddy135 Mon Jan 13 02:54:22 2014
[0]PETSC ERROR: Libraries linked from /home1/00924/Reddy135/LocalApps/petsc/mpi_
rScalar_Debug/lib
[0]PETSC ERROR: Configure run at Tue Nov 12 17:14:19 2013
[0]PETSC ERROR: Configure options --with-debugging=1 --with-petsc-arch=mpi_rScal
ar_Debug --with-petsc-dir=/home1/00924/Reddy135/LocalApps/petsc --with-mpi-dir=/
opt/apps/intel13/mvapich2/1.9/ --with-scalar-type=real --with-blas-lapack-dir=/o
pt/apps/intel/13/composer_xe_2013.2.146/mkl/lib/intel64/ --with-clanguage=c --do
wnload-boost=1 --download-umfpack=1 --download-superlu_dist=1 --download-metis=1
 --download-parmetis=1 --download-mumps=1 --download-scalapack=1 --download-blac
s=1 --download-triangle=1 --download-ctetgen=1
[0]PETSC ERROR: ----------------------------------------------------------------
--------
[0]PETSC ERROR: DMPlexMatSetClosure() line 5640 in src/dm/impls/plex/plex.c

Breakpoint 2.1, FEMLIB_M::computejacobian_vp (this=0x230c5b0, dm=36335024, x=464
92464, jac=46384176, pcmat=46384176, flag=0, ierr=85) at /home1/00924/Reddy135/p
rojects/utgds/src/fem/VariationalProblemBoundProcedures.F90:779
779         call DMPlexMatSetClosure(dm, PETSC_NULL_OBJECT, PETSC_NULL_OBJECT,Ja
c, &

How do i fix this ?

Reddy
-- 
-----------------------------------------------------
Dharmendar Reddy Palle
Graduate Student
Microelectronics Research center,
University of Texas at Austin,
10100 Burnet Road, Bldg. 160
MER 2.608F, TX 78758-4445
e-mail: dharmareddy84 at gmail.com
Phone: +1-512-350-9082
United States of America.
Homepage: https://webspace.utexas.edu/~dpr342