Karen Lee 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. > > > > > ohh, sorry for misunderstanding, you may write a interpolation function subroutine to calculate the polynominals of the exact function, but do not assemble the exact_function because it will be very a lot of terms. Newton or Lagrange method should be able to do this linear interpolation, I did 1D interpolation with 20 nodes but have no idea on 3D interpolation. then just bring the x,y,z into your interpolation subroutine and return the exact variable on each node.
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
