We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 6d05a96 commit e0fd7acCopy full SHA for e0fd7ac
1 file changed
src/libpsc/psc_push_fields/marder_impl.hxx
@@ -43,20 +43,19 @@ inline void correct(const Grid_t& grid, E1& efield, const Int3& efield_ib,
43
continue;
44
}
45
46
- Int3 l = grid.ibn;
47
- Int3 r = grid.ldims + grid.ibn;
48
auto e_comp = efield.view(_all, _all, _all, d, p);
49
50
- gt::gslice s1x = _s(l[0], r[0]);
51
- gt::gslice s1y = _s(l[1], r[1]);
52
- gt::gslice s1z = _s(l[2], r[2]);
53
-
+ Int3 l = grid.ibn;
+ Int3 r = grid.ldims + grid.ibn;
+ gt::gslice s1[3] = {_s(l[0], r[0]), _s(l[1], r[1]), _s(l[2], r[2])};
54
gt::gslice s2[3] = {_s(l[0], r[0]), _s(l[1], r[1]), _s(l[2], r[2])};
55
- s2[d] = _s(l[d] + 1, r[d] + 1);
+ s2[d].start += 1;
+ s2[d].stop += 1;
56
57
- e_comp.view(s1x, s1y, s1z) =
58
- e_comp.view(s1x, s1y, s1z) +
59
- (res.view(s2[0], s2[1], s2[2]) - res.view(s1x, s1y, s1z)) * fac[d];
+ e_comp.view(s1[0], s1[1], s1[2]) =
+ e_comp.view(s1[0], s1[1], s1[2]) +
+ (res.view(s2[0], s2[1], s2[2]) - res.view(s1[0], s1[1], s1[2])) *
+ fac[d];
60
61
62
0 commit comments