By the Monte-Carlo method, magnetic properties have been investigated for Eu-doped GaN involving internal nanostructures induced by nanoscale spinodal decomposition, where these nanostructures are spontaneously or artificially formed. In the present simulations, hysteretic and nonhysteretic magnetization curves are observed in the systems with the large-sized and small-sized nanostructures, respectively. These nanostructures affect the blocking temperatures as well. Furthermore, they influence temperature-dependent energy barriers of spin flipping; therefore, the simulations suggest that the magnetization is thermally stable. However, we observe that the blocking temperatures are smaller than the experimental values, which may be due to atomic vacancies.