static void updateLineX(HOOD_NEW& hoodNew, int indexEnd, HOOD_OLD& hoodOld, unsigned /* nanoStep */)
{
for (int i = hoodOld.index(); i < indexEnd; ++i, ++hoodOld) {
ShortVec tmp;
tmp.load_aligned(&hoodNew->sum() + i * C);
for (const auto& j: hoodOld.weights(0)) {
ShortVec weights;
ShortVec values;
weights.load_aligned(j.second());
values.gather(&hoodOld->value(), j.first());
tmp += values * weights;
}
tmp.store_aligned(&hoodNew->sum() + i * C);
}
}
static void updateLineX(HOOD_NEW& hoodNew, int indexEnd, HOOD_OLD& hoodOld, unsigned /* nanoStep */)
{
for (int i = hoodOld.index(); i < indexEnd; ++i, ++hoodOld) {
ShortVec tmp;
tmp.load_aligned(&hoodNew->sum() + i * C);
for (const auto& j: hoodOld.weights(0)) {
ShortVec weights;
ShortVec values;
weights.load_aligned(j.second());
// fixme: is this gahter actually correct? shouldn't we use offset 0 for the gather? see also hpxperformancetests/main.cpp
values.gather(&hoodOld->value(), j.first());
tmp += values * weights;
}
tmp.store_aligned(&hoodNew->sum() + i * C);
}
}
inline static void updateLineX(HOOD_NEW& hoodNew, int indexEnd, HOOD_OLD& hoodOld, unsigned /* nanoStep */)
{
typedef LibFlatArray::short_vec<double, C> ShortVec;
for (; hoodNew.index() < indexEnd; hoodNew += C, ++hoodOld) {
ShortVec x = &hoodOld->x();
ShortVec y = &hoodOld->y();
ShortVec cReal = &hoodOld->cReal();
ShortVec cImag = &hoodOld->cImag();
for (int i = 0; i < ITERATIONS; ++i) {
ShortVec cRealOld = cReal;
cReal = cReal * cReal - cImag * cImag;
cImag = ShortVec(2.0) * cImag * cRealOld;
}
for (const auto& j: hoodOld.weights(0)) {
ShortVec weights;
ShortVec otherX;
ShortVec otherY;
weights.load_aligned(j.second());
otherX.gather(&hoodOld->x(), j.first());
otherY.gather(&hoodOld->y(), j.first());
cReal += otherX * weights;
cImag += otherY * weights;
}
&hoodNew->x() << x;
&hoodNew->y() << y;
&hoodNew->cReal() << cReal;
&hoodNew->cImag() << cImag;
}
}
inline static void updateLineX(HOOD_NEW& hoodNew, int indexEnd, HOOD_OLD& hoodOld, unsigned /* nanoStep */)
{
for (; hoodOld.index() < indexEnd; ++hoodNew, ++hoodOld) {
UnstructuredBusyworkCellWithUpdateLineX& self = *hoodNew;
self = hoodOld[hoodOld.index()];
for (int i = 0; i < ITERATIONS; ++i) {
double cRealOld = self.cReal;
self.cReal = self.cReal * self.cReal - self.cImag * self.cImag;
self.cImag = 2 * self.cImag * cRealOld;
}
for (auto i = hoodOld.begin(); i != hoodOld.end(); ++i) {
self.cReal += hoodOld[i.first()].x * i.second();
self.cImag += hoodOld[i.first()].y * i.second();
}
}
}
