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