← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/property/hnd_hyperfine_zora.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/property/hnd_hyperfine_zora.F
1
1
      subroutine hnd_hyperfine_zora(rtdb,basis,geom)
2
 
c $Id: hnd_hyperfine_zora.F 21478 2011-11-09 17:02:28Z niri $
 
2
c $Id: hnd_hyperfine_zora.F 23558 2013-02-11 03:42:43Z niri $
3
3
      implicit none
4
4
#include "errquit.fh"
5
5
#include "global.fh"
251
251
        write(*,*) 'after-creating g_rhs ...'
252
252
        write(*,102)npol,nocc(1),nocc(2),
253
253
     &              nclosed(1),nclosed(2),
254
 
     &              nvirt(1),nvirt(2),scftyp
 
254
     &              nvirt(1),nvirt(2),scftyp,ndens
255
255
 102   format('BEF pre-fock::npol=',i3,' nocc=(',i3,',',i3,') ',
256
256
     &       'nclos=(',i3,',',i3,') ',
257
 
     &       'nvirt=(',i3,',',i3,') scftyp=',a,')')
 
257
     &       'nvirt=(',i3,',',i3,') scftyp=',a,',ndens=',i3,')')
258
258
       endif
259
259
       call add_H10_hfine(g_rhs, ! out: ga_rhs(a,i)=ga_rhs(a,i)+H10(a,i)
260
260
     &                    ga_Fji_hfine,vectors, 
307
307
     &                 mt_log,1,skip_cphf_ev_hyp))        
308
308
     &  skip_cphf_ev_hyp = .false.       
309
309
      if (.not.(skip_cphf_ev_hyp)) then ! ----- if-skip-cphf--START
310
 
       if (ga_nodeid().eq.0)
311
 
     &  write(*,*) 'COMPUTE cphf hyp data ...'
 
310
c       if (ga_nodeid().eq.0)
 
311
c     &  write(*,*) 'COMPUTE cphf hyp data ...'
312
312
c     Write ga_rhs to disk 
313
313
       call util_file_name('cphf_rhs',.true.,.true.,cphf_rhs)
314
314
       call util_file_name('cphf_sol',.true.,.true.,cphf_sol)
444
444
       enddo ! end-loop-ix
445
445
      enddo ! end-loop-iatom
446
446
c +++++++++ print-total-dia-transferred +++ END
447
 
       do ispin=1,npol
 
447
       do ispin=1,ndens 
448
448
        if (.not.ga_destroy(g_dens(ispin))) call 
449
449
     &    errquit('hnd_hfine_zora: ga_destroy failed g_dens',
450
450
     &    0,GA_ERR)
652
652
     &    errquit('hnd_hyperfine_zora: ga_destroy failed vectors',
653
653
     &             0,GA_ERR)
654
654
      enddo
 
655
        if (.not. ga_destroy(g_AtNr1)) call errquit(
 
656
     &    'hnd_giaox_zora: ga_destroy failed ',0, GA_ERR)   
655
657
      if (.not.ma_pop_stack(l_zan)) call
656
658
     &    errquit('hnd_hyperfine_zora: ma_pop_stack failed k_zan',
657
659
     &             0,MA_ERR)
1388
1390
     &                   basis,nbf,iat,rtdb)
1389
1391
c +++++++++ NOW: do ddot product to get diagonalized tensor +++ END
1390
1392
c ++++++++++ CHECK diagonalization in para hyperfine ++++++++ END
 
1393
      if (.not. ga_destroy(g_p10)) call errquit(
 
1394
     &  'create_munu4nbo: ga_destroy failed ',0, GA_ERR)
1391
1395
      enddo ! end-loop-iat
1392
1396
        if (ga_nodeid().eq.0)
1393
1397
     &   write(*,*) 'CHCooooooooooooo',
1427
1431
     &  'create_munu4nbo: ga_destroy failed ',0, GA_ERR)
1428
1432
      if (.not. ga_destroy(g_sdens1)) call errquit(
1429
1433
     &  'create_munu4nbo: ga_destroy failed ',0, GA_ERR)
1430
 
      if (.not. ga_destroy(g_p10)) call errquit(
1431
 
     &  'create_munu4nbo: ga_destroy failed ',0, GA_ERR)
1432
1434
      if (.not. ga_destroy(g_c1)) call errquit(
1433
1435
     &  'create_munu4nbo: ga_destroy failed ',0, GA_ERR)
1434
1436
      do i=1,npol
1532
1534
       ahi(1) = nbf
1533
1535
       ahi(2) = nbf
1534
1536
       ahi(3) = 3
1535
 
       if (.not.nga_create(MT_DBL,3,ahi,'D10 matrix',alo,g_d1)) call 
 
1537
       if (.not.nga_create(MT_DBL,3,ahi,'get_P10_rot: g_d1 matrix',
 
1538
     &                    alo,g_d1)) call 
1536
1539
     &     errquit('g_d1: nga_create failed g_d1',0,GA_ERR)
1537
 
       if (.not.nga_create(MT_DBL,3,ahi,'D10 matrix',alo,g_p1)) call 
 
1540
       if (.not.nga_create(MT_DBL,3,ahi,'get_P10_rot: g_p1 matrix',
 
1541
     &                     alo,g_p1)) call 
1538
1542
     &     errquit('g_p1: nga_create failed g_p1',0,GA_ERR)
1539
1543
       call ga_zero(g_p1)
1540
1544
      do ispin=1,npol