include file fortran

[Tahar Amari]

Hello,

I finally succeed building Petsc and the fortran  code thanks to your  
advice .

However when I run the  code (piece of it is below)

      	if(pe==0) write(*,*) "before MatCopy"
       CALL MatCopy(dvoeta_petsc,lhs_petsc,SAME_NONZERO_PATTERN,
      & ierror)
      	if(pe==0) write(*,*) "after MatCopy"


I got the following message (in between my printing for debugging) .

  before MatCopy

[0]PETSC ERROR: --------------------- Error Message  
------------------------------------
[0]PETSC ERROR: Corrupt argument: see http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Corrupt 
!
[0]PETSC ERROR: Invalid Pointer to Object: Parameter # 1!
[0]PETSC ERROR:  
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.0.0, Patch 5, Mon Apr 13  
09:15:37 CDT 2009
[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: mh4d on a macx named imac-de-tahar-amari.local by  
amari Mon May 18 13:47:42 2009
[0]PETSC ERROR: Libraries linked from /Data/Poub1/petsc-3.0.0-p5/macx/ 
lib
[0]PETSC ERROR: Configure run at Sat May 16 13:49:21 2009
[0]PETSC ERROR: Configure options -with-cc=gcc --download-mpich=1 -- 
download-f-bl

after MatCopy


I looked at the piece of code which define the first argument  
"dvoeta_petsc"

which is

          Mat, SAVE :: dvoeta_petsc
           CALL init_petsc_mat_vector(dvoeta_petsc,diagonal=dvoeta)


This code was working with Petsc version 2.xxx
I wonder if some changes to version 3 might  be the reason ?


May it be the call to

         CALL MatSetOption(matrix,MAT_NO_NEW_NONZERO_LOCATIONS,ierr)

in this function please ?

c     subprogram 8. init_petsc_mat_vector
c     Make a petsc matrix out of an operator
c-----------------------------------------------------------------------
       SUBROUTINE init_petsc_mat_vector(matrix,v2v_operator_vv,diagonal,
      &bctype)
c-----------------------------------------------------------------------
       USE operator_mod
#include "include/finclude/petsc.h"
#include "include/finclude/petscvec.h"
#include "include/finclude/petscmat.h"
#include "include/finclude/petscao.h"

c-----------------------------------------------------------------------
       Mat, INTENT(INOUT) :: matrix
       INTERFACE
         SUBROUTINE v2v_operator_vv(x)
           USE operator_mod
           TYPE (v2v_operators), INTENT(INOUT) :: x(:)
         END SUBROUTINE v2v_operator_vv
       END INTERFACE
       OPTIONAL :: v2v_operator_vv
       TYPE (vertex_scalar), INTENT(IN), OPTIONAL :: diagonal
       INTEGER(i4), INTENT(IN), OPTIONAL :: bctype
       INTEGER(i4) :: m,n,ierr,i,n1,m1,nv_local_inside,
      &  iv_local,iv_couple,d_nz,o_nz
       TYPE(v2v_operators),  ALLOCATABLE :: v2voperator(:)
       REAL(r8) :: v(3,3)
       LOGICAL, SAVE :: first=.true.
c-----------------------------------------------------------------------

       CALL get_grid(thegrid,nv_local_inside=nv_local_inside)
       CALL get_nonzero(thegrid,d_nz,o_nz)

       CALL MatCreateMPIAIJ(comm_grid,3*nv_local_inside,
      & 3*nv_local_inside,PETSC_DECIDE,PETSC_DECIDE,3*d_nz,
      & PETSC_NULL_INTEGER,3*o_nz,PETSC_NULL_INTEGER,matrix,ierr)

       ALLOCATE(v2voperator(nv_local_inside))
c-----------------------------------------------------------------------
c See if it is a diagonal operator or not.
c-----------------------------------------------------------------------
       IF (PRESENT(diagonal)) THEN
         IF(first) CALL terminator("Error in init_petsc_mat_vector: "//
      &  "You shouldn't call it the first time with a diagonal  
operator")
         IF(PRESENT(bctype)) THEN
           CALL v2v_diagonal_vv(v2voperator,diagonal,bctype)
         ELSE
           CALL v2v_diagonal_vv(v2voperator,diagonal,1)
         ENDIF
       ELSE
         CALL v2v_operator_vv(v2voperator)
       ENDIF

       DO iv_local=1,nv_local_inside
         CALL index_vv2petsc(m,iv_local,1)
         DO i=1,SIZE(v2voperator(iv_local)%neighbors)
           iv_couple=v2voperator(iv_local)%neighbors(i)
           CALL index_vv2petsc(n,iv_couple,1)
           v=v2voperator(iv_local)%coupling_with(i)%a
           DO m1=0,2
             DO n1=0,2
               CALL MatSetValues(matrix,1,m+m1,1,n+n1,v(1+m1,1+n1),
      &        INSERT_VALUES,ierr)
             ENDDO
           ENDDO
         ENDDO
       ENDDO

       CALL MatAssemblyBegin(matrix,MAT_FINAL_ASSEMBLY,ierr)
       CALL MatAssemblyEnd(matrix,MAT_FINAL_ASSEMBLY,ierr)

       IF(first) THEN
         CALL MatSetOption(matrix,MAT_NO_NEW_NONZERO_LOCATIONS,ierr)
         first=.false.
       ENDIF

       CALL destroy(v2voperator)
       DEALLOCATE(v2voperator)

       RETURN