On Tue, Jun 24, 2014 at 10:40 AM, Sahai, Amal <sah...@illinois.edu> wrote:
> I have a few questions regarding importing ExodusII meshes created in Gridgen
> to Libmesh. I have created a simple structured mesh in Gridgen with the .exo
> format. Also, I made sure that the boundary conditions were defined as side
> sets. I face no problems reading the mesh using the ExodusII_IO object:
>
> ExodusII_IO exo(mesh);
> exo.read("BUMP100x100.exo");
>
> I get the following mesh message when I use mesh.print_info()
>
> Mesh Information:
> mesh_dimension()=2
> spatial_dimension()=3
> n_nodes()=10000
> n_local_nodes()=0
> n_elem()=9801
> n_local_elem()=0
> n_active_elem()=9801
> n_subdomains()=1
> n_partitions()=1
> n_processors()=1
> n_threads()=1
> processor_id()=0
>
> The issue is, when I try accessing the elements or the nodes of the mesh
> through the element iterators, I am unable to do so. A case in point is the
> following loop:
>
> MeshBase::const_element_iterator el =
> mesh.active_local_elements_begin();
> const MeshBase::const_element_iterator end_el =
> mesh.active_local_elements_end();
> for ( ; el != end_el; ++el)
> std::cout<<"AAAAAAAAAAAAAA\n";
>
> This loop is never entered. Is there something wrong with how I am importing
> the mesh? Or is there something that I need to correct while making the mesh
> in Gridgen? Please advice.
This is certainly a strange error. I originally suspected that your
mesh might contain element types not supported by libmesh, but that
would have generated an error message during read... is there any way
you can share this mesh file (google doc or something similar) with me
directly so that I can take a closer look?
--
John
------------------------------------------------------------------------------
Open source business process management suite built on Java and Eclipse
Turn processes into business applications with Bonita BPM Community Edition
Quickly connect people, data, and systems into organized workflows
Winner of BOSSIE, CODIE, OW2 and Gartner awards
http://p.sf.net/sfu/Bonitasoft
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
