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

2022-05-23 Thread BBah
Hello Mr. Bangerth, thank you for the response. As mentioned in my last reply I was able to solve the issue. Thank you again for being helpful to the deal ii community and responding to our questions. Best regards Wolfgang Bangerth schrieb am Freitag, 20. Mai 2022 um 18:29:35 UTC+2: > On 5/20/2

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

2022-05-20 Thread Wolfgang Bangerth
On 5/20/22 01:20, BBah wrote: Hello Mr.Bangerth. Thank you very much for your input and your help. I did implement the same loop actually in my code but i am getting a segmentation fault but i guess this is a minor issue that i cannot find. Maybe because it is my first code in c++ and using dea

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

2022-05-20 Thread BBah
I was able to solve this issue by using this in the output_results void: FEValues fe_values(fe, qf_cell, update_JxW_values); *Vector norm_of_stress(triangulation.n_active_cells());* *{* *for (auto &cell : triangulation.active_cell_iterators()){* * fe_values.reinit(cel

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

2022-05-20 Thread BBah
Hello Mr.Bangerth. Thank you very much for your input and your help. I did implement the same loop actually in my code but i am getting a segmentation fault but i guess this is a minor issue that i cannot find. Maybe because it is my first code in c++ and using deal ii. My debugger seems to no

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

2022-05-20 Thread BBah
My debugger seems to not work properly i do not know why. Maybe you can find the issue if i post some of the code: This is my PointHistory class with old_stress included: * template * * class PointHistory* * {* * public:* *PointHistory()* *: old_stress(SymmetricTensor<2, dim,NumberTyp

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

2022-05-20 Thread BBah
Like u mentioned i should be able to use the same loop without the workstream::run function but i do not know why this is not working. I guess it is a really small issue but i cannot find it. BBah schrieb am Freitag, 20. Mai 2022 um 09:03:00 UTC+2: > Hello Mr.Bangerth. I already implemented the

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

2022-05-20 Thread BBah
Hello Mr.Bangerth. I already implemented the old stress in my PointHistory class to get the stresses form material->get_tau() (like in step 44). This is my loop right now which actually the same as in step 18: Vector norm_of_stress(triangulation.n_active_cells()); for (auto &cell : tr

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

2022-05-19 Thread Wolfgang Bangerth
On 5/19/22 00:47, BBah wrote: 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

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

2022-05-19 Thread BBah
One more hint maybe. I have just modified the code " Quasi-Static Finite-Strain Compressible Elasticity " for my purposes. If anybody know how to put out stresses from this code i would be able to implement it into my code too. BBah schrieb am Donnerstag, 19. Mai 2022 um 12:17:07 UTC+2: > I ju

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

2022-05-19 Thread BBah
I just can't find any sample or example how to put out my stress tensor. It should be pretty easy because my stress is stored for each cell in my PointHistory class but how to access these stresses and put them out ? I hope somebody can help me BBah schrieb am Donnerstag, 19. Mai 2022 um 08:47:

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