From: Carlos A. <bd...@my...> - 2013-02-27 23:25:24
|
Thank you very much. That was actually very helpful. I had to save the mesh as quad4 before importing it to libmesh. One thing that I found though is that libmesh is not reading the number of local nodes "n_local_nodes()=0". there for the number of degrees of freedom is zero. I am printing mesh info reading the .e file vs using MeshTools::Generation::build_square (...) and this is what I get. *****EXODUS OUTPUT******* Mesh Information: mesh_dimension()=2 spatial_dimension()=3 n_nodes()=121 n_local_nodes()=0 n_elem()=100 n_local_elem()=0 n_active_elem()=100 n_subdomains()=1 n_partitions()=1 n_processors()=1 n_threads()=1 processor_id()=0 EquationSystems n_systems()=1 System #0, "Poisson" Type "LinearImplicit" Variables="u" Finite Element Types="LAGRANGE" Approximation Orders="FIRST" n_dofs()=0 n_local_dofs()=0 n_constrained_dofs()=0 n_local_constrained_dofs()=0 n_vectors()=1 n_matrices()=1 DofMap Sparsity Average On-Processor Bandwidth <= 0 Average Off-Processor Bandwidth <= 0 Maximum On-Processor Bandwidth <= 0 Maximum Off-Processor Bandwidth <= 0 DofMap Constraints Number of DoF Constraints = 0 ******************* vs. ****************mesh generation output****** Mesh Information: mesh_dimension()=2 spatial_dimension()=3 n_nodes()=256 n_local_nodes()=256 n_elem()=225 n_local_elem()=225 n_active_elem()=225 n_subdomains()=1 n_partitions()=1 n_processors()=1 n_threads()=1 processor_id()=0 EquationSystems n_systems()=1 System #0, "Poisson" Type "LinearImplicit" Variables="u" Finite Element Types="LAGRANGE" Approximation Orders="FIRST" n_dofs()=256 n_local_dofs()=256 n_constrained_dofs()=0 n_local_constrained_dofs()=0 n_vectors()=1 n_matrices()=1 DofMap Sparsity Average On-Processor Bandwidth <= 8.26562 Average Off-Processor Bandwidth <= 0 Maximum On-Processor Bandwidth <= 9 Maximum Off-Processor Bandwidth <= 0 DofMap Constraints Number of DoF Constraints = 0 ************************************************ Is there an extra step I need to do to read the exodus mesh??? Thank you in advance. On Wed, Feb 27, 2013 at 1:06 PM, John Peterson <jwp...@gm...> wrote: > > > > > On Feb 27, 2013, at 12:23 PM, "Paul T. Bauman" <ptb...@gm...> wrote: > > > On Wed, Feb 27, 2013 at 11:51 AM, Carlos Acosta <bd...@my...> > wrote: > > > >> but when I import an exodus mesh I get errors > >> recognizing the exodus format. > >> It says: > >> > >> > >> ERROR! Unrecognized element type_str: SHELL4 > > > > Hmm, this is strange because Exodus reader should understand SHELL4. > > But it doesn't quote work right: shell4's actually have six "sides" and so > some of our side id mapping stuff doesn't work. Definitely stick to > Quad4's. I think I will patch this to go back to making shell4 throw an > error, maybe something like "please export quad4 from cubit instead"... > > > > > Could > > you rerun with METHOD=dbg and see if an assert is tripped? How big is > your > > mesh? > > It's possible he's using an old version, maybe a sourceforge SVN version? > > > > > ------------------------------------------------------------------------------ > > Everyone hates slow websites. So do we. > > Make your web apps faster with AppDynamics > > Download AppDynamics Lite for free today: > > http://p.sf.net/sfu/appdyn_d2d_feb > > _______________________________________________ > > Libmesh-users mailing list > > Lib...@li... > > https://lists.sourceforge.net/lists/listinfo/libmesh-users > -- Carlos Acosta Berlinghieri Graduate Research Assistant Department of Biomedical Engineering Computational Center University of Texas at San Antonio (210)458-4680 (lab phone) |