search

ARudik / feelpp

Automatically exported from code.google.com/p/feelpp
0 stars 0 forks source link

trace or createSubmesh in parallel #13

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago 
I meet a problem with trace and createSubmesh in parallel. 
First we create our mesh M_mesh=createGMSHMesh(.., proc id) in parallel. Then

 M_Xh = space_type::New( M_mesh, index ); // index is the proc id

auto sol = M_Xh->element();
auto trace_mesh = M_mesh->trace( markedfaces(M_mesh,1) );
auto Th = trace_space_type::New( trace_mesh );

auto t = vf::project( _space=Th, _range=elements(trace_mesh), 
_expr=vf::idv(sol) ); 

This calculation does not execute in parallel.

Original issue reported on code.google.com by samakea...@gmail.com on 28 Mar 2012 at 4:49