Hi all,
I am trying to reproduce the code in Step 38 that extracts the boundary
mesh, but I'd love to use a parallel::shared::Triangulation rather than a
Triangulation
<http://www.dealii.org/8.4.1/doxygen/deal.II/classTriangulation.html>
class. In a nutshell, I did this:
template <int dim>
class SmallStrainBeltramiDiffusionAndMechanicalProblem
{
public:
[ ... ]
private:
[ ... ]
// Surface manifold dimension
static const unsigned int manifold_dim = dim-1;
// FEM member variables
parallel::shared::Triangulation<dim> triangulation;
parallel::shared::Triangulation<manifold_dim,dim>
manifold_triangulation;
DoFHandler<dim> dof_handler;
DoFHandler<manifold_dim,dim> manifold_dof_handler;
MappingQ<manifold_dim,dim> manifold_mapping;
[ ... ]
}
template <int dim>
void SmallStrainBeltramiDiffusionAndMechanicalProblem<dim>::create_coarse_grid
( bool neumann )
{
const Point<dim> center;
GridGenerator::hyper_ball (triangulation, center, 1.);
triangulation.refine_global ( 6 );
GridGenerator::extract_boundary_mesh (triangulation,
manifold_triangulation); //, boundary_ids);
[ ...]
}
I got this error message:
*--------------------------------------------------------*
*An error occurred in line <220> of file
<../source/distributed/tria_base.cc> in function*
* virtual types::subdomain_id dealii::parallel::Triangulation<1,
2>::locally_owned_subdomain() const*
*The violated condition was: *
* dim > 1*
*The name and call sequence of the exception was:*
* ExcNotImplemented()*
*Additional Information: *
*You are trying to use functionality in deal.II that is currently not
implemented. In many cases, this indicates that there simply didn't appear
much of a need for it, or that the author of the original code did not have
the time to implement a particular case. If you hit this exception, it is
therefore worth the time to look into the code to find out whether you may
be able to implement the missing functionality. If you do, please consider
providing a patch to the deal.II development sources (see the deal.II
website on how to contribute).*
*Stacktrace:*
*-----------*
*#0 2 libdeal_II.g.8.4.1.dylib 0x0000000108b2c3f7
_ZNK6dealii8parallel13TriangulationILi1ELi2EE23locally_owned_subdomainEv +
231: 2 libdeal_II.g.8.4.1.dylib 0x0000000108b2c3f7
_ZNK6dealii8parallel13TriangulationILi1ELi2EE23locally_owned_subdomainEv *
*#1 3 libdeal_II.g.8.4.1.dylib 0x00000001081ce7f1
_ZNK6dealii12CellAccessorILi1ELi2EE8is_ghostEv + 577: 3
libdeal_II.g.8.4.1.dylib 0x00000001081ce7f1
_ZNK6dealii12CellAccessorILi1ELi2EE8is_ghostEv *
*#2 4 libdeal_II.g.8.4.1.dylib 0x0000000108b2cb8f
_ZN6dealii8parallel13TriangulationILi1ELi2EE19update_number_cacheEv + 1455:
4 libdeal_II.g.8.4.1.dylib 0x0000000108b2cb8f
_ZN6dealii8parallel13TriangulationILi1ELi2EE19update_number_cacheEv *
*#3 5 libdeal_II.g.8.4.1.dylib 0x0000000108b33f6d
_ZN6dealii8parallel6shared13TriangulationILi1ELi2EE20create_triangulationERKNSt3__16vectorINS_5PointILi2EdEENS4_9allocatorIS7_EEEERKNS5_INS_8CellDataILi1EEENS8_ISE_EEEERKNS_11SubCellDataE
+ 413: 5 libdeal_II.g.8.4.1.dylib 0x0000000108b33f6d
_ZN6dealii8parallel6shared13TriangulationILi1ELi2EE20create_triangulationERKNSt3__16vectorINS_5PointILi2EdEENS4_9allocatorIS7_EEEERKNS5_INS_8CellDataILi1EEENS8_ISE_EEEERKNS_11SubCellDataE
*
*#4 6 libdeal_II.g.8.4.1.dylib 0x0000000107c383b9
_ZN6dealii13GridGenerator21extract_boundary_meshINS_8parallel6shared13TriangulationELi2ELi2EEENSt3__13mapINT_IXmiT0_Li1EEXT1_EE13cell_iteratorENS7_IXT0_EXT1_EE13face_iteratorENS5_4lessIS9_EENS5_9allocatorINS5_4pairIKS9_SB_EEEEEERKSA_RS8_RKNS5_3setIhNSC_IhEENSE_IhEEEE
+ 4537: 6 libdeal_II.g.8.4.1.dylib 0x0000000107c383b9
_ZN6dealii13GridGenerator21extract_boundary_meshINS_8parallel6shared13TriangulationELi2ELi2EEENSt3__13mapINT_IXmiT0_Li1EEXT1_EE13cell_iteratorENS7_IXT0_EXT1_EE13face_iteratorENS5_4lessIS9_EENS5_9allocatorINS5_4pairIKS9_SB_EEEEEERKSA_RS8_RKNS5_3setIhNSC_IhEENSE_IhEEEE
*
*#5 7 heat-eq 0x0000000100058754
_ZN20SmallStrainMechanics48SmallStrainBeltramiDiffusionAndMechanicalProblemILi2EE18create_coarse_gridEb
+ 196: 7 heat-eq 0x0000000100058754
_ZN20SmallStrainMechanics48SmallStrainBeltramiDiffusionAndMechanicalProblemILi2EE18create_coarse_gridEb
*
*#6 8 heat-eq 0x000000010005432f
_ZN20SmallStrainMechanics48SmallStrainBeltramiDiffusionAndMechanicalProblemILi2EE16do_timestep_zeroEjPK26TimeIntegrationDataManagerb
+ 191: 8 heat-eq 0x000000010005432f
_ZN20SmallStrainMechanics48SmallStrainBeltramiDiffusionAndMechanicalProblemILi2EE16do_timestep_zeroEjPK26TimeIntegrationDataManagerb
*
*#7 9 heat-eq 0x00000001000343cd
_ZN20SmallStrainMechanics48SmallStrainBeltramiDiffusionAndMechanicalProblemILi2EE3runEjb
+ 717: 9 heat-eq 0x00000001000343cd
_ZN20SmallStrainMechanics48SmallStrainBeltramiDiffusionAndMechanicalProblemILi2EE3runEjb
*
*#8 10 heat-eq 0x0000000100033b17 main + 391:
10 heat-eq 0x0000000100033b17 main *
*#9 11 libdyld.dylib 0x00007fffa8348235 start + 1:
11 libdyld.dylib 0x00007fffa8348235 start *
*--------------------------------------------------------*
Am I missing something or attempting at doing something that has not been
implemented yet?
Thanks
Alberto
--
Informativa sulla Privacy: http://www.unibs.it/node/8155
--
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
