Dear Prof.  Bangerth 

The definition of  class PointHistory is as follows:

````
 24 class PointHistory : public TransferableQuadraturePointData

 23 {

 22 public:

 21   void

 20   update(std::string name, double solution, std::string scheme = 
"latest");

 19   void

 18   update_independent(std::string name,

 17                      double      solution,

 16                      std::string scheme = "latest");

````

It is used to store my hist variable while assembling the stiffness matrix. 
Thanks for your suggestion, I'll take a look.

Best regards,

Junxiang.


On Friday, January 17, 2025 at 1:06:12 PM UTC+8 Wolfgang Bangerth wrote:

> On 1/16/25 18:40, Junxiang Wang wrote:
> > 
> > const std::vector<std::shared_ptr<PointHistory>> lqph = 
> > quadrature_point_history.get_data(cell);
> > 
> > to record some history variables (e.g., max and min) at quadrature_point.
> > 
> > Is there a way to output the nodal value of the history value? I am 
> currently 
> > using DataPostprocessorVector to handle my output.
>
> I don't know what the PointHistory class does (did you mis-spell the name 
> of a 
> deal.II class?), but if I understand you correctly, then this is the value 
> (or 
> multiple values) of the solution at a specific point. DataOut does not 
> have 
> facilities to output point values itself, but you can work around this by 
> creating a ParticleHandler that has particles at quadrature points (I 
> think 
> either step-70 or step-83) shows this, then transfer your history variable 
> from quadrature points to the particles, and then you can output this 
> information via particles::DataOut. This information can then be 
> visualized 
> alongside the field data you output separately via DataOut.
>
> Best
> W.
>
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: bang...@colostate.edu
> 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 dealii+unsubscr...@googlegroups.com.
To view this discussion visit 
https://groups.google.com/d/msgid/dealii/80b93c98-9bed-4182-93f2-38552da57ab8n%40googlegroups.com.

Reply via email to