docs: def _getScalesRand(self):
          """
          Return a vector of random scales
          """
          if self.P>1:
              scales = []
              for term_i in range(self.n_terms):
                  _scales = SP.randn(self.diag[term_i].shape[0])
                  if self.offset[term_i]>0:
                      _scales = SP.concatenate((_scales,SP.array([SP.sqrt(self.offset[term_i])])))
                  scales.append(_scales)
              scales = SP.concatenate(scales)
          else:
              scales=SP.randn(self.vd.getNumberScales())
  
          return scales