NeoPZ is an object oriented programming environment which implements hp adaptive finite element algorithms. The PZ classes can be used standalone but can also be linked with blas, boost, pthread, log4cxx and metis. The advanced finite element technologies implemented are : h and hp adaptivity for one, two and three dimensional meshes the variational formulation of the conservation law is separated from the generation of the approximation space. Therefore NeoPZ can be used for virtually any system of differential equation Various storage patterns for the global stiffness matrix : full, banded, skyline, sparse, frontal Iterative solvers with configurable preconditioners Post processing interface to opendx and vtk