[Dune] [Dune-Commit] dune-grid r7996 - trunk/dune/grid/common

[Dedner, Andreas]

Hi Benedikt.

I am a bit surprised since from earlier (a lot earlier) emails
I got the impression that this removal was due to a discussion between Peter, yourself and (I think) Oliver.
I can not find that emails any more where this discussion was started.
All of this is a few years old and the boundarySegmentIndex was the result of this.

At the moment it is not possible to easily replace the original boundaryId with the boundarySegmentIndex
and if you have to do quite a bit of work yourself.
To get the information out of the dgf parser there is a method
 template< class GridImp, template< class > class IntersectionImp >
 const DGFBoundaryParameter::type &
 GridPtr::parameters ( const Intersection< GridImp, IntersectionImp > & intersection ) const
which returns the boundaryId for a given intersection (assuming that it lies on the boundary). Then you could store these ids in a vector
using intersection.boundarySegmentIndex().

Parallel is (as Oli mentioned) still a problem.
If gmsh is better suited then DGF I can not say. No idea how
boundary identifiers are transported through the gmsh reader
but the problem with parallel are the same.

I should add that I have never testes this since I also still use
the boundaryId....

Best
Andreas

________________________________
From: dune-bounces+a.s.dedner=warwick.ac.uk at dune-project.org [dune-bounces+a.s.dedner=warwick.ac.uk at dune-project.org] on behalf of Benedikt Oswald [benedikt.oswald at psi.ch]
Sent: Wednesday, April 18, 2012 1:02 PM
To: dune
Subject: Re: [Dune] [Dune-Commit] dune-grid r7996 - trunk/dune/grid/common

Dear Dune,


I was woken up somewhat by the recent discussion on 'boundaryId'
and the fact that it is marked at 'deprecated'. Now, since we need the boundarId
quite much, what would be the 'clean Dune style' of using id's on boundary faces ?

Should we use

template<class GridImp, template< class > class IntersectionImp>
size_t Dune::Intersection<http://www.dune-project.org/doc/doxygen/html/classDune_1_1Intersection.html>< GridImp, IntersectionImp >::boundarySegmentIndex        (                )      const [inline]

instead of:


template<class GridImp, template< class > class IntersectionImp>
int Dune::Intersection<http://www.dune-project.org/doc/doxygen/html/classDune_1_1Intersection.html>< GridImp, IntersectionImp >::boundaryId     (                )      const [inline]



How does boundarySegmentIndex connect to boundary tags in the DGF file ?

Also, how would we use this in the parallel case ? This is most relevant for us.

Would it be wiser to not used DGF but read directly from gmsh ?


Thanks a lot for your enlightenment and have a great day!

Benedikt




--------------------------------------------------------------------------------------------------------------------------------------------------------------------
Benedikt Oswald, Dr. sc. techn., dipl. El. Ing. ETH, www.psi.ch<http://www.psi.ch>, GFA
Paul Scherrer  Institute (PSI), CH-5232 Villigen, Suisse, benedikt.oswald at psi.ch<mailto:benedikt.oswald at psi.ch>, +41(0)56 310 32 12

"Passion is required for any great work, and for the revolution passion and audacity are required in big doses.",
 Ernesto 'Che' Guevara, Letter to his parents.
--------------------------------------------------------------------------------------------------------------------------------------------------------------------

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20120418/04b9e8dd/attachment.htm>