Skip to content

Commit

Permalink
AsterX: don't sync Bs
Browse files Browse the repository at this point in the history
  • Loading branch information
lwJi committed Oct 3, 2024
1 parent 985e0a3 commit f814710
Show file tree
Hide file tree
Showing 2 changed files with 9 additions and 14 deletions.
17 changes: 6 additions & 11 deletions AsterX/schedule.ccl
Original file line number Diff line number Diff line change
Expand Up @@ -44,24 +44,21 @@ SCHEDULE AsterX_ComputedBstagFromA IN AsterX_InitialGroup
{
LANG: C
READS: Avec_x(everywhere) Avec_y(everywhere) Avec_z(everywhere)
WRITES: dBx_stag(interior) dBy_stag(interior) dBz_stag(interior)
SYNC: dBx_stag dBy_stag dBz_stag
WRITES: dBx_stag(everywhere) dBy_stag(everywhere) dBz_stag(everywhere)
} "Calculate dBstag from curl of A"

SCHEDULE AsterX_ComputedBFromdBstag IN AsterX_InitialGroup AFTER AsterX_ComputedBstagFromA
{
LANG: C
READS: dBx_stag(everywhere) dBy_stag(everywhere) dBz_stag(everywhere)
WRITES: dB(interior)
SYNC: dB
WRITES: dB(everywhere)
} "Calculate centered dB from dBstag"

SCHEDULE AsterX_ComputeBFromdB IN AsterX_InitialGroup AFTER AsterX_ComputedBFromdBstag
{
LANG: C
READS: dB(interior) ADMBaseX::metric(interior)
WRITES: HydroBaseX::Bvec(interior)
SYNC: HydroBaseX::Bvec
READS: dB(everywhere) ADMBaseX::metric(everywhere)
WRITES: HydroBaseX::Bvec(everywhere)
} "Calculate centered B from densitized B"

SCHEDULE AsterX_Prim2Con_Initial IN AsterX_InitialGroup AFTER AsterX_ComputeBFromdB
Expand Down Expand Up @@ -107,16 +104,14 @@ SCHEDULE AsterX_ComputedBstagFromA IN AsterX_Con2PrimGroup
{
LANG: C
READS: Avec_x(everywhere) Avec_y(everywhere) Avec_z(everywhere)
WRITES: dBx_stag(interior) dBy_stag(interior) dBz_stag(interior)
SYNC: dBx_stag dBy_stag dBz_stag
WRITES: dBx_stag(everywhere) dBy_stag(everywhere) dBz_stag(everywhere)
} "Calculate dBstag from curl of A"

SCHEDULE AsterX_ComputedBFromdBstag IN AsterX_Con2PrimGroup AFTER AsterX_ComputedBstagFromA
{
LANG: C
READS: dBx_stag(everywhere) dBy_stag(everywhere) dBz_stag(everywhere)
WRITES: dB(interior)
SYNC: dB
WRITES: dB(everywhere)
} "Calculate centered dB from dBstag"

SCHEDULE AsterX_Con2Prim IN AsterX_Con2PrimGroup AFTER AsterX_ComputedBFromdBstag
Expand Down
6 changes: 3 additions & 3 deletions AsterX/src/computeBfromA.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@ template <int dir> void ComputeStaggeredB(CCTK_ARGUMENTS) {

constexpr array<int, dim> face_centred = {!(dir == 0), !(dir == 1),
!(dir == 2)};
grid.loop_int_device<face_centred[0], face_centred[1], face_centred[2]>(
grid.loop_all_device<face_centred[0], face_centred[1], face_centred[2]>(
grid.nghostzones,
[=] CCTK_DEVICE(const PointDesc &p) CCTK_ATTRIBUTE_ALWAYS_INLINE {
// Neighbouring "plus" and "minus" cell indices
Expand Down Expand Up @@ -68,7 +68,7 @@ extern "C" void AsterX_ComputedBFromdBstag(CCTK_ARGUMENTS) {
DECLARE_CCTK_ARGUMENTSX_AsterX_ComputedBFromdBstag;
DECLARE_CCTK_PARAMETERS;

grid.loop_int_device<1, 1, 1>(
grid.loop_all_device<1, 1, 1>(
grid.nghostzones,
[=] CCTK_DEVICE(const PointDesc &p) CCTK_ATTRIBUTE_ALWAYS_INLINE {
// Neighbouring "plus" and "minus" cell indices
Expand All @@ -88,7 +88,7 @@ extern "C" void AsterX_ComputeBFromdB(CCTK_ARGUMENTS) {
DECLARE_CCTK_ARGUMENTSX_AsterX_ComputeBFromdB;
DECLARE_CCTK_PARAMETERS;

grid.loop_int_device<1, 1, 1>(
grid.loop_all_device<1, 1, 1>(
grid.nghostzones,
[=] CCTK_DEVICE(const PointDesc &p) CCTK_ATTRIBUTE_ALWAYS_INLINE {
/* Interpolate metric terms from vertices to center */
Expand Down

0 comments on commit f814710

Please sign in to comment.