Stiffness matrix of manipulators with passive joints: computational aspects