Attempt formatting

    // Maps the field indices onto the range of -d/2 < x < d/2, where d is the scale of the field in coordinate x.

    // This means, {x,y,z} = (0,0,0) is in the center of the field.

    template <typename T, size_t N>

    bool DetectorField<T, N>::get_grid_index(size_t& index,

                                               const double x,

                                               const double y,

                                               const double z,

                                               const bool extrapolate_z) const noexcept {

    bool DetectorField<T, N>::get_grid_index(

        size_t& index, const double x, const double y, const double z, const bool extrapolate_z) const noexcept {

        // Compute indices

        // If the number of bins in x or y is 1, the field is assumed to be 2-dimensional and the respective index