<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Hi Andrea<br>
<br>
try this :<br>
<br>
const SubEntityIteratorType endFaceIt = pointSet.template
endSubEntity<1>( localFaceNumber );<br>
for( SubEntityIteratorType faceDofIt = pointSet.template
beginSubEntity<1>( localFaceNumber );<br>
faceDofIt != endFaceIt; ++faceDofIt )<br>
{<br>
for( int i = 0; i < RangeType::dimension; ++i )<br>
{<br>
const int index = (*faceDofIt ) * RangeType::dimension + i;<br>
local[ index ] = ...<br>
}<br>
}<br>
<br>
where local is a local Lagrange discrete function.<br>
In the same way you can modify the local stiffness matrix in the
index-th row/column.<br>
<br>
<br>
<br>
Best Tobias<br>
<br>
<br>
On 04/11/2013 07:12 PM, Sacconi, Andrea wrote:<br>
</div>
<blockquote
cite="mid:EC7DA3C0917A914BAFB620BB60189EBD1362F5F2@icexch-m2.ic.ac.uk"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<style id="owaParaStyle" type="text/css">
<!--
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
-->
P {margin-top:0;margin-bottom:0;}</style>
<div style="direction: ltr;font-family: Tahoma;color:
#000000;font-size: 10pt;">
<div style="direction:ltr; font-family:Tahoma; color:#000000;
font-size:10pt">
<div style="direction:ltr; font-family:Tahoma; color:#000000;
font-size:10pt">Hi all,<br>
<br>
I have a quick question for you.<br>
<br>
Suppose you need to find the global DoFs associated to the
Dirichlet nodes. The purpose is, of course, to modify the
rows of the stiffness matrix in the standard way, as well as
the corresponding components on the rhs.<br>
When the unknown is a scalar variable, I don't have
problems: I loop over the entities, extract the
LagrangePointSetType set for each entity, use
FaceDofIteratorType, ...<br>
<br>
But when the unknown is a vector-valued function, how do I
proceed? I have tried to do exactly the same, and I noticed
that the same node can refer to two different DoFs.<br>
The problem is, how do I prescribe the value for the
function? Since it's a vector-valued function, where do I
prescribe the first component and where the second?
<br>
In other words, how is the DoFs set created?<br>
<br>
Thanks in advance for your help!<br>
<br>
Cheers,<br>
Andrea<br>
<div><br>
<div style="font-family:Tahoma; font-size:13px"><br>
__________________________________________________________<br>
<br>
Andrea Sacconi<br>
PhD student, Applied Mathematics<br>
AMMP Section, Department of Mathematics, Imperial
College London,<br>
London SW7 2AZ, UK<br>
<a class="moz-txt-link-abbreviated" href="mailto:a.sacconi11@imperial.ac.uk">a.sacconi11@imperial.ac.uk</a><br>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
dune-fem mailing list
<a class="moz-txt-link-abbreviated" href="mailto:dune-fem@dune-project.org">dune-fem@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune-fem">http://lists.dune-project.org/mailman/listinfo/dune-fem</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Tobias Malkmus <tomalk at mathematik.uni-freiburg.de>
Mathematisches Institut Tel: +49 (0) 761 203 5651
Abt. für Angewandte Mathematik Universität Freiburg
Hermann-Herder-Str. 10
79104 Freiburg
</pre>
</body>
</html>