Yep,<br><br>I have been looking at it and have used it already quite extensively. Up to now, for stability analysis, I have been using snapshots taken from the Linearized DNS and then used them to form a Krylov subspace for an Arnoldi algorithm. However, what I would like to do know is to use Nek as blackbox for matrix-vector product only. In the perturbation mode, eventhough the Jacobian matrix is not explicitely formed, the equations still read du/dt = J * u<br>

<br>What I am interested in is, given u(0), to output only:<br><br>u(1) = J * u(0)<br>u(2) = J * u(1)<br>...<br>u(N) = J * u(N-1)<br><br>Not sure though if that's clear enough?<br><br><br><br><div class="gmail_quote">

On 21 October 2011 12:37,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<br>
Jean-Christophe,<br>
<br>
Have you already looked at the perturbation code in nek?<br>
It already does these things.<br>
<br>
Also, it's not clear that you want everything in strong<br>
form, particularly the 2nd-order spatial terms, since<br>
the solution is only C0?<br>
<br>
Paul<div class="im"><br>
<br>
<br>
On Fri, 21 Oct 2011, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Hi Nek's,<br>
<br>
Given a random initial vector u and the Jacobian matrix of the linearized<br>
Navier-Stokes, I would like to calculate the matrix-vector product J * u. In<br>
Nek, that would consists in calculating only the right hand side of the<br>
equations, that is:<br>
<br>
RHS = ( -U.grad(u) - u.grad(U) - grad(p) + 1/Re Lap(u)    ;    div(u) )<br>
<br>
I have a few questions on how to compute those different terms in userchk<br>
(Nek is used in post-process mode):<br>
<br></div>
  - I would tend to use convop() for the convection terms, however I'm not<div class="im"><br>
  sure whereas it would output the strong or weak form (and I need the strong<br>
  one obviously)?<br></div>
  - Same question regarding grad(p) and div(u) if I use opgradt() and<div class="im"><br>
  opdiv() to calculate them?<br>
<br>
Last but not least, for the Laplacian I've read a previous mail where it is<br>
said that it can be computed using a sequence of wgradm1() and vec_dssum()<br>
so I guess I kown how to compute at the strong form for it.<br>
<br>
Sincerely,<br>
-- <br>
Jean-Christophe<br>
<br>
</div></blockquote>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Jean-Christophe<br>