qq = self._compute_qq(p, qdq, R)