Heteroscedasticity testing has a long history and is still an important matter in the linear model. There exist many types of tests, but they are limited in use to their own specific cases and sensitive to normality. Here, we propose a dimension test approach to heteroscedasticity. The proposed test overcomes the shortcomings of the existing methods, so that it is robust to normality and is unified in sense that it is applicable in the linear model with multi-dimensional response. Numerical studies confirm that the proposed test is favorable over the existing tests with moderate sample sizes, and real data analysis is presented.