forum@abinit.org
Subject: The ABINIT Users Mailing List ( CLOSED )
List archive
- From: Xavier Gonze <gonze@pcpm.ucl.ac.be>
- To: forum@abinit.org
- Subject: Re: [abinit-forum] gfortran 4.2.1 bug: conjg(transpose(z(:,:))
- Date: Fri, 31 Aug 2007 08:24:03 +0200
Dear Takeshi Nishimatsu,
Thanks ! Likely the very last modif in v5.4.3 .
Xavier
On 31 Aug 2007, at 03:01, Takeshi Nishimatsu wrote:
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
- gfortran 4.2.1 bug: conjg(transpose(z(:,:)), Takeshi Nishimatsu, 08/31/2007
- Re: [abinit-forum] gfortran 4.2.1 bug: conjg(transpose(z(:,:)), Anglade Pierre-Matthieu, 08/31/2007
- Re: [abinit-forum] gfortran 4.2.1 bug: conjg(transpose(z(:,:)), Xavier Gonze, 08/31/2007
Archive powered by MHonArc 2.6.16.