<div dir="ltr"><div dir="ltr">On Thu, Oct 27, 2022 at 1:20 PM Yann Jobic <<a href="mailto:yann.jobic@univ-amu.fr">yann.jobic@univ-amu.fr</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">There is this factor also in front of the rectangle, which should 1/4 <br>
for a unit rectangle of dim 2, but it's also 0.5.<br>
</blockquote><div><br></div><div>Hi Yann,</div><div><br></div><div>The reference element in Petsc is [-1, -1] to [1, 1]. This matches FEniCS/Firedrake and also makes</div><div>it easier to use orthogonal polynomials.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thanks,<br>
Yann<br>
<br>
<br>
Le 10/27/2022 à 6:23 PM, Mark Adams a écrit :<br>
> Area of a unit triangle is (1x1)/2 I would guess.<br>
> <br>
> On Thu, Oct 27, 2022 at 11:32 AM Yann Jobic <<a href="mailto:yann.jobic@univ-amu.fr" target="_blank">yann.jobic@univ-amu.fr</a> <br>
> <mailto:<a href="mailto:yann.jobic@univ-amu.fr" target="_blank">yann.jobic@univ-amu.fr</a>>> wrote:<br>
> <br>
>     Hello,<br>
> <br>
>     I'm trying to understand how to use the lowlevel Petsc Finit Element<br>
>     Framework. I've got few questions about it.<br>
> <br>
>     1) I'm checking the determinant of the Jacobian transformation from<br>
>     real<br>
>     space to the parametric one. The source code is in :<br>
>     <a href="https://petsc.org/main/src/dm/impls/plex/plexgeometry.c.html#DMPlexComputeTriangleGeometry_Internal" rel="noreferrer" target="_blank">https://petsc.org/main/src/dm/impls/plex/plexgeometry.c.html#DMPlexComputeTriangleGeometry_Internal</a> <<a href="https://petsc.org/main/src/dm/impls/plex/plexgeometry.c.html#DMPlexComputeTriangleGeometry_Internal" rel="noreferrer" target="_blank">https://petsc.org/main/src/dm/impls/plex/plexgeometry.c.html#DMPlexComputeTriangleGeometry_Internal</a>><br>
>     I don't quite understand the 0.5 factor in front of it.<br>
>     Geometrically, it should be the ratio of the area/volume of the reel<br>
>     element over the referenced one. It's not the case here no ?<br>
>     It seems that, for other elements, this factor still applies.<br>
>     What this factor represent ?<br>
> <br>
>     2) I only checked triangle elements here. The ordering of the test<br>
>     functions and its derivatives is different it seems, with also the 0.5<br>
>     factor. Am i mistaken ?<br>
> <br>
>     I could not find any informations about it in the user-doc. Maybe it's<br>
>     written somewhere else ?<br>
> <br>
>     Many Thanks,<br>
> <br>
>     Best regards,<br>
> <br>
>     Yann<br>
> <br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>