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))