"I *assume* that your goal is to visualize the data at the actual quadrature points, not on quarters of the elements. In that case, may I suggest that you create a ParticleHandler with particles located at quadrature points (for which there is already a generator function). Each particle would then carry one "property" and you'd just have to assign that property the value of your solution at that point (or whatever other quantity you want visualized at this point). The final step is to use the particles::DataOut class to output the information at the particle locations."
My idea was to actually plot the data at quarters of the elements (this includes the quadrature point, of course, but looks nicer that just a scatter plot -- in my opinion). I have to check how the output looks like when using particles::DataOut class. I just scanned the documentation of that class, which says: 'This class currently only supports witing the particle position and their ID and does not allow to write the properties attached to the particles'. Given that, I can not write the properties (thus the qp values) to file, can I? Best Simon Wolfgang Bangerth <[email protected]> schrieb am Fr., 26. Aug. 2022, 22:09: > On 8/26/22 10:31, Simon wrote: > > > > > > To this end, I think I need a finite element with piecewise constant > > shape functions. Each shape function is '1' at one quarter of the > > reference domain and '0' at the remaining three quarters. > > Or is there an easier way to do this? > > If not, I probably would have to write my own finite element class, > > given that such an element does not exist in dealii? > > It is possible to implement such an element, but like Daniel already > said, it's maybe more work than you want to do. > > I *assume* that your goal is to visualize the data at the actual > quadrature points, not on quarters of the elements. In that case, may I > suggest that you create a ParticleHandler with particles located at > quadrature points (for which there is already a generator function). > Each particle would then carry one "property" and you'd just have to > assign that property the value of your solution at that point (or > whatever other quantity you want visualized at this point). The final > step is to use the particles::DataOut class to output the information at > the particle locations. > > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > -- > 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/2ede856e-6144-caba-4858-7525b03d55cc%40colostate.edu > . > -- 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/CAM50jEvi1C3Ykzm2B41zFn76UJ51LmgbSY6xKkwXtQpNHxGSxg%40mail.gmail.com.
