code: def vinet_v(p, v0, k0, k0p, min_strain=0.01):
      
      if isuncertainties([p, v0, k0, k0p]):
          f_u = np.vectorize(uct.wrap(vinet_v_single), excluded=[1, 2, 3, 4])
          return f_u(p, v0, k0, k0p, min_strain=min_strain)
      else:
          f_v = np.vectorize(vinet_v_single, excluded=[1, 2, 3, 4])
          return f_v(p, v0, k0, k0p, min_strain=min_strain)
docs: find volume at given pressure
  
      :param p: pressure in GPa
      :param v0: unit-cell volume in A^3 at 1 bar
      :param k0: bulk modulus at reference conditions
      :param k0p: pressure derivative of bulk modulus at reference conditions
      :param min_strain: defining minimum v/v0 value to search volume for
      :return: unit cell volume at high pressure in A^3
      :note: wrapper function vetorizing vinet_v_single