def harmonizer_tester(n,m,harms,bw,Rsphere,mode="NONE"): """ This function shows the data produced by the harmonizer function at a number of data points given a certain set of harmonic components see 'anothertester' """ meas_points=[]#SHT.measurement_points(bw) for theta,phi in SHT.measurement_points(bw): meas_points.append(Vector([Rsphere,theta,phi],type="SPHERICAL")) d,rdata=SHT.harmonizer(meas_points,harms,mode=mode,pointType="VECTOR"); harmplotterMAT(bw,harms,mode=mode) l=SHT.SHT(bw,realdata=rdata,mode=mode,Rref=Rsphere) pprint(l) return
def multitester(): """ This function generates some data and then when the SHT is taken it should have the same harmonic components as were origininally generated """ bw=8; #harmtype='_{n,m}' harmtype='_{n}^{m}' Rsphere=0.08 meas_points=[] for theta,phi in SHT.measurement_points(bw): meas_points.append(Vector([Rsphere,theta,phi],type="SPHERICAL")) mode = "COS_SIN" harms=[((0,0),(1.0,0.0)), ((1,0),(1.0,0.0)), ((2,0),(1.0,0.0)), ((3,0),(1.0,0.0)), ((4,0),(2.0,0.0)), ((1,1),(1.0,1.0)), ((2,1),(1.0,1.0)), ((2,2),(0.0,1.0)), ((7,6),(0.0,20))] #harms=[((1,1),(1.0,0.0))] d,rdata=SHT.harmonizer(meas_points,harms,mode=mode,pointType="VECTOR",harmtype=harmtype) SHT.harmplotterMAT(bw,harms,Rsphere,mode=mode,harmtype=harmtype,type="3D") print "I am about to SHT" l=SHT.SHT(bw,realdata=rdata,mode=mode,Rref=Rsphere,harmtype=harmtype) print "I SHTed" pprint(l) return
def val(self,theta,phi): d,rdata=SHT.harmonizer([(theta,phi)],self.harms,mode=self.mode) result=float(rdata[0]) return result