[dune-pdelab] dune-pdelab Digest, Vol 77, Issue 1
Dodwell, Timothy
T.Dodwell at exeter.ac.uk
Fri May 13 19:43:35 CEST 2016
Thanks Steffan, I’ll give it a go!
That seems a bit odd, I assumed you must be able to do it as I have used ‘alpha-skeleton’ in a localOperator with UG before!
Cheers
T
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 <http://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>
On 13/05/2016, 18:22, "dune-pdelab-bounces+t.dodwell=exeter.ac.uk at dune-project.org on behalf of dune-pdelab-request at dune-project.org" <dune-pdelab-bounces+t.dodwell=exeter.ac.uk at dune-project.org on behalf of dune-pdelab-request at dune-project.org> wrote:
>Send dune-pdelab mailing list submissions to
> dune-pdelab at dune-project.org
>
>To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab
>or, via email, send a message with subject or body 'help' to
> dune-pdelab-request at dune-project.org
>
>You can reach the person managing the list at
> dune-pdelab-owner at dune-project.org
>
>When replying, please edit your Subject line so it is more specific
>than "Re: Contents of dune-pdelab digest..."
>
>
>Today's Topics:
>
> 1. Problem with codim 1 iterator and UG Grid (Dodwell, Timothy)
> 2. Re: Problem with codim 1 iterator and UG Grid (Steffen M?thing)
> 3. Re: Problem with codim 1 iterator and UG Grid (Ansgar Burchardt)
>
>
>----------------------------------------------------------------------
>
>Message: 1
>Date: Fri, 13 May 2016 14:40:51 +0000
>From: "Dodwell, Timothy" <T.Dodwell at exeter.ac.uk>
>To: "dune-pdelab at dune-project.org" <dune-pdelab at dune-project.org>
>Subject: [dune-pdelab] Problem with codim 1 iterator and UG Grid
>Message-ID: <FB04C448-EC48-444B-B3F1-F5E46718EECD at exeter.ac.uk>
>Content-Type: text/plain; charset="utf-8"
>
>Hi all,
>
>I am implementing implicit error estimator, which requires me to loop over all entities of codim 1
>
>When I do this, something like
>
> // Make grid using UGGrid
> typedef Dune::UGGrid<dim> GRID;
> GRID grid;
> Dune::GridFactory<GRID> factory(&grid);
> Dune::GmshReader<GRID>::read(factory,gridName);
> factory.createGrid();
> //grid.loadBalance(newcomm);
>
> typedef typename GRID::LeafGridView GV;
> GV gv = grid.leafGridView(); // Get finest grid
>
> typedef typename GV::template Codim< 1 >::Iterator CodimIterator;
>
> const CodimIterator codimEnd = gv.template end< 1 >();
>
> int size = 0;
> for( CodimIterator it = gv.template begin< 1 >(); it != codimEnd; ++it )
> {
> ++size;
> }
>
>I get the following error
>
>
>/home/tjd20/git_workspace/dune/dune/dune-grid/dune/grid/uggrid/uggridleafiterator.hh:141:56: error: invalid conversion from void*' to 'Dune::UG_NS<2>::Entity<1>::T* {aka UG::D2::edge*}' [-fpermissive]
>
> this->virtualEntity_.setToTarget(UG_NS<dim>::succ(this->virtualEntity_.getTarget()), grid_);
>
>Works fine for codim 0 and 2
>
>Anyway ideas?
>
>Cheers
>
>Tim
>
>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<http://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>
>-------------- next part --------------
>An HTML attachment was scrubbed...
>URL: <http://lists.dune-project.org/pipermail/dune-pdelab/attachments/20160513/81761546/attachment-0001.html>
>
>------------------------------
>
>Message: 2
>Date: Fri, 13 May 2016 17:07:15 +0200
>From: Steffen M?thing <steffen.muething at iwr.uni-heidelberg.de>
>To: "Dodwell, Timothy" <T.Dodwell at exeter.ac.uk>
>Cc: dune-pdelab mailing list <dune-pdelab at dune-project.org>
>Subject: Re: [dune-pdelab] Problem with codim 1 iterator and UG Grid
>Message-ID:
> <88F63FF3-8728-48DB-AE15-034F34F681FF at iwr.uni-heidelberg.de>
>Content-Type: text/plain; charset="utf-8"
>
>Hi Tim,
>
>> Am 13.05.2016 um 16:40 schrieb Dodwell, Timothy <T.Dodwell at exeter.ac.uk>:
>>
>> Hi all,
>>
>> I am implementing implicit error estimator, which requires me to loop over all entities of codim 1
>>
>> When I do this, something like
>>
>> // Make grid using UGGrid
>> typedef Dune::UGGrid<dim> GRID;
>> GRID grid;
>> Dune::GridFactory<GRID> factory(&grid);
>> Dune::GmshReader<GRID>::read(factory,gridName);
>> factory.createGrid();
>> //grid.loadBalance(newcomm);
>>
>> typedef typename GRID::LeafGridView GV;
>> GV gv = grid.leafGridView(); // Get finest grid
>>
>> typedef typename GV::template Codim< 1 >::Iterator CodimIterator;
>>
>> const CodimIterator codimEnd = gv.template end< 1 >();
>>
>> int size = 0;
>> for( CodimIterator it = gv.template begin< 1 >(); it != codimEnd; ++it )
>> {
>> ++size;
>> }
>>
>> I get the following error
>>
>> /home/tjd20/git_workspace/dune/dune/dune-grid/dune/grid/uggrid/uggridleafiterator.hh:141:56: error: invalid conversion from void*' to 'Dune::UG_NS<2>::Entity<1>::T* {aka UG::D2::edge*}' [-fpermissive]
>> this->virtualEntity_.setToTarget(UG_NS<dim>::succ(this->virtualEntity_.getTarget()), grid_);
>>
>> Works fine for codim 0 and 2
>
>AFAIK UG doesn?t support codim 1 iterators, it only has codim 1 entities (which you can get by calling subEntity() on
>cells). But I?m not entirely sure?
>
>Cheers
>Steffen
>
>>
>> Anyway ideas?
>>
>> Cheers
>>
>> Tim
>>
>> 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
>> Citations: @Google-Scholar
>> _______________________________________________
>> dune-pdelab mailing list
>> dune-pdelab at dune-project.org
>> http://lists.dune-project.org/mailman/listinfo/dune-pdelab
>
>-------------- next part --------------
>A non-text attachment was scrubbed...
>Name: signature.asc
>Type: application/pgp-signature
>Size: 455 bytes
>Desc: Message signed with OpenPGP using GPGMail
>URL: <http://lists.dune-project.org/pipermail/dune-pdelab/attachments/20160513/f8bdcb6f/attachment-0001.pgp>
>
>------------------------------
>
>Message: 3
>Date: Fri, 13 May 2016 18:06:56 +0200
>From: Ansgar Burchardt <Ansgar.Burchardt at tu-dresden.de>
>To: dune-pdelab at dune-project.org
>Subject: Re: [dune-pdelab] Problem with codim 1 iterator and UG Grid
>Message-ID: <1463155616.18274.1.camel at tu-dresden.de>
>Content-Type: text/plain; charset="utf-8"
>
>Hi,
>
>On Fri, 2016-05-13 at 17:07 +0200, Steffen M?thing wrote:
>> AFAIK UG doesn?t support codim 1 iterators, it only has codim 1
>> entities (which you can get by calling subEntity() on
>> cells). But I?m not entirely sure?
>
>UG only has iterators for elements and vertices.??However you can use
>GeometryGrid to get iterators over all codimensions (it iterates over
>elements and their subentities); at least once [1] gets merged.
>
>Ansgar
>
>? [1] https://gitlab.dune-project.org/core/dune-grid/merge_requests/54
>-------------- next part --------------
>A non-text attachment was scrubbed...
>Name: signature.asc
>Type: application/pgp-signature
>Size: 819 bytes
>Desc: This is a digitally signed message part
>URL: <http://lists.dune-project.org/pipermail/dune-pdelab/attachments/20160513/99ee590e/attachment.pgp>
>
>------------------------------
>
>_______________________________________________
>dune-pdelab mailing list
>dune-pdelab at dune-project.org
>http://lists.dune-project.org/mailman/listinfo/dune-pdelab
>
>
>End of dune-pdelab Digest, Vol 77, Issue 1
>******************************************
More information about the dune-pdelab
mailing list