diff --git a/demos/getting_started_5_mf.py b/demos/getting_started_5_mf.py
index e6bf3c2a67dc688a4bc78394ccbe65aecc042922..2d6db9fa1f3689b48945ee9e673d52f2ff018f62 100644
--- a/demos/getting_started_5_mf.py
+++ b/demos/getting_started_5_mf.py
@@ -81,7 +81,7 @@ def main():
     cfmaker.set_amplitude_total_offset(0., (1e-2, 1e-6))
     correlated_field = cfmaker.finalize()
 
-    normalized_amp = cfmaker.normalized_amplitudes()
+    normalized_amp = cfmaker.get_normalized_amplitudes()
     pspec1 = normalized_amp[0]**2
     pspec2 = normalized_amp[1]**2
     DC = SingleDomain(correlated_field.target, position_space)
diff --git a/src/library/correlated_fields.py b/src/library/correlated_fields.py
index 3183525be5843799850216627370f06b451dc63e..b037f0a29862058b44189eddf76da4098cc35a2f 100644
--- a/src/library/correlated_fields.py
+++ b/src/library/correlated_fields.py
@@ -707,7 +707,7 @@ class CorrelatedFieldMaker:
             ht = HarmonicTransformOperator(ht.target,
                                            self._target_subdomains[i][amp_space],
                                            space=spaces[i]) @ ht
-        a = list(self.normalized_amplitudes())
+        a = list(self.get_normalized_amplitudes())
         for ii in range(n_amplitudes):
             co = ContractionOperator(hspace, spaces[:ii] + spaces[ii + 1:])
             pp = a[ii].target[amp_space]
@@ -758,7 +758,7 @@ class CorrelatedFieldMaker:
         """Returns the added fluctuations operators used in the model"""
         return self._a
 
-    def normalized_amplitudes(self):
+    def get_normalized_amplitudes(self):
         """Returns the normalized amplitude operators used in the final model
 
         The amplitude operators are corrected for the otherwise degenerate
@@ -798,7 +798,7 @@ class CorrelatedFieldMaker:
                  ' no unique set of amplitudes exist because only the',
                  ' relative scale is determined.')
             raise NotImplementedError(s)
-        normal_amp = self.normalized_amplitudes()[0]
+        normal_amp = self.get_normalized_amplitudes()[0]
 
         expand = ContractionOperator(
             normal_amp.target, len(normal_amp.target) - 1
diff --git a/test/test_operators/test_correlated_fields.py b/test/test_operators/test_correlated_fields.py
index ba4c1f350ab4eecd9e74734a1200692966c1b6d2..c748e79f7c072ce05b2a416e5e417d17cb306a71 100644
--- a/test/test_operators/test_correlated_fields.py
+++ b/test/test_operators/test_correlated_fields.py
@@ -142,7 +142,7 @@ def testAmplitudesInvariants(sspace, N):
     assert op.target[-2] == sspace
     assert op.target[-1] == fsspace
 
-    for ampl in fa.normalized_amplitudes():
+    for ampl in fa.get_normalized_amplitudes():
         ift.extra.check_operator(ampl, 0.1*ift.from_random(ampl.domain), ntries=10)
     ift.extra.check_operator(op, 0.1*ift.from_random(op.domain), ntries=10)