def METHOD_NAME(self):
    # Bin edges
    be_mea = get_uneven_bins_edges([-3, -2, 1, 2], [80, 80, 10])
    be_sig = np.arange(0.50, 0.70, 0.01)
    # Set params
    imls1 = self.oqp.imtls[self.imts[0]]
    imls2 = self.oqp.imtls[self.imts[1]]
    len1 = len(imls1)-1
    # Compute the MRD: indirect
    imt1, imt2 = self.imts
    mrdi = calc_mean_rate_dist(self.ctx, 1, self.cmaker, self.crosscorr,
                               imt1, imt2, be_mea, be_sig)
    # Compute the MRD: direct
    mrdd = np.zeros((len1, len1, len(self.cmaker.gsims)))
    update_mrd(self.ctx, self.cmaker, self.crosscorr, mrdd)
    np.testing.assert_almost_equal(mrdi[:, :, 0], mrdd)
    if PLOT:
        imlc1 = np.diff(imls1) / 2 + imls1[:-1]
        imlc2 = np.diff(imls2) / 2 + imls2[:-1]
        fig, axs = plt.subplots(1, 1)
        fig.set_size_inches(9, 6)
        plt1 = plt.contourf(np.log(imlc1), np.log(imlc2), mrdd[:, :, 0])
        _ = plt.contour(np.log(imlc1), np.log(imlc2), mrdi[:, :, 0, 0],
                        colors='orange', linestyles='dashed')
        _ = plt.colorbar(mappable=plt1)
        _ = plt.title('MRD')
        axs.set_xlabel(f'{self.imts[0]}')
        axs.set_ylabel(f'{self.imts[1]}')
        plt.savefig('./figs/mrd_test_mrd_comparison.png')
        plt.show()