Lazily evaluated arrays in Python
Modes and nodes for high-order discretizations
Choose your favorite numpy-workalike!