template <int dim, class InVector, class OutVector, int spacedim>
void interpolate (const DoFHandler<dim,spacedim> &dof_1,
const DoFHandler<dim,spacedim> &dof_2,
-- const FullMatrix<double> &transfer,
++ const FullMatrix<double> &transfer,
const InVector &data_1,
OutVector &data_2);
for (; h != endh; ++h, ++l)
{
- // @whattodo Some of the problems in petsc_vector_base.h are instantiated here.
- /* h->get_dof_values(data_1, cell_data_1); // <- complaint line 655 */
- h->get_dof_values(data_1, cell_data_1);
++ // @whattodo Some of the problems in petsc_vector_base.h are
++ // instantiated here. Should look into DoFCellAccessor too. It
++ // looks like *maybe* vector_tools needs to be templated (run
++ // grep -R "double" on that file). Additionally, *this* file
++ // needs to be templated with "number".
++ /* h->get_dof_values(data_1, cell_data_1); */
+ Assert ((false), ExcMessage ("This function is corrupt: @whattodo"));
+
transfer.vmult(cell_data_2, cell_data_1);
l->get_dof_indices (local_dof_indices);