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