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

2016-12-20 Thread Denis Davydov
On Tuesday, December 20, 2016 at 5:22:50 PM UTC+1, Martin Kronbichler wrote: > > Hi Denis, > > VectorizedArray r0 = op.do_vectorized(index); > VectorizedArray r1 = op.do_vectorized(index+nvecs); > VectorizedArray r2 = op.do_vectorized(index+2*nvecs); > VectorizedAr

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 sh

[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 { VectorizedArray