Hello, Roy, Thank you for your reply.
Because I want to use Jacobian_free Newton_Krylov method to solve the problem, if I setup two EquationSystems for each domain, I will get two rhs and connect them, with a matrix_free preconditioner I can get the full implicit solver, is that advisable? If it is OK, the only thing I need to do is write a suitable send_list for parellel like the unconformed mesh? By the way, because my implicit discontinuous code need so much memory, so I want to use Jacobian_free Newton_Krylov method and write a matrix free preconditioner, can you give me some suggestions about this? Do I need to write a JFNK_systems inherit from the class "explicit_system"? Thank you very much! Best regards! Luyi 在 2008-07-31四的 08:09 -0500，Roy Stogner写道： > moment I think the answer is "not efficiently". Although our > meshes can be flagged with subdomain ids, our DofMap still assumes > that every subdomain in the mesh has the same variables defined on it. > What we need is an API for the user to say "this variable only exists > on subdomains 1, 5, and 13", as well as the DofMap code to respect > those requests. Unless you want to dig into the library source code > and add that feature (for which a patch would be much appreciated, > naturally), I think the best you can do for a fully coupled solve is > define all variables everywhere, then manually restrict T to be 0 on > domain A and all the other variables to be 0 on domain B. -- Gas Turbine Research Center email: [EMAIL PROTECTED] ------------------------------------------------------------------------- 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
