On the computational aspects of a symmetric multidomain BEM for elastoplastic analysis