Dear Michal,

Could you please be a bit more specific? What was the error that you obtained? What vector did you pass into FEEvaluation::read_dof_values/distribute_local_to_global? I guess that you handed in a simple vector, e.g. LinearAlgebra::distributed::Vector and got a segmentation fault. But to add the assertion I need more details. The issues that in principle, we do support using multiple components out of a scalar FE - but in that case the user must hand in a BlockVector or std::vector<VectorType>. We could catch that case in the machinery, though.

Thanks!

Best,
Martin


On 07.10.2017 23:00, Michał Wichrowski wrote:
Dear deal.II developers,
I've found out that in Matrix-free framework there is no assertion that checks if number of components of finite element matches number components used in cell operations. By mistake I used scalar FE instead vector and it took me quite a long time to figure it out. The program worked and even covered in reasonable number of iterations. I think FEEvaluation is right place to put this assertion.

Michał Wichrowski
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com <mailto:dealii+unsubscr...@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to