On Mon, Jul 26, 2010 at 11:41 AM, David Van Wie <[email protected]> wrote: > Hi all, > > I'm using libMesh as part of a tool that solves some E+M problems. One > task we have is to find the current distribution through different > conductors. Overall, this task is pretty similar to example 14 from the > website. > > Many of the conductors we're studying are small -- on the order of a 10s > of microns wide. We had been performing our calculations in units of > microns, but recently changed to meters and started to see crashes in > equation_systems.reinit() when running in parallel. These crashes happen > in both our code and in example 14, though in slightly different places > (tracebacks follow). The attached input file will cause example 14 to > crash.
I'm just guessing but if you're doing refinement, you're probably failing some of our built-in geometric tests which run when new points are inserted. These tests use a default absolute tolerance that may be too small for the Mesh you are working with. > Any thoughts on the best approach to take here? Solve a suitably non-dimensionalized version of your problem instead. This may make other parameters, say, conductivity, much larger or smaller but will keep the geometric tests normalized to 1. -- John ------------------------------------------------------------------------------ The Palm PDK Hot Apps Program offers developers who use the Plug-In Development Kit to bring their C/C++ apps to Palm for a share of $1 Million in cash or HP Products. Visit us here for more details: http://ad.doubleclick.net/clk;226879339;13503038;l? http://clk.atdmt.com/CRS/go/247765532/direct/01/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
