  It is understood that strong magnetic fields affect the structure of neutron
stars. Nevertheless, many calculations for magnetized neutron stars are still
being performed using symmetric solutions of Einstein's equations. In this
conference proceeding, we review why this is not the correct procedure and we
also discuss the effects of magnetic fields on the stellar population and
temperature profiles.
