Inconsistent results with gfortran

In [test1](https://github.com/gha3mi/pinverse/blob/main/test/test1.f90), the pinv function in MATLAB was compared with this implementation. The results were consistent with the Intel compilers (ifort (IFORT) 2021.8.0 20221119 and ifx (IFORT) 2023.0.0 20221201) and NVIDIA (nvfortran 23.1-0). However, there were differences when using gfortran (GNU Fortran (GCC) 13.1.0).