<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Nov 17, 2013 at 1:42 PM, Subramanya Sadasiva <span dir="ltr"><<a href="mailto:potaman@outlook.com" target="_blank">potaman@outlook.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
I am trying to solve a degenerate cahn hilliard inequality using petsc_DM and virs through libmesh. The problem is that the jacobian for this system has a nearly singular block owing to the \delta t term multiplying the K<br>

<br>
[  K + M             M  ] [ \delta phi]  =  -R_phi<br>
[    M            K deltat] [\delta mu ]  =  -R_mu<br>
<br>
This leads to the unfortunate situation that the time stepping actually fails for very small timesteps.  a solution that seems possible is to add a constant mass matrix to the K delta t, making the system<br>
<br>
[  K + M             M          ] [ \delta phi]  =  -R_phi<br>
[    M            K deltat +M] [\delta mu ]  =  -R_mu + M \delta mu<br>
<br>
And using the \delta mu from the previous timestep. I was just wondering if there is any way to get the pre-line search value of the \delta mu?  Is it okay to get the \alpha from the linesearch and divide the previous change by the \alpha?<br>
</blockquote><div><br></div><div>Why wouldn't you just scale the problem (Jacobi)?</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">

I hope I am making sense here.<br>
Thanks,<br>
Subramanya<br>
<br>
<br>
</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>