On Tue, Mar 5, 2013 at 10:23 AM, Robert <[email protected]> wrote:

> Hello,
>
> I am currently catching errors due to element inversion (det F < 0)
> during assembly. If such an error occurs I want to reduce the timestep
> and try again.
> This works well for a single process. When trying to run the whole
> thing in parallel I am running into trouble, because usually just one of
> the processes has an inverted element and the whole thing deadlocks.
> Does anyone know how to handle this situations?
>

In general you have to catch the exception, do a global communication step
(e.g. MPI_Reduce) to figure out if any processor has an exception, and if
so, raise the exception on all processors.

This is something that should probably be handled at the library level.  No
doubt Roy has some ideas how to do it (or has already implemented it?).

-- 
John
------------------------------------------------------------------------------
Everyone hates slow websites. So do we.
Make your web apps faster with AppDynamics
Download AppDynamics Lite for free today:
http://p.sf.net/sfu/appdyn_d2d_feb
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to