[petsc-users] Modified Taylor-Hood elements with piece-wise constant pressure for Stokes equation

[Jed Brown]

Justin Chang <jychang48 at gmail.com> writes:

> I originally solve that example problem using LU. 

I'd like to learn why LU didn't notice that the system is singular.
(The checks are not reliable, but this case should be pretty obviously
bad.)

> But when I solve this one:
>
> http://fenicsproject.org/documentation/dolfin/1.5.0/python/demo/documented/stokes-iterative/python/documentation.html
>
> By simply running their code as is for TH and adding the one like I
> mentioned for MTH, I get the following outputs when I pass in -ksp_monitor
> -ksp_view and -log_summary

Thanks, Justin.  Could you please run with

 -ksp_type gmres -ksp_gmres_restart 1000 -ksp_monitor_true_residual -ksp_monitor_singular_value

Also, can you run a tiny problem size (like less than 1000 dofs) with

  -pc_type svd -pc_svd_monitor

Thanks.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150602/261081af/attachment.pgp>