If you are working only on Cartesian grids, these matrices will be same for all
cells, except maybe for scalar factors. In this case, you can first compute
them on a reference cell. For assembly, you loop over cells and perform a small
matrix-vector product on each cell. Similar things can be do
Hello all,
I am a beginner in dealii. I want to solve a linear, transient advection
equation explicitly in two dimensions using DG. The resulting discrete
equation will have a mass matrix as the system matrix and a sum of terms
which depend on previous solution (multiplied by mass, differentia