#endif
double precision g(maxvar),f1
integer nvarx
+ double precision energia(0:n_ene)
#ifdef LBFGS
maxiter=maxmin
coordtype='CARTESIAN'
enddo
nvarx=k
write (iout,*) "Variables set up nvarx",nvarx
+ write (iout,*) "Before energy minimization"
+ call etotal(energia(0))
+ call enerprint(energia(0))
#ifdef LBFGS
c
c From tinker