Add tested features (Kronecker, isotropic) and their tests

.gitignore

@@ -1,3 +1,4 @@
+
 # Dev environment
 local_*
 
@@ -6,9 +7,24 @@ __pycache__/
 *.py[cod]
 *$py.class
 
+# IDEs and editors
+.vscode/
+.DS_Store
+
+# Local files
+examples/
+benchmarks/
+helper/
+TODO.md
+
+
 # C extensions
 *.so
 
+# Images
+
+*.png
+
 # Distribution / packaging
 .Python
 build/

linox/__init__.py

@@ -1,3 +1,5 @@
+# __init__.py
+
 r"""`linox`: Linear operators in JAX.
 
 This package provides a collection of linear operators for JAX, including:
@@ -23,19 +25,37 @@
 
 # Import functions from _arithmetic module
 from ._arithmetic import (
+    AddLinearOperator,
+    InverseLinearOperator,
+    ProductLinearOperator,
+    PseudoInverseLinearOperator,
+    ScaledLinearOperator,
+    TransposedLinearOperator,
     congruence_transform,
     diagonal,
     is_square,
     is_symmetric,
+    kron,
+    lcholesky,
+    ldet,
+    leigh,
     linverse,
+    lpinverse,
+    lpsolve,
+    lqr,
+    lsolve,
     lsqrt,
+    slogdet,
+    svd,
     symmetrize,
     transpose,
 )
 
 # Import classes from other modules
 from ._block import BlockDiagonal, BlockMatrix, BlockMatrix2x2
 from ._eigen import EigenD
+from ._isotropicadd import IsotropicAdditiveLinearOperator
+from ._kernel import ArrayKernel
 from ._kronecker import Kronecker
 from ._linear_operator import LinearOperator
 from ._low_rank import (
@@ -46,16 +66,21 @@
 )
 from ._matrix import Diagonal, Identity, Matrix, Ones, Scalar, Zero
 from ._permutation import Permutation
-from .utils import todense
+from ._toeplitz import Toeplitz
+from .utils import allclose, todense
 
 # Explicitly declare public API
 __all__ = [
+    "AddLinearOperator",
+    "ArrayKernel",
     "BlockDiagonal",
     "BlockMatrix",
     "BlockMatrix2x2",
     "Diagonal",
     "EigenD",
     "Identity",
+    "InverseLinearOperator",
+    "IsotropicAdditiveLinearOperator",
     "IsotropicScalingPlusSymmetricLowRank",
     "Kronecker",
     "LinearOperator",
@@ -64,15 +89,31 @@
     "Ones",
     "Permutation",
     "PositiveDiagonalPlusSymmetricLowRank",
+    "ProductLinearOperator",
+    "PseudoInverseLinearOperator",
     "Scalar",
+    "ScaledLinearOperator",
     "SymmetricLowRank",
+    "Toeplitz",
+    "TransposedLinearOperator",
     "Zero",
+    "allclose",
     "congruence_transform",
     "diagonal",
     "is_square",
     "is_symmetric",
+    "kron",
+    "lcholesky",
+    "ldet",
+    "leigh",
     "linverse",
+    "lpinverse",
+    "lpsolve",
+    "lqr",
+    "lsolve",
     "lsqrt",
+    "slogdet",
+    "svd",
     "symmetrize",
     "todense",
     "transpose",