code: def compute_dual_rmetric(self,Ynew=None):
          
          usedY = self.Y if Ynew is None else Ynew
          rieman_metric = RiemannMetric(usedY, self.laplacian_matrix)
          return rieman_metric.get_dual_rmetric()
docs: Helper function to calculate the