docs: def getVariances(self):
          """
          get variances
          """
          var = []
          var.append(self.Cr.K().diagonal())
          if self.bgRE:
              var.append(self.Cg.K().diagonal())
          var.append(self.Cn.K().diagonal())
          var = sp.array(var)
          return var