[Dune] Question about the GenericReferenceElements

Oliver Sander sander at mi.fu-berlin.de
Sat Oct 23 11:03:26 CEST 2010


Hi Martin!
Thanks for the information.  I think discussion in Münster is a good
idea.  There are actually more issues with the reference elements
that need discussion (a FlySpray entry will appear in a few days).

best,
Oliver

Am 22.10.2010 20:22, schrieb Martin Nolte:
> Hi Oliver,
>
> when implementing the generic reference element, I implemented all methods of
> the original ReferenceElement class. Especially, this included the method
>
> global ( const FieldVector<  ctype, dim-codim>  &local, int i, int c ) const.
>
> When I saw the redundancy, I simply added another method, dropping the redundant
> argument (for convenience).
>
> Back then, I somehow assumed that the interface for the ReferenceElement has
> been thoroughly discussed. Over time, this impression has changed. Since the old
> ReferenceElement is now gone, I would be very happy to discuss what methods
> should be present on the GenericReferenceElement. In particular, I would very
> much like to get rid of the (i,c,ii,cc) methods. In my mind it is (c,i,cc,ii)
> and, since the types don't differ, I frequently get weird errors from this.
>
> I would very much appreciate, if an "interface" for the reference element could
> be discussed on the developer meeting in Münster.
>
> Best,
>
> Martin
>
>
> On 10/22/2010 02:43 PM, Oliver Sander wrote:
>    
>> Dear Dune!
>> I just had a look at the GenericReferenceElement class.
>> I learned that said class has methods
>>
>> template<  int codim>
>> FieldVector<  ctype, dim>
>> global( const FieldVector<  ctype, dim-codim>  &local, int i )
>>
>> and
>>
>> template<  int codim>
>> FieldVector<  ctype, dim>
>> global( const FieldVector<  ctype, dim-codim>  &local, int i, int c ) const
>>
>> Their documentations appear to be identical, with the exception of a
>> note saying "The runtime argument c is redundant and must equal codim"
>> for the second method.  And indeed the method throws an exception
>> if c!=codim.
>>
>>
>> Why do we have this seemingly pointless method as part
>> of our grid interface?
>>
>> Best,
>> Oliver
>>
>> -- 
>>
>> ************************************************************************
>> * Oliver Sander                ** email: sander at mi.fu-berlin.de        *
>> * Freie Universität Berlin     ** phone: + 49 (30) 838 75348           *
>> * Institut für Mathematik      ** URL  : page.mi.fu-berlin.de/~sander  *
>> * Arnimallee 6                 ** -------------------------------------*
>> * 14195 Berlin, Germany        ** Member of MATHEON (www.matheon.de)   *
>> ************************************************************************
>>
>>
>>
>> _______________________________________________
>> Dune mailing list
>> Dune at dune-project.org
>> http://lists.dune-project.org/mailman/listinfo/dune
>>      
>    


-- 
************************************************************************
* Oliver Sander                ** email: sander at mi.fu-berlin.de        *
* Freie Universität Berlin     ** phone: + 49 (30) 838 75348           *
* Institut für Mathematik      ** URL  : page.mi.fu-berlin.de/~sander  *
* Arnimallee 6                 ** -------------------------------------*
* 14195 Berlin, Germany        ** Member of MATHEON (www.matheon.de)   *
************************************************************************





More information about the Dune mailing list