Hi all
I am trying to solve the following time-dependent variant of step-20
(here p is a vector and u a scalar):
u_t + div p = f
p + a grad u = 0
u(0) = u_initial
I define my finite element as
<file:///C:\Users\Hitor%20Pina\Desktop\dealII700-doc\doxygen\deal.II\classFE
System.html> FESystem<dim> finite_element (
<file:///C:\Users\Hitor%20Pina\Desktop\dealII700-doc\doxygen\deal.II\classFE
__Q.html> FE_Q<dim>(1), dim+1), the first dim components pertaining to the
vector p and the last component to the scalar u.
The difficulty I am facing is how to interpolate the initial condition: if I
simply use VectorTools::interpolate (dof_handler, u_initial, u0) deal.ii
complains saying that there is a dimension mismatch, which is true. How do I
tell deal.ii to use in the interpolate function only the scalar part of the
finite element?
--
TIA
Heitor
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
