@@ -166,10 +166,10 @@ TaskStatus CoolingFunctionCalculateFourForce(MeshBlockData<Real> *rc, const doub
166
166
167
167
// Light Bulb with Liebendorfer model
168
168
const bool is_monopole_cart =
169
- (typeid (PHOEBUS_GEOMETRY) == typeid (Geometry::MonopoleCart));
169
+ (typeid (PHOEBUS_GEOMETRY) == typeid (Geometry::MonopoleCart));
170
170
const bool is_monopole_sph =
171
- (typeid (PHOEBUS_GEOMETRY) == typeid (Geometry::MonopoleSph));
172
-
171
+ (typeid (PHOEBUS_GEOMETRY) == typeid (Geometry::MonopoleSph));
172
+
173
173
auto &coords = pmb->coords ;
174
174
using Transformation_t = Geometry::SphericalToCartesian;
175
175
auto const &geom_pkg = pmb->packages .Get (" geometry" );
@@ -189,14 +189,14 @@ TaskStatus CoolingFunctionCalculateFourForce(MeshBlockData<Real> *rc, const doub
189
189
DEFAULT_LOOP_PATTERN, " CoolingFunctionCalculateFourForce" , DevExecSpace (), kb.s ,
190
190
kb.e , jb.s , jb.e , ib.s , ib.e ,
191
191
KOKKOS_LAMBDA (const int k, const int j, const int i) {
192
- const Real x1 = coords.Xc <1 >(k, j, i);
192
+ const Real x1 = coords.Xc <1 >(k, j, i);
193
193
const Real x2 = coords.Xc <2 >(k, j, i);
194
194
const Real x3 = coords.Xc <3 >(k, j, i);
195
195
Real Cart[3 ];
196
196
Real s2c[3 ][3 ], c2s[3 ][3 ];
197
197
Real r;
198
198
199
- if (is_monopole_sph){
199
+ if (is_monopole_sph) {
200
200
r = std::abs (x1);
201
201
} else {
202
202
r = std::sqrt (x1 * x1 + x2 * x2 + x3 * x3);
0 commit comments