> 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?
The function does not currently seem to have an argument that denotes a component mask. In that case, I'd say simply interpolate all parts of the solution, i.e. scalar + vectors, from a function with dim+1 components. You can always set the vector components in your function to zero, or something else artificial. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
