scico
scico copied to clipboard
lanl
Scientific Computational Imaging COde
PR #154 will add dependency `cached_property` to `examples_requirements.txt`, required due to a bug in `xdesign` packaging. This issue is a reminder to undo this change once tomography/xdesign#77 has been addressed.
Arithmetic operations (+, -, *, ...) between a SCICO BlockArray and a NumPy ndarray are defined in the [BlockArray docs](https://scico.readthedocs.io/en/latest/_autosummary/scico.blockarray.html#module-scico.blockarray) to work in two cases: (1) when the ndarray is...
While not documented, the `ADMM` class `__init__` supports setting `f=None`. This capability is used in several example scripts, e.g. https://github.com/lanl/scico/blob/029a595e562420cb4e4559995d0b7808448d1c0d/examples/scripts/denoise_l1tv_admm.py#L62-L64 This support should either be documented or removed. The latter...
Current type annotation usage is in need of some attention: 1. Many functions and methods still do not have type annotations, e.g. `jvp`, `jhvp`, and `vjp` in `scico.operator`. 2. Figure...
It would be useful to have a general linop that can implement a stacked version of another linop using `jax.vmap`. An example application would be in constructing a 3D Abel...
We should decide what `input_dtype` and `output_dtype` mean, then document it. Proposed `input_dtype` and `output_dtype` policy: - All instantiated `Operator`s must have `input_dtype` and `output_dtype` properties - Possible data types...
In _real scientific imaging code_ one of our team members has wanted to sum each block of a block array along 2 of their 3 dimensions. Can we extend the...
Function `_gradient_filters` in `scico/linop/_circconv.py` does not appear to be used by any other code in the package or examples. It appears to be a vestige of some previous functionality that...
JAX works in single precision by default, and won't even let you create double precision arrays unless an environment variable is set (`JAX_ENABLE_X64=True`) or a special command is run when...
The example scripts could do with improvement on a number of fronts - [x] Where isotropic TV is used, change regularization term notation from `||Cx||_1` to `||Cx||_{2,1}`. - [x] Script...