const unsigned int n1 = dof1.get_fe().dofs_per_cell;
const unsigned int n2 = dof2.get_fe().dofs_per_cell;
- Vector<double> u1_local(n1); // <- Vector<std::complex>
- Vector<double> u2_local(n2); // <- Vector<std::complex>
+ Vector<double> u1_local(n1);
+ Vector<double> u2_local(n2);
std::vector<types::global_dof_index> dofs(n2);
- FullMatrix<double> matrix(n2,n1); // ?
+ FullMatrix<double> matrix(n2,n1);
get_projection_matrix(dof1.get_fe(), dof2.get_fe(), matrix);
while (cell2 != end)
{
-
- // @whattodo Some of the problems in petsc_vector_base.h are instantiated here.
- // cell1->get_dof_values(u1, u1_local); // <- complaint line 655
- Assert ((false), ExcMessage ("This function is corrupt: @whattodo"));
+ cell1->get_dof_values(u1, u1_local);
matrix.vmult(u2_local, u1_local);
cell2->get_dof_indices(dofs);
-for (deal_II_dimension : DIMENSIONS; deal_II_space_dimension : SPACE_DIMENSIONS; VEC : SERIAL_VECTORS)
+// TODO[TY]: allow these for COMPLEX_SERIAL_VECTORS where applicable
+for (deal_II_dimension : DIMENSIONS; deal_II_space_dimension : SPACE_DIMENSIONS; VEC : REAL_SERIAL_VECTORS)
{
namespace FETools
\{
#if deal_II_dimension == deal_II_space_dimension
-#ifndef PETSC_USE_COMPLEX
template
void back_interpolate<deal_II_dimension>
(const DoFHandler<deal_II_dimension> &, const VEC &,
const VEC &,
const DoFHandler<deal_II_dimension> &, const ConstraintMatrix &,
VEC &);
-#endif
template
void project_dg<deal_II_dimension>
(const DoFHandler<deal_II_dimension> &, const VEC &,