LiScM2 is a MATLAB package for the numerical solutions of some partial differential evolution equations u_t+Lu=N(u,\nabla u), with boundary and initial conditions, where L is a 2D linear elliptic operator (Laplace operator for this version) and N is a nonlinear part.
The algorithm is based on the eigenfunctions expansion together with the Lyapunov-Schmidt reduction procedure. The eigenfunction basis of the linear part L of the problem is used to represent the solution at every time level and these eigenfunctions are calculated in a preprocessing stage. A boundary-fitted computational grid is generated in order to allow a complex geometry of the physical domain in the 2D case. |