From: Wolfgang Bangerth <[email protected]>
To: Ted Kord <[email protected]>
Cc: [email protected]
Sent: Sun, 24 January, 2010 21:23:16
Subject: Re: [deal.II] How to Obtain the Coordinates of all Degrees of Freedom
> Is the mapping from the support points one-to-one relative to the
> entries in the solution vector ?
That depends on the finite element you use. For example, if you solve the
elasticity equation, you would have displacements in x, y, and z direction
defined in the vertices of the cells, so there will be three degrees of
freedom that are all defined at the same support point.
> For example, can I reassemble the
> solution_vector 'manually' in exactly the same order as given in the
> support_points vector and pass it back as the initial values for a
> parabolic pde?
I have no idea what exactly you want to do, but I can also not see any
fault with what you write :-)
Best
W.
===
One degree of freedom per node is sufficient for the moment. So,
I want to solve an ode at every node.
I then want to assemble each solution from each ode/node into a vector, i.e.,
(old_solution_u).
Then pass the vector as the initial condition for a parabolic pde.
Solve the pde.
Then, for each component of the pde solution_vector, use that as an initial
condition for the odes again.
Then repeat.
So, assembling the vector is the tricky part. I'd like to do it in exactly the
same way the the pde writes out each nodal solution. that's why I was wondering
if the DoFTools::map_dofs_to_support_points gives these in order of evaluation.
Does it? Or, is there a better way?
Regards
Ted
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
