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

« back to all changes in this revision

Viewing changes to src/util/ga_extra.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
      subroutine ga_scale_lh(g_a,vector)
2
 
c $Id: ga_extra.F 19707 2010-10-29 17:59:36Z d3y133 $
 
2
c $Id: ga_extra.F 23018 2012-10-29 17:39:42Z d3y133 $
3
3
      implicit none
4
4
#include "global.fh"
5
5
#include "mafdecls.fh"
20
20
      integer ma_type, dim1, dim2, n
21
21
      integer i, j, ilo, ihi, jlo, jhi
22
22
      double precision number
 
23
      double complex znumber
23
24
c
24
25
      call ga_sync()
25
26
      call ga_inquire(g_a, ma_type, dim1, dim2)
26
 
      if (ma_type.ne.mt_dbl) then
27
 
         write(0,*) ' ma_type ',ma_type,' mt_dbl ',mt_dbl
 
27
      if (ma_type.eq.mt_dbl) then
 
28
      else if (ma_type.eq.mt_dcpl) then
 
29
      else
 
30
        write(0,*) ' ma_type ',ma_type,' mt_dbl  ',mt_dbl
 
31
        write(0,*) ' ma_type ',ma_type,' mt_dcpl ',mt_dcpl
28
32
        call errquit
29
 
     $     ('ga_trace_diagonal: array is not real', g_a,0)
 
33
     $     ('ga_trace_diagonal: array is not real nor complex ', g_a,0)
30
34
      endif
31
35
      n = dim1
32
36
c
34
38
c
35
39
      call ga_distribution(g_a, ga_nodeid(), ilo, ihi, jlo, jhi)
36
40
      if (ilo.gt.0 .and. jlo.gt.0) then
37
 
         do i = ilo,ihi
38
 
            do j = jlo,jhi
39
 
               call ga_get(g_a, i, i, j, j, number, 1)
40
 
               number = number*vector(i)
41
 
               call ga_put(g_a,i, i, j, j, number, 1)
42
 
            enddo
43
 
         enddo
 
41
         if (ma_type.eq.MT_DBL) then
 
42
            do i = ilo,ihi
 
43
               do j = jlo,jhi
 
44
                  call ga_get(g_a, i, i, j, j, number, 1)
 
45
                  number = number*vector(i)
 
46
                  call ga_put(g_a,i, i, j, j, number, 1)
 
47
               enddo
 
48
            enddo
 
49
         else if (ma_type.eq.MT_DCPL) then
 
50
            do i = ilo,ihi
 
51
               do j = jlo,jhi
 
52
                  call ga_get(g_a, i, i, j, j, znumber, 1)
 
53
                  znumber = znumber*vector(i)
 
54
                  call ga_put(g_a,i, i, j, j, znumber, 1)
 
55
               enddo
 
56
            enddo
 
57
         endif
44
58
      endif
45
59
c
46
60
      call ga_sync()