The parallel performance and efficiency of boundary element methods in the solution of inelastic problems are investigated. A parallel algorithm based on the so-called initial stress scheme for inelastic analysis is developed and implemented in BEMGEO (Boundary Element Method program for GEOmechanics applications). Numerical examples solved on an SP2 of 34 thin nodes are presented which highlight the advantage of boundary element methods in high performance computing.