Numerical methods such as boundary element and finite element methods are widely used for the stress analysis in solid mechanics. This study presents boundary element method based on the displacement discontinuity formulation to solve general problems of interaction between hydraulic fracturing and discontinuities. The crack tip element and a higher order boundary displacement collocation technique are used to study the hydraulic fracture propagation and its interaction with the pre-existing cracks and discontinuities in an elastic rock mass. The maximum tangential stress criterion (or -criterion) and the strain energy density criterion (SED) are used to obtain the fracture path and the results of both criteria are compared with each other. The comparison of numerical method with the results brought in the literature shows a good performance of the method in the case of interacting cracks.