[petsc-users] Problem with GAMG elasticity case:'createdefaultdata' not set(?) need to support NULL data!

[Einar Sørheim]

Hi
I am testing out the petsc ksp solvers in our own FEM code, in particluar
we are interseted in the GAMG preconditioner.
Applying it to our mechanical problem we get some error messages from
petsc, which I suspect is due to some necessary init calls that are
missing. First the error message:
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Petsc has generated inconsistent data!
[0]PETSC ERROR: 'createdefaultdata' not set(?) need to support NULL data!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013
[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: D:\Programming\gits\Alsim\SharedModules -
Copy\Source\Release\Alsim_nypetsc.exe on a arch-mswin-c-opt named CMP3 by
admeinar Wed Oct 16 18:23:32 2013
[0]PETSC ERROR: Libraries linked from
/cygdrive/d/Programming/petsc-3.4.3/arch-mswin-c-opt/lib
[0]PETSC ERROR: Configure run at Wed Oct 16 14:46:01 2013
[0]PETSC ERROR: Configure options --with-debugging=0 --with-openmp=yes
--with-pthread=no --with-cc="win32fe icl" --with-fc="win32fe ifort"
--with-cxx="win32fe icl" --with-blas-lapack-dir="C:/Program Files
(x86)/Intel/Composer XE 2013/mkl/lib/intel64" --download-superlu
--with-sowing=0 --with-c2html=0 --with-mpi-dir="C:/Program Files
(x86)/Intel/MPI/4.1.0.028/em64t" --useThreads=0 --useThreads=0
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: PCSetUp_GAMG() line 609 in
src/ksp/pc/impls/gamg/D:\PROGRA~1\PETSC-~1.3\src\ksp\pc\impls\gamg\gamg.c
[0]PETSC ERROR: PCSetUp() line 890 in
src/ksp/pc/interface/D:\PROGRA~1\PETSC-~1.3\src\ksp\pc\INTERF~1\precon.c
[0]PETSC ERROR: KSPSetUp() line 278 in
src/ksp/ksp/interface/D:\PROGRA~1\PETSC-~1.3\src\ksp\ksp\INTERF~1\itfunc.c

In short our code looks like the following:


      call KSPGetPC(ksp,pc,ierr)
      tol = 1.e-20
      call
KSPSetTolerances(ksp,tol,PETSC_DEFAULT_DOUBLE_PRECISION,PETSC_DEFAULT_DOUBLE_PRECISION,imaxiter,ierr)
      select case (i_Precond)

      case ( 3 )
         call PCSetType(pc,PCGAMG,ierr)
!         call PCGAMGInitializePackage()
!         call PCCreate_GAMG(pc,ierr)
      case DEFAULT
         call PCSetType(pc,PCJACOBI,ierr)
      end select

!      call PCSetType(pc,PCJACOBI,ierr)
      select case (i_solver)
      case ( 0 )
         call KSPSetType(ksp,KSPCG,ierr)
      case DEFAULT
         call KSPSetType(ksp,KSPCG,ierr)
      end select
      call KSPCGSetType(ksp,KSP_CG_SYMMETRIC,ierr)
      if (i_nodedof==3) then
         write(*,*) "Set coord, number of nodes is:", i_nodes/i_nodedof
         call
VecCreateSeqWithArray(MPI_COMM_SELF,3,i_nodes,coords_,vec_coords,ierr)
         call MatNullSpaceCreateRigidBody(vec_coords,Matnull,ierr)
         call MatSetNearNullSpace(A,Matnull,ierr)
         call MatNullSpaceDestroy(Matnull,ierr)
         call PetscViewerASCIIOpen(PETSC_COMM_WORLD, "Coordvec.m",
viewer,ierr)
         call PetscViewerSetFormat( viewer, PETSC_VIEWER_ASCII_MATLAB,ierr)
         call VecView(vec_coords,viewer,ierr)
         call PetscViewerDestroy( viewer,ierr )
         call VecDestroy(vec_coords,ierr)
         write(*,*) "Finish setting null space ierr =", ierr
!         call PCSetCoordinates( pc, 3,i_nodes/i_nodedof, coords_, ierr )
!         CHKERRQ(ierr)
      end if


!      call KSPSetFromOptions(ksp,ierr)
      call KSPSetInitialGuessNonzero(ksp,PETSC_TRUE,ierr)
      call KSPSetUp(ksp,ierr)
      call PetscTime(t0,ierr)
      write(*,*) "Time setup", t0-t1
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
!                      Solve the linear system
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

      call KSPSolve(ksp,b,x,ierr)
      call PetscTime(t1,ierr)

-- 
Einar Sørheim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131016/2180d993/attachment.html>