Coordinate-space

operations

pfc::DiscreteField<double, 3> field({32, 32, 32}, {0,0,0},
                                     {0.0,0.0,0.0}, {0.5,0.5,0.5});
 
// Map index to physical coordinate
auto coords = field.map_indices_to_coordinates({10, 10, 10});
// coords = {5.0, 5.0, 5.0}  (10 * 0.5 spacing)
 
// Map coordinate to nearest index
auto indices = field.map_coordinates_to_indices({7.3, 8.1, 9.8});
// indices = {15, 16, 20}  (rounded to nearest)
 
// Check if coordinate is in bounds
bool valid = field.inbounds({12.0, 8.0, 4.0});