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

Reply via email to