code: def get_covaried_params(mus, evecsCV):
      
      mus = numpy.array(mus, copy=False)
      # If original inputs were floats we need to make this a 2D array
      if len(mus.shape) == 1:
          resize_needed = True
          mus = mus[:,None]
      else:
          resize_needed = False
  
      xis = ((mus.T).dot(evecsCV)).T
  
      if resize_needed:
          xis = numpy.ndarray.flatten(xis)
  
      return xis
docs: Function to rotate from position(s) in the mu_i coordinate system into the
      position(s) in the xi_i coordinate system
  
      Parameters
      -----------
      mus : list of floats or numpy.arrays
          Position of the system(s) in the mu coordinate system
      evecsCV : numpy.matrix
          This matrix is used to perform the rotation to the xi_i
          coordinate system.
  
      Returns
      --------
      xis : list of floats or numpy.arrays
          Position of the system(s) in the xi coordinate system