Thank you for your help.
In view of system::currrent_solution(), whether may I use it after adaptive
mesh refinement and/or parallel mode (update() should be needed)?
If the system has several variables, which will be dealt with by
system::current_solution(). After all, it can't select the variables
by the parameter "global_dof_number".
In high order shaper function (such as second order Lagrange), there are one
point (likely the midpoint of one edge) in the mesh corresponding to the
variable, How to get the match between the variable and its space location.
Another question is about renumbering. In libmesh, renumbering of point in
the mesh is used in adaptive mesh refinement and parallel mode. Actually,
point renumbering affects the computerational
performance in parallel mode. whether can I not change the index of
points after renumbering in libmesh.
Renumbering in libmesh is optimized in view of the computerational
Thanks a lot.
On 9/10/07, Roy Stogner <roystgnr@...> wrote:
> On Mon, 10 Sep 2007, Yujie wrote:
> > I know the degree of freedom numbering is done element-by-element. I
> > know if one uses high order shape functions (p>1), he will obtain more
> > solution number than the nodes in the mesh. However, now I use
> > Lagrange shape functions, In one tetrahedron element, to my knowledge,
> > varialbe is relevant to one vertex of tetrahedron. Is it right?
> In this case, you can get the solution corresponding to a single Node
> Unlike the two methods which I previously suggested to you, this
> method will fail on non-vertex nodes on non-Lagrange elements, on
> interior nodes on quadratic Lagrange elements, and on non-vertex nodes
> on non-isoparametric Lagrange elements.