Hello,
I have configured libmesh with --enable-parmesh and I am
reading a partitioned mesh with Nemesis_IO. My mesh is too
big to refine it before partitioning - so I tried the
following:
Mesh new_mesh;
Nemesis_IO(new_mesh).read("meshfile");
if (!new_mesh. is_prepared ())
new_mesh.prepare_for_use();
new_mesh.print_info();
MeshRefinement meshrefine (new_mesh);
meshrefine.uniformly_refine(1);
new_mesh.print_info();
EquationSystems new_equation_systems (new_mesh);
//Add some systems
...
//
new_equation_systems.print_info();
So, I uniformly refine the mesh, then I build the
EquationSystems object.
For the *.print_info() 's I get the following output:
serial: 0
Mesh Information:
mesh_dimension()=3
spatial_dimension()=3
n_nodes()=869738
n_local_nodes()=224297
n_elem()=647844
n_local_elem()=161961
n_active_elem()=647844
n_subdomains()=4
n_partitions()=1
n_processors()=4
n_threads()=1
processor_id()=0
refine mesh
Mesh Information:
mesh_dimension()=3
spatial_dimension()=3
n_nodes()=6953271
n_local_nodes()=1764933
n_elem()=5830596
n_local_elem()=1457649
n_active_elem()=5182752
n_subdomains()=4
n_partitions()=1
n_processors()=4
n_threads()=1
processor_id()=0
EquationSystems
n_systems()=1
System #0, "Diffusion"
Type "TransientLinearImplicit"
Variables="T"
Finite Element Types="LAGRANGE"
Approximation Orders="FIRST"
n_dofs()=869738 // !!!!!!!!!!!!!!!!!!!!!!!
n_local_dofs()=224297
n_constrained_dofs()=0
n_vectors()=3
DofMap Sparsity
Average On-Processor Bandwidth <= 14.7421
Average Off-Processor Bandwidth <= 0.0864746
Maximum On-Processor Bandwidth <= 65
Maximum Off-Processor Bandwidth <= 16
DofMap Constraints
Number of DoF Constraints = 0
Although I refine the mesh before I build the
EquationSystems object, the n_dofs() correspond to the
coarse mesh, not to the refined one.
Any suggestions?
Robert
------------------------------------------------------------------------------
RSA(R) Conference 2012
Save $700 by Nov 18
Register now
http://p.sf.net/sfu/rsa-sfdev2dev1
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
