On Tue, Apr 28, 2015 at 6:09 PM, Xujun Zhao <[email protected]> wrote:
> Hi folks, > > In my problem, I would like to assemble my global matrix and rhs vector > separately, because at each time step, my global matrix doesn't change. > What I need to do is to simply update my rhs vector and resolve it with my > solver. Is there a way to do this in libMesh? I noticed when we > attach_assemble_function, both the matrix and rhs vector are reconstructed. > Have a look at System::assemble_before_solve. If you set that to false then you can call your own assembly functions whenever you like, and the assembly function won't automatically get called when you call solve. This way you could implement an assembly function that has bool arguments assemble_matrix and assemble_rhs and call this whenever you need to with the bools set appropriately. David ------------------------------------------------------------------------------ One dashboard for servers and applications across Physical-Virtual-Cloud Widest out-of-the-box monitoring support with 50+ applications Performance metrics, stats and reports that give you Actionable Insights Deep dive visibility with transaction tracing using APM Insight. http://ad.doubleclick.net/ddm/clk/290420510;117567292;y _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
