Field3DParallel + more FCI changes #3197

dschwoerer · 2025-11-06T08:06:22Z

Here a summary of the changes:

Field3DParallel

A Field3DParallel is essentially a Field3D, that preserves the parallel slices.
This allows to avoid unnecessary computation of of parallel fields, that are not needed.
It also ensures parallel slices are present (for FCI) when later parallel derivatives are computed.
So any function that will call DDY later on, can change the signature to Field3DParallel to declare that parallel slices need to be present for FCI.

New Y-Boundary Operators

Custom sheath boundary conditions can now be implemented using an abstraction. The documentation is on RTD
That allows having one code path for FCI and FA, for upper and lower sheath BC, all in one place.

Various small fixups

loadParallelMetrics

For FCI the parallel slices of the metric components are loaded, as they are sometimes needed.
This makes FCI essentially 3D only.
There is now also the option to say a field is not allow to compute the parallel slices, to avoid overwriting them: allowCalcParallelSlices.

Tracking on failure

Only if the simulation fails to evolve (Currently euler and pvode only):
The different components of timederivate are dumped to a BOUT.debug.*.nc which can be handy for figuring out which term is causing the instability. No more having to re-run the simulation with some terms disabled to debug such issues. This causes (an unkown) memory overhead - but ONLY if the solver has failed already.

More names of parallel boundary region + cleanup of code

Support things like bndry_par_xin or bndry_par_yup.

Add monotonichermitespline for all cases

Merged into hermitespline and works also in parallel.

Support several parallel slices for FCI

MYG=2 is now supported.

Previously only the first boundary point was set

For FCI we need to be able to access "random" data from the adjacent slices. If they are split in x-direction, this requires some tricky communication pattern. It can be used like this: ``` // Create object GlobalField3DAccess fci_comm(thismesh); // let it know what data points will be required: // where IndG3D is an index in the global field, which would be the // normal Ind3D if there would be only one proc. fci_comm.get(IndG3D(i, ny, nz)); // If all index have been added, the communication pattern will be // established. This has to be called by all processors in parallel fci_comm.setup() // Once the data for a given field is needed, it needs to be // communicated: GlobalField3DAccessInstance global_data = fci_comm.communicate(f3d); // and can be accessed like this BoutReal data = global_data[IndG3D(i, ny, nz)]; // ny and nz in the IndG3D are always optional. ```

If they are two instances of the same template, this allows to have an if in the inner loop that can be optimised out.

…igher-order

This reverts commit d7d7c91.

We have the info, so might as well store it

Taken from hermes-3, adopted for higher order

lower_bound takes into account the data is sorted

Ensures we all ways check for monotonicity

Tags were different for sender and receiver

Otherwise mpi might wait for the wrong request.

github-actions

clang-tidy made some suggestions

There were too many comments to post at once. Showing the first 25 out of 219. Check the log or trigger a new build to see more.

github-actions · 2025-11-10T09:14:03Z

externalpackages/PVODE/include/pvode/band.h

 ******************************************************************/

-#define PVODE_BAND_ELEM(A,i,j) ((A->data)[j][i-j+(A->smu)])
+#define PVODE_BAND_ELEM(A, i, j) (((A)->data)[j][(i) - (j) + ((A)->smu)])


warning: function-like macro 'PVODE_BAND_ELEM' used; consider a 'constexpr' template function [cppcoreguidelines-macro-usage]

#define PVODE_BAND_ELEM(A, i, j) (((A)->data)[j][(i) - (j) + ((A)->smu)]) ^

github-actions · 2025-11-10T09:14:03Z

externalpackages/PVODE/include/pvode/band.h

 ******************************************************************/

-#define PVODE_BAND_COL(A,j) (((A->data)[j])+(A->smu))
+#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu))


warning: function-like macro 'PVODE_BAND_COL' used; consider a 'constexpr' template function [cppcoreguidelines-macro-usage]

#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu)) ^

github-actions · 2025-11-10T09:14:03Z

externalpackages/PVODE/precon/band.h

 ******************************************************************/

-#define BAND_COL(A, j) (((A->data)[j]) + (A->smu))
+#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu))


warning: function-like macro 'PVODE_BAND_COL' used; consider a 'constexpr' template function [cppcoreguidelines-macro-usage]

#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu)) ^

github-actions · 2025-11-10T09:14:03Z

include/bout/boundary_iterator.hxx

+#include <algorithm>
+#include <functional>
+
+class BoundaryRegionIter {


warning: class 'BoundaryRegionIter' defines a default destructor but does not define a copy constructor, a copy assignment operator, a move constructor or a move assignment operator [cppcoreguidelines-special-member-functions]

class BoundaryRegionIter { ^

github-actions · 2025-11-10T09:14:03Z

include/bout/field.hxx

 inline T toFieldAligned(const T& f, const std::string& region = "RGN_ALL") {
  static_assert(bout::utils::is_Field_v<T>, "toFieldAligned only works on Fields");
+  ASSERT3(f.getCoordinates() != nullptr);
  return f.getCoordinates()->getParallelTransform().toFieldAligned(f, region);


warning: member access into incomplete type 'Coordinates' [clang-diagnostic-error]

return f.getCoordinates()->getParallelTransform().toFieldAligned(f, region); ^

Additional context

include/bout/index_derivs_interface.hxx:218: in instantiation of function template specialization 'toFieldAligned' requested here

const T f_aligned = is_unaligned ? toFieldAligned(f, "RGN_NOX") : f; ^

include/bout/index_derivs_interface.hxx:227: in instantiation of function template specialization 'bout::derivatives::index::DDY' requested here

return DDY(f.asField3D(), outloc, method, region); ^

include/bout/field_data.hxx:42: forward declaration of 'Coordinates'

class Coordinates; ^

github-actions · 2025-11-10T09:14:06Z

include/bout/field3d.hxx

+
+  int tracking_state{0};
+  Options* tracking{nullptr};
+  std::string selfname;


warning: member variable 'selfname' has protected visibility [cppcoreguidelines-non-private-member-variables-in-classes]

std::string selfname; ^

github-actions · 2025-11-10T09:14:06Z

include/bout/field3d.hxx

+  explicit Field3DParallel(Types... args) : Field3D(std::move(args)...) {
+    ensureFieldAligned();
+  }
+  Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }


warning: std::move of the const variable 'f' has no effect; remove std::move() or make the variable non-const [performance-move-const-arg]

Suggested change

Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }

Field3DParallel(const Field3D& f) : Field3D(f) { ensureFieldAligned(); }

github-actions · 2025-11-10T09:14:08Z

include/bout/field3d.hxx

+    ensureFieldAligned();
+  }
+  Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }
+  Field3DParallel(const Field2D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }


warning: std::move of the const variable 'f' has no effect; remove std::move() or make the variable non-const [performance-move-const-arg]

Suggested change

Field3DParallel(const Field2D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }

Field3DParallel(const Field2D& f) : Field3D(f) { ensureFieldAligned(); }

github-actions · 2025-11-10T09:14:08Z

include/bout/field3d.hxx

+};
+
+Field3DParallel Field3D::asField3DParallel() { return Field3DParallel(*this); }
+const Field3DParallel Field3D::asField3DParallel() const {


warning: return type 'const Field3DParallel' is 'const'-qualified at the top level, which may reduce code readability without improving const correctness [readability-const-return-type]

include/bout/field3d.hxx:533:

- inline const Field3DParallel asField3DParallel() const; + inline Field3DParallel asField3DParallel() const;

Suggested change

const Field3DParallel Field3D::asField3DParallel() const {

Field3DParallel Field3D::asField3DParallel() const {

github-actions · 2025-11-10T09:14:08Z

include/bout/field3d.hxx

+
+inline Field3DParallel
+filledFrom(const Field3DParallel& f,
+           std::function<BoutReal(int yoffset, Ind3D index)> func) {


warning: the parameter 'func' is copied for each invocation but only used as a const reference; consider making it a const reference [performance-unnecessary-value-param]

Suggested change

std::function<BoutReal(int yoffset, Ind3D index)> func) {

const std::function<BoutReal(int yoffset, Ind3D index)>& func) {

github-actions

clang-tidy made some suggestions

There were too many comments to post at once. Showing the first 25 out of 219. Check the log or trigger a new build to see more.

github-actions · 2025-11-10T09:14:13Z

externalpackages/PVODE/include/pvode/band.h

 ******************************************************************/

-#define PVODE_BAND_ELEM(A,i,j) ((A->data)[j][i-j+(A->smu)])
+#define PVODE_BAND_ELEM(A, i, j) (((A)->data)[j][(i) - (j) + ((A)->smu)])


warning: function-like macro 'PVODE_BAND_ELEM' used; consider a 'constexpr' template function [cppcoreguidelines-macro-usage]

#define PVODE_BAND_ELEM(A, i, j) (((A)->data)[j][(i) - (j) + ((A)->smu)]) ^

github-actions · 2025-11-10T09:14:13Z

externalpackages/PVODE/include/pvode/band.h

 ******************************************************************/

-#define PVODE_BAND_COL(A,j) (((A->data)[j])+(A->smu))
+#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu))


warning: function-like macro 'PVODE_BAND_COL' used; consider a 'constexpr' template function [cppcoreguidelines-macro-usage]

#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu)) ^

github-actions · 2025-11-10T09:14:13Z

externalpackages/PVODE/precon/band.h

 ******************************************************************/

-#define BAND_COL(A, j) (((A->data)[j]) + (A->smu))
+#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu))


warning: function-like macro 'PVODE_BAND_COL' used; consider a 'constexpr' template function [cppcoreguidelines-macro-usage]

#define PVODE_BAND_COL(A, j) ((((A)->data)[j]) + ((A)->smu)) ^

github-actions · 2025-11-10T09:14:13Z

include/bout/boundary_iterator.hxx

+#include <algorithm>
+#include <functional>
+
+class BoundaryRegionIter {


warning: class 'BoundaryRegionIter' defines a default destructor but does not define a copy constructor, a copy assignment operator, a move constructor or a move assignment operator [cppcoreguidelines-special-member-functions]

class BoundaryRegionIter { ^

github-actions · 2025-11-10T09:14:14Z

include/bout/field.hxx

 inline T toFieldAligned(const T& f, const std::string& region = "RGN_ALL") {
  static_assert(bout::utils::is_Field_v<T>, "toFieldAligned only works on Fields");
+  ASSERT3(f.getCoordinates() != nullptr);
  return f.getCoordinates()->getParallelTransform().toFieldAligned(f, region);


warning: member access into incomplete type 'Coordinates' [clang-diagnostic-error]

return f.getCoordinates()->getParallelTransform().toFieldAligned(f, region); ^

Additional context

include/bout/index_derivs_interface.hxx:218: in instantiation of function template specialization 'toFieldAligned' requested here

const T f_aligned = is_unaligned ? toFieldAligned(f, "RGN_NOX") : f; ^

include/bout/index_derivs_interface.hxx:227: in instantiation of function template specialization 'bout::derivatives::index::DDY' requested here

return DDY(f.asField3D(), outloc, method, region); ^

include/bout/field_data.hxx:42: forward declaration of 'Coordinates'

class Coordinates; ^

github-actions · 2025-11-10T09:14:16Z

include/bout/field3d.hxx

+
+  int tracking_state{0};
+  Options* tracking{nullptr};
+  std::string selfname;


warning: member variable 'selfname' has protected visibility [cppcoreguidelines-non-private-member-variables-in-classes]

std::string selfname; ^

github-actions · 2025-11-10T09:14:16Z

include/bout/field3d.hxx

+  explicit Field3DParallel(Types... args) : Field3D(std::move(args)...) {
+    ensureFieldAligned();
+  }
+  Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }


warning: std::move of the const variable 'f' has no effect; remove std::move() or make the variable non-const [performance-move-const-arg]

Suggested change

Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }

Field3DParallel(const Field3D& f) : Field3D(f) { ensureFieldAligned(); }

github-actions · 2025-11-10T09:14:16Z

include/bout/field3d.hxx

+    ensureFieldAligned();
+  }
+  Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }
+  Field3DParallel(const Field2D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }


warning: std::move of the const variable 'f' has no effect; remove std::move() or make the variable non-const [performance-move-const-arg]

Suggested change

Field3DParallel(const Field2D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }

Field3DParallel(const Field2D& f) : Field3D(f) { ensureFieldAligned(); }

github-actions · 2025-11-10T09:14:16Z

include/bout/field3d.hxx

+};
+
+Field3DParallel Field3D::asField3DParallel() { return Field3DParallel(*this); }
+const Field3DParallel Field3D::asField3DParallel() const {


warning: return type 'const Field3DParallel' is 'const'-qualified at the top level, which may reduce code readability without improving const correctness [readability-const-return-type]

include/bout/field3d.hxx:533:

- inline const Field3DParallel asField3DParallel() const; + inline Field3DParallel asField3DParallel() const;

Suggested change

const Field3DParallel Field3D::asField3DParallel() const {

Field3DParallel Field3D::asField3DParallel() const {

github-actions · 2025-11-10T09:14:16Z

include/bout/field3d.hxx

+
+inline Field3DParallel
+filledFrom(const Field3DParallel& f,
+           std::function<BoutReal(int yoffset, Ind3D index)> func) {


warning: the parameter 'func' is copied for each invocation but only used as a const reference; consider making it a const reference [performance-unnecessary-value-param]

Suggested change

std::function<BoutReal(int yoffset, Ind3D index)> func) {

const std::function<BoutReal(int yoffset, Ind3D index)>& func) {

ZedThree · 2025-11-17T16:14:51Z

Please can you split this into one PR per change? It's 100 files and 5k lines of changes, it's going to be very difficult to review. GitHub will only show a single file at once.

dschwoerer · 2025-11-18T12:42:40Z

Sure!
Is ready for re-review / merge: #2889
I will keep working on the other bits ...

Information like perpendicular J, Bxyz and g_22 are only known to the grid generator, but are needed for more correct parallel derivatives

J on the parallel slices is not accurate.

This avoids mostly the need for annotations using BOUT_USE_TRACK. setName could thus be dropped.

Add unit test as well as fix

github-actions

clang-tidy made some suggestions

There were too many comments to post at once. Showing the first 25 out of 288. Check the log or trigger a new build to see more.

github-actions · 2026-01-29T08:59:22Z

include/bout/coordinates.hxx

+  const FieldMetric& Jxz() const;
+
+private:
+  mutable std::optional<FieldMetric> _g_22_ylow, _g_22_yhigh;


warning: no header providing "std::optional" is directly included [misc-include-cleaner]

include/bout/coordinates.hxx:40:

+ #include <optional>

github-actions · 2026-01-29T08:59:22Z

include/bout/coordinates.hxx

+
+private:
+  mutable std::optional<FieldMetric> _g_22_ylow, _g_22_yhigh;
+  mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;


warning: declaration uses identifier '_Jxz_centre', which is a reserved identifier [bugprone-reserved-identifier]

Suggested change

mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;

mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, Jxz_centre;

github-actions · 2026-01-29T08:59:22Z

include/bout/coordinates.hxx

+
+private:
+  mutable std::optional<FieldMetric> _g_22_ylow, _g_22_yhigh;
+  mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;


warning: declaration uses identifier '_Jxz_yhigh', which is a reserved identifier [bugprone-reserved-identifier]

Suggested change

mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;

mutable std::optional<FieldMetric> _Jxz_ylow, Jxz_yhigh, _Jxz_centre;

github-actions · 2026-01-29T08:59:22Z

include/bout/coordinates.hxx

+
+private:
+  mutable std::optional<FieldMetric> _g_22_ylow, _g_22_yhigh;
+  mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;


warning: declaration uses identifier '_Jxz_ylow', which is a reserved identifier [bugprone-reserved-identifier]

Suggested change

mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;

mutable std::optional<FieldMetric> Jxz_ylow, _Jxz_yhigh, _Jxz_centre;

github-actions · 2026-01-29T08:59:22Z

include/bout/deriv_store.hxx


+template <typename T>
+auto& getStore() {
+  if constexpr (std::is_same<T, Field3DParallel>::value) {


warning: no header providing "Field3DParallel" is directly included [misc-include-cleaner]

include/bout/deriv_store.hxx:37:

- #include <bout/scorepwrapper.hxx> + #include "bout/field3d.hxx" + #include <bout/scorepwrapper.hxx>

github-actions · 2026-01-29T08:59:23Z

include/bout/fv_ops.hxx

 #include <bout/mesh.hxx>
+#if CHECK > 0
+#include <bout/output_bout_types.hxx>
+#endif


warning: included header output_bout_types.hxx is not used directly [misc-include-cleaner]

Suggested change

#endif

#endif

github-actions · 2026-01-29T08:59:23Z

include/bout/fv_ops.hxx

-  // Left and right cell face values
-  BoutReal L, R;
+  /// Cell centre values
+  BoutReal c;


warning: no header providing "BoutReal" is directly included [misc-include-cleaner]

include/bout/fv_ops.hxx:7:

- #include "bout/build_defines.hxx" + #include "bout/bout_types.hxx" + #include "bout/build_defines.hxx"

github-actions · 2026-01-29T08:59:23Z

include/bout/fv_ops.hxx

+  BoutReal m;
+  BoutReal p;
+  BoutReal mm = BoutNaN;
+  BoutReal pp = BoutNaN;


warning: no header providing "BoutNaN" is directly included [misc-include-cleaner]

BoutNaN; ^

github-actions · 2026-01-29T08:59:23Z

include/bout/fv_ops.hxx

+        BoutReal flux = NAN;

-        if (mesh->lastY(i) && (j == mesh->yend) && !mesh->periodicY(i)) {
+        if (is_last_y && (j == mesh->yend) && !is_periodic_y) {


warning: no header providing "NAN" is directly included [misc-include-cleaner]

Real flux = NAN; ^

github-actions · 2026-01-29T08:59:24Z

include/bout/fv_ops.hxx

+    const auto& v_down = v_in.ydown();
+
+    Field3D result{emptyFrom(f_in)};
+    BOUT_FOR(i, f_in.getRegion("RGN_NOBNDRY")) {


warning: no header providing "emptyFrom" is directly included [misc-include-cleaner]

emptyFrom(f_in)}; ^

dschwoerer added 30 commits January 13, 2025 11:15

Fix extrapolaton / interpolation

7e1067a

Fix parallel boundary to interpolate into the boundary

d416b5d

Previously only the first boundary point was set

Ensure data is sorted

6cca6a4

fix typo

62b62bf

Calculate valid for the general case

628a6ce

Unify XZMonotonicHermiteSpline and XZMonotonicHermiteSpline

40dac59

If they are two instances of the same template, this allows to have an if in the inner loop that can be optimised out.

Use x-splitting for monotonichermitespline test

9a98a8b

Fix position of maybe_unused for old gcc

cb068ea

Ensure PETSC_SUCCESS is defined

d7d7c91

Merge remote-tracking branch 'origin/next' into fci-auto-with-debug-h…

3cc13b9

…igher-order

Revert "Ensure PETSC_SUCCESS is defined"

a05201b

This reverts commit d7d7c91.

Define PETSC_SUCCESS for old petsc versions

7d73c5c

Set region for lagrange4pt

ebf2d07

Fix split parallel slices

ae1fcfe

Add option to split + allocate

8b7d97d

fix parallel neumann BC

ad56d8e

set name for parallel component

2a48341

We have the info, so might as well store it

Add limitFree

f73813b

Taken from hermes-3, adopted for higher order

add setAll to parallel BC

ff7525d

Add monotonic check also to other code branches

1c44f98

use lower_bound instead of find

42c0958

lower_bound takes into account the data is sorted

Do not shadow mapping

21d1985

Ensure setup has been called

4f39c1d

Use pointer to data

cabdc4c

Call setup before communicator is used

aa7938d

Deduplicate code

d1dee59

Ensures we all ways check for monotonicity

Fix tags for comm

5818327

Tags were different for sender and receiver

Use pointer instead of std::vector

aa2d8b6

Do not reuse requests if the array is still in use

007fed0

Otherwise mpi might wait for the wrong request.

github-actions bot reviewed Nov 10, 2025

View reviewed changes

Fix merge

f46e2fd

dschwoerer mentioned this pull request Nov 10, 2025

Port some Hermes-3 finite volume operators #3200

Open

dschwoerer added 3 commits November 10, 2025 13:13

Move parallel region creation to defaultRegions

35e86c5

Add Bxy to reference grid

a4473ee

Move check up

5a04011

dschwoerer added 17 commits January 20, 2026 09:25

Add metric info on cell faces for FCI

b524bba

Information like perpendicular J, Bxyz and g_22 are only known to the grid generator, but are needed for more correct parallel derivatives

ensure parallel fields are allocated

96ad527

Ensure FCI path is always used for FCI

c132458

FCI: Use cell face values for derivatives

a3206db

Do not load parallel J

38233c3

J on the parallel slices is not accurate.

Avoid warning for unused var

9e5f448

Merge remote-tracking branch 'origin/next' into f3dwy-port-fv-ops

c5e3669

Add backtrace to tracked operations

ef3683e

This avoids mostly the need for annotations using BOUT_USE_TRACK. setName could thus be dropped.

Only call iterateBndryLowerY() for non-FCI

b3cc11c

Ensure f3d.asField3DParallel() *= x updates f3d

a799155

Add unit test as well as fix

ensure also parallel fields are updated in-place

b959cc7

Clarify exception message

234940c

Only compute parallel fields for 3D Metrics

f4f8362

Always check whether main component is constant or doZero is set

1031ccf

Remove stray ,

78296ce

Add missing header

32af728

Load J and set J to sane value

f5d72eb

dschwoerer force-pushed the f3dwy branch from 4ee2443 to f5d72eb Compare January 29, 2026 08:45

dschwoerer added 2 commits January 29, 2026 08:45

Apply black changes

655d9a9

Apply clang-format changes

c0ac486

github-actions bot reviewed Jan 29, 2026

View reviewed changes

	Field3DParallel(const Field3D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }
	Field3DParallel(const Field3D& f) : Field3D(f) { ensureFieldAligned(); }

	Field3DParallel(const Field2D& f) : Field3D(std::move(f)) { ensureFieldAligned(); }
	Field3DParallel(const Field2D& f) : Field3D(f) { ensureFieldAligned(); }

	const Field3DParallel Field3D::asField3DParallel() const {
	Field3DParallel Field3D::asField3DParallel() const {

	std::function<BoutReal(int yoffset, Ind3D index)> func) {
	const std::function<BoutReal(int yoffset, Ind3D index)>& func) {

	mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;
	mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, Jxz_centre;

	mutable std::optional<FieldMetric> _Jxz_ylow, _Jxz_yhigh, _Jxz_centre;
	mutable std::optional<FieldMetric> Jxz_ylow, _Jxz_yhigh, _Jxz_centre;

Field3DParallel + more FCI changes #3197

Are you sure you want to change the base?

Field3DParallel + more FCI changes #3197

Uh oh!

Conversation

dschwoerer commented Nov 6, 2025 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!