[dune-pdelab] adaptivity with a composite grid function space
Dodwell, Timothy
T.Dodwell at exeter.ac.uk
Mon Jan 18 17:33:44 CET 2016
Dear all,
On the latest Dune PDElab release (on Mac OS X 10.10) I am trying to implement adaptive FEM for 3D elasticity (vector valued problem)
I get a problem when calling adapt_grid ( below is the compiler error )
Looks like adapt_grid does not work for composite grid function spaces, is that right? Is there a work around?
Thanks in advance.
Tim
Tims-MacBook-Pro:src timdodwell$ make dune-uq
Scanning dependencies of target dune-uq
[ 50%] Building CXX object src/CMakeFiles/dune-uq.dir/dune-uq.cc.o
In file included from /Users/timdodwell/git_workspace/dune/dune-uq/src/dune-uq.cc:83:
/Users/timdodwell/git_workspace/dune/dune-pdelab/dune/pdelab/adaptivity/adaptivity.hh:387:51: error: no member
named 'FiniteElementMap' in
'Dune::PDELab::PowerCompositeGridFunctionSpaceTraits<Dune::PDELab::PartitionViewEntitySet<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, Dune::PartitionSet<31> >,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LexicographicOrderingTag, 3>'
using FiniteElement = typename GFS::Traits::FiniteElementMap::Traits::FiniteElement;
~~~~~~~~~~~~~^
/Users/timdodwell/git_workspace/dune/dune-pdelab/dune/pdelab/adaptivity/adaptivity.hh:565:17: note: in
instantiation of template class
'Dune::PDELab::replay_visitor<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::PDELab::istl::BlockVector<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::BlockVector<Dune::FieldVector<double, 1>, std::__1::allocator<Dune::FieldVector<double, 1> > > >,
Dune::PDELab::istl::BlockVector<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >, Dune::BlockVector<Dune::FieldVector<int,
1>, std::__1::allocator<Dune::FieldVector<int, 1> > > > >' requested here
Visitor visitor(gfsu,u,uc,_leaf_offset_cache);
^
/Users/timdodwell/git_workspace/dune/dune-pdelab/dune/pdelab/adaptivity/adaptivity.hh:638:20: note: in
instantiation of member function 'Dune::PDELab::GridAdaptor<Dune::UGGrid<3>,
Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::PDELab::istl::BlockVector<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::BlockVector<Dune::FieldVector<double, 1>, std::__1::allocator<Dune::FieldVector<double, 1> > > >,
Dune::PDELab::L2Projection<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::PDELab::istl::BlockVector<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::BlockVector<Dune::FieldVector<double, 1>, std::__1::allocator<Dune::FieldVector<double, 1> > > > >
>::replayData' requested here
grid_adaptor.replayData(grid,gfs,projection,x1,transferMap1);
^
/Users/timdodwell/git_workspace/dune/dune-uq/src/user_inputs/../UQ_functions/../user_inputs/user_model_3Delasticity.hh:571:17: note:
in instantiation of function template specialization 'Dune::PDELab::adapt_grid<Dune::UGGrid<3>,
Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::PDELab::istl::BlockVector<Dune::PDELab::CompositeGridFunctionSpace<Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none,
1>, Dune::PDELab::LexicographicOrderingTag,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> >,
Dune::PDELab::GridFunctionSpace<Dune::GridView<Dune::UGGridLevelGridViewTraits<const Dune::UGGrid<3> >
>, Dune::PDELab::PkLocalFiniteElementMap<Dune::GridView<Dune::UGGridLevelGridViewTraits<const
Dune::UGGrid<3> > >, double, double, 1>, Dune::PDELab::ConformingDirichletConstraints,
Dune::PDELab::istl::VectorBackend<Dune::PDELab::istl::Blocking::none, 1>,
Dune::PDELab::LeafOrderingTag<Dune::PDELab::EmptyParams> > >,
Dune::BlockVector<Dune::FieldVector<double, 1>, std::__1::allocator<Dune::FieldVector<double, 1> > > >
>' requested here
adapt_grid( grid_, gfs, xf,1);
^
/Users/timdodwell/git_workspace/dune/dune-uq/src/UQ_functions/STD_MLMC.hh:278:29: note: in instantiation of
member function 'MODEL<Dune::UGGrid<3> >::user_adaptive_fem_driver' requested here
Qf = model_.user_adaptive_fem_driver(l,z);
^
/Users/timdodwell/git_workspace/dune/dune-uq/src/UQ_functions/STD_MLMC.hh:102:61: note: in instantiation of
function template specialization 'STD_MLMC<MODEL<Dune::UGGrid<3> > >::getSample<COEFF>' requested here
for (int i = 0; i < numSamples; i++){ Ytmp[i] = getSample(L,z); } // Compute Samples on ea...
^
/Users/timdodwell/git_workspace/dune/dune-uq/src/dune-uq.cc:140:12: note: in instantiation of member function
'STD_MLMC<MODEL<Dune::UGGrid<3> > >::apply' requested here
mlmc.apply();
^
1 error generated.
make[3]: *** [src/CMakeFiles/dune-uq.dir/dune-uq.cc.o] Error 1
make[2]: *** [src/CMakeFiles/dune-uq.dir/all] Error 2
make[1]: *** [src/CMakeFiles/dune-uq.dir/rule] Error 2
make: *** [src/CMakeFiles/dune-uq.dir/rule] Error 2
Dr Tim Dodwell
Lecturer in Engineering Mathematics
Rm 221 - Harrison Building
College of Engineering, Mathematics & Physical Sciences
University of Exeter
Exeter
Devon
EX4 4PY
mail: t.dodwell(at)exeter.ac.uk
tel: +44 (0)1392 725899
mob: +44 (0)7745 622870
web: http://emps.exeter.ac.uk/engineering/staff/td336
Papers and Pre-prints: @Research-Gate<https://www.researchgate.net/profile/Timothy_Dodwell>
Citations: @Google-Scholar<https://scholar.google.co.uk/citations?user=lPpjRfUAAAAJ&hl=en>
[http://www.exeter.ac.uk/media/universityofexeter/communicationservices/design/corporate_sig.jpg]
This email and any attachment may contain information that is confidential, privileged, or subject to copyright, and which may be exempt from disclosure under applicable legislation. It is intended for the addressee only. If you received this message in error, please let me know and delete the email and any attachments immediately. The University will not accept responsibility for the accuracy/completeness of this email and its attachments.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20160118/d30525c8/attachment.htm>
More information about the dune-pdelab
mailing list