Thanks you Peter, The error message is quite long and tricky (beyond my terminal can show) when I set VectorizedArray<double, width> with width=1. Knowing that deal.II do support this function is great. I‘ll try debug the code first and came back with further progress or questions :)
Best, Wayne 在2023年2月23日星期四 UTC+8 16:00:43<peterr...@gmail.com> 写道： > Hi Wayne, > > What is you error message? > > > I replace all FEEvaluation, MatrixFree and VectorizedArray construction > with a VectorizedArray<number, width>, but make returns some intricate > error. Is specify vectorization length possible? I think I've seen one in > dealii/test ... > > Yes, it is possible and the to specify the additional template argument is > the right approach. > > I have introduced this template parameter for the hyper.deal project: > https://github.com/hyperdeal/hyperdeal, where we create a tensor product > of MatrixFree instances (where one is vectorized and the other is not). > > Take a look at the test: > https://github.com/dealii/dealii/blob/3b785178375a8ea233fada96cf987aed6d53583e/tests/matrix_free/ecl.h#L52 > > The important aspect is that the vector length has to match the ones > supported by given ISA. See table 1 in the release paper of deal.II 9.2: > https://dealii.org/deal92-preprint.pdf > > Peter > > On Thursday, February 23, 2023 at 7:09:07 AM UTC+1 yy.wayne wrote: > >> Hi, >> is it possible to specify vectorization length(<deafult)? >> I'm running same code on local computer and server, whose cpu supports >> SSE2+4.10GHz >> and AVX512+2.60GHz respectively. The speed up for iteration time is only >> 2 on server, which should be 8*2.6 / 2*4.1 > 2.5 times. Therefore I'd like >> to run code with different >> vectorization length on same computer to test Vectorization's power. >> >> I replace all FEEvaluation, MatrixFree and VectorizedArray construction >> with a VectorizedArray<number, width>, but make returns some intricate >> error. Is specify vectorization length possible? I think I've seen one in >> dealii/test ... >> > -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/880f9fc6-763d-4d25-826f-2c23dc7dee39n%40googlegroups.com.
