[MOAB-dev] iMesh_setAdjTable() causes error in iMesh_load()

[Jed Brown]

On Thu, Jun 9, 2011 at 18:17, Iulian Grindeanu <iulian at mcs.anl.gov> wrote:

> Can you send me the portion of the code that
> sets the table
> loads the mesh
>
> I should be able to reproduce it here, and I will let you know
>
>
The table is set here:

https://github.com/jedbrown/dohp/blob/master/src/fs/mesh/interface/mesh.c#L938

The actual load call is here:

https://github.com/jedbrown/dohp/blob/master/src/fs/mesh/impls/serial/serial.c#L17

no "inoptions" are set.


You can get my code with

git clone --recursive git://github.com/jedbrown/dohp.git
cd dohp
mkdir build
cd build
cmake .. -DCMAKE_BUILD_TYPE=Debug
make -j4 ellip
src/fs/tests/ellip -dmesh_in ~/dl/he6.h5m -dirichlet 1,2,3
-dmesh_intermediate_adjacencies -snes_monitor -ksp_monitor -exact 4
-const_bdeg 1 -pc_type ilu -viewdhm

You need to have petsc-dev installed for this to go through.

This should write a file called ellip.dhm. It fails on read with recent
MOAB, and fails on write (after the solver converges) with my stack trace (I
use PETSc's error handling).


[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: iMesh(22) iBase_NOT_SUPPORTED: iMesh_save:ERROR saving a
mesh.  (MOAB Error Code: MB_FAILURE)!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Development HG revision:
4fd025355c9048e04aaa020799f60fbb1a911361  HG Date: Tue Jun 07 23:33:45 2011
+0200
[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: src/fs/tests/ellip on a ompi named kunyang by jed Thu Jun  9
18:55:52 2011
[0]PETSC ERROR: Libraries linked from /home/jed/petsc/ompi/lib
[0]PETSC ERROR: Configure run at Wed Jun  8 20:31:36 2011
[0]PETSC ERROR: Configure options --download-blacs --download-chaco
--download-hypre=externalpackages/hypre-2.7.0b.tar.gz --download-ml
--download-mumps --download-parms --download-scalapack --download-scotch
--download-spai --download-spooles --download-sundials --download-superlu
--download-superlu_dist --with-blas-lapack-dir=/usr --with-c2html
--with-cholmod-dir=/usr --with-clib-autodetect=0 --with-cxxlib-autodetect=0
--with-fortranlib-autodetect=0 --with-hdf5-dir=/usr --with-java-dir=/usr
--with-lgrind --with-mpi-dir=/usr --with-parmetis-dir=/usr
--with-shared-libraries --with-single-library=0 --with-sowing
--with-umfpack-dir=/usr -LIBS="-L/usr/lib/openmpi -lmpi_f90 -lmpi_f77
-lmpi_cxx -lstdc++ -lgfortran" PETSC_ARCH=ompi
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: dFSView_Cont_DHM() line 98 in
/home/jed/dohp/src/fs/impls/cont/contviewdhm.c
[0]PETSC ERROR: VecView_Dohp_FSCont_DHM() line 199 in
/home/jed/dohp/src/fs/impls/cont/contviewdhm.c
[0]PETSC ERROR: VecView_Dohp_FSCont() line 263 in
/home/jed/dohp/src/fs/impls/cont/contviewdhm.c
[0]PETSC ERROR: VecView() line 685 in
/home/jed/petsc/src/vec/vec/interface/vector.c
[0]PETSC ERROR: main() line 725 in /home/jed/dohp/src/fs/tests/ellip.c
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20110609/45a6d98b/attachment.htm>