Array-style

indexing

pfc::DiscreteField<double, 3> field({64, 64, 64}, {0,0,0},
                                     {0.0,0.0,0.0}, {1.0,1.0,1.0});
 
// Access by 3D index
field[{10, 20, 30}] = 1.0;
 
// Access by linear index
field[0] = 0.5;  // First element
 
// Iterate over all elements
for (size_t i = 0; i < field.get_data().size(); i++) {
    field[i] *= 2.0;
}