diff --git a/doc/index.rst b/doc/index.rst
index 85fbd27078b46a530c8d851d19ee5230599f7342..2040ebafe8a3da3944ba30779a1e13b3ee598fb5 100644
--- a/doc/index.rst
+++ b/doc/index.rst
@@ -24,7 +24,7 @@ TetraX Documentation
 
 .. automodule:: tetrax
 
-**Version**: 1.3.1
+**Version**: 1.3.3
 
 TetraX is a package for finite-element-method (FEM) micromagnetic modeling with the aim to provide user friendly and versatile micromagnetic workflows. Apart from energy minimizers and an LLG solver, it aims to provide implementations of several FEM dynamic-matrix approaches to numerically calculate the normal modes and associated frequencies for magnetic specimen of different geometries such as confined samples, infinitely long waveguides, or infinitely extended multilayers. Apart from ferromagnets, the package also supports antiferromagnets as an experimental feature.
 
diff --git a/doc/release_notes.rst b/doc/release_notes.rst
index 6a95a4750cdc4ff9faa0e4e29f0cc4ac1c86f4c2..142fad93d87e7b7fb618d63f145a3854c28fa489 100644
--- a/doc/release_notes.rst
+++ b/doc/release_notes.rst
@@ -1,6 +1,15 @@
 Release Notes
 =============
 
+Version 1.3.3
+-------------
+`2023-12-20`
+
+Hotfixes
+^^^^^^^^
+- fixed wrong scaling factor in iDMI tensor (We thank Joo-Von Kim for highlighting this issue)
+
+
 Version 1.3.2
 -------------
 `2023-10-17`
diff --git a/setup.py b/setup.py
index 99af59e92e3678bee8a77e6f25591fea40357989..ff17b3bd7fb84840c7d21466057f6cab40037e46 100644
--- a/setup.py
+++ b/setup.py
@@ -37,7 +37,7 @@ extensions = [
 
 setup(
     name='TetraX',
-    version='1.3.2',
+    version='1.3.3',
     entry_points={'console_scripts': ['TetraX=fem_core:main']},
     packages=find_packages(),
     ext_modules=cythonize(extensions, language_level=3),
diff --git a/tetrax/core/operators/dmi_interfacial.py b/tetrax/core/operators/dmi_interfacial.py
index 064a3973311994d2011e63e25271a21ef7c2b263..8f89b4ff2b2c27cee1f371ee52e9f979a4cf161a 100644
--- a/tetrax/core/operators/dmi_interfacial.py
+++ b/tetrax/core/operators/dmi_interfacial.py
@@ -38,21 +38,17 @@ def _construct_diff_and_pi_matrices(sample):
 class InterfacialDMIOperator(LinearOperator):
 
     def __init__(self, sample):
-#        self.fac = 2 * sample.Didmi * (sample.Msat ** 2 * mu_0 * sample.scale)
-#        self.fac = 2 * sample.Didmi * (mu_0 * sample.scale)
-#        self.Msat = sample.Msat
-#        self.k = 0
+        self.fac = 2 * sample.Didmi / (mu_0 * sample.scale)
+        self.Msat = sample.Msat
         self.diff_mat, self.Pi_mat = _construct_diff_and_pi_matrices(sample)
-#        self.sparse_mat = self.fac * (self.diff_mat + self.k * self.Pi_mat)
-#        self.shape = self.sparse_mat.shape
         self.update(sample)
 
     def update(self,sample):
-        self.fac = 2 * sample.Didmi * (mu_0 * sample.scale)
+        self.fac = 2 * sample.Didmi / (mu_0 * sample.scale)
         self.Msat = sample.Msat
-#        self.fac = 2 * sample.Didmi * (sample.Msat ** 2 * mu_0 * sample.scale)
         self.k = 0
-        self.sparse_mat = self.fac * (self.diff_mat + self.k * self.Pi_mat)
+        self.set_k(self.k)
+#        self.sparse_mat = self.fac * (self.diff_mat + self.k * self.Pi_mat)
         self.shape = self.sparse_mat.shape
 
     def set_k(self, k):