net rate coverage

Author: mrbuche

statmechcrack/isometric.py

@@ -337,7 +337,7 @@ def k_net_isometric(self, v):
                 ...     v = model.v(Np/N, ensemble='isometric')
                 ...     _ = plt.semilogy(
                 ...         3*model.kappa*(v - 1)/N**2,
-                ...         model.k_net_isometric(v),
+                ...         model.k_net(v, ensemble='isometric'),
                 ...         label=r'$N=$'+str(N))
                 >>> _ = plt.xlabel(r'$3\kappa\Delta v/N^2$')
                 >>> _ = plt.ylabel(r'$k^\mathrm{net}/k_\mathrm{ref}$')

statmechcrack/isotensional.py

@@ -325,13 +325,13 @@ def k_net_isotensional(self, p):
 
                 >>> import numpy as np
                 >>> import matplotlib.pyplot as plt
-                >>> from statmechcrack import CrackIsotensional
+                >>> from statmechcrack import Crack
                 >>> Np = np.logspace(-1, np.log(6)/np.log(10), 33)
                 >>> _ = plt.figure()
                 >>> for N in [4, 8, 16, 64]:
-                ...     model = CrackIsotensional(N=N)
+                ...     model = Crack(N=N)
                 ...     _ = plt.semilogy(
-                ...         Np, model.k_net_isotensional(Np/N),
+                ...         Np, model.k_net(Np/N, ensemble='isotensional'),
                 ...         label=r'$N=$'+str(N))
                 >>> _ = plt.xlabel(r'$Np$')
                 >>> _ = plt.ylabel(r'$k^\mathrm{net}/k_\mathrm{ref}$')