Jonathan and Nicola,<div><br></div><div>Thanks for the suggestions. I tried different libraries on the machine I was using and was not able to avoid the problem. However, I found out that if I use pw2wannier and get the .amn, .mmn, and .eig files, take those to my local machine, and run Wannier as a standalone program rather than as a library, that everything is easily wannierized. </div>
<div><br></div><div>So I suppose this suggests that either the difference is the math libraries I'm using on my local machine, some subtlety between the Wannier standalone program and running as a library, or another issue possibly related to the code which calls Wannier as a library.</div>
<div><br></div><div>Anyway, I had the Wannier standalone program write out the two Wannier rotation U matrices, multiplied them together, and then had my original code use that rotation matrix and continue on with its work. Tests on pressures for which both methods work yield the same answer so I am content with this. </div>
<div><br></div><div>Thanks again,</div><div>Brad<br><br><div class="gmail_quote">On Thu, Dec 8, 2011 at 12:40 AM, <span dir="ltr"><<a href="mailto:wannier-request@quantum-espresso.org">wannier-request@quantum-espresso.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Send Wannier mailing list submissions to<br>
<a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:wannier-request@quantum-espresso.org">wannier-request@quantum-espresso.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:wannier-owner@quantum-espresso.org">wannier-owner@quantum-espresso.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Wannier digest..."<br>
<br>
<br>
Today's Topics:<br><br> 1. Re: dis_project ZGESVD error (Jonathan Yates)<br>
2. Re: dis_project ZGESVD error (Nicola Marzari)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 7 Dec 2011 11:11:03 +0000<br>
From: Jonathan Yates <<a href="mailto:jonathan.yates@materials.ox.ac.uk">jonathan.yates@materials.ox.ac.uk</a>><br>
Subject: Re: [Wannier] dis_project ZGESVD error<br>
To: "<a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a>" <<a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a>><br>
Message-ID: <<a href="mailto:7F921E85-8DAB-4ED3-BB90-E396772FCFA7@materials.ox.ac.uk">7F921E85-8DAB-4ED3-BB90-E396772FCFA7@materials.ox.ac.uk</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
<br>
On 6 Dec 2011, at 19:25, Brad Malone wrote:<br>
<br>
> Hi, I am using wannier in library mode for calculations on a system as a function of pressure. I am able to successfully wannierize the system for 20,25,30, and 40 GPa with essentially no change in the p.win file except for a slight change to the frozen window range (to account for the increase in the bandwidth as a function of pressure).<br>
><br>
> However, when I try to do it for 50 GPa, for the same system and parameters, I get the error:<br>
><br>
><br>
> Unitarised projection of Wannier functions<br>
> ------------------------------------------<br>
> A_mn = <psi_m|g_n> --> S = A.A^+ --> U = S^-1/2.A<br>
> In dis_project... ERROR: IN ZGESVD IN dis_project<br>
> K-POINT NKP= 1 INFO= 7<br>
> Exiting.......<br>
> dis_project: problem in ZGESVD 1<br>
><br>
> I was originally using 'random' projections, since that is what has worked for me for the previous pressures. When I saw this message I tried fiddling with the projections, the number of wannier functions, and the number of input bands, but nothing changed (except changing the number of wannier functions changes the value for INFO).<br>
<br>
Brad,<br>
<br>
The error means that that singular value decomposition has failed to converge. It's being used to compute the initial guess for the unitary rotation.<br>
<br>
I've never encountered this error before, and from your description it is not at all obvious what might cause it. Perhaps it is worth trying a different set of math libraries. (eg if you are using mkl, try the netlib lapack).<br>
<br>
Otherwise, if no-one else has a solution, you could send the files (win, eig, amn, mmn) to me.<br>
<br>
Jonathan<br>
<br>
<br>
<br>
<br>
--<br>
Department of Materials, University of Oxford, Parks Road, Oxford, OX1 3PH, UK<br>
tel: <a href="tel:%2B44%20%280%291865%20612797" value="+441865612797">+44 (0)1865 612797</a> <a href="http://users.ox.ac.uk/~oums0549/" target="_blank">http://users.ox.ac.uk/~oums0549/</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 07 Dec 2011 13:19:14 +0100<br>
From: Nicola Marzari <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>><br>
Subject: Re: [Wannier] dis_project ZGESVD error<br>
Cc: "<a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a>" <<a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a>><br>
Message-ID: <<a href="mailto:4EDF59C2.6040606@epfl.ch">4EDF59C2.6040606@epfl.ch</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
<br>
<br>
Dear Brad,<br>
<br>
<br>
a guess is that your random projections, projected on the kohn-sham<br>
manifold, end up having a couple of them that are almost linearly<br>
dependent.<br>
<br>
Can you choose instead a regular mesh of centers for the projections,<br>
or maybe something slightly more physical?<br>
<br>
If your system is large, you could try the gamma sampling only -<br>
that typically does not need projections, and you could use the results<br>
as centers for the next set of projections?<br>
<br>
nic<br>
<br>
<br>
On 12/7/11 12:11 PM, Jonathan Yates wrote:<br>
><br>
> On 6 Dec 2011, at 19:25, Brad Malone wrote:<br>
><br>
>> Hi, I am using wannier in library mode for calculations on a system as a function of pressure. I am able to successfully wannierize the system for 20,25,30, and 40 GPa with essentially no change in the p.win file except for a slight change to the frozen window range (to account for the increase in the bandwidth as a function of pressure).<br>
>><br>
>> However, when I try to do it for 50 GPa, for the same system and parameters, I get the error:<br>
>><br>
>><br>
>> Unitarised projection of Wannier functions<br>
>> ------------------------------------------<br>
>> A_mn =<psi_m|g_n> --> S = A.A^+ --> U = S^-1/2.A<br>
>> In dis_project... ERROR: IN ZGESVD IN dis_project<br>
>> K-POINT NKP= 1 INFO= 7<br>
>> Exiting.......<br>
>> dis_project: problem in ZGESVD 1<br>
>><br>
>> I was originally using 'random' projections, since that is what has worked for me for the previous pressures. When I saw this message I tried fiddling with the projections, the number of wannier functions, and the number of input bands, but nothing changed (except changing the number of wannier functions changes the value for INFO).<br>
><br>
> Brad,<br>
><br>
> The error means that that singular value decomposition has failed to converge. It's being used to compute the initial guess for the unitary rotation.<br>
><br>
> I've never encountered this error before, and from your description it is not at all obvious what might cause it. Perhaps it is worth trying a different set of math libraries. (eg if you are using mkl, try the netlib lapack).<br>
><br>
> Otherwise, if no-one else has a solution, you could send the files (win, eig, amn, mmn) to me.<br>
><br>
> Jonathan<br>
><br>
><br>
><br>
><br>
<br>
<br>
--<br>
<br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
<br>
<br>
End of Wannier Digest, Vol 49, Issue 5<br>
**************************************<br>
</blockquote></div><br></div>