On Mon, Dec 8, 2014 at 5:10 PM, David Knezevic <david.kneze...@akselos.com>
wrote:
>
> John suggested that I put this code into my fork on github so that it's
>> easier for others to view and comment on. I'll take care of that shortly.
>>
>
>
> OK, it's up on my fork on github, see here
> <https://github.com/dknez/libmesh/tree/systems_of_equations_ex7/examples/systems_of_equations/systems_of_equations_ex7>
> .
>
I believe there was an error in the formulation I was using yesterday (I
had an extra det(F) which wasn't supposed to be there). I've updated the
code on my fork, and it looks like it's working well now. (Note that it
converges well as long as I use the PETSc command line option
"-snes_linesearch_type basic".)
David
------------------------------------------------------------------------------
Download BIRT iHub F-Type - The Free Enterprise-Grade BIRT Server
from Actuate! Instantly Supercharge Your Business Reports and Dashboards
with Interactivity, Sharing, Native Excel Exports, App Integration & more
Get technology previously reserved for billion-dollar corporations, FREE
http://pubads.g.doubleclick.net/gampad/clk?id=164703151&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
