HvD: More docs.

src/selci/anal.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_anal(irange, roots, ept, ptnorm, enew, ncinew,
      &           ptnorm_mp, ept_mp)
 *
@@ -68,3 +70,4 @@ c
      $ 20(:f15.8,2x,f15.8,2x,f15.8/))
 c
       end
+C> @}

src/selci/axb.f

@@ -1,3 +1,6 @@
+C> \ingroup selci
+C> @{
+C>
       subroutine selci_axb(a,mrowa,b,mrowb,c,mrowc,ncol,nlink,nrow)
 *
 * $Id$
@@ -87,3 +90,5 @@ c
       enddo
 c
       end
+C>
+C> @}

src/selci/fcidim.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       double precision function selci_fcidim(nelec, multi, issss,
      $     norbs, isym)
 *
@@ -104,3 +106,4 @@ c
       selci_fcidim = r1(0,0,0)
 c
       end
+C> @}

src/selci/getij.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_getij(i,j,ioconf,joconf,nintpo,nbitpi,iocc)
 *
 * $Id$
@@ -48,3 +50,4 @@ c
       endif
 c
       end
+C> @}

src/selci/gettim.f

@@ -1,3 +1,6 @@
+C> \ingroup selci
+C> @{
+C>
       subroutine selci_gettim(cpud,elapsd)
 *
 * $Id$
@@ -18,3 +21,5 @@
       timediff = util_cpusec() - elast
       cpud = dble(timediff)
       end
+C>
+C> @}

src/selci/gtijkl.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_gtijkl(ind,ioconf,joconf,nintpo,nbitpi,
      $     iocc,icase)
 *
@@ -76,3 +78,4 @@ c        endif
       endif
 c
       end
+C> @}

src/selci/idmin.f

@@ -1,3 +1,6 @@
+C> \ingroup selci
+C> @{
+C>
       integer function selci_idmin(n,a,ia)
 *
 * $Id$
@@ -19,3 +22,5 @@ c
 c
       selci_idmin = ind
       end
+C>
+C> @}

src/selci/inicij.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_inicij(w1,w2)
 *
 * $Id$
@@ -82,3 +84,4 @@ c
  50   continue
 c
       end
+C> @}

src/selci/initci.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_initci(title, type)
 *
 * $Id$
@@ -82,3 +84,4 @@ c
      &       ' no. of orbitals per sym   ',8i5/)
  5    format(' symmetry of the orbitals'//10(26i3/:))
       end
+C> @}

src/selci/input.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine setinp(ir,iw)
 *
 * $Id$
@@ -363,3 +365,4 @@ c     in jbuf and transferred to junke
  2    locatc=i
       return
       end
+C> @}

src/selci/iprtri.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_iprtri(id,n)
 *
 * $Id$
@@ -31,3 +33,4 @@ c
  8028 format(6x,8(2x,i5,2x))
  8048 format(i5,1x,8i9)
       end
+C> @}

src/selci/isum.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       integer function selci_isum(n,m,im)
 *
 * $Id$
@@ -14,3 +16,4 @@ c
       selci_isum = is
 c
       end
+C> @}

src/selci/load.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_load(h, g, int12, int34)
 *
 * $Id$
@@ -64,3 +66,4 @@ c 30   continue
 c
 c      write(6,*) ' ngot in load ',ngot
       end
+C> @}

src/selci/loada.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_loada(h, g, int12, int34)
 *
 * $Id$
@@ -70,3 +72,4 @@ c
       close(iflint,status='keep')
 c
       end
+C> @}

src/selci/loadab.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_loadab(h, g, int12, int34)
 *
 * $Id$
@@ -68,3 +70,4 @@ c
       close(iflint,status='keep')
 c
       end
+C> @}

src/selci/loadg.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_loadg(g,int12,int34)
 *
 * $Id$
@@ -77,3 +79,4 @@ c
 c      write(nlist,'(1x,1p5e15.7)')sum
       call util_flush(6)
       end
+C> @}

src/selci/loadw.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_loadw(iflwmt, w1, w2, nfmax, nfmax2, nsmax)
 *
 * $Id$
@@ -25,3 +27,4 @@ c         call selci_output(w2(1,1,i),1,nfmax2,1,nfmax,nfmax2,nfmax,1)
 c 20   continue
 c
       end
+C> @}

src/selci/mkdele.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mkdele(numf, nroot, deltae, vc, roots,
      $     diag, nfmax, dnorm)
 *
@@ -23,3 +25,4 @@ c
  10   continue
 c
       end
+C> @}

src/selci/mkindx.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mkindx(ioconf,indxci,ncinew,nnsmax)
 *
 * $Id$
@@ -45,3 +47,4 @@ c$doit asis
       indxci(noconf+1) = ncinew
 c
       end
+C> @}

src/selci/mkinta.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mkinta(int12, int34)
 *
 * $Id$
@@ -37,3 +39,4 @@ c
       numint = ipt + int34(kl)
 c
       end
+C> @}

src/selci/mknsneed.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mknsneed(ioconf,noconf,nintpo,nbitpi,
      $     norbs,nsneed)
 *
@@ -51,4 +53,4 @@ c
       enddo
 c
       end
-
+C> @}

src/selci/mkvc.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mkvc(q, h,g,int12,int34,w1,w2,ioconf,indxci,
      $     iexv,iconf,ncold,ci,vc,iocc,ipack,numf,work1)
 *
@@ -88,3 +90,4 @@ c
       junk = selci_mfree(iwork2)
 c
       end
+C> @}

src/selci/mpdele.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mpdele(numf, nroot, deltae_mp, vc,
      $     roots_mp, nfmax, dnorm_mp,iocc)
 *
@@ -28,3 +30,4 @@ c:old            deltae_mp(ir) = deltae_mp(ir) - vc(i,ir)*coeff
  10   continue
 c
       end
+C> @}

src/selci/mxma.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mxma(a,iac,iar, b,ibc,ibr, c,icc,icr,
      $     nar,nac,nbc)
 *
@@ -52,3 +54,4 @@ c
 c
       return
       end
+C> @}

src/selci/mxmb.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_mxmb(a,mcola,mrowa,b,mcolb,mrowb,
      1     r,mcolr,mrowr, ncol,nlink,nrow)
 *
@@ -29,3 +31,4 @@ c
  50   continue
 c
       end
+C> @}

src/selci/oeq.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       logical function selci_oeq(noconf,ioconf,joconf,nintpo)
 *
 * $Id$
@@ -46,3 +48,4 @@ c
       selci_oeq = .false.
 c     
       end
+C> @}

src/selci/ointer.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       logical function selci_ointer(noconf,ioconf,joconf,nintpo)
 *
 * $Id$
@@ -95,3 +97,4 @@ c
       selci_ointer = .false.
 c     
       end
+C> @}

src/selci/orbanal.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_orbanal(noconf, ioconf, nintpo, nbitpi, norbs,
      $     map)
       implicit none
@@ -114,3 +116,4 @@ c
       enddo
 c
       end
+C> @}

src/selci/pconls.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_pconls(ioconf)
 *
 * $Id$
@@ -13,3 +15,4 @@ c
  10   continue
 c
       end
+C> @}

src/selci/popcnt32.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       integer function selci_popcnt32(i)
 *
 * $Id$
@@ -391,3 +393,4 @@ c
 *c
 *      end
 *mdc*endif
+C> @}

src/selci/prtcon.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_prtcon(ifllog, norbs, ioconf, nintpo, nbitpi)
 *
 * $Id$
@@ -8,3 +10,4 @@ c
       call selci_wrtcon(ifllog, iocc, norbs)
 c
       end
+C> @}

src/selci/prtri.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_prtri(d,n)
 *
 * $Id$
@@ -31,3 +33,4 @@ c
  8028 format(6x,6(4x,i4,4x))
  8048 format(i5,1x,6f12.6)
       end
+C> @}

src/selci/ptidy.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_ptidy(ioconf, nintpo, noconf, ncold, nconmx,
      $     nroot, ept, enew, irange, ptnorm, ptnorm_mp, ept_mp)
 *
@@ -98,3 +100,4 @@ c        call util_flush(6)
       endif
 c
       end
+C> @}

src/selci/putinb.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_putinb(iflham,ibase,numf,hbuf,ninhb,jbuf,ninjb)
 *
 * $Id$
@@ -75,3 +77,4 @@ c
       ninjb = 0
 c
       end
+C> @}

src/selci/rd1.F

@@ -1,9 +1,11 @@
 *
 * $Id$
 *
-c
-c     routines to read columbus integrals
-c
+C> \ingroup selci
+C> @{
+C>
+C> \brief Routines to read columbus integrals
+C>
       subroutine selci_rd1(iflint,h,norbs)
       implicit real*8 (a-h,o-z)
       parameter(lenbuf=2047, n1eabf=1636, n2eabf=1364)
@@ -282,3 +284,4 @@ c        ...big-endian.
 *mdc*endif
       return
       end
+C> @}

src/selci/rdconf.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_rdconf(iflcon,ioconf,indxci,nintpo,noconf)
 *
 * $Id$
@@ -10,3 +12,4 @@ c
       read (iflcon) indxci
 c
       end
+C> @}

src/selci/rdhcon.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_rdhcon(iflcon, title, multi, nelec, issss, norbs,
      &     nnsmax, nci, noconf, nintpo, nbitpi, nbpsy, isym, nsym,
      &     inttyp, nsneed)
@@ -14,3 +16,4 @@ c
      &     nsneed
 c
       end
+C> @}

src/selci/rdhead.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_rdhead(title,type)
 *
 * $Id$
@@ -56,3 +58,4 @@ c
       endif
 c
       end
+C> @}

src/selci/rdhint.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_rdhint(iflint,ntitle,dtitle,repnuc,nnorbs,
      $     inttyp, moints)
 *
@@ -79,3 +81,4 @@ c
       endif
 c
       end
+C> @}

src/selci/rdhwmt.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_rdhwmt(iflwmt,mmulti, nsmax, nf, nfmax, nfmax2)
 *
 * $Id$
@@ -11,3 +13,4 @@ c
       nfmax2 = nf(nsmax-2)
 c
       end
+C> @}

src/selci/readh.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_readh(iflham,hh,nnci)
 *
 * $Id$
@@ -8,3 +10,4 @@ c
       read (iflham) hh
 c
       end
+C> @}

src/selci/restin.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_restin(iconf,ept,ept_mp,enew,ptnorm,ptnorm_mp,
      &                  nroot,irange,nrange,iwpt,
      &                  noconf,ioconf,nintpo,nconmx,ncold)
@@ -84,3 +86,4 @@ c
       close(unit=irest,status='keep')
       goto 90000
       end
+C> @}

src/selci/restk.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine restk
 *
 * $Id$
@@ -39,3 +41,4 @@ c
 c
       return
       end
+C> @}

src/selci/slect.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_slect(q,h,g,int12,int34,w1,w2,ioconf,indxci,
      $     roots,ci,nconmx,ncold,thresh,ept,enew,irange, ptnorm,
      &     ptnorm_mp,ept_mp,roots_mp,min2c)
@@ -154,3 +156,4 @@ c
            junk = selci_mfree(ivc)
 c     
            end
+C> @}

src/selci/stool.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_stool(iconf,ept,ept_mp,enew,ptnorm,ptnorm_mp,
      &     nroot,irange,nrange,iwpt,
      &     noconfi,ioconf,nintpo,nconmx,ncold)
@@ -164,3 +166,4 @@ c
 90912 call errquit('stool: error opening restart file',0, DISK_ERR)
 c     
       end
+C> @}

src/selci/tester.F

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_tester(q,h, g, int12, int34, w1, w2, ioconf,
      $     indxci,
      $     roots, ci, nconmx, thresh, ept, enew, iconf, ncold,
@@ -125,3 +127,4 @@ c
       junk = selci_mfree(iwork1)
 c
       end
+C> @}

src/selci/tqlrat.f

@@ -1,10 +1,51 @@
+C> \ingroup selci
+C> @{
+C>
+C> \brief Find the eigenvalues of a symmetric tridiagonal matrix
+C>
+C>
+C> This subroutine is a translation of the algol procedure
+C> tqlrat [1].
+C>
+C> This subroutine finds the eigenvalues of a symmetric
+C> tridiagonal matrix by the rational ql method.
+C>
+C> Calls pythag for  dsqrt(a*a + b*b) .
+C>
+C> Questions and comments should be directed to Burton S. Garbow,
+C> mathematics and computer science div, argonne national laboratory
+C>
+C> This version dated august 1983.
+C>
+C>
+C> ### References ###
+C>
+C> [1] C.H. Reinsch, "Algorithm 464: eigenvalues of a real, symmetric,
+C>     tridiagonal matrix [F2]", Comm. ACM. <b>16</b> (1973) 689,
+C>     doi: <a href="http://dx.doi.org/10.1145/355611.362539">
+C>     10.1145/355611.362539</a>.
+C>
       subroutine selci_tqlrat(n,d,e2,ierr)
 *
 * $Id$
 *
 c
-      integer i,j,l,m,n,ii,l1,mml,ierr
-      double precision d(n),e2(n)
+      integer n !< [Input] The order of the matrix
+      double precision d(n) !< [In/Output] The diagonal of the matrix
+                            !< - On input: the diagonal of the 
+                            !<   tridiagonal matrix
+                            !< - On output: the eigenvalues of the 
+                            !<   matrix
+      double precision e2(n) !< [In/Output] The off diagonal elements
+                             !< - On input: the squares of the 
+                             !<   subdiagonal elements of the input 
+                             !<   matrix in the positions 2:n,
+                             !<   element e2(1) is arbitrary
+                             !< - On output: destroyed
+      integer ierr !< [Output] Error flag
+                   !< - 0: routine completed successfully
+                   !< - j: Eigenvalue j and further are not converged
+      integer i,j,l,m,ii,l1,mml
       double precision b,c,f,g,h,p,r,s,t,selci_epslon,selci_pythag
 c
 c     this subroutine is a translation of the algol procedure tqlrat,
@@ -187,4 +228,5 @@ c
    20 selci_pythag = p
       return
       end
-
+C>
+C> @}

src/selci/tred1.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_tred1(nm,n,a,d,e,e2)
 *
 * $Id$
@@ -136,4 +138,4 @@ c
 c
       return
       end
-
+C> @}

src/selci/vadd.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_vadd(n,a,ia,b,ib,c,ic)
 *
 * $Id$
@@ -11,3 +13,4 @@ C
  10   continue
 c
       end
+C> @}

src/selci/wrtcon.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_wrtcon(ifllog,iocc,norbs)
 *
 * $Id$
@@ -29,6 +31,7 @@ c
       endif
 c
       end
+C> @}
 
             
       

src/selci/wtconf.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_wtconf(ifl,irefo,nn,indxci,mm)
 *
 * $Id$
@@ -8,3 +10,4 @@ c
       write(1) indxci
 c
       end
+C> @}

src/selci/wthcon.f

@@ -1,3 +1,5 @@
+C> \ingroup selci
+C> @{
       subroutine selci_wthcon(iflcon, title, multi, nelec, issss,
      $     norbs, 
      &     nnsmax, nci, noconf, nintpo, nbitpi, nbpsy, isym, nsym,
@@ -20,3 +22,4 @@ c      write(6,*) title
 c      write(6,*) multi, nelec, issss, norbs, nnsmax, nci, noconf,
 c     &     nintpo, nbitpi, nbpsy, isym, inttyp, nsneed
       end
+C> @}