Improve interface for passing matrices and/or network reduction kwargs

[m-bossart]

Currently the kwargs for applying network reductions must match the ptdf matrix if it is passed directly.
For example, this will fail because ptdf was computed without additional network reductions which contradicts the kwargs. We should have a better interface so users don't need to make changes in both places.

ptdf = PNM.PTDF(sys)
template = PSI.ProblemTemplate(
                    PSI.NetworkModel(
                    PSI.PTDFPowerModel;
                    PTDF_matrix = ptdf,
                    use_slacks = true, 
                    reduce_radial_branches = true,     
                    reduce_degree_two_branches = true,
                                )
                            )

[jd-lara]

This is resolved with the common core for VirtualPTDF and MODF