[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