def testf_CigaxgKnownValues(self):
"""Cigaxg should give known result with known input"""
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))
xgout = (cxg1.value,cxg2.value,cxg3.value,cxg4.value)
self.assertAlmostEqual(xgout[0],xg[0],1,name+'[0]'+repr(xgout))
self.assertAlmostEqual(xgout[1],xg[1],1,name+'[1]'+repr(xgout))
self.assertAlmostEqual(xgout[2],xg[2],1,name+'[2]'+repr(xgout))
self.assertAlmostEqual(xgout[3],xg[3],1,name+'[3]'+repr(xgout))
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))