If you want to form a RHS by integrating f_i Phi_i Phi_j, You could form a Mass matrix and then multiply with your vector of nodal values.
Rahul On Thu, Mar 25, 2010 at 11:40 PM, Karen Lee <[email protected]> wrote: > I'm afraid you misunderstood. I don't have the function that when given x, > y, z values gives me the function value. What I do have is just the values > at the nodes of the mesh, which need to be linearly interpolated such that I > will have something like exact_function. which gives me the value when > supplied with any x, y, z. > > > > On Thu, Mar 25, 2010 at 10:54 PM, Liang <[email protected]> wrote: > >> Karen Lee wrote: >> >>> I guess I'm not clear how to do this: Load data as a solution into that, >>> and >>> query >>> it when you're integrating your real system. >>> >>> I have: >>> Mesh mesh(3); >>> MeshData mesh_data(mesh); >>> mesh_data.activate(); >>> mesh.read (mesh_file, &mesh_data); >>> mesh_data.read(mesh_file); >>> EquationSystems equation_systems (mesh); >>> >>> >>> equation_systems.add_system<ExplicitSystem> ("RHS"); >>> equation_systems.get_system("RHS").add_variable("R", FIRST); >>> >>> After that, I'm not clear how exactly to load data as a solution in the >>> code. My goal is to get a linearly interpolated function of my data on the >>> nodes (in the form of exact_solution, such that I get the function value >>> out >>> when supplying x, y and z). >>> >>> Hope that clarifies things, and sorry for the multiple emails... >>> >>> Karen >>> >>> ------------------------------------------------------------------------------ >>> Download Intel® Parallel Studio Eval >>> Try the new software tools for yourself. Speed compiling, find bugs >>> proactively, and fine-tune applications for parallel performance. >>> See why Intel Parallel Studio got high marks during beta. >>> http://p.sf.net/sfu/intel-sw-dev >>> _______________________________________________ >>> Libmesh-users mailing list >>> [email protected] >>> https://lists.sourceforge.net/lists/listinfo/libmesh-users >>> >>> >>> >> so you already have the function, which is obtained from your discreted >> data? >> then just put the function as the exact_function. >> I think you are trying the 3D case, start from a 2d will be easier. >> >> Liang >> > ------------------------------------------------------------------------------ > Download Intel® Parallel Studio Eval > Try the new software tools for yourself. Speed compiling, find bugs > proactively, and fine-tune applications for parallel performance. > See why Intel Parallel Studio got high marks during beta. > http://p.sf.net/sfu/intel-sw-dev > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users > ------------------------------------------------------------------------------ Download Intel® Parallel Studio Eval Try the new software tools for yourself. Speed compiling, find bugs proactively, and fine-tune applications for parallel performance. See why Intel Parallel Studio got high marks during beta. http://p.sf.net/sfu/intel-sw-dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
