[BUG] MR Deconvolve raises a segmentation fault

[sfarrens]

Running mr_deconv with the following option (i.e. those used in PySAP tests) raises a segmentation fault.

mr_deconv -d 5 -u 3 -n 2 -g 5 image.fits psf.fits res.mr

⬇️

[1]    31617 segmentation fault  mr_deconv -d 5 -u 3 -n 2 -g 5 image.fits psf.fits res.mr

The example image and PSF data can be obtained from PySAP.

import pysap
from pysap.data import get_sample_data
image = get_sample_data(dataset_name="astro-galaxy").data
psf = get_sample_data(dataset_name="astro-psf").data
pysap.io.save(image, "image.fits")
pysap.io.save(psf, "psf.fits")

This needs to be investigated and solved.

[jstarck]

Option -u is only valid in combination with -t24 (undecimated Wavelet Transform). So a fix could be to have the program returning a user error when -u is used without -t24.