On Fri, Nov 11, 2011 at 2:58 PM, Juha Jäykkä <span dir="ltr"><<a href="mailto:juhaj@iki.fi">juhaj@iki.fi</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
> We pass in the time of the stage when your RHSFunction is called. If you<br>
<br>
Right, how can I have missed that?!? Thanks.<br>
<br>
> haven't eliminated boundary conditions, your system is likely<br>
> differential-algebraic, in which case you can use the IFunction interface<br>
<br>
Eliminated? I am not sure what you mean. My boundaries are "periodic", but<br>
thanks to the nature of my problem, they are periodic only up to a gauge<br>
transformation. I need to perform that gauge transform at each time step.<br>
<br>
I would be more than happy not to have to do that, but I cannot really see a<br>
way out of it, so, if you have any ideas, I am very interested in hearing<br>
them.<br>
<br>
> improved in petsc-dev, so consider using it if you want the latest<br>
<br>
Thanks for the tip, but I am stuck at 3.1 because I need TAO, too and it still<br>
does not support 3.2.<br></blockquote><div><br></div><div>What exactly are you using TAO to do?</div><div><br></div><div> Thanks,</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;">
Cheers,<br>
<font color="#888888">Juha<br>
</font></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<br>