Skip to content

Commit 342924b

Browse files
JamesMcClungJamesMcClung
committed
current_deposition: format to be obviously correct
1 parent a8717d1 commit 342924b

1 file changed

Lines changed: 39 additions & 53 deletions

File tree

‎src/include/psc/current_deposition.hxx‎

Lines changed: 39 additions & 53 deletions
Original file line numberDiff line numberDiff line change
@@ -21,32 +21,22 @@ public:
2121
real_t h = (1.f / real_t(12.f)) * dx.prod();
2222
Real3 fnq = qni_wni * fnqs_;
2323

24-
curr.add(0, i[0], i[1], i[2],
25-
fnq[0] * (dx[0] * (1.f - xa[1]) * (1.f - xa[2]) + h));
26-
curr.add(0, i[0], i[1] + 1, i[2],
27-
fnq[0] * (dx[0] * (xa[1]) * (1.f - xa[2]) - h));
28-
curr.add(0, i[0], i[1], i[2] + 1,
29-
fnq[0] * (dx[0] * (1.f - xa[1]) * (xa[2]) - h));
30-
curr.add(0, i[0], i[1] + 1, i[2] + 1,
31-
fnq[0] * (dx[0] * (xa[1]) * (xa[2]) + h));
32-
33-
curr.add(1, i[0], i[1], i[2],
34-
fnq[1] * (dx[1] * (1.f - xa[0]) * (1.f - xa[2]) + h));
35-
curr.add(1, i[0] + 1, i[1], i[2],
36-
fnq[1] * (dx[1] * (xa[0]) * (1.f - xa[2]) - h));
37-
curr.add(1, i[0], i[1], i[2] + 1,
38-
fnq[1] * (dx[1] * (1.f - xa[0]) * (xa[2]) - h));
39-
curr.add(1, i[0] + 1, i[1], i[2] + 1,
40-
fnq[1] * (dx[1] * (xa[0]) * (xa[2]) + h));
41-
42-
curr.add(2, i[0], i[1], i[2],
43-
fnq[2] * (dx[2] * (1.f - xa[0]) * (1.f - xa[1]) + h));
44-
curr.add(2, i[0] + 1, i[1], i[2],
45-
fnq[2] * (dx[2] * (xa[0]) * (1.f - xa[1]) - h));
46-
curr.add(2, i[0], i[1] + 1, i[2],
47-
fnq[2] * (dx[2] * (1.f - xa[0]) * (xa[1]) - h));
48-
curr.add(2, i[0] + 1, i[1] + 1, i[2],
49-
fnq[2] * (dx[2] * (xa[0]) * (xa[1]) + h));
24+
// clang-format off
25+
curr.add(0, i[0], i[1], i[2], fnq[0] * (dx[0] * (1.f - xa[1]) * (1.f - xa[2]) + h));
26+
curr.add(0, i[0], i[1] + 1, i[2], fnq[0] * (dx[0] * ( xa[1]) * (1.f - xa[2]) - h));
27+
curr.add(0, i[0], i[1], i[2] + 1, fnq[0] * (dx[0] * (1.f - xa[1]) * ( xa[2]) - h));
28+
curr.add(0, i[0], i[1] + 1, i[2] + 1, fnq[0] * (dx[0] * ( xa[1]) * ( xa[2]) + h));
29+
30+
curr.add(1, i[0], i[1], i[2], fnq[1] * (dx[1] * (1.f - xa[2]) * (1.f - xa[0]) + h));
31+
curr.add(1, i[0], i[1], i[2] + 1, fnq[1] * (dx[1] * ( xa[2]) * (1.f - xa[0]) - h));
32+
curr.add(1, i[0] + 1, i[1], i[2], fnq[1] * (dx[1] * (1.f - xa[2]) * ( xa[0]) - h));
33+
curr.add(1, i[0] + 1, i[1], i[2] + 1, fnq[1] * (dx[1] * ( xa[2]) * ( xa[0]) + h));
34+
35+
curr.add(2, i[0], i[1], i[2], fnq[2] * (dx[2] * (1.f - xa[0]) * (1.f - xa[1]) + h));
36+
curr.add(2, i[0] + 1, i[1], i[2], fnq[2] * (dx[2] * ( xa[0]) * (1.f - xa[1]) - h));
37+
curr.add(2, i[0], i[1] + 1, i[2], fnq[2] * (dx[2] * (1.f - xa[0]) * ( xa[1]) - h));
38+
curr.add(2, i[0] + 1, i[1] + 1, i[2], fnq[2] * (dx[2] * ( xa[0]) * ( xa[1]) + h));
39+
// clang-format on
5040
};
5141

5242
GT_INLINE void operator()(Curr& curr, const Int3& i, real_t qni_wni,
@@ -56,43 +46,39 @@ public:
5646
real_t h = (1.f / real_t(12.f)) * dx.prod();
5747
Real3 fnq = qni_wni * fnqs_;
5848

59-
curr.add(0, i[0], i[1], i[2], fnq[0] * (dx[0] * (1.f - xa[2])));
60-
curr.add(0, i[0], i[1], i[2] + 1, fnq[0] * (dx[0] * (xa[2])));
49+
// clang-format off
50+
curr.add(0, i[0], i[1], i[2], fnq[0] * (dx[0] * (1.f - xa[2]) ));
51+
curr.add(0, i[0], i[1], i[2] + 1, fnq[0] * (dx[0] * ( xa[2]) ));
6152

62-
curr.add(1, i[0], i[1], i[2],
63-
fnq[1] * (dx[1] * (1.f - xa[0]) * (1.f - xa[2]) + h));
64-
curr.add(1, i[0] + 1, i[1], i[2],
65-
fnq[1] * (dx[1] * (xa[0]) * (1.f - xa[2]) - h));
66-
curr.add(1, i[0], i[1], i[2] + 1,
67-
fnq[1] * (dx[1] * (1.f - xa[0]) * (xa[2]) - h));
68-
curr.add(1, i[0] + 1, i[1], i[2] + 1,
69-
fnq[1] * (dx[1] * (xa[0]) * (xa[2]) + h));
53+
curr.add(1, i[0], i[1], i[2], fnq[1] * (dx[1] * (1.f - xa[2]) * (1.f - xa[0]) + h));
54+
curr.add(1, i[0], i[1], i[2] + 1, fnq[1] * (dx[1] * ( xa[2]) * (1.f - xa[0]) - h));
55+
curr.add(1, i[0] + 1, i[1], i[2], fnq[1] * (dx[1] * (1.f - xa[2]) * ( xa[0]) - h));
56+
curr.add(1, i[0] + 1, i[1], i[2] + 1, fnq[1] * (dx[1] * ( xa[2]) * ( xa[0]) + h));
7057

71-
curr.add(2, i[0], i[1], i[2], fnq[2] * (dx[2] * (1.f - xa[0])));
72-
curr.add(2, i[0] + 1, i[1], i[2], fnq[2] * (dx[2] * (xa[0])));
58+
curr.add(2, i[0], i[1], i[2], fnq[2] * (dx[2] * (1.f - xa[0]) ));
59+
curr.add(2, i[0] + 1, i[1], i[2], fnq[2] * (dx[2] * ( xa[0]) ));
60+
// clang-format on
7361
}
7462

7563
GT_INLINE void operator()(Curr& curr, const Int3& i, real_t qni_wni,
7664
const Real3& dx, const Real3& xa,
7765
dim_yz tag_dim) const
7866
{
7967
real_t h = (1.f / real_t(12.f)) * dx.prod();
80-
8168
Real3 fnq = qni_wni * fnqs_;
82-
curr.add(0, i[0], i[1], i[2],
83-
fnq[0] * (dx[0] * (1.f - xa[1]) * (1.f - xa[2]) + h));
84-
curr.add(0, i[0], i[1] + 1, i[2],
85-
fnq[0] * (dx[0] * (xa[1]) * (1.f - xa[2]) - h));
86-
curr.add(0, i[0], i[1], i[2] + 1,
87-
fnq[0] * (dx[0] * (1.f - xa[1]) * (xa[2]) - h));
88-
curr.add(0, i[0], i[1] + 1, i[2] + 1,
89-
fnq[0] * (dx[0] * (xa[1]) * (xa[2]) + h));
90-
91-
curr.add(1, i[0], i[1], i[2], fnq[1] * (dx[1] * (1.f - xa[2])));
92-
curr.add(1, i[0], i[1], i[2] + 1, fnq[1] * (dx[1] * (xa[2])));
93-
94-
curr.add(2, i[0], i[1], i[2], fnq[2] * (dx[2] * (1.f - xa[1])));
95-
curr.add(2, i[0], i[1] + 1, i[2], fnq[2] * (dx[2] * (xa[1])));
69+
70+
// clang-format off
71+
curr.add(0, i[0], i[1], i[2], fnq[0] * (dx[0] * (1.f - xa[1]) * (1.f - xa[2]) + h));
72+
curr.add(0, i[0], i[1] + 1, i[2], fnq[0] * (dx[0] * ( xa[1]) * (1.f - xa[2]) - h));
73+
curr.add(0, i[0], i[1], i[2] + 1, fnq[0] * (dx[0] * (1.f - xa[1]) * ( xa[2]) - h));
74+
curr.add(0, i[0], i[1] + 1, i[2] + 1, fnq[0] * (dx[0] * ( xa[1]) * ( xa[2]) + h));
75+
76+
curr.add(1, i[0], i[1], i[2], fnq[1] * (dx[1] * (1.f - xa[2]) ));
77+
curr.add(1, i[0], i[1], i[2] + 1, fnq[1] * (dx[1] * ( xa[2]) ));
78+
79+
curr.add(2, i[0], i[1], i[2], fnq[2] * (dx[2] * (1.f - xa[1]) ));
80+
curr.add(2, i[0], i[1] + 1, i[2], fnq[2] * (dx[2] * ( xa[1]) ));
81+
// clang-format on
9682
}
9783

9884
private:

0 commit comments

Comments
 (0)