On Wed, Aug 28, 2013 at 8:31 PM, John Peterson <[email protected]> wrote:
> Are you planning to do adaptivity ever? The indexing into this
> std::vector may get a bit tricky in that scenario...
>
> Also, when running in parallel, you probably don't have to store mass
> matrices for non-local elements, but if you use a vector you'll waste a lot
> of space storing empty mass matrices.
>
Yes I want to do adaptive and parallel also.
When I loop over mesh elements like this
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)
{
this loops only over the elements local to the processor, right ? In that
case I would assemble and store the mass matrix only for the locally held
elements.
Thanks
praveen
------------------------------------------------------------------------------
Learn the latest--Visual Studio 2012, SharePoint 2013, SQL 2012, more!
Discover the easy way to master current and previous Microsoft technologies
and advance your career. Get an incredible 1,500+ hours of step-by-step
tutorial videos with LearnDevNow. Subscribe today and save!
http://pubads.g.doubleclick.net/gampad/clk?id=58040911&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
