docs: def cpower(a,b):
      '''
      cpower(a,b) is equivalent to a**b except that it also operates over sparse arrays; though it
      must reify them to do so.
      '''
      if sps.issparse(a): a = a.toarray()
      if sps.issparse(b): b = b.toarray()
      return a ** b