Vectorization of astrophysical code for massively parallel supercomputers