Thanks for the help. The actual systems are slightly different so I am willing to pay for triplicate matrix storage. However, I was hoping to save the DofMap storage if possible. That does not seem easy, so I will have to see how close I come to running out of memory. The information on the initial guess will be useful. Adam
________________________________ From: Kirk, Benjamin (JSC-EG) [mailto:[EMAIL PROTECTED] Sent: Tue 9/30/2008 8:45 PM To: Adam Arbree; libmesh-users@lists.sourceforge.net Subject: Re: [Libmesh-users] (no subject) By default if you declare three systems you will triple the storage, even if the fe types and everything are identical. It would be pretty straightforward to force the matrix storage to be shared by all three systems, but you will still have three times the DofMap storage (and restriction/prolongation overhead in the case of amr). But yes, if you do nothing special you will get triplicate maps which should stay consistent. One of the more hackish ways to handle the situation would be to only declare one system and use a static counter inside your assemble() function to actually assemble (and thus solve) the different systems. Roy or John?? As for the other issue, system.solution is uised as the initial guess to the iterative siolver, so if you copy the contents between systems you should achieve the desired effect. -Ben ----- Original Message ----- From: Adam Arbree <[EMAIL PROTECTED]> To: libmesh-users@lists.sourceforge.net <libmesh-users@lists.sourceforge.net> Sent: Tue Sep 30 18:01:40 2008 Subject: [Libmesh-users] (no subject) I have an application where I need to solve three nearly identical systems. For each system, the matrix entries and the rhs will be slightly different, but the mesh and finite element type and degree will be identical. I have three questions: 1) I would like to reuse the DofMap between each of these systems to save memory (and any other memory that is possible, I am not sure what the overhead per system is); 2) if it is not convenient to avoid duplication of the map memory, I would like to construct the system such that at least the three maps are the same (will this happen automatically?); 3) can I specify an initial guess for the solution when using an iterative solver so that I can use the solution of one system as the initial guess for the another. Thanks you very much Adam Arbree ------------------------------------------------------------------------- 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 Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------- 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 Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
