Computer simulation of ultrahigh-density perpendicular magnetic recording