An iterative multi-fidelity approach for model order reduction of multi-dimensional input parametric PDE systems