In this paper the problem of computing bifurcation diagrams for large-scale nonlinear parameter-dependent steady state systems which arise following the spatial discretization of semilinear PDEs is ...