[Pw_forum] GIPAW: error in output
Lorenzo Paulatto
lorenzo.paulatto at impmc.upmc.fr
Wed Oct 26 09:52:13 CEST 2011
In data 25 ottobre 2011 alle ore 19:05:45, Carlo Nervi
<carlo.nervi at unito.it> ha scritto:
> It should be the line with dclamda (see below) - no idea how to furter
> help the developers ;-)
This *could* originate at line 137 where "t" is computed by h_psi (which
is really ch_psi_all); t is then used to compute c.
In another contest I found cg_solve_all call to ch_psi_all to be extremely
easy to get wrong, as it secretly uses a few global variable. The gipaw
code uses its own version of ch_psi_all, but how to be sure that is is
that one being called an not the one from phonon? In principle the linking
order would guarantee it, but better safe then sorry.
The two subroutines have the same name, but use different global
variables. I would try this, just to be sure:
* in GIPAW/ch_psi_all.f90 and GIPAW/greenfunction.f90 change "ch_psi_all"
to "ch_psi_all_gipaw" (2 changes for file, 4 total)
* recompile
* run
cheers
--
Lorenzo Paulatto IdR @ IMPMC/UPMC CNRS & Université Paris 6
phone: +33 (0)1 44275 084 / skype: paulatz
www: http://www.impmc.upmc.fr/~paulatto/
mail: 23-24/4é14 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05
More information about the Pw_forum
mailing list