  A flexible spectral code for the study of general relativistic
magnetohydrodynamics is presented. Aiming at investigating the physics of
slowly rotating magnetized compact stars, this new code makes use of various
physically motivated approximations. Among them, the relativistic anelastic
approximation is a key ingredient of the current version of the code. In this
article, we mainly outline the method, putting emphasis on algorithmic
techniques that enable to benefit as much as possible of the non-dissipative
character of spectral methods, showing also a potential astrophysical
application and providing a few illustrative tests.
