That's fair. To be more specific, I'm essentially trying to incorporate the
multigrid code from step 16 into an HDG solver that's very similar to step
51.
My code fails when I try to call the build(DoFHandler<dim>) function on my
MGTransferPrebuilt object as shown below
(similar to what's done in the void LaplaceProblem<dim>::solve() function
in step 16)
MGTransferPrebuilt<Vector,double>> mg_transfer(mg_constrained_dofs);
mg_transfer.build(this->dofh_uhat);
Note: dofh_uhat is a DoFHandler<dim> object corresponding to the
FE_FaceQ<dim,spacedim> finite element space.
*ERROR MESSAGE*
An error occurred in line <374> of file
</home/--------/software/dealii-9.4.0/source/fe/fe.cc> in function
const dealii::FullMatrix<double>& dealii::FiniteElement<<anonymous>,
<anonymous> >::get_prolongation_matrix(unsigned int, const
dealii::RefinementCase<dim>&) const [with int dim = 2; int spacedim = 2]
The violated condition was:
prolongation[refinement_case - 1][child].n() == this->n_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.
Stacktrace:
-----------
#0 /home/--------/software/dealii_current_build/lib/libdeal_II.g.so.9.4.0:
dealii::FiniteElement<2, 2>::get_prolongation_matrix(unsigned int,
dealii::RefinementCase<2> const&) const
#1 /home/--------/software/dealii_current_build/lib/libdeal_II.g.so.9.4.0:
void dealii::MGTransferPrebuilt<dealii::Vector<double> >::build<2,
2>(dealii::DoFHandler<2, 2> const&)
On Thursday, May 4, 2023 at 10:44:23 PM UTC-4 Wolfgang Bangerth wrote:
> On 5/4/23 19:26, Aditya Karthik Saravanakumar wrote:
> > I was able to set up all the MGLevelObjects but the MGTransferPrebuilt
> class
> > isn't able to generate the prolongation matrices.
>
> Aditya,
> I don't know enough about multigrid to help, but you might get better
> answers
> if you are specific about how your code looks and what actually happens.
> Saying "isn't able to" is too non-specific for anyone to venture a guess
> what
> the problem might be.
>
> Best
> W.
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: bang...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
>
>
--
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/b90d30e4-e5ee-479a-8dd9-bea9ea7cb48en%40googlegroups.com.
