Open LLVM bug for

[gonzalobg]

constexpr T dist2(vec<T, N> a, vec<T, N> b) {
  T tmp = T(0.);
  for (dim_t i = 0; i < N; ++i) {
    T di = std::abs(a[i] - b[i]); // Improves perf to change it to:     T di = a[i] - b[i];
    tmp += di * di;
  }
  return tmp;