def testf_CxgaigKnownValues(self):
"""Cxgaig should give known result with known input"""
for name,proj,dims,xg,ig in self.knownValues:
(cxg1,cxg2,cxg3,cxg4) = (ct.c_float(xg[0]),ct.c_float(xg[1]),ct.c_float(xg[2]),ct.c_float(xg[3]))
(cig1,cig2,cig3,cig4) = (ct.c_int(0),ct.c_int(0),ct.c_int(0),ct.c_int(0))
istat = rmn.f_cxgaig(proj,
ct.byref(cig1),ct.byref(cig2),
ct.byref(cig3),ct.byref(cig4),
ct.byref(cxg1),ct.byref(cxg2),
ct.byref(cxg3),ct.byref(cxg4))
igout = (cig1.value,cig2.value,cig3.value,cig4.value)
self.assertEqual(igout,ig,name+repr(ig)+' != '+repr(igout))
def testf_Sanity(self):
"""cigaxg(cxgaig(n))==n for all n"""
for name, proj, dims, xg, ig in self.knownValues:
(cxg1, cxg2, cxg3, cxg4) = (ct.c_float(0.), ct.c_float(0.),
ct.c_float(0.), ct.c_float(0.))
(cig1, cig2, cig3, cig4) = (ct.c_int(ig[0]), ct.c_int(ig[1]),
ct.c_int(ig[2]), ct.c_int(ig[3]))
istat = rmn.f_cigaxg(
proj, ct.byref(cxg1), ct.byref(cxg2), ct.byref(cxg3),
ct.byref(cxg4), ct.byref(cig1), ct.byref(cig2), ct.byref(cig3),
ct.byref(cig4))
istat = rmn.f_cxgaig(
proj, ct.byref(cig1), ct.byref(cig2), ct.byref(cig3),
ct.byref(cig4), ct.byref(cxg1), ct.byref(cxg2), ct.byref(cxg3),
ct.byref(cxg4))
igout = (cig1.value, cig2.value, cig3.value, cig4.value)
self.assertEqual(igout, ig, name + repr(ig) + ' != ' + repr(igout))
