This paper, as a continuation of our previous investigation [Nucl. Phys. B 1005 (2024) 116582] aims to study the glassy random matrices with quenched Wigner disorder. In this previous work, we have constructed a renormalization group based on the effective deterministic kinetic spectrum emerging from large N limit, and we extended approximate solutions using standard vertex expansion, at the leading order of the derivative expansion. Now in the following work, by introducing the non-trivial Ward identities which come from the (U(N))×2 symmetry broken of the effective kinetic action, we provide in the deep IR the explicit solution of the functional renormalization group for a model with quartic coupling by solving the Hierarchy to all orders in the local sector, which in particular imply the vanishing of the anomalous dimension. The numerical investigations confirm the first-order phase transition discovered in the vertex expansion framework, both in the active and passive schemes. Finally, we extend the discussion to hermitian matrices.