<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body ocsi="0" fpstyle="1" bgcolor="#FFFFFF">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Tobias,<br>
<br>
it works perfectly!<br>
Thank you very much!<br>
<div><br>
Andrea<br>
<div style="font-family:Tahoma; font-size:13px">__________________________________________________________<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.sacconi11@imperial.ac.uk<br>
</div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF574474"><font color="#000000" face="Tahoma" size="2"><b>From:</b> Tobias Malkmus [tomalk@mathematik.uni-freiburg.de]<br>
<b>Sent:</b> 11 April 2013 19:25<br>
<b>To:</b> Sacconi, Andrea<br>
<b>Cc:</b> dune-fem@dune-project.org<br>
<b>Subject:</b> Re: [dune-fem] How to get DoFs for Dirichlet nodes with vector-valued unknowns<br>
</font><br>
</div>
<div></div>
<div>
<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 type="cite"><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}
-->
BODY {direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;}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" target="_blank">a.sacconi11@imperial.ac.uk</a><br>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader" target="_blank"></fieldset> <br>
<pre>_______________________________________________
dune-fem mailing list
<a class="moz-txt-link-abbreviated" href="mailto:dune-fem@dune-project.org" target="_blank">dune-fem@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune-fem" target="_blank">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>
</div>
</div>
</div>
</body>
</html>