Hi Paras, Thanks for the reply!
``` > dealii::FEValuesExtractorScalar pressureExtractor(2); //assuming you have > three scalar fields and pressure is the last one > dealii::Vector<double> solVecPerCell(nDOFsPerCell); > feValues[pressureExtractor].get_function_values(totalSolution, > solVecPerCell); // here only the values corresponding to pressure are > filled. > ``` > Suppose I have 2 DG elements of order 3, i.e., 16 DoFs per element. U, v, p -> 2*(16)*3 = 96 DoF in totalSolution. It's my understanding that this will extract the pressure values from totalSolution into the *cell-local* Vector<double> solVecPerCell of length 16 (assuming FEValues was instantiated with a quadrature rule that contained the nodal locations, rather than the usual quadrature points). However, this isn't what I want to do. I want to take all 32 of the p nodal values and add them to another field with a DoFHandler instantiated using a single FE_DGQ element, rather than a system--therefore with a solution Vector<double> containing 32 DoF. Maybe the way to go is a BlockVector object as you alternatively suggest, but my problem doesn't inherently have blocks. I'm simply looking to add, for example, the pressure field to another scalar DG field unrelated to the system; I'd rather keep the FESystem organized by (u,v,p) on each cell. Best, Corbin -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/49866be9-6017-4dde-8f10-8b421863e6a4n%40googlegroups.com.
