The Schrödinger equation is at heart of Bose Einstein Condensates with the celebrated Gross Pitaevskii equation. Such dispersive partial differential equations have many preserved quantities such as the mass, the energy or the momentum. It is therefore crucial to build numerical schemes that preserve these invariants. We will present in this talk various way to take this property into account and the pros and cons of the various methods.