Re: [Wien] maximum number of inequivalent k-points in case.outputkgen
Dear Dr. Kroeker I increased the NKP parameter in the param.inc file (NKP = 2), then rebuild wn_readbands and still the same error in program Thanks in advance, Arena ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] L2main - QTL-B Error
Dear Wien2k User, I initialized the 2D monolayer at different K points 4*4*1, 5*5*1, 6*6*1, and 7*7*1 and also varied the energy parameter 0.3 to 1.1 but when i gave the k points 8*8*1 then it gives the *L2main - QTL-B Error. *So i changed the E parameter from 0.3 to 1.1 and also the RKmax value varied from 7 to 6 but i could not initialize. so please give me idea why *L2main - QTL-B Error *is occuring by increasing the number of Kpoints how we can resolve this problem. -- Thanks & Regards Rajneesh Chaurasiya Research Scholar IIT Jodhpur, India Mob. No. +91-9584499697 +91-7610950803 ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] Compilation error with libxc
On 12/24/2016 03:24 PM, Dr. K. C. Bhamu wrote: Dear Users, I compiled Wien2k_16.1 without any error on latest 64 bit Ubuntu OS on Lenovo laptop with ifort+cc (2015). If I install with LIBXC, I am getting below error(compile.msg file is attached here): Compile time errors (if any) were: SRC_lapw0/compile.msg:inputpars.F(6): error #7013: This module file was not generated by any release of this compiler. [XC_F03_LIB_M] When you install the libxc you probably have to tell him, that your fortran compiler is ifort (and not gfortran). SRC_lapw0/compile.msg:inputpars.F(17): error #6457: This derived type name has not been declared. [XC_F03_FUNC_T] SRC_lapw0/compile.msg:inputpars.F(18): error #6457: This derived type name has not been declared. [XC_F03_FUNC_INFO_T] SRC_lapw0/compile.msg:inputpars.F(161): error #6404: This name does not have a type, and must have an explicit type. [XC_FUNC] SRC_lapw0/compile.msg:inputpars.F(161): error #6404: This name does not have a type, and must have an explicit type. [XC_UNPOLARIZED] SRC_lapw0/compile.msg:inputpars.F(162): error #6404: This name does not have a type, and must have an explicit type. [XC_INFO] SRC_lapw0/compile.msg:inputpars.F(162): error #6404: This name does not have a type, and must have an explicit type. [XC_F03_FUNC_GET_INFO] SRC_lapw0/compile.msg:inputpars.F(163): error #6404: This name does not have a type, and must have an explicit type. [XC_F03_FUNC_INFO_GET_FAMILY] SRC_lapw0/compile.msg:inputpars.F(163): error #6608: In a CASE statement, the case-expr must be of type INTEGER, CHARACTER, or LOGICAL. [XC_F03_FUNC_INFO_GET_FAMILY] SRC_lapw0/compile.msg:inputpars.F(164): error #6404: This name does not have a type, and must have an explicit type. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(164): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(164): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(164): error #6404: This name does not have a type, and must have an explicit type. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(164): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(164): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(173): error #6608: In a CASE statement, the case-expr must be of type INTEGER, CHARACTER, or LOGICAL. [XC_F03_FUNC_INFO_GET_FAMILY] SRC_lapw0/compile.msg:inputpars.F(174): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(174): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(174): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(174): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(189): error #6608: In a CASE statement, the case-expr must be of type INTEGER, CHARACTER, or LOGICAL. [XC_F03_FUNC_INFO_GET_FAMILY] SRC_lapw0/compile.msg:inputpars.F(190): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(190): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(190): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(190): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(195): error #6608: In a CASE statement, the case-expr must be of type INTEGER, CHARACTER, or LOGICAL. [XC_F03_FUNC_INFO_GET_FAMILY] SRC_lapw0/compile.msg:inputpars.F(196): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(196): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_MGGA] SRC_lapw0/compile.msg:inputpars.F(196): error #6601: In a CASE statement, the case-value must be a constant expression. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:inputpars.F(196): error #6612: In a CASE statement, the case-value must be of type INTEGER, CHARACTER, or LOGICAL. [XC_FAMILY_HYB_MGGA] SRC_lapw0/compile.msg:/tmp/ifort150RO5.i(813): catastrophic error: Too many errors, exiting SRC_lapw0/compile.msg:make[1]: *** [inputpars.o] Error 1 SRC_lapw0/compile.msg:make: *** [seq]
Re: [Wien] L2main - QTL-B Error
In most cases where QL-B occurs, the case.struct file is wrong. There are some cases (samll RMTs), where the default input can give problems. I have described how to analyse qtl-b errors several times in the mailing list. Just the message: qtl-b errors occur, please help; gives us no information. And yes, the problems can occur only on certain k-meshes. Am 27.12.2016 um 11:51 schrieb Rajneesh Chaurasiya: Dear Wien2k User, I initialized the 2D monolayer at different K points 4*4*1, 5*5*1, 6*6*1, and 7*7*1 and also varied the energy parameter 0.3 to 1.1 but when i gave the k points 8*8*1 then it gives the *L2main - QTL-B Error. *So i changed the E parameter from 0.3 to 1.1 and also the RKmax value varied from 7 to 6 but i could not initialize. so please give me idea why *L2main - QTL-B Error *is occuring by increasing the number of Kpoints how we can resolve this problem. -- Thanks & Regards Rajneesh Chaurasiya Research Scholar IIT Jodhpur, India Mob. No. +91-9584499697 +91-7610950803 ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- Peter Blaha Inst.Materials Chemistry TU Vienna Getreidemarkt 9 A-1060 Vienna Austria +43-1-5880115671 ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] mpi compilation problem
The attached hmsec.F for lapwso contains the old and new Scalapack routines. Add -Dold_scalapack to the parallel compiler options. Please note: There are cases, where the old Scalapack diagonalization fails. On 12/23/2016 03:47 PM, cesar wrote: Hi, I'm having a problem to get WIEN2k_16 installed. I can compile wien2k_14.2 perfectly but wien2k_16 is impossible (LIBXC and ELPA will not be included for now). The problem seems related with mpi versions for lapw1 and lapwso : seclr4.o: In function `seclr4_': seclr4_tmp_.F:(.text+0x981): undefined reference to `pdsyevr_' seclr4_tmp_.F:(.text+0x1330): undefined reference to `pdsyevr_' make[1]: *** [lapw1_mpi] Error 1 make[1]: Leaving directory `/home/wien2k/wien2k_16/SRC_lapw1' make: *** [rp] Error 2 hmsec.o: In function `hmsec_': hmsec.F:(.text+0x332c): undefined reference to `pzheevr_' hmsec.F:(.text+0x3aa5): undefined reference to `pzheevr_' make[1]: *** [lapwso_mpi] Error 1 make[1]: Leaving directory `/home/wien2k/wien2k_16/SRC_lapwso' make: *** [para] Error 2 I'm using the old : openmpi 1.6 and composerxe-2011.2.137. For now, I don't like so much to the idea of updating ... Are there any other way to get the wien2k 16 ? Thanks in advance. César ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- P.Blaha -- Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna Phone: +43-1-58801-165300 FAX: +43-1-58801-165982 Email: bl...@theochem.tuwien.ac.atWIEN2k: http://www.wien2k.at WWW: http://www.imc.tuwien.ac.at/TC_Blaha -- hmsec.F.gz Description: GNU Zip compressed data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Optic error
Dear all I am trying to simulate Cu2ZnSnS4 structure. I have done the initialization and scf run. Now in property calculation i calculated the DOS. However during optical property calculation i am getting error of "Cannot allocate array - overflow on array size calculation". I have used K points equal to 1000 for optic calculation. I also have tried with lower no of k points i.e 8, 64 etc. but the problem remained same. What i am observing is that Kmin, Kmax in the error is showing 1310722. which i think is not reasonable. For optic calculation i have done x kgen x lapw1 x lapw2 -fermi x optic after this i am getting this problem [ambeshst@IITJ_HPC kczgs]$ x optic emin,emax,nbvalmax -5.003.00 creating ust kmin,kmax-1310722-1310722-1310722 1310722 1310722 1310722 forrtl: severe (179): Cannot allocate array - overflow on array size calculation. Image PCRoutineLine Source opticc 004E123A Unknown Unknown Unknown opticc 004DFDB5 Unknown Unknown Unknown opticc 00491D26 Unknown Unknown Unknown opticc 00457215 Unknown Unknown Unknown opticc 0047AA3A Unknown Unknown Unknown opticc 0042617A planew_86 planew_tmp.f opticc 004380FC mom_mat_ 588 sph-UP_tmp.f opticc 0042064D MAIN__453 opmain.f opticc 0040337C Unknown Unknown Unknown libc.so.6 00344221D9C4 Unknown Unknown Unknown opticc 00403289 Unknown Unknown Unknown 27.105u 0.040s 0:28.25 96.0%0+0k 0+0io 17pf+0w error: command /home/IITJHOME/ambeshst/win2k13mpi/opticc optic.def failed I have checked the memory issue but it should not be the case since i am simulating the structure on HPC which has very high memory. Thanks -- *Goutam Kumar Gupta* *Mob:8561995547* ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] maximum number of inequivalent k-points in case.outputkgen
Dear Prof. Blaha and Dr. Kroeker Thank you for your prompt responses and attached file. As suggested Dr. Kroeker, I increased the numbers in arbmsh.f. The stars disappeared in outputkgen files. I also used the wien2k files attached by prof. Blaha and it works properly. However, I still have problems with xcrysden. When I increased both MAX_NKPT=99 and MAX_IRNKP= 2 in original wn_readbakgen.f file and I get ERROR: while executing wn_readbands program: Error Message: STOP NKP TOO SMALL. It seems this error I get when the Wien2k No. Of INEQIVALEN K-points is above (in my case, for 200 000 k-points, I get 13776 No. Of INEQUIVALENT K-POINTS and the xcrysden ERROR, but for 140 000 kpoints with 9324 INEQUIVALENT k-points there everything is working properly). For the file attached by prof. Blaha, I get errors during compilation (without any changes inside): gfortran -O2 -c wn_readbakgen.f wn_readbakgen.f:110.33: read(line,'(43x,i7,err=111)') nkpt 1 Error: Positive width required in format specifier E at (1) wn_readbakgen.f:111.72: goto 112 1 Error: Label 112 referenced at (1) is never defined make[1]: *** [wn_readbakgen.o] Error 1 Thanks in advance, Arena, Moscow Re: [Wien] maximum number of inequivalent k-points in case.outputkgen Peter Blaha Wed, 21 Dec 2016 03:34:38 -0800 I attach a few subroutines for SRC_kgen, which have been changed for larger formats. They should now not produce "***" up to 100 k-points. The problem is, that for 1 000 000 k-points, you also have to change wn_readbakgen.f in xcrysden. (Or you do not change the "6000 format line" in reduz.f ( and stay at max 99 k-points.) When it works, please let me know. On 12/20/2016 06:52 PM, Arena Konta wrote: Dear win2k-user, I am running wien version 14.2 and have similar problem to that observed in thread http://wien.zeus.theochem.tuwien.ac.narkive.com/pXOQN5tX/maximum-number-of-k-points - the answer does not work for me. I also try to do calculation a Fermi surface with a very dense k-mesh (above 10, Spacegroup: 129_P4/nmm), e.g.: x kgen 12 or using xcrysden In my case.outputkgen file, the maximum number of inequivalent k-points can only be lower than . When the number of k points is larger, it gives * in the rows, as shown below: case.outputkgen 9997( 0.234479 0.456618 0.00)9998 ( 0.234479 0.456618 0.012504) ( 0.234479 0.456618 0.025008) ( 0.234479 0.456618 0.037512) ( 0.234479 0.456618 0.050016) ( 0.234479 0.456618 0.062520) ( 0.234479 0.456618 0.075024) ( 0.234479 0.456618 0.087528) ( 0.234479 0.456618 0.100032) ( 0.234479 0.456618 0.112536) ( 0.234479 0.456618 0.125040) ( 0.234479 0.456618 0.137544) ( 0.234479 0.456618 0.150048) ( 0.234479 0.456618 0.162553) in the reduz.f file i have 100 format(i6,4x,3i4,i8,2f10.5), so everything seems to be ok. I checked this out using gfortran and ifort compilation on two different machine type and still the issue persists. Can you help me? Regards, Arena, Moscow ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html -- P.Blaha -- Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna Phone: +43-1-58801-165300 FAX: +43-1-58801-165982 Email: bl...@theochem.tuwien.ac.atWIEN2k: http://www.wien2k.at WWW: http://www.imc.tuwien.ac.at/TC_Blaha -- Attachment: kgen-update.tar.gz Description: GNU Zip compressed data ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] maximum number of inequivalent k-points in case.outputkgen
To get around the error in xcrysden, you need to increase the NKP parameter in the param.inc file that comes with the parts of SRC_spaghetti included with xcrysden (xcrysden-1.5.60/F/SRC_spaghetti/param.inc) and rebuild wn_readbands -- Dr. Martin Kroekermar...@ruby.chemie.uni-freiburg.de c/o Prof.Dr. Caroline Roehr Institut fuer Anorganische und Analytische Chemie der Universitaet Freiburg ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] mpi compilation problem
Hi, Thank you in advance. I will check your solución as soon as i could but I am pretty sure it will be fine. Seems like the new mkl does not installs cluster libs under a non-commertial license. For now, I ve got no errors during the compilation by using the static libscalapak.a from netlib scalapac (2.0.2) and called from the same directory as the old mkl scalapack. But, I m affraid this lib is not as highly optimized as the mkl one . Enviado desde mi iPhone > El 27 dic 2016, a las 15:17, Peter Blaha> escribió: > > The attached hmsec.F for lapwso contains the old and new Scalapack routines. > > Add -Dold_scalapack to the parallel compiler options. > > Please note: There are cases, where the old Scalapack diagonalization fails. > > > >> On 12/23/2016 03:47 PM, cesar wrote: >> Hi, >> >> I'm having a problem to get WIEN2k_16 installed. >> >> I can compile wien2k_14.2 perfectly but wien2k_16 is impossible (LIBXC >> and ELPA will not be included for now). >> >> The problem seems related with mpi versions for lapw1 and lapwso : >> >> seclr4.o: In function `seclr4_': >> seclr4_tmp_.F:(.text+0x981): undefined reference to `pdsyevr_' >> seclr4_tmp_.F:(.text+0x1330): undefined reference to `pdsyevr_' >> make[1]: *** [lapw1_mpi] Error 1 >> make[1]: Leaving directory `/home/wien2k/wien2k_16/SRC_lapw1' >> make: *** [rp] Error 2 >> >> >> hmsec.o: In function `hmsec_': >> hmsec.F:(.text+0x332c): undefined reference to `pzheevr_' >> hmsec.F:(.text+0x3aa5): undefined reference to `pzheevr_' >> make[1]: *** [lapwso_mpi] Error 1 >> make[1]: Leaving directory `/home/wien2k/wien2k_16/SRC_lapwso' >> make: *** [para] Error 2 >> >> >> I'm using the old : openmpi 1.6 and composerxe-2011.2.137. >> >> For now, I don't like so much to the idea of updating ... >> Are there any other way to get the wien2k 16 ? >> >> Thanks in advance. >> César >> ___ >> Wien mailing list >> Wien@zeus.theochem.tuwien.ac.at >> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien >> SEARCH the MAILING-LIST at: >> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html > > -- > > P.Blaha > -- > Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna > Phone: +43-1-58801-165300 FAX: +43-1-58801-165982 > Email: bl...@theochem.tuwien.ac.atWIEN2k: http://www.wien2k.at > WWW: http://www.imc.tuwien.ac.at/TC_Blaha > -- > > ___ > Wien mailing list > Wien@zeus.theochem.tuwien.ac.at > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > SEARCH the MAILING-LIST at: > http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html