[SyneRBI/SIRF] Subsets : SIRF algebra (Issue #1111)

Vaggelis Papoutsellis

unread,

May 23, 2022, 9:09:55 AM5/23/22

to SyneRBI/SIRF, Subscribed

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"

—
Reply to this email directly, view it on GitHub, or unsubscribe.
You are receiving this because you are subscribed to this thread.

Evgueni Ovtchinnikov

unread,

Jun 6, 2022, 12:23:13 PM6/6/22

to SyneRBI/SIRF, Subscribed

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

—
Reply to this email directly, view it on GitHub, or unsubscribe.

You are receiving this because you are subscribed to this thread.

Vaggelis Papoutsellis

unread,

Jun 23, 2022, 7:53:59 AM6/23/22

to SyneRBI/SIRF, Subscribed

Notebook using the subsets functionality is here.

—
Reply to this email directly, view it on GitHub, or unsubscribe.

You are receiving this because you are subscribed to this thread.

Kris Thielemans

unread,

Oct 23, 2025, 5:26:26 AM10/23/25

to SyneRBI/SIRF, Subscribed

KrisThielemans left a comment (SyneRBI/SIRF#1111)

would be good to add a test for this as per https://github.com/SyneRBI/SIRF/blob/82aecc5574360c9486484cb2722d2befddf4ee9d/src/xSTIR/pSTIR/tests/test_algebra.py#L72-L73
adding get_subset([0])

—
Reply to this email directly, view it on GitHub, or unsubscribe.

You are receiving this because you are subscribed to this thread.

Edoardo Pasca

unread,

Oct 24, 2025, 8:20:26 AM10/24/25

to SyneRBI/SIRF, Subscribed

paskino left a comment (SyneRBI/SIRF#1111)

#1342 does not close this issue. I would suggest to keep this open

—
Reply to this email directly, view it on GitHub, or unsubscribe.

You are receiving this because you are subscribed to this thread.

Edoardo Pasca

unread,

Oct 27, 2025, 6:18:01 AM10/27/25

to SyneRBI/SIRF, Subscribed

paskino left a comment (SyneRBI/SIRF#1111)

With #1342 unit tests for subset algebra are added for both storage scheme file and memory. For the first one the unit tests expect that they raise an exception.