Well, thinking about it, it seems that spacedim=3 is not relevant here. I
can do
const FE_DGQArbitraryNodes<dim,dim> fe1(QGauss<1>(fe.degree+1));
const FE_DGQArbitraryNodes<dim,dim> fe2(QGauss<1>(fe.degree));
FETools::get_back_interpolation_matrix(fe1, fe2, filter_matrix);
This should give me the same matrix as in the case dim=2, spacedim=3, right
???
Best
praveen
On Wed, Oct 19, 2016 at 7:22 PM, Praveen C <cprav...@gmail.com> wrote:
> Dear all
>
> I am doing this
>
> const FE_DGQArbitraryNodes<dim,spacedim> fe(QGauss<1
> >(degree+1));
>
> filter_matrix.reinit(fe.dofs_per_cell, fe.dofs_per_cell);
>
> const FE_DGQArbitraryNodes<dim,spacedim> fe1(QGauss<1>(degree));
>
> FETools::get_back_interpolation_matrix(fe, fe1, filter_matrix);
>
>
> but this gives me compilation error
>
> $ make
>
> [ 20%] *Linking CXX executable swm*
>
> Undefined symbols for architecture x86_64:
>
> "void dealii::FETools::get_back_interpolation_matrix<2, double,
> 3>(dealii::FiniteElement<2, 3> const&, dealii::FiniteElement<2, 3> const&,
> dealii::FullMatrix<double>&)", referenced from:
>
> SWM::SWProblem<3>::make_grid_and_dofs() in swm.cc.o
>
> ld: symbol(s) not found for architecture x86_64
>
>
> Best
> praveen
>
--
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.
For more options, visit https://groups.google.com/d/optout.
