Hello I found that TransientExplicitSystem does not have any mass matrix. I wrote my explicit dg code with TransientLinearImplicitSystem and it works. However I end up computing mass matrix every time step and I am not able to solve the matrices at element level.
I want to use the TransientExplicitSystem. I want to store the mass matrix as std::vector< DenseMatrix<Real> > with one small mass matrix per element. I can then solve at element level using cholesky_solve function available in DenseMatrix. But where do I store this matrix ? Is it good to store it in the parameters and then access it in the assemble function ? Thanks praveen On Sun, Aug 25, 2013 at 1:38 PM, Praveen C <[email protected]> wrote: > Hello > I am writing an explicit RKDG code using libmesh. Having read some > examples, I am thinking of using a TransientExplicitSystem. Is this the > correct approach ? > > Since it is explicit, I only need to assemble mass matrix once. In the > time stepping only the right hand side needs to be assembled. > > Being DG, the mass matrix is block diagonal and can be solved per element. > Its not clear how to do this. Can you point to some existing examples ? > > Also I want to use a multi-stage RK like the 3-stage SSP RK scheme of > Shu-Osher. > > Thanks > praveen > ------------------------------------------------------------------------------ Learn the latest--Visual Studio 2012, SharePoint 2013, SQL 2012, more! Discover the easy way to master current and previous Microsoft technologies and advance your career. Get an incredible 1,500+ hours of step-by-step tutorial videos with LearnDevNow. Subscribe today and save! http://pubads.g.doubleclick.net/gampad/clk?id=58040911&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
