Hi all, I'm having a problem with trying to initialise a vector solution. Background: I am solving two systems of equations back to back repeatedly. I need the solution for the 2nd system to solve the first one at the next 'timestep', say. But initially/at the first 'timestep', I solve the first system assuming the solution to the 2nd system is known, and therefore I want an interpolation of some initial function onto the solution of the 2nd system to solve the first system.
The second system is a darcy-type, so I am using a code pretty much analogous to step-20, with the same FEs. Now, VectorTools::interpolate(dof_handler, exact_solution, solution); gives me a dimension mismatch error - I'm assuming that this is due to the fact that RT elements can't distinguish between the x and y components. I then tried adding a component mask in various ways to put the velocity and pressure solutions in separately, but I still get the same error. Doing a component mask with solution.block(0) separately also gives me a number of dofs mismatch error. Is there any way around this? The above interpolation works if i have Fe_Q elements (like in step-22) so I've narrowed it down to the fact that it is because of the function space, but I need the function spaces as they are. I can't seem to find another way around this using the interpolation function, and hope that there is actually a simple way around this? Many thanks in advance. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
