def tst_integral_triple_zernike2d(nmax):
Bnmk_obj = Bnmk(nmax)
#Bnmk_obj.print_bnmk()
coef_table = []
for m in range(nmax+1):
C_m_3n = math.nmk_array(nmax)
for n1 in range(m,nmax+1,2):
for n2 in range(m,n1+1,2):
for n3 in range(m,n2+1,2):
value = integrate_triple_zernike2d(n1,n2,n3,m,Bnmk_obj)
C_m_3n.set_coef(n1,n2,n3,value)
# print m,n1,n2,n3,value.real
coef_table.append( C_m_3n )
return coef_table
def __init__(self, nmax): self.nmax=nmax self.Bnmk=math.nmk_array(nmax) self.initialize_bnmk()
