Dear dealii-users, I am having an issue interpolating functions (defined via expressions) into a PETScWrappers::MPI::BlockVector
The expressions are ordered according to component. But when setting up the system I change to DoFRenumbering::subdomain_wise (dof_handler); to be able to use contiguous rows and columns. Then later executing the calls: initial_conditions.initialize(FunctionParser<dim>::default_variable_names(), initial_conditions_strings, std::map<std::string, double>()); VectorTools::interpolate (dof_handler, initial_conditions, solution_old); mixes up the values for each of the solution blocks. The above problem does not exist if I renumber by component, i.e. DoFRenumbering::component_wise(dof_handler); but that is not how I want to renumber. It seems a silly problem but I cannot find any examples or documentation on how to overcome this issue. Any help will be appreciated. Regards, -- 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.
