[Dune] question about JacobianInverseTransposed

Martin Nolte nolte at mathematik.uni-freiburg.de
Mon Aug 4 09:59:34 CEST 2014


PS: Maybe the AffineGeometry is a better place to learn the implementation of a 
DUNE geometry, because it is slimmer. The MultiLinearGeometry is a quite complex 
construct with a few tweaks to gain performance. Some of them might not even be 
necessary with recent compilers. For example, some of the integral_constant 
magic could be removed if we rely on constant propagation. However, this would 
rely on people recognizing this opportunity and not "improving" it without thought.

Best,

Martin

On 08/04/2014 09:48 AM, Aleksejs Fomins wrote:
> Dear Dune,
>
> I would be happy if someone could explain me the initializer of the
> MultilinearGeometry:
>
>
> ---------------------------------------------------------------
>
>    typedef FieldMatrix< ctype, coorddimension, mydimension > Base;
>
>     public:
>       void setup ( const JacobianTransposed &jt )
>       {
>         detInv_ = MatrixHelper::template rightInvA< mydimension,
> coorddimension >( jt, static_cast< Base & >( *this ) );
>       }
>
> ---------------------------------------------------------------
>
> Namely, I am not sure what does "static_cast< Base & >( *this )" do, and
> thus where the inverse of jt is stored?
>
> Thank you,
> Aleksejs
>
>
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
>

-- 
Dr. Martin Nolte <nolte at mathematik.uni-freiburg.de>

Universität Freiburg                                   phone: +49-761-203-5630
Abteilung für angewandte Mathematik                    fax:   +49-761-203-5632
Hermann-Herder-Straße 10
79104 Freiburg, Germany



More information about the Dune mailing list