ModeCalc.py File Reference

Go to the source code of this file.

Namespaces
namespace	ModeCalc

Functions
	ModeCalc.arctan2_2pi (yy, xx)
	ModeCalc.sph_design (func, t=5, author='Hardin', **kwargs)
	ModeCalc.gauss_legendre (func, N=40, x_low=0, x_up=pi, **kwargs)
	Gauss-Legendre 1D integration.
	ModeCalc.trapezoidal (func, N=40, x_low=0, x_up=pi, **kwargs)
	Gauss 1D integration.
	ModeCalc.prod_quad (func, N=20, M=40, **kwargs)
	Gaussian product quadrature using Trapezoidal for azimuthal direction and Gauss-Legendre for polar angle.
	ModeCalc.file_len (fname)

Variables
	ModeCalc.pi = np.pi
str	ModeCalc.scheme = "t-design"
	main code scheme="Gauss_Legendre"
str	ModeCalc.file_name = "zz.csv"
	ModeCalc.f = open(file_name)
	ModeCalc.dat = f.readlines()[1:]
int	ModeCalc.tstep = 0
int	ModeCalc.told = 1e6
	ModeCalc.tem = list(map(float,line.split(', ')))
	ModeCalc.tnew = tem[0]
int	ModeCalc.N_tracer = 0
int	ModeCalc.ID1 = -100
list	ModeCalc.coors = []
list	ModeCalc.coor = [tem[3],tem[4],tem[5]]
	ModeCalc.data_set = np.zeros((N_tracer,tstep,7))
int	ModeCalc.t_count = 0
int	ModeCalc.N_count = 0
	ModeCalc.thetas
	ModeCalc.ws
float	ModeCalc.zs = np.cos(np.pi*thetas/2 + np.pi/2)*0.4
	ModeCalc.zip_iterator = zip(zs,ws)
	ModeCalc.w_dict = dict(zip_iterator)
int	ModeCalc.l_mode = 2;
	ready for mode decomposition remember to add sin(theta) for GL scheme
	ModeCalc.x = data_set[n][t][3];
	ModeCalc.p4re = data_set[n][t][5]
tuple	ModeCalc.sintheta = z/(x**2+y**2+z**2)**0.5;
float	ModeCalc.sinphi = x/(x**2+y**2)**0.5
int	ModeCalc.cos2phi = 2*sinphi*cosphi
	ModeCalc.ModeRe = np.zeros(tstep)
	start real surface integral here
	ModeCalc.ModeIm = np.zeros(tstep)
tuple	ModeCalc.w = (1.0/40)*(2*np.pi)*data_set[n][t][6]*(np.pi/2.0)