Subsets : SIRF algebra

[epapoutsellis]

For a subset acquisition data acq_subset the following

tmp = acq_subset*0.0

raises this error

ERROR: write_basic_interfile_PDFS_header: Error casting the projdata to one of its geometries: Cylindrical/BlocksOnCylindrical/Genreic
---------------------------------------------------------------------------
error                                     Traceback (most recent call last)
Input In [24], in <cell line: 1>()
----> 1 tmp = acq_subset*0.0

File ~/Edo/build/INSTALL/python/sirf/SIRF.py:397, in DataContainer.__mul__(self, other)
    394     z.handle = pysirf.cSIRF_axpby \
    395         (a.ctypes.data, self.handle, zero.ctypes.data, self.handle)
    396     z.src = 'mult'
--> 397     check_status(z.handle)
    398     return z
    400 return NotImplemented

File ~/Edo/build/INSTALL/python/sirf/Utilities.py:459, in check_status(handle, stack)
    454 line = pyiutil.executionErrorLine(handle)
    455 errorMsg = \
    456     repr(msg) + ' exception caught at line ' + \
    457     repr(line) + ' of ' + file + '; ' + \
    458     'the reconstruction engine output may provide more information'
--> 459 raise error(errorMsg)

error: ??? "'write_basic_interfile_PDFS_header: Error casting the projdata to one of its geometries: Cylindrical/BlocksOnCylindrical/Genreic' exception caught at line 162 of /home/vaggelis/Edo/build/sources/SIRF/src/common/csirf.cpp; the reconstruction engine output may provide more information"

@evgueni-ovtchinnikov @KrisThielemans @paskino

related to #1028

[evgueni-ovtchinnikov]

@epapoutsellis the issue is being investigated, meanwhile you may like to try memory storage scheme, which seems to work.

[epapoutsellis]

Notebook using the subsets functionality is here.