[petsc-users] Fieldsplit with LSC for constrained elasticity/poroelasticity?

[Tabrez Ali]

Jed

You're right and it is precisely what I am trying to understand as well. 
I am attaching the comparison to the analytical solution. As you can see 
the agreement is almost perfect, except at the discontinuity (y=0 in the 
analytical solution).

So why am I getting the correct solution when according to true 
residuals I should not?

Tabrez

On 10/28/2014 11:18 AM, Jed Brown wrote:
> Tabrez Ali <stali at geology.wisc.edu> writes:
>
>> Matt and Jed
>>
>> But it does appear to work well in practice (as in I get the correct
>> solution).
> Your "correct" solution has a residual seven orders of magnitude bigger
> than the initial guess.  Does that not give you pause?

-------------- next part --------------
A non-text attachment was scrubbed...
Name: petsc_vs_analytical.pdf
Type: application/pdf
Size: 19410 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141028/fa7c72a7/attachment-0001.pdf>