Michael: For beta.F and handle_gradient simply change the line 779 and 125 to  if (if_scf)  and if (do_summary_print) respectively. <br><br>I am looking at the run time error.<br><br>Ajith<br><br><div class="gmail_quote">On Wed, Jul 18, 2012 at 12:13 PM, Michael Banck <span dir="ltr"><<a href="mailto:mbanck@debian.org" target="_blank">mbanck@debian.org</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Ajith,<br>
<div><br>
On Wed, Jul 18, 2012 at 10:37:48AM -0400, Ajith perera wrote:<br>
> I have the geop.tar.gz in the following location<br>
> <a href="https://dl.dropbox.com/u/60183859/geopt.tar.gz" target="_blank">https://dl.dropbox.com/u/60183859/geopt.tar.gz</a>. Please unzip/untar it in<br>
> ACESIII directory.<br>
<br>
</div>Thanks for that!<br>
<br>
I now encountered another two errors:<br>
<br>
beta.F:779.15:<br>
<br>
           if (if_scf .eq. .true.) then<br>
               1<br>
Error: Logicals at (1) must be compared with .eqv. instead of .eq.<br>
<br>
and a similar one in handle_gradient.f around line 125.  I am attaching<br>
my current path and our GNUmakefile.chssi for reference as well.<br>
<br>
Also, I now got a better backtrace for the runtime issue:<br>
<br>
(gdb) bt<br>
#0  dlamc3_ (a=1, b=-0) at dlamc3.f:52<br>
#1  0x0837072e in dlamc1_ (beta=0, t=0, rnd=.FALSE., ieee1=.FALSE.) at<br>
dlamc1.f:186<br>
#2  0x08368c4f in dlamc2_ (beta=1146246233, t=0, rnd=<a href="tel:3218793352" value="+13218793352" target="_blank">3218793352</a>, eps=0,<br>
emin=0, rmin=0, emax=0, rmax=0) at dlamc2.f:124<br>
#3  0x08360566 in dlamch_ (cmach=..., _cmach=12) at dlamch.f:90<br>
#4  0x0835ec3c in dsyev_ (jobz=..., uplo=..., narg=3, a=..., ldaarg=3,<br>
w=..., work=..., lworkarg=9, info=0, _jobz=1,<br>
    _uplo=1) at dsyev.F:177<br>
#5  0x082d0a09 in eig_ (a=..., b=..., junk=3, n=3, sort=0) at eig.F:79<br>
#6  0x082e3599 in symmetry_ (scratch=..., qtmp=..., newq=...,<br>
nosilent=.TRUE.) at symmetry.F:120<br>
#7  0x082d7129 in geopt_ (z=..., memreq=12180600) at geopt.F:516<br>
#8  0x082d2df7 in aces2_joda_main_ () at aces2_joda_main.F:107<br>
#9  0x08067efd in scf_init_ () at scf_init.F:43<br>
#10 0x08062c67 in aces3 () at beta.F:160<br>
#11 0x0805e1a4 in main (argc=1, argv=0xbfdb3928 './xaces3\000') at<br>
beta.F:1055<br>
#12 0x40d2de16 in ?? ()<br>
#13 0x0805e1c5 in _start ()<br>
<br>
As far as I can tell, DLAMCH are trying to figure out some machine<br>
constants (see comment in dup/dsyev.F) and running in some tight loops.<br>
Possibly it is some issue with the blas/lapack interface.<br>
<br>
<br>
Best regards,<br>
<br>
Michael<br>
</blockquote></div><br>