Add __all__ lists to slinalg and nlinalg

aesara/tensor/nlinalg.py

@@ -821,3 +821,21 @@ def tensorsolve(a, b, axes=None):
     """
 
     return TensorSolve(axes)(a, b)
+
+
+__all__ = [
+    "pinv",
+    "inv",
+    "trace",
+    "matrix_dot",
+    "det",
+    "eig",
+    "eigh",
+    "qr",
+    "svd",
+    "lstsq",
+    "matrix_power",
+    "norm",
+    "tensorinv",
+    "tensorsolve",
+]

aesara/tensor/slinalg.py

@@ -597,3 +597,14 @@ class ExpmGrad(Op):
 
 
 expm = Expm()
+
+__all__ = [
+    "cholesky",
+    "solve",
+    "solve_lower_triangular",
+    "solve_upper_triangular",
+    "solve_symmetric",
+    "eigvalsh",
+    "kron",
+    "expm",
+]