pygeoinf 1.2.6__tar.gz → 1.2.8__tar.gz

{pygeoinf-1.2.6 → pygeoinf-1.2.8}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: pygeoinf
-Version: 1.2.6
+Version: 1.2.8
 Summary: A package for solving geophysical inference and inverse problems
 License: BSD-3-Clause
 Author: David Al-Attar and Dan Heathcote
@@ -64,22 +64,15 @@ git clone https://github.com/da380/pygeoinf.git
 cd pygeoinf
 poetry install
 ```
-
-You can install the optional dependencies for running tests, building documentation, or running the tutorials by using the --with flag.
+You can install all optional dependencies for development—including tools for running the test suite, 
+building the documentation, and running the Jupyter tutorials—by using the ```--with``` flag and specifying the ```dev``` group.
 
 ```bash
-# To install dependencies for running the test suite
-poetry install --with tests
-
-# To install dependencies for building the documentation
-poetry install --with docs
+# Install all development dependencies (for tests, docs, and tutorials)
+poetry install --with dev
+```
 
-# To install dependencies for running the Jupyter tutorials
-poetry install --with tutorials
 
-# You can also combine them
-poetry install --with tests,docs,tutorials
-```
 
 ## Documentation
 
@@ -245,9 +238,9 @@ The output of the above script will look similar to the following figure:
 
 ## Future Plans
 
-`pygeoinf` is under active development. Future work will focus on expanding the library's capabilities to address a broader range of geophysical problems. Key areas for development include:
+`pygeoinf` is under active development. Current work is focused on expanding the library's capabilities to address a broader range of geophysical problems. Key areas for development include:
 
-* **Generalised Backus-Gilbert Methods**: Implementation of a generalised Backus-Gilbert framework for linear inference problems, building on the work of [Al-Attar(2021)](https://arxiv.org/abs/2104.12256). The focus will be on constructing direct estimates of specific properties of interest (i.e., linear functionals of the model) from data, without needing to first solve for the full model itself.
+* **Generalised Backus-Gilbert Methods**: Implementation of a generalised Backus-Gilbert framework for linear inference problems. The focus will be on constructing direct estimates of specific properties of interest (i.e., linear functionals of the model) from data, without needing to first solve for the full model itself.
 
 * **Non-linear Optimisation**: Extension of the current optimisation framework to handle non-linear inverse problems. This will involve creating a general interface where users can provide their own non-linear forward mapping, a misfit functional, and methods for computing gradients (and optionally Hessians) for use in gradient-based optimisation algorithms.
 

{pygeoinf-1.2.6 → pygeoinf-1.2.8}/README.md

@@ -43,22 +43,15 @@ git clone https://github.com/da380/pygeoinf.git
 cd pygeoinf
 poetry install
 ```
-
-You can install the optional dependencies for running tests, building documentation, or running the tutorials by using the --with flag.
+You can install all optional dependencies for development—including tools for running the test suite, 
+building the documentation, and running the Jupyter tutorials—by using the ```--with``` flag and specifying the ```dev``` group.
 
 ```bash
-# To install dependencies for running the test suite
-poetry install --with tests
-
-# To install dependencies for building the documentation
-poetry install --with docs
+# Install all development dependencies (for tests, docs, and tutorials)
+poetry install --with dev
+```
 
-# To install dependencies for running the Jupyter tutorials
-poetry install --with tutorials
 
-# You can also combine them
-poetry install --with tests,docs,tutorials
-```
 
 ## Documentation
 
@@ -224,9 +217,9 @@ The output of the above script will look similar to the following figure:
 
 ## Future Plans
 
-`pygeoinf` is under active development. Future work will focus on expanding the library's capabilities to address a broader range of geophysical problems. Key areas for development include:
+`pygeoinf` is under active development. Current work is focused on expanding the library's capabilities to address a broader range of geophysical problems. Key areas for development include:
 
-* **Generalised Backus-Gilbert Methods**: Implementation of a generalised Backus-Gilbert framework for linear inference problems, building on the work of [Al-Attar(2021)](https://arxiv.org/abs/2104.12256). The focus will be on constructing direct estimates of specific properties of interest (i.e., linear functionals of the model) from data, without needing to first solve for the full model itself.
+* **Generalised Backus-Gilbert Methods**: Implementation of a generalised Backus-Gilbert framework for linear inference problems. The focus will be on constructing direct estimates of specific properties of interest (i.e., linear functionals of the model) from data, without needing to first solve for the full model itself.
 
 * **Non-linear Optimisation**: Extension of the current optimisation framework to handle non-linear inverse problems. This will involve creating a general interface where users can provide their own non-linear forward mapping, a misfit functional, and methods for computing gradients (and optionally Hessians) for use in gradient-based optimisation algorithms.
 

{pygeoinf-1.2.6 → pygeoinf-1.2.8}/pygeoinf/__init__.py

@@ -9,6 +9,7 @@ from .random_matrix import (
     random_svd,
     random_eig,
     random_cholesky,
+    random_diagonal,
 )
 
 from .hilbert_space import (
@@ -34,7 +35,11 @@ from .nonlinear_operators import NonLinearOperator
 
 from .linear_operators import (
     LinearOperator,
-    DiagonalLinearOperator,
+    MatrixLinearOperator,
+    DenseMatrixLinearOperator,
+    SparseMatrixLinearOperator,
+    DiagonalSparseMatrixLinearOperator,
+    NormalSumOperator,
 )
 
 

{pygeoinf-1.2.6 → pygeoinf-1.2.8}/pygeoinf/gaussian_measure.py

@@ -28,11 +28,11 @@ from scipy.sparse import diags
 from scipy.stats import multivariate_normal
 
 
-from .hilbert_space import EuclideanSpace, HilbertModule
+from .hilbert_space import EuclideanSpace, HilbertModule, Vector
 
 from .linear_operators import (
     LinearOperator,
-    DiagonalLinearOperator,
+    DiagonalSparseMatrixLinearOperator,
 )
 
 from .direct_sum import (
@@ -40,17 +40,6 @@ from .direct_sum import (
 )
 
 
-# This block only runs for type checkers, not at runtime, to prevent
-# circular import errors while still allowing type hints.
-if TYPE_CHECKING:
-    from .hilbert_space import HilbertSpace, EuclideanSpace
-    from .operators import LinearOperator, DiagonalLinearOperator
-    from .direct_sum import BlockDiagonalLinearOperator
-
-# Define a generic type for vectors in a Hilbert space
-Vector = TypeVar("Vector")
-
-
 class GaussianMeasure:
     """
     Represents a Gaussian measure on a Hilbert space.
@@ -76,29 +65,29 @@ class GaussianMeasure:
         inverse_covariance_factor: LinearOperator = None,
     ) -> None:
         """
-        Initializes the GaussianMeasure.
+         Initializes the GaussianMeasure.
 
         The measure can be defined in several ways, primarily by providing
-        either a covariance operator or a covariance factor.
-
-        Args:
-            covariance (LinearOperator, optional): A self-adjoint and positive
-                semi-definite linear operator on the domain.
-            covariance_factor (LinearOperator, optional): A linear operator L
-                such that the covariance C = L @ L*.
-            expectation (vector, optional): The expectation (mean) of the
-                measure. Defaults to the zero vector of the space.
-            sample (callable, optional): A function that returns a random
-                sample from the measure. If a `covariance_factor` is given,
-                a default sampler is created.
-            inverse_covariance (LinearOperator, optional): The inverse of the
-                covariance operator (the precision operator).
-            inverse_covariance_factor (LinearOperator, optional): A factor Li
-                of the inverse covariance, such that C_inv = Li.T @ Li.
-
-        Raises:
-            ValueError: If neither `covariance` nor `covariance_factor`
-                is provided.
+         either a covariance operator or a covariance factor.
+
+         Args:
+             covariance (LinearOperator, optional): A self-adjoint and positive
+                 semi-definite linear operator on the domain.
+             covariance_factor (LinearOperator, optional): A linear operator L
+                 such that the covariance C = L @ L*.
+             expectation (vector, optional): The expectation (mean) of the
+                 measure. Defaults to the zero vector of the space.
+             sample (callable, optional): A function that returns a random
+                 sample from the measure. If a `covariance_factor` is given,
+                 a default sampler is created.
+             inverse_covariance (LinearOperator, optional): The inverse of the
+                 covariance operator (the precision operator).
+             inverse_covariance_factor (LinearOperator, optional): A factor Li
+                 of the inverse covariance, such that C_inv = Li.T @ Li.
+
+         Raises:
+             ValueError: If neither `covariance` nor `covariance_factor`
+                 is provided.
         """
         if covariance is None and covariance_factor is None:
             raise ValueError(
@@ -186,8 +175,8 @@ class GaussianMeasure:
                 "Standard deviation vector does not have the correct length"
             )
         euclidean = EuclideanSpace(domain.dim)
-        covariance_factor = DiagonalLinearOperator(
-            euclidean, domain, standard_deviations, galerkin=True
+        covariance_factor = DiagonalSparseMatrixLinearOperator.from_diagonal_values(
+            euclidean, domain, standard_deviations
         )
         return GaussianMeasure(
             covariance_factor=covariance_factor,