Re: [dealii-developers] do_vectorized() in struct Dot and alike for shared memory parallel reduction

2016-12-20 Thread Martin Kronbichler
Hi Denis, | VectorizedArrayr0 =op.do_vectorized(index); VectorizedArrayr1 =op.do_vectorized(index+nvecs); VectorizedArrayr2 =op.do_vectorized(index+2*nvecs); VectorizedArrayr3 =op.do_vectorized(index+3*nvecs); | where | constunsignedintnvecs =VectorizedArray::n_array_elements; | I wonder

[dealii-developers] do_vectorized() in struct Dot and alike for shared memory parallel reduction

2016-12-20 Thread Denis Davydov
Hi All, I was digging around lac/vectorized_operations_internal.h and noticed that objects that are used in SIMD to do reduction are of the form: template struct Dot {