def METHOD_NAME(self, theta, init_state, tol=1e-14, nshots=100000):
    """
    Args:
        theta: list or numpy.array with the angles to be used in the circuit
        init_state: numpy.array with the quantum state to be Schmidt-decomposed
        nshots: int number of runs of the circuit during the sampling process (default=10000)
    Returns:
        numpy.float64 with the value of the Von Neumann entropy for the given bipartition
    """
    Schmidt = self.Schmidt_coeff(theta, init_state, nshots=nshots)
    Schmidt = Schmidt**2
    non_zero_coeff = np.array([coeff for coeff in Schmidt if coeff > tol])
    return -np.sum(non_zero_coeff * np.log2(non_zero_coeff))