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

« back to all changes in this revision

Viewing changes to src/nwpw/nwpwlib/ion/shake/molecular.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: molecular.F 20171 2011-04-01 18:22:27Z bylaska $
 
2
* $Id: molecular.F 23425 2013-01-15 16:51:25Z bylaska $
3
3
*
4
4
 
5
5
*     ***************************
115
115
         call pspw_molecule_msize(msize)
116
116
         do m=1,msize
117
117
            call pspw_molecule_data(m,na,indxa,blist,cyclic)
 
118
            if (cyclic.eq.1)  then
 
119
               call ion_add_constraint(na)
 
120
            else
 
121
               call ion_add_constraint(na-1)
 
122
            end if
 
123
            
118
124
            if (taskid.eq.MASTER) then
119
125
            if (cyclic.eq.1) then
120
126
               write(*,'(A,A,10(A,I3,I3,A,F6.3,A))') 
162
168
      return
163
169
      end 
164
170
 
 
171
*     *******************************************
 
172
*     *                                         *
 
173
*     *      molecular_number_constraints       *
 
174
*     *                                         *
 
175
*     *******************************************
 
176
      integer function molecular_number_constraints()
 
177
      implicit none
 
178
 
 
179
      integer m,msize,na,nb,nc,cyclic,indxa(20)
 
180
      double precision blist(20)
 
181
 
 
182
      logical value
 
183
      common /pspw_molecular/ value
 
184
 
 
185
      nc = 0
 
186
      if (value) then
 
187
         call pspw_molecule_msize(msize)
 
188
         do m=1,msize
 
189
            call pspw_molecule_data(m,na,indxa,blist,cyclic)
 
190
            nb = na-1
 
191
            if (cyclic.eq.1) nb = nb+1
 
192
            nc = nc + nb
 
193
         end do
 
194
      end if
 
195
 
 
196
      molecular_number_constraints = nc
 
197
      return
 
198
      end
 
199
 
 
200
 
165
201
 
166
202
*     ***************************
167
203
*     *                         *
354
390
              write(*,120) indxa(nb),indxa(1),blist(nb),
355
391
     >                     dbl_mb(indx2),dbl_mb(indx),
356
392
     >                       dbl_mb(indx3)-dbl_mb(indx)**2
357
 
                indx  = indx + 1
358
393
              indx  = indx + 1
359
394
              indx2 = indx2 + 1
360
395
              indx3 = indx3 + 1
364
399
     >                       dbl_mb(indx2),dbl_mb(indx),
365
400
     >                       dbl_mb(indx3)-dbl_mb(indx)**2
366
401
                indx  = indx + 1
367
 
                indx  = indx + 1
368
402
                indx2 = indx2 + 1
369
403
                indx3 = indx3 + 1
370
404
              end do
390
424
      return
391
425
      end
392
426
 
 
427
 
393
428
      subroutine molecular_sub_force_average(rion,fion,af,af2,ac)
394
429
      implicit none
395
430
      real*8 rion(3,*)