Re: [deal.II] Re: Outputting stress tensor at each quadrature point

2022-05-18 Thread BBah
Hello Mr. Bangerth, i am just trying to put out the norm of the stresses like in step 18. Is there a way to save my stress tensor for each cell and put it out ? In step 18 there is a loop through each cell and the locally owned are used to put out the stresses. Is there a way to do the same wit

Re: [deal.II] Re: Outputting stress tensor at each quadrature point

2022-05-18 Thread Wolfgang Bangerth
BBah, I have to admit that I don't understand what precisely your question is -- your first email does not contain enough information what you have already tried, or what exactly you want to do. But regarding what you have below: segmentation faults are the best kind of error because you ca

[deal.II] Re: Outputting stress tensor at each quadrature point

2022-05-18 Thread BBah
My PointHistory class does include old_stress but running this code leads to a segmentation fault: template class PointHistory { public: PointHistory() : old_stress(SymmetricTensor<2, dim,NumberType>()) {} virtual ~PointHistory() {} // The first function is use

[deal.II] Outputting stress tensor at each quadrature point

2022-05-18 Thread BBah
Hello everybody, i am trying to output my stresses by using a similar approach like in Step 18: *Vector norm_of_stress(triangulation