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
