Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Use a block sparse format for the observation operator #2

Open
wants to merge 39 commits into
base: master
Choose a base branch
from
Open

Use a block sparse format for the observation operator #2

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
39 commits
Select commit Hold shift + click to select a range
6903568
First stab at writing functions to put obs op in BSR form.
DWesl Sep 20, 2018
b55e4bc
Mark tests for new features as expected failures.
DWesl Nov 19, 2018
e795464
Return the influence functions as bsr_matrix directly, rather than wr…
DWesl Mar 4, 2019
e9fbff8
Test that the bsr matrix represents things in the right places.
DWesl Mar 4, 2019
d9dab0b
Remove code for DataArray, since support for that was dropped.
DWesl Mar 4, 2019
bc90450
Ignore additional testing artifacts.
DWesl Mar 4, 2019
0a471c6
Add a function to calculate the product of a separable covariance and…
DWesl Mar 5, 2019
e044a63
Test that the new quadratic form function is doing the right thing.
DWesl Mar 5, 2019
5527380
Set up the framework to get this used.
DWesl Mar 5, 2019
a57cde5
Get tests working on python2.
DWesl Mar 5, 2019
1c9f21f
Test that the new quadratic form fails properly when it won't work.
DWesl Mar 5, 2019
2ef6b90
Fix a few style issues in the inversion code.
DWesl Mar 5, 2019
8a52598
Simplify the handling of reduced-resolution variants in method_common.
DWesl Mar 5, 2019
65ce042
Allow requesting a shape for the influence function and start using t…
DWesl Mar 6, 2019
f8dfccc
Add the remappers and influence function helpers to the documentation.
DWesl Mar 6, 2019
2c57bbc
A first attempt to resample BSR-format data.
DWesl Mar 7, 2019
fcdec41
Get the BSR observation operators used in the inversion.
DWesl Mar 8, 2019
97a279a
Use a lower spatial and higher temporal resolution for uncertainties.
DWesl Mar 8, 2019
516bde6
Add more print statements to try to narrow down where the slowdown oc…
DWesl Mar 10, 2019
33b500d
Explain reasoning for ordering decisions in temporal remapping code.
DWesl Mar 10, 2019
03c073f
Be more explicit about the shape tests.
DWesl Mar 10, 2019
bbae686
Add a few more checks to the code style test run.
DWesl Mar 10, 2019
794d353
Include a description for 'Helpers for real-data inversions'.
DWesl Mar 10, 2019
34183fd
Wrap a few overlong lines in the pseudo-data inversion.
DWesl Mar 10, 2019
6498a93
Merge branch 'master' into bsr_obs_op
DWesl Mar 15, 2019
455bf64
Use sparse matrices for the influence function resolution coarsening.
DWesl Mar 20, 2019
3bb9526
Remove several now-unneeded print statements.
DWesl Mar 20, 2019
a05d9b1
Update the timings in the run script.
DWesl Mar 20, 2019
ef3140b
Wrap some long lines and otherwise make flake8 happy.
DWesl Mar 20, 2019
f0970c6
Merge the style fixups into bsr_obs_op on server.
DWesl Mar 20, 2019
127293e
Check that the returned matrix stores flux times in increasing order.
DWesl Mar 20, 2019
9247958
Raise an error on unsupported inputs to cholesky.
DWesl Mar 22, 2019
6bb0fe6
Store the variational descent methods in constants.
DWesl Mar 22, 2019
ca58b84
Test the output with a BSR matrix against that with a dense matrix.
DWesl Mar 22, 2019
1979922
Decrease default planner effort to speed tests.
DWesl Mar 25, 2019
766938f
Stop reversing the shape of a square matrix.
DWesl Apr 4, 2019
af312fa
Test a few more postconditions of the covariance.
DWesl Apr 4, 2019
a8dc1c6
Include a transpose for a _SumLinearOperator.
DWesl Apr 8, 2019
a4b5800
Fix a syntax error in the tests.
DWesl Apr 8, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion .gitignore
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ build/

# Test artifacts
.tox/
.coverage
.coverage*
*_flymake.py
.mypy_cache

Expand Down
7 changes: 7 additions & 0 deletions doc/source/inversion.observation_operator.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
inversion\.observation\_operator module
=======================================

.. automodule:: inversion.observation_operator
:members:
:undoc-members:
:show-inheritance:
7 changes: 7 additions & 0 deletions doc/source/inversion.remapper.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
inversion\.remapper module
==========================

.. automodule:: inversion.remapper
:members:
:undoc-members:
:show-inheritance:
18 changes: 15 additions & 3 deletions doc/source/inversion.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -19,9 +19,9 @@ Inversion functions

All inversion functions have the same signature and give similar
answers: the difference is how they get there. PSAS and Variational
methods use iterative solvers. Optimal Interpolation uses a
Gauss-Jordan solver. Variational methods use a different but
equivalent formulation of the problem.
methods use iterative solvers for linear systems, where Optimal
Interpolation uses Gauss-Jordan elimination. Variational methods use
a different but equivalent formulation of the problem.

.. toctree::
inversion.optimal_interpolation
Expand All @@ -34,6 +34,18 @@ High-level wrappers
.. toctree::
inversion.wrapper

Helpers for real-data inversions
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Helper functions to calculate the prior error covariance and
observation operator at reduced resolution so the inversion routines
can calculate the posterior error covariance in reasonable time and in
a reasonable amount of storage.

.. toctree::
inversion.observation_operator
inversion.remapper

Other utilities
~~~~~~~~~~~~~~~

Expand Down
Loading