def METHOD_NAME(self, value):
    if self._validate_args:
        self._validate_sample(value)
    # Eq. 2.1 in [1]
    skew_prob = torch.log1p(
        (self.skewness * torch.sin((value - self.base_dist.mean) % (2 * pi))).sum(
            -1
        )
    )
    return self.base_dist.METHOD_NAME(value) + skew_prob