ode_solve.py File Reference

Classes
class	ode_solve.Grid1D
class	ode_solve.FEM_coefficient_matrix_generator
class	ode_solve.FEM_load_vector_generator
class	ode_solve.ode_solve

Namespaces
	ode_solve

Functions
def	ode_solve.solve (A, rhs)
def	ode_solve.load (f)
def	ode_solve.assemble (d, k, ikt1, ikt2, M1, M2, D1, D2, D0)
def	ode_solve.mult (A, x)
def	ode_solve.mult_transpose (A, x)
def	ode_solve.quadrature_points_local (order)
def	ode_solve.quadrature_weights_local (order)

Variables
int	ode_solve.L = lambda z : z + 1
int	ode_solve.dLdz = lambda z : z**0
int	ode_solve.t = 0
int	ode_solve.dof = 2**10
float	ode_solve.domain_height = 0.9
tuple	ode_solve.grid = (np.linspace(0,domain_height,dof)/domain_height)**2 *domain_height
int	ode_solve.d = 1
int	ode_solve.k1 = 1
int	ode_solve.k2 = 1
int	ode_solve.n = 2
	manufactured solution and derivatives (u=0 at z=0 dudz=0 at z=domain_height) More...
tuple	ode_solve.const = (n + 1/2)*np.pi/domain_height
	ode_solve.psi = lambda z : np.sin(const*z)
tuple	ode_solve.dpsidz = lambda z : const*np.cos(const*z)
tuple	ode_solve.d2psidz2 = lambda z : -const**2*np.sin(const*z)
	ode_solve.f
	ode_solve.kappa = np.vectorize(lambda z: L(z)**2)
	ode_solve.ode_solve_inst = ode_solve(dof, kappa, domain_height=domain_height, grid=grid)
	ode_solve.psi_approx = ode_solve_inst(d, f, k1, k2, t=t)
	ode_solve.z = np.linspace(0,domain_height,len(psi_approx))