This package takes in a 2D nodal grid and produces the laplacian with Dirichlet's, Neumann's, or Robin boundary condition (BC). You can mix these BCs per cells on the border. This means you can have one type of BC on one edge and another type on another edge, or even on a single edge one cell can have one type and other cells can have other types of BC. The grid must be a structural grid but it can be general curvilinear.
You can generate the operator using 2nd, 4th, or 6th order Castillo-Grone's Mimetic Difference Operators. For interpolation you can use any order of polynomial that you want.
Once the operator is generated, you can use any solver of your choice to solve the equations. As long as your grid is not changed the operator stays the same and you can solve multiple problems.
Start with the one of the Test*.m functions to familiarize yourself with the code.

If you use this package or a modified version of it, please cite the followings:
- M. Abouali, "Investigating Castillo-Grone's Mimetic Difference Operators in Development of Geophysical Fluid Dynamics Models Implemented on GPGPUs", Ph.D. Thesis, San Diego State University, May 2014.
- Abouali M., Castillo J.E., Solving Poisson equation with Robin boundary condition on a curvilinear mesh using high order mimetic discretization methods, Mathematics and Computers in Simulation, (2014), http://dx.doi.org/10.1016/j.matcom.2014.10.004

NOTE: To run this package you need to have "Interpolant/Extrapolant 2D/3D" downloaded and within your MATLAB path. This package can be found on http://www.mathworks.com/matlabcentral/fileexchange/41669-interpolantextrapolant-2d3d-data