← Back to branch summary

~ubuntu-branches/ubuntu/utopic/nwchem/utopic

« back to all changes in this revision

Viewing changes to src/nwpw/nwpwlib/D3dB/D3dB-new.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
http://bugs.debian.org/696361
* 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:
src/nwpw/nwpwlib/D3dB/D3dB-new.F
2
2
#define NBLOCKS 2
3
3
 
4
4
*
5
 
* $Id: D3dB-new.F 20420 2011-06-13 20:19:47Z bylaska $
 
5
* $Id: D3dB-new.F 23927 2013-03-30 00:21:12Z edo $
6
6
*
7
7
 
8
8
*     ***********************************************************
2345
2345
 
2346
2346
      indx = 1
2347
2347
      do j=1,(ny*nq*ne)
 
2348
CDIR$ NOVECTOR
2348
2349
         do i=nx,2,-1
2349
2350
            B(indx+i) = A(indx+i-1)
2350
2351
         end do
2718
2719
      end
2719
2720
 
2720
2721
 
 
2722
*     ***********************************
 
2723
*     *                                 *
 
2724
*     *        D3dB_rr_SqrAdd           *
 
2725
*     *                                 *
 
2726
*     ***********************************
 
2727
 
 
2728
      subroutine D3dB_rr_SqrAdd(nb,A,C)
 
2729
      implicit none
 
2730
      integer nb
 
2731
      real*8 A(*)
 
2732
      real*8 C(*)
 
2733
 
 
2734
#include "D3dB.fh"
 
2735
 
 
2736
 
 
2737
      integer i
 
2738
 
 
2739
      do i=1,n2ft3d_map(nb)
 
2740
         C(i) = C(i) + A(i)**2
 
2741
      end do
 
2742
      return
 
2743
      end
 
2744
 
 
2745
 
 
2746
 
2721
2747
 
2722
2748
*     ***********************************
2723
2749
*     *                                 *
3832
3858
      integer i
3833
3859
 
3834
3860
      do i=1,n2ft3d_map(nb)
3835
 
         B(i) = B(i) + alpha* B(i)
 
3861
         B(i) = B(i) + alpha* A(i)
3836
3862
      end do
3837
3863
 
3838
3864
      return
3861
3887
      !do q=1,nq(nb)
3862
3888
      !do j=1,ny(nb)
3863
3889
      !do i=1,nx(nb)
 
3890
CDIR$ NOVECTOR
3864
3891
      do index = 1,n2ft3d_map(nb)
3865
3892
         !index = i + (j-1)*(nx(nb)+2) + (q-1)*(nx(nb)+2)*ny(nb)
3866
3893
         if (dabs(B(index)) .le. eta) then
3900
3927
      !do q=1,nq(nb)
3901
3928
      !do j=1,ny(nb)
3902
3929
      !do i=1,nx(nb)
 
3930
CDIR$ NOVECTOR
3903
3931
      do index = 1,n2ft3d_map(nb)
3904
3932
         !index = i + (j-1)*(nx(nb)+2) + (q-1)*(nx(nb)+2)*ny(nb)
3905
3933
         if (dabs(A(index)) .le. eta) then