New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
DataOutResample for multiple components #14258
DataOutResample for multiple components #14258
Conversation
8bf608e
to
29d31b1
Compare
for (unsigned int comp = 0; comp < dh.get_fe_collection().n_components(); | ||
++comp) | ||
{ | ||
const auto values = VectorTools::point_values<1>( | ||
rpe, dh, data_ptr->vector, VectorTools::EvaluationFlags::avg, comp); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This is the only relevant change. Loop over all components and call VectorTools::point_values
for each component.
/rebuild |
|
||
Assert(n_components == 1, | ||
Assert(n_components == 1 && first_selected_component == 0, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
That's the only place where you use first_selected_component
, right? Did I overlook something?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
It is used here:
dealii/include/deal.II/numerics/vector_tools_evaluate.h
Lines 446 to 453 in 29d31b1
std::make_unique<FEPointEvaluation<n_components, | |
dim, | |
spacedim, | |
typename VectorType::value_type>>( | |
cache.get_mapping(), | |
cell->get_fe(), | |
update_flags, | |
first_selected_component); |
Would you document the new parameter in the header file? I would also like to have a test if that's alright :) |
29d31b1
to
d9230d5
Compare
@marcfehling Done! |
FYI @mschreter @jh66637