[deal.II] Some step tutorials break in the current documentation

Hi deal.II developers, It seems that some tutorials, say step-1 and step-47 (may be incomplete), breaks in the current version of documentation. I'm not fimilar with the d

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

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

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

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

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

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

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