Skip to Content.
Sympa Menu

forum - gfortran 4.2.1 bug: conjg(transpose(z(:,:))

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

gfortran 4.2.1 bug: conjg(transpose(z(:,:))


Chronological Thread 
  • From: Takeshi Nishimatsu <takeshi@physics.rutgers.edu>
  • To: forum@abinit.org
  • Subject: gfortran 4.2.1 bug: conjg(transpose(z(:,:))
  • Date: Thu, 30 Aug 2007 21:01:23 -0400

Hi,

I encountered this bug in another program.
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=32962
(This bug is fixed on trunk of gcc.)
So I did grep in ABINIT src.

$ pwd ; grep -i 'conjg(transpose' */*.F90
abinit-5.3.5/src
15gw/setup_hamilt.F90:
h(:,:,ik,is)=matmul(conjg(transpose(c(:,:,ik,is))),ctmp(:,:))
15gw/setup_hamilt.F90: htotal(:,:,is)=0.5 * ( htotal(:,:,is) +
conjg(transpose(htotal(:,:,is))) )

$ pwd ; grep -i 'conjg(transpose' */*.F90
abinit-5.4.2/src
15gw/csigme.F90: call
calc_sig_ppm(sp%npwc,nomega,rhotwgp(1:sp%npwc),conjg(transpose(botsq(:,:))),&
15gw/setup_hamilt.F90: h(:,:,ik,is)=
matmul(conjg(transpose(c(:,:,ik,is))),ctmp(:,:))
15gw/setup_hamilt.F90: htotal(:,:,is)= 0.5*( htotal(:,:,is) +
conjg(transpose(htotal(:,:,is))) )
15gw/symf12.F90: ! herm(:,:,io) =
half*(eps(:,:,io)+conjg(transpose(eps(:,:,io))))
15gw/symf12.F90: ! aherm(:,:,io) =
-(0.,1)*half*(eps(:,:,io)-conjg(transpose(eps(:,:,io))))

What I did is only grep.
So it might be false alarm for ABINIT. Sorry,
--
Takeshi Nishimatsu





Archive powered by MHonArc 2.6.16.

Top of Page