I get rounded edges also with RegionFunction:
RegionFunction
ListDensityPlot[S, ColorFunction -> "SolarColors", DataRange -> rstar {{-1, 1}, {-1, 1}}, RegionFunction -> Function[{x, y, z}, Norm[{x, y}] < .95 rstar]]