<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hi.<br>
You should check out the definition of the 3 different geometries
of the intersection in detail.<br>
<br>
Only intersection.geometry() maps into real space. <br>
geometryInInside and geometryInOutsdide map between reference
elements. They map from the<br>
reference element of the intersection into the reference element
of the inside/ouside entity. So<br>
geometryInInside.global(x) is in fact a local coordinate within
the reference element of the inside entity.<br>
<br>
If the global geometry is changed like in GeometryGrid then only
the geometries into real space need to be changed<br>
the embeddings between different reference elements i.e.
geometryInInside do not have to be touched.<br>
<br>
Best<br>
Andreas<br>
<br>
On 17/09/14 16:38, Aleksejs Fomins wrote:<br>
</div>
<blockquote cite="mid:5419AB0E.7070206@lspr.ch" type="cite">
<pre wrap="">Dear Duners,
I started dissecting the geometrygrid/intersection.hh today, and I am a
bit puzzled how it works. Below I will provide my questions which will
correspond to certain lines in the code
Line 131:
FieldVector< ctype, dimension > x( geometryInInside().global( local ) );
Where can I find the implementation of
hostIntersection().geometryInInside() ?
--------------------------------------------------------
Line 132:
const typename ElementGeometryImpl::JacobianInverseTransposed &jit =
insideGeo_.jacobianInverseTransposed( x );
We have just defined x to be a global coordinate in the line 131,
however in the element geometry implementations known to me (Affine and
Multilinear) the jacobianInverseTransposed method is only defined for
local coordinates. How does this work?
--------------------------------------------------------
Line 136/137:
jit.mv( refNormal, normal );
normal *= jit.detInv();
We take a unit vector refNormal, multiply it by a matrix, and then
multiply it by the value of the integration element. Why do we not
normalize the vector normal after matrix multiplication?
Regards,
Aleksejs
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Dune mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org">Dune@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a>
</pre>
</blockquote>
<br>
</body>
</html>