Another question I have is, is there any say to speed up the reading of a large mesh (about 55M elements)? It seemed that a lot of time was spent looking for neighbors... (when renumbering is turned off. When it's not, much time is spent renumbering stuff...)
Thanks, Karen On Thu, Apr 7, 2011 at 3:39 PM, Karen Lee <[email protected]> wrote: > So my code with MeshData worked fine, but is not parallelizable. So > I'm trying to get some help of understanding how to use > ExplicitSystem. > > The only example I found from the examples that contain ExplicitSystem > is Example 2, and being a newbie to C++ does not help me understand > how I can give values to and access the variables from ExplicitSystem. > > 1. What is the format of the input to > > void read_parallel_data (Xdr &io, const bool read_additional_data) > > ? I would like to put my data for all the nodes in .xdr (or is it > .xda?) form. I have .xda files of meshes from previous runs where I > output the .xda files, but I'm not sure what the format of .xdr is. If > I have a single field of 5's (5 nodes), do I just do the following? > > 5 > 1 > 1 > 1 > 1 > 1 > > I have this code within a loop that goes through all the elements > from my previous program that uses MeshData, which I now wanna get rid > of using: > > for (unsigned int i=0; i<phi.size(); i++) { > const std::vector<Number>& data = > mesh_data.get_data(elem->get_node(i)); > eps(qp) += data[1] * phi[i][qp]; > rho(qp) += data[0] * phi[i][qp]; > } > > I guess I just need a function that takes the place of > mesh_data.get_data(elem->get_node(i)) and takes the > "elem->get_node(i)"th item from the explicit system.(I would like to > have 3 of such variables, and just be able to access the ith element > of each of them.) > > Any help will be deeply appreciated. > > Thanks, > Karen > > > > > > ---------- Forwarded message ---------- > From: Roy Stogner <[email protected]> > Date: Wed, Mar 30, 2011 at 1:17 AM > Subject: Re: [Libmesh-users] Adaptive Mesh: how to update data values? > To: Karen Lee <[email protected]> > Cc: libmesh-users <[email protected]> > > > > On Wed, 30 Mar 2011, Karen Lee wrote: > >> HOWEVER, now that I'm trying to do adaptive refinement, it doesn't >> look like these data values are refined. > > Nope. The short answer: don't use MeshData. In this case, the > primary reason is that there's no good way for the library to know how > you want it interpreted and interpolated. Use a separate > ExplicitSystem with a finite element data field instead, and then the > data will be projected from mesh to mesh in as consistent a way as > possible. > --- > Roy > ------------------------------------------------------------------------------ Xperia(TM) PLAY It's a major breakthrough. An authentic gaming smartphone on the nation's most reliable network. And it wants your games. http://p.sf.net/sfu/verizon-sfdev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
