1 if { [ catch {exec which qchem } ] } then {

2   set msg "Unsupported test: Requires qchem set up in path)"

3   puts stdout $msg

4   return $msg

5 }

6

7 #

8 # Molpro example using ChemShell basis library

9 # Molpro uses harmonic basis by default, switch off for comparison

10 #

11 z_create zmatrix=z {

12 zmatrix angstrom

13 o

14 h 1 1

15 h 1 1 2 100

16 }

17 z_to_c zmatrix=z coords=c

18 matrix e

19 energy coords=c theory=qchem : { harmonic=no basisspec = { {sto-3g *} {dzp o} } } energy=e

20 set final_energy  [ get_matrix_element matrix=e indices= {0 0 } ]

21 delete_object e

22

23 # Test result

24 set etest -76.01244105

25

26 if { abs( $final_energy - $etest) > 1.0e-10 } {

27    chemerr "Failed .. E= $final_energy, should be $etest"

28 } else {

29    return "OK"

30 }