On Fri, Aug 12, 2011 at 4:21 AM, robert <robert.bod...@unil.ch> wrote:
> Hello,
>
> again, I have a question as far as reading meshes with tetgen is
> concered.
> I read a mesh by
>
> Mesh new_mesh;
> EquationSystems new_equation_systems (new_mesh);
> TetGenIO TETGEN(new_mesh);
>
> OStringStream inmesh;
> inmesh<<"geometry_wings/export_pov/out/"<<ReadMesh[actual_pulse]<<"_box.1.ele";
> TETGEN.read(inmesh.str().c_str());
>
> Afterwards, if I look at the equation system I get n_dofs() = 0;
This makes sense, you haven't added any Systems to your
EquationSystems object. All the n_dofs() function does is loop over
the Systems stored in the EquationSystems object and add up the dofs
of each...
> Up to now I have looped through the nodes and
>
> MeshBase::node_iterator nd = mesh.nodes_begin();
> const MeshBase::node_iterator end_nd = mesh.nodes_end();
>
> int proc = 0;
> for ( ; nd != end_nd; ++nd){
> Node* NOD = *nd;
> NOD->processor_id() = proc;
>
> }
>
> My question now is, if the code still works in parallel?
Probably, but the real question is why you'd want to do this. If you
are running in parallel it will almost certainly screw up the domain
decomposition which has been determined by the Partitioner classes...
--
John
------------------------------------------------------------------------------
uberSVN's rich system and user administration capabilities and model
configuration take the hassle out of deploying and managing Subversion and
the tools developers use with it. Learn more about uberSVN and get a free
download at: http://p.sf.net/sfu/wandisco-dev2dev
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
