On Oct 15, 2008, at 2:33 AM, Norbert Stoop wrote: Parallel is where the issue is... and is why I haven't enabled that capability yet. I haven't yet done the work to make EpetraMatrix work in parallel like I did for EpetraVector.
But if you have it working in serial... by all means send a patch. > I understand that Trilinos support is unofficial and work in > progress. I > am still new to libmesh, and Trilinos in particular, but will invest > some more hours into this problem. Not sure if that will be of any > help, > though... Any help is appreciated... I'm currently being pulled in about 7 different directions.... so even if small progress is being made in this direction... that will be better than none. > What puzzles me is that the error is already fully present after the > first nonlinear iteration, ie. after solving the linearized system. Yes... if there is a problem in parallel then you probably aren't getting the correct summations of off-processor contributions when you do a residual evaluation.... so you would see it _immediately_. Your very first residual evaluation would be wrong. Like I say though.... that shouldn't be broken anymore (for pure matrix free). Are you sure you've updated your libMesh? There was a small bug where closing the EpetraVector was not calling GlobalAssemble... but I think I committed the fix for that a little bit ago. Derek ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
