Hello all,
I am encountered with this error while using adaptive mesh refinement in parallel.
virtual const FullMatrix<double> &dealii::FiniteElement<2, 2>::get_prolongation_matrix(const unsigned int, const RefinementCase<dim> &) const
The violated condition was:
prolongation[refinement_case - 1][child].n() == this->dofs_per_cell
Additional information:
You are trying to access the matrices that describe how to embed a finite element function on one cell into the finite element space on one of its children (i.e., the 'embedding' or 'prolongation' matrices). However, the current finite element can either not define this sort of operation, or it has not yet been implemented.
A little background of my problem.
Deal ii version: 9.1.1
I am using Weak Galerkin FE so the FESystem is (FE_DGQ<dim>(degree), 1, FE_FaceQ<dim>(degree), 1).
After mesh is refined, I need to interpolate solution from old to new mesh using parallel::distributed::SolutionTransfer and this is where this error pops out. I believe there is some conflict using solutiontransfer with the FE_FaceQ element.
Does anyone know how to solve this?
Thanks,
Yu