Evan, > I'm been working on trying to formulate a local DG scheme in deal.ii for a > reaction-diffusion system. The general formulations in the examples for > the time dependent cases though -- for example in step-33 -- are sort of CG > based implicit solvers, while I want to build a local DG based explicit > solver. Moving from the continuous to discontinuous elements seems fairly > straight forward, but what has been causing me some severe problems is > constructing a system explicitly (and locally) using the usual > system_matrix formulation, or otherwise. My feeling is this should be > fairly easy, but a few words of initial advice might go a long way at this > point.
I'm not entirely sure what exactly you are trying, but when you say "explicit" I assume you want to work matrix-free, right? In that case, you may want to consider step-37, which you can get by checking out the current development version of deal.II and which is also here: http://www.dealii.org/developer/doxygen/deal.II/step_37.html It shows how to solve linear systems without assembling a global matrix. If that isn't what you were looking for, can you explain a bit more what exactly the problem is? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
