<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Nov 15, 2013 at 5:31 PM, Geoffrey Irving <span dir="ltr"><<a href="mailto:irving@naml.us" target="_blank">irving@naml.us</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On Fri, Nov 15, 2013 at 3:21 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>

> On Fri, Nov 15, 2013 at 5:09 PM, Geoffrey Irving <<a href="mailto:irving@naml.us">irving@naml.us</a>> wrote:<br>
>><br>
>> Is the second equation the weak form of incompressibility?  If so, why<br>
>> isn't q u?  It's used as q several places throughout the code, so it<br>
>> seems more likely to be my misunderstanding than a typo.<br>
>><br>
>>   < \nabla v, \nabla u + {\nabla u}^T > - < \nabla\cdot v, p > + < v, f > = 0<br>
>>   < q, \nabla\cdot v >                                                     = 0<br>
><br>
> u, v \in V and p, q \in P, and hopefully \div V = P for stability.<br>
<br>
</div>Are you using the same variables as<br>
<br>
    <a href="http://fenicsproject.org/documentation/dolfin/dev/python/demo/pde/stokes-iterative/python/documentation.html" target="_blank">http://fenicsproject.org/documentation/dolfin/dev/python/demo/pde/stokes-iterative/python/documentation.html</a><br>

<br>
If so, maybe that line should be<br>
<br>
    < q, \nabla\cdot u ><br>
<br>
I.e., "u" instead of "v"?</blockquote><div><br></div><div>Damn, you are right.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
Geoffrey<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>