I wonder if it is possible to simply use the all_close method from ArrayBase provided by the ndarray crate itself to simplify the implementation of the trait of this library. It also offers an eps value for approximate equality.
Have you considered using it?