Closes: https://github.com/dealii/dealii/issues/10146
// evaluate the geometry in double precision; thus, for floats we need
// to do things in two sweeps and convert the final result.
constexpr unsigned int n_lanes = VectorizedArrayType::size();
+
using VectorizedDouble =
VectorizedArray<double,
- ((std::is_same<Number, float>::value && n_lanes > 1) ?
- n_lanes / 2 :
- n_lanes)>;
+ ((std::is_same<Number, float>::value &&
+ VectorizedArrayType::size() > 1) ?
+ VectorizedArrayType::size() / 2 :
+ VectorizedArrayType::size())>;
// Create a ShapeInfo object to provide the necessary interpolators to
// the various quadrature points. Note that it is initialized with the
In the beginning the Universe was created. This has made a lot of
people very angry and has been widely regarded as a bad move.
Douglas Adams