new_contents: import unittest
  import numpy as np
  from pyquante2.dft.functionals import xs,cvwn
  from pyquante2.dft.reference import data
  
  class test_dft(unittest.TestCase):
      def test_xs(self):
          na = data['xs'][:,0]
          nb = data['xs'][:,1]
          fa,dfa = xs(na)
          fb,dfb = xs(nb)
          max_f = np.max(fa+fb-data['xs'][:,5])
          max_dfa = np.max(dfa-data['xs'][:,6])
          max_dfb = np.max(dfb-data['xs'][:,7])
          self.assertAlmostEqual(max_f,0,5) ## Fix this!
          self.assertAlmostEqual(max_dfa,0)
          self.assertAlmostEqual(max_dfb,0)
  
      def test_cvwn(self):
          na = data['xs'][:,0]
          nb = data['xs'][:,1]
          f,dfa,dfb = cvwn(na,nb)
          max_f = np.max(f-data['cvwn'][:,5])
          max_dfa = np.max(dfa-data['cvwn'][:,6])
          max_dfb = np.max(dfb-data['cvwn'][:,7])
          self.assertAlmostEqual(max_f,0)
          self.assertAlmostEqual(max_dfa,0)
          #self.assertAlmostEqual(max_dfb,0,6) ## Fix this!
  
  
  def runsuite(verbose=True):
      if verbose: verbosity=2
      else: verbosity=1
      suite = unittest.TestLoader().loadTestsFromTestCase(test_dft)
      unittest.TextTestRunner(verbosity=verbosity).run(suite)
      return
  
  def debugsuite():
      import cProfile,pstats
      cProfile.run('runsuite()','prof')
      prof = pstats.Stats('prof')
      prof.strip_dirs().sort_stats('time').print_stats(15)
  
  if __name__ == '__main__':
      import sys
      if "-d" in sys.argv:
          debugsuite()
      else:
          runsuite()
      
      
          
  
subject: Test suite for lda functionals