[deal.II] How to use component_is_part_of_tensor

2020-11-24 Thread Audrey Collard-Daigneault
Hi everyone! I'm trying to export data as tensors to vtu format (for ParaView) with the add_data_vector() function with DataComponentInterpretation:: component_is_part_of_vector. Unfortunately, I don't know how to use it properly and I don't find any example or clue about it in the deal.II doc

[deal.II] Re: Post-processing velocities with Trilinos vectors during simulation

2020-11-10 Thread Audrey Collard-Daigneault
Hi Bruno, Thanks for your help! The problem was actually my "sum_reynolds_stress_dt += reynolds_stress_dt;" line. It needed to be outside of the for loop to be effective. I didn't know about LinearAlgebra, but it should propably be useful! Same for the way to work with BlockVectors! About the "

[deal.II] Post-processing velocities with Trilinos vectors during simulation

2020-11-04 Thread Audrey Collard-Daigneault
Hi everyone! I'm working on post-processing velocities with Trilinos solution vectors during the simulation on Lethe. Calculating average velocities and pressures (, , and ) works well using Trilinos vectors with no ghost cells and .equ(...) function. The calculation of average Reynolds stress