Function almost_equal

Tests/Particles/NeighborParticles/MDParticleContainer.cpp:98–104

    };

    bool almost_equal (ParticleReal lhs, ParticleReal rhs)
    {
        auto const scale = std::max({ParticleReal(1.0), std::abs(lhs), std::abs(rhs)});
        auto const tol = std::max(ParticleReal(1.0e-12),
                                  ParticleReal(64.0) * std::numeric_limits<ParticleReal>::epsilon());
        return std::abs(lhs - rhs) <= tol * scale;
    }

    void get_position_unit_cell(Real* r, const IntVect& nppc, int i_part)
    {