svmbir
svmbir copied to clipboard
Fast code for parallel or fan beam tomographic reconstruction
The current C code reduces the number of parallel threads to stabilize convergence when the positivity constraint is turned off. If we implement the relaxation factor for voxel update, then...
Current code sets auto_sigma_p to be `1.0 * auto_sigma_prior(sino, delta_channel, sharpness)` (1/5 of auto_sigma_x). In the cone beam code we found out that `2.0 * auto_sigma_prior(sino, delta_channel, sharpness)` is generally...
I pushed a branch `oo_interface` that includes a definition of a new callable class `Recon` that wraps `svmbir.recon()`. What it currently does is very simple. Instances of Recon have the...
Right now the `_sino_indicator()` function does not use the weight values. This could be a problem if there are large erroneous values in the sinogram with weights of 0. Should...
It seems very reasonable to iterate over parameter arrays such as ```python num_threads = np.array([2, 4]) # num_threads[n] is of type numpy.int64 for n in range(2): result[n] = svmbir.recon( ...,...
Add functions to "list" and "clear" sysmatrix files from cache directory. We should use "touch" or something equivalent to update the time stamp of the sysmatrix file each time it's...
Hi! Currently svmbir prints to console all it's messages. This makes ist somewhat hacky to use inside other projects, which have a propper logging interface or a CLI. There basically...