Dans le domaine des hyperfréquences, le comportement magnétique des matériaux ferri ou ferro magnétiques est modélisé par le couplage des équations d’évolution de Landau et Lifchitz (micromagnétisme) et des équations de la magnétostatique. Le terme démagnétisant (magnétostatique) est non local, ce qui rend la simulation dynamique coûteuse pour des systèmes de taille suffisamment importante pour pouvoir observer les phénomènes clef (déplacement de parois, configuration en domaines). Afin de résoudre ce problème, on propose une méthode d’intégration des équations de la magnétostatique s’appuyant sur l’utilisation de la forme particulière du problème discret dans le cas d’un maillage cubique régulier: les matrices Toeplitz multi-niveaux.