Fast computing of the Moore-Penrose inverse matrix