The proposed methodology, combining the Newton-Raphson method based on the consistent nonlocal stiffness matrix with direct and indirect solvers, has been successfully implemented into an existing object-oriented finite element environment [9]. The numerical performance is evaluated using three examples. The computations are performed on Pentium III workstations with 1GB of memory, running at 800 MHz (for 2D tests) and 933 MHz (for the 3D test) under the Linux operating system. The egcs-2.91.66 C++ compiler is used. All the reported times are measured as user times. The nonlocal isotropic damage model with a Rankine-type equivalent strain measure is used in all the examples.