old_contents: from BigDFT.Systems import System
  from BigDFT.Fragments import Fragment
  from BigDFT.Atoms import Atom
  
  at = Atom({"He": [0, 0, 0]})
  frag = Fragment([at])
  sys = System({"FRA:0": frag})
  
  from BigDFT.Calculators import SystemCalculator
  code = SystemCalculator()
  
  from BigDFT.Inputfiles import Inputfile
  inp = Inputfile()
  inp.set_xc("PBE")
  inp.set_hgrid(0.4)
  
  log = code.run(sys=sys, input=inp)
  print(log.energy)
  
  
new_contents: # Test that we can import everything
  import futile
  import gi
  import BigDFT
  
  # Test A Full Calculation with PyBigDFT
  from BigDFT.Systems import System
  from BigDFT.Fragments import Fragment
  from BigDFT.Atoms import Atom
  
  at = Atom({"He": [0, 0, 0]})
  frag = Fragment([at])
  sys = System({"FRA:0": frag})
  
  from BigDFT.Calculators import SystemCalculator
  code = SystemCalculator()
  
  from BigDFT.Inputfiles import Inputfile
  inp = Inputfile()
  inp.set_xc("PBE")
  inp.set_hgrid(0.4)
  
  log = code.run(sys=sys, input=inp)
  print(log.energy)
  
  # Test a calculation using GIBinding
  from BigDFT.Calculators import GIBinding
  code = GIBinding()
  inp["posinp"] = sys.get_posinp()
  code.set(inp)
  out = code.run()
  print(out.eKS)
  
  
  
subject: Test more imports