MPI

subdomain with offset

// Subdomain starting at global index (10, 20, 0)
pfc::DiscreteField<double, 3> subdomain(
    {32, 32, 64},           // local size
    {10, 20, 0},            // global offset
    {0.0, 0.0, 0.0},        // global origin
    {0.5, 0.5, 0.5}         // spacing
);
// Physical subdomain: [5,21) x [10,26) x [0,32)

Note

Bounds are calculated as: [origin + offset*dx, origin + (offset+size)*dx)

Use consistent origin/discretization across all subdomains for correctness