Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry