> I did not know that resize can take a second argument. I was instead > searching for a resize function for FullMatrix.
You'd use FullMatrix::reinit > I am using subdivided_hyper_cube to make the grid. Then I know the order of > the face indices and cell numbering. I dont plan to do any refinement. So I > first compute flux using some Riemann solver at all cell faces and store > them. Then while assembling the rhs, I am using these fluxes. Basically I > dont want to compute flux at a face twice. That's the classical tradeoff: you want to use memory to save CPU time. That's sometimes possible, but not always. You may have a cases where you can't. If you want to see how to avoid computing it twice, take a look at step-12 (I believe). W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
