Numerical prediction of crack growth is an important problem in computational mechanics. The difficulty in this problem arises from the basic incompatibility of cracks with the partial differential equations that are used in the classical theory of solid mechanics. The spatial derivatives needed for these partial differential equations to make sense do not exist on a crack tip or surface. Any numerical method derived from these equations inherits this difficulty in modeling cracks. In spite of the progress that has been made in developing meshfree methods aimed at modelling fracture such meshfree techniques generally require some special method of evaluating the spatial derivatives on each crack surface.