~ubuntu-branches/ubuntu/trusty/nwchem/trusty-proposed

« back to all changes in this revision

Viewing changes to src/nwpw/band/minimizer/band_minimizer.F

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Daniel Leidert, Andreas Tille, Michael Banck
  • Date: 2013-07-04 12:14:55 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130704121455-5tvsx2qabor3nrui
Tags: 6.3-1
* New upstream release.
* Fixes anisotropic properties (Closes: #696361).
* New features include:
  + Multi-reference coupled cluster (MRCC) approaches
  + Hybrid DFT calculations with short-range HF 
  + New density-functionals including Minnesota (M08, M11) and HSE hybrid
    functionals
  + X-ray absorption spectroscopy (XAS) with TDDFT
  + Analytical gradients for the COSMO solvation model
  + Transition densities from TDDFT 
  + DFT+U and Electron-Transfer (ET) methods for plane wave calculations
  + Exploitation of space group symmetry in plane wave geometry optimizations
  + Local density of states (LDOS) collective variable added to Metadynamics
  + Various new XC functionals added for plane wave calculations, including
    hybrid and range-corrected ones
  + Electric field gradients with relativistic corrections 
  + Nudged Elastic Band optimization method
  + Updated basis sets and ECPs 

[ Daniel Leidert ]
* debian/watch: Fixed.

[ Andreas Tille ]
* debian/upstream: References

[ Michael Banck ]
* debian/upstream (Name): New field.
* debian/patches/02_makefile_flags.patch: Refreshed.
* debian/patches/06_statfs_kfreebsd.patch: Likewise.
* debian/patches/07_ga_target_force_linux.patch: Likewise.
* debian/patches/05_avoid_inline_assembler.patch: Removed, no longer needed.
* debian/patches/09_backported_6.1.1_fixes.patch: Likewise.
* debian/control (Build-Depends): Added gfortran-4.7 and gcc-4.7.
* debian/patches/10_force_gcc-4.7.patch: New patch, explicitly sets
  gfortran-4.7 and gcc-4.7, fixes test suite hang with gcc-4.8 (Closes:
  #701328, #713262).
* debian/testsuite: Added tests for COSMO analytical gradients and MRCC.
* debian/rules (MRCC_METHODS): New variable, required to enable MRCC methods.

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
1
*
2
 
* $Id: band_minimizer.F 21108 2011-09-18 00:35:31Z bylaska $
 
2
* $Id: band_minimizer.F 23444 2013-01-15 19:08:17Z bylaska $
3
3
*
4
4
***********************************************************************
5
5
*                     band_minimizer                                                              *
88
88
      character*12 control_boundry
89
89
      external     control_boundry
90
90
 
91
 
      logical  brillioun_print
 
91
      logical  brillioun_print,control_Mulliken
92
92
      integer  brillioun_nbrillioun,brillioun_nbrillq
93
93
      real*8   brillioun_weight_brdcst,brillioun_ks_brdcst
94
94
      real*8   brillioun_k_brdcst,brillioun_weight
95
 
      external brillioun_print
 
95
      external brillioun_print,control_Mulliken
96
96
      external brillioun_nbrillioun,brillioun_nbrillq
97
97
      external brillioun_weight_brdcst,brillioun_ks_brdcst 
98
98
      external brillioun_k_brdcst,brillioun_weight
148
148
c      character*255 cpsp_comment,comment
149
149
c      external      cpsp_comment
150
150
 
 
151
      integer  ion_nconstraints,ion_ndof
 
152
      external ion_nconstraints,ion_ndof
 
153
 
151
154
    
152
155
 
153
156
 
343
346
         else
344
347
           write(luout,1130) "spin orbit"
345
348
         end if
 
349
 
 
350
         call v_bwexc_print(luout,control_gga())
346
351
         
347
 
         IF (control_gga().eq.0) THEN
348
 
            write(luout,1131) 'Vosko et al parameterization'
349
 
         ELSE IF (control_gga().eq.10) THEN
350
 
            write(luout,1131)
351
 
     >      'PBE96 (White and Bird) parameterization'
352
 
         ELSE IF (control_gga().eq.11) THEN
353
 
            write(luout,1131)
354
 
     >      'BLYP (White and Bird) parameterization'
355
 
         ELSE IF (control_gga().eq.12) THEN
356
 
            write(luout,1131)
357
 
     >      'revPBE (White and Bird) parameterization'
358
 
         ELSE IF (control_gga().eq.13) THEN
359
 
            write(luout,1131)
360
 
     >      'PBEsol (White and Bird) parameterization'
361
 
 
362
 
         ELSE IF (control_gga().eq.110) THEN
363
 
            write(luout,1131) 
364
 
     >      'PBE0 (White and Bird) parameterization'
365
 
         ELSE IF (control_gga().eq.112) THEN
366
 
            write(luout,1131)
367
 
     >      'revPBE0 (White and Bird) parameterization'
368
 
         ELSE IF (control_gga().eq.113) THEN
369
 
            write(luout,1131) 
370
 
     >      'BNL (White and Bird) parameterization'
371
 
         ELSE IF (control_gga().eq.114) THEN
372
 
            write(luout,1131) 
373
 
     >      'HSE (White and Bird) parameterization'
374
 
         ELSE IF (control_gga().eq.200) THEN
375
 
            write(luout,1131) 'Hartree-Fock'
376
 
         ELSE
377
 
            write(luout,1131) 'unknown parameterization'
378
 
            call errquit('bad exchange_correlation',0, INPUT_ERR)
379
 
         END IF
380
352
         call band_print_HFX(luout)
381
353
         write(luout,1140)
382
354
         do ia = 1,ion_nkatm()
414
386
         end do
415
387
         write(luout,1200) cx,cy,cz
416
388
         write(luout,1210) gx,gy,gz
 
389
         write(luout,1211) ion_nconstraints(),ion_ndof()
417
390
         endif
418
391
         write(luout,1220) en(1),en(cpsi_ispin()),
419
392
     >                 ' (fourier space)'
621
594
cc      call band_kp_structure(rtdb)
622
595
 
623
596
*     **** calculate Mulliken Populations ***
 
597
      if (control_Mulliken()) then
 
598
         call cpsi_Mulliken(rtdb)
 
599
      end if
 
600
 
624
601
 
625
602
*     **** write geometry to rtdb ****
626
603
      call ion_write(rtdb)
740
717
     >       ' ) - atomic mass= ',F7.3,' - fixed')
741
718
 1200 FORMAT(5X,'   G.C.  ',' (',3F11.5,' )')
742
719
 1210 FORMAT(5X,'   C.O.M.',' (',3F11.5,' )')
 
720
 1211 FORMAT(5X,'   number of constraints = ', I6,' ( DOF = ',I6,' )' )
743
721
 1220 FORMAT(/' number of electrons: spin up=',F8.2,
744
722
     >                          '  spin down=',F8.2,A)
745
723
 1221 FORMAT( ' number of orbitals:  spin up=',I8,