Commit 01d3cc4f authored by Uliana Alekseeva's avatar Uliana Alekseeva

ylm4.f changed to ylm4.f90

parent f391706f
...@@ -51,12 +51,12 @@ inpgen/lapw_input.f inpgen/struct_input.f inpgen/write_struct.f ...@@ -51,12 +51,12 @@ inpgen/lapw_input.f inpgen/struct_input.f inpgen/write_struct.f
io/calculator.f global/ss_sym.f global/soc_sym.f math/inv3.f io/rw_symfile.f io/calculator.f global/ss_sym.f global/soc_sym.f math/inv3.f io/rw_symfile.f
global/sort.f kpoints/kptgen_hybrid.f kpoints/od_kptsgen.f kpoints/bravais.f kpoints/divi.f kpoints/brzone.f global/sort.f kpoints/kptgen_hybrid.f kpoints/od_kptsgen.f kpoints/bravais.f kpoints/divi.f kpoints/brzone.f
kpoints/kptmop.f kpoints/kpttet.f init/bandstr1.F kpoints/ordstar.f kpoints/fulstar.f kpoints/kprep.f kpoints/kptmop.f kpoints/kpttet.f init/bandstr1.F kpoints/ordstar.f kpoints/fulstar.f kpoints/kprep.f
kpoints/tetcon.f kpoints/kvecon.f init/boxdim.f math/ylm4.f global/radsra.f math/intgr.F global/differ.f math/inwint.f kpoints/tetcon.f kpoints/kvecon.f init/boxdim.f global/radsra.f math/intgr.F global/differ.f math/inwint.f
math/outint.f xc-pot/gaunt.f math/grule.f math/outint.f xc-pot/gaunt.f math/grule.f
) )
set(inpgen_F90 ${inpgen_F90} global/constants.f90 io/xsf_io.f90 set(inpgen_F90 ${inpgen_F90} global/constants.f90 io/xsf_io.f90
eigen/vec_for_lo.f90 eigen/orthoglo.F90 juDFT/usage_data.F90 eigen/vec_for_lo.f90 eigen/orthoglo.F90 juDFT/usage_data.F90 math/ylm4.f90
global/enpara.f90 global/chkmt.f90 inpgen/inpgen.f90 inpgen/set_inp.f90 inpgen/inpgen_help.f90 io/rw_inp.f90 juDFT/juDFT.F90 global/find_enpara.f90 global/enpara.f90 global/chkmt.f90 inpgen/inpgen.f90 inpgen/set_inp.f90 inpgen/inpgen_help.f90 io/rw_inp.f90 juDFT/juDFT.F90 global/find_enpara.f90
inpgen/closure.f90 inpgen/inpgen_arguments.F90 inpgen/closure.f90 inpgen/inpgen_arguments.F90
juDFT/info.F90 juDFT/stop.F90 juDFT/args.F90 juDFT/time.F90 juDFT/init.F90 juDFT/sysinfo.F90 io/w_inpXML.f90 kpoints/julia.f90 global/utility.F90 juDFT/info.F90 juDFT/stop.F90 juDFT/args.F90 juDFT/time.F90 juDFT/init.F90 juDFT/sysinfo.F90 io/w_inpXML.f90 kpoints/julia.f90 global/utility.F90
......
...@@ -21,7 +21,7 @@ math/rfft.F ...@@ -21,7 +21,7 @@ math/rfft.F
math/sphbes.f math/sphbes.f
math/sphpts.f math/sphpts.f
math/util.F math/util.F
math/ylm4.f math/ylm4.f90
math/difcub.f math/difcub.f
) )
set(fleur_F90 ${fleur_F90} set(fleur_F90 ${fleur_F90}
......
...@@ -6,9 +6,7 @@ ...@@ -6,9 +6,7 @@
INTEGER, SAVE :: lmaxd = -1 ! initial value INTEGER, SAVE :: lmaxd = -1 ! initial value
public ylm4,ylmnorm_init public ylm4,ylmnorm_init
CONTAINS CONTAINS
SUBROUTINE ylm4( SUBROUTINE ylm4(lmax,v,ylm)
> lmax,v,
< ylm)
!************************************************************ !************************************************************
! generate the spherical harmonics for the vector v ! generate the spherical harmonics for the vector v
! using a stable upward recursion in l. (see notes ! using a stable upward recursion in l. (see notes
...@@ -121,8 +119,8 @@ ...@@ -121,8 +119,8 @@
lmaxd = lmax lmaxd = lmax
fpi = 4.0*pimach() fpi = 4.0*pimach()
!$ if (omp_in_parallel()) call juDFT_error( !$ if (omp_in_parallel()) call juDFT_error( &
!$ + "ylmnorm should not called in parallel",calledby="ylm4.f") !$ "ylmnorm should not called in parallel",calledby="ylm4.f")
DO l=0,lmax DO l=0,lmax
lm0 = l*(l+1) + 1 lm0 = l*(l+1) + 1
cd=1.0 cd=1.0
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment