sagemathgh-36594: Replace relative imports by absolute ones in `sage.{matrix,matroids,modules,stats}`
    
<!-- ^^^^^
Please provide a concise, informative and self-explanatory title.
Don't put issue numbers in there, do this in the PR body below.
For example, instead of "Fixes sagemath#1234" use "Introduce new method to
calculate 1+1"
-->
<!-- Describe your changes here in detail -->

<!-- Why is this change required? What problem does it solve? -->
<!-- If this PR resolves an open issue, please link to it here. For
example "Fixes sagemath#12345". -->
- Cherry-picked from sagemath#35095
- This is for sagemath#36228
<!-- If your change requires a documentation PR, please link it
appropriately. -->

### 📝 Checklist

<!-- Put an `x` in all the boxes that apply. -->
<!-- If your change requires a documentation PR, please link it
appropriately -->
<!-- If you're unsure about any of these, don't hesitate to ask. We're
here to help! -->
<!-- Feel free to remove irrelevant items. -->

- [x] The title is concise, informative, and self-explanatory.
- [ ] The description explains in detail what this PR is about.
- [ ] I have linked a relevant issue or discussion.
- [ ] I have created tests covering the changes.
- [ ] I have updated the documentation accordingly.

### ⌛ Dependencies

<!-- List all open PRs that this PR logically depends on
- sagemath#12345: short description why this is a dependency
- sagemath#34567: ...
-->

<!-- If you're unsure about any of these, don't hesitate to ask. We're
here to help! -->
    
URL: sagemath#36594
Reported by: Matthias Köppe
Reviewer(s): David Coudert, Matthias Köppe

Author: Release Manager

src/sage/matrix/action.pyx

@@ -62,7 +62,7 @@ AUTHOR:
 
 import operator
 
-from .matrix_space import MatrixSpace, is_MatrixSpace
+from sage.matrix.matrix_space import MatrixSpace, is_MatrixSpace
 from sage.modules.free_module import FreeModule, is_FreeModule
 from sage.structure.coerce cimport coercion_model
 from sage.categories.homset import Hom, End

src/sage/matrix/args.pyx

@@ -983,7 +983,7 @@ cdef class MatrixArgs:
         if self.space is None:
             global MatrixSpace
             if MatrixSpace is None:
-                from .matrix_space import MatrixSpace
+                from sage.matrix.matrix_space import MatrixSpace
             self.space = MatrixSpace(self.base, self.nrows, self.ncols,
                     sparse=self.sparse, **self.kwds)
 
@@ -1114,7 +1114,7 @@ cdef class MatrixArgs:
         if not (e.flags.c_contiguous is True or e.flags.f_contiguous is True):
             raise TypeError('numpy matrix must be either c_contiguous or f_contiguous')
 
-        from .constructor import matrix
+        from sage.matrix.constructor import matrix
         from sage.rings.real_double import RDF
         from sage.rings.complex_double import CDF
 

src/sage/matrix/change_ring.pyx

@@ -2,9 +2,9 @@
 Functions for changing the base ring of matrices quickly
 """
 
-from .matrix_space import MatrixSpace
-from .matrix_real_double_dense cimport Matrix_real_double_dense
-from .matrix_integer_dense cimport Matrix_integer_dense
+from sage.matrix.matrix_space import MatrixSpace
+from sage.matrix.matrix_real_double_dense cimport Matrix_real_double_dense
+from sage.matrix.matrix_integer_dense cimport Matrix_integer_dense
 
 from sage.rings.real_double import RDF
 

src/sage/matrix/constructor.pyx

@@ -14,7 +14,7 @@ General matrix Constructor and display options
 #                  http://www.gnu.org/licenses/
 #*****************************************************************************
 
-from .args cimport MatrixArgs
+from sage.matrix.args cimport MatrixArgs
 from sage.structure.global_options import GlobalOptions
 
 
@@ -654,7 +654,7 @@ def matrix(*args, **kwds):
 Matrix = matrix
 
 
-from .special import *
+from sage.matrix.special import *
 
 
 @matrix_method

src/sage/matrix/matrix.pxd

@@ -1 +1 @@
-from .matrix2 cimport Matrix
+from sage.matrix.matrix2 cimport Matrix

src/sage/matrix/matrix0.pyx

@@ -47,7 +47,7 @@ from sage.rings.integer_ring import is_IntegerRing
 
 import sage.modules.free_module
 
-from .matrix_misc import row_iterator
+from sage.matrix.matrix_misc import row_iterator
 
 _Fields = Fields()
 _IntegralDomains = IntegralDomains()
@@ -1750,7 +1750,7 @@ cdef class Matrix(sage.structure.element.Matrix):
              100 x 100 dense matrix over Integer Ring]
 
         """
-        from .constructor import options
+        from sage.matrix.constructor import options
         if self._nrows <= options.max_rows() and self._ncols <= options.max_cols():
             return self.str()
         if self.is_sparse():
@@ -2051,7 +2051,7 @@ cdef class Matrix(sage.structure.element.Matrix):
 
         # only use floating point formatting for inexact types that have
         # custom implementation of __format__
-        from .constructor import options
+        from sage.matrix.constructor import options
         prec = options.precision()
         if prec is None or callable(rep_mapping) or not entries \
                 or type(entries[0]).__format__ is Element.__format__ \
@@ -2135,7 +2135,7 @@ cdef class Matrix(sage.structure.element.Matrix):
               -0.35104242112828943    0.5084492941557279]
                -0.9541798283979341   -0.8948790563276592]
         """
-        from .constructor import options
+        from sage.matrix.constructor import options
         if self._nrows <= options.max_rows() and self._ncols <= options.max_cols():
             return self.str(character_art=True)
         else:
@@ -2164,7 +2164,7 @@ cdef class Matrix(sage.structure.element.Matrix):
             sage: unicode_art(A)
             100 x 100 dense matrix over Integer Ring
         """
-        from .constructor import options
+        from sage.matrix.constructor import options
         if self._nrows <= options.max_rows() and self._ncols <= options.max_cols():
             return self.str(unicode=True, character_art=True)
         else:
@@ -2375,7 +2375,7 @@ cdef class Matrix(sage.structure.element.Matrix):
             [  5  25]
             [125 625]
         """
-        from .constructor import matrix
+        from sage.matrix.constructor import matrix
         return matrix(self.nrows(), self.ncols(), [e(*args, **kwargs) for e in self.list()])
 
     ###################################################
@@ -3850,7 +3850,7 @@ cdef class Matrix(sage.structure.element.Matrix):
             raise ValueError("length of v must be at most the number of rows of self")
         if not self._nrows:
             return self.parent().row_space().zero_vector()
-        from .constructor import matrix
+        from sage.matrix.constructor import matrix
         v = matrix(list(v)+[0]*(self._nrows-len(v)))
         return (v * self)[0]
 
@@ -3927,7 +3927,7 @@ cdef class Matrix(sage.structure.element.Matrix):
             raise ValueError("length of v must be at most the number of columns of self")
         if not self._ncols:
             return self.parent().column_space().zero_vector()
-        from .constructor import matrix
+        from sage.matrix.constructor import matrix
         v = matrix(self._ncols, 1, list(v)+[0]*(self._ncols-len(v)))
         return (self * v).column(0)
 
@@ -6178,7 +6178,7 @@ def set_max_rows(n):
     """
     from sage.misc.superseded import deprecation
     deprecation(30552, "'set_max_rows' is replaced by 'matrix.options.max_rows'")
-    from .constructor import options
+    from sage.matrix.constructor import options
     options.max_rows = n-1
 
 def set_max_cols(n):
@@ -6195,5 +6195,5 @@ def set_max_cols(n):
     """
     from sage.misc.superseded import deprecation
     deprecation(30552, "'set_max_cols' is replaced by 'matrix.options.max_cols'")
-    from .constructor import options
+    from sage.matrix.constructor import options
     options.max_cols = n-1

src/sage/matrix/matrix1.pxd

@@ -1,4 +1,4 @@
-from .matrix0 cimport Matrix as Matrix0
+from sage.matrix.matrix0 cimport Matrix as Matrix0
 
 cdef class Matrix(Matrix0):
     cdef _stack_impl(self, bottom) noexcept

src/sage/matrix/matrix2.pxd

@@ -12,7 +12,7 @@ Generic matrices
 #                  http://www.gnu.org/licenses/
 #*****************************************************************************
 
-from .matrix1 cimport Matrix as Matrix1
+from sage.matrix.matrix1 cimport Matrix as Matrix1
 
 cdef class Matrix(Matrix1):
     cdef _det_by_minors(self, Py_ssize_t level) noexcept

src/sage/matrix/matrix2.pyx

@@ -99,7 +99,7 @@ from sage.arith.numerical_approx cimport digits_to_bits
 from copy import copy
 
 import sage.modules.free_module
-from . import berlekamp_massey
+from sage.matrix import berlekamp_massey
 from sage.modules.free_module_element import is_FreeModuleElement
 from sage.matrix.matrix_misc import permanental_minor_polynomial
 
@@ -8832,7 +8832,7 @@ cdef class Matrix(Matrix1):
         output_window = output.matrix_window()
 
 
-        from . import strassen
+        from sage.matrix import strassen
         strassen.strassen_window_multiply(output_window, self_window, right_window, cutoff)
         return output
 
@@ -8870,7 +8870,7 @@ cdef class Matrix(Matrix1):
             self._echelon_in_place_classical()
             return
 
-        from . import strassen
+        from sage.matrix import strassen
         pivots = strassen.strassen_echelon(self.matrix_window(), cutoff)
         self.cache('pivots', pivots)
         verbose('done with strassen', tm)
@@ -8909,7 +8909,7 @@ cdef class Matrix(Matrix1):
             ...
             IndexError: matrix window index out of range
         """
-        from . import matrix_window
+        from sage.matrix import matrix_window
         if nrows == -1:
             nrows = self._nrows - row
         if ncols == -1:
@@ -15117,7 +15117,7 @@ cdef class Matrix(Matrix1):
         except (OverflowError, TypeError):
             from sage.rings.real_mpfr import RealField
             # Try using MPFR, which handles large numbers much better, but is slower.
-            from .misc_mpfr import hadamard_row_bound_mpfr
+            from sage.matrix.misc_mpfr import hadamard_row_bound_mpfr
             R = RealField(53, rnd='RNDU')
             A = self.change_ring(R)
             m1 = hadamard_row_bound_mpfr(A)

src/sage/matrix/matrix_complex_ball_dense.pxd

@@ -1,5 +1,5 @@
 from sage.libs.arb.types cimport acb_mat_t
-from .matrix_dense cimport Matrix_dense
+from sage.matrix.matrix_dense cimport Matrix_dense
 from sage.matrix.matrix_generic_dense cimport Matrix_generic_dense
 from sage.structure.parent cimport Parent