docs: def multiplyC(self, alpha):
          """multiply C with a scalar and update all related internal variables (dC, D,...)"""
          self.C *= alpha
          if self.dC is not self.C:
              self.dC *= alpha
          self.D *= alpha**0.5