| /dports/math/py-ducc0/ducc0-0.22.0/python/test/ |
| H A D | test_sht.py | 25 def nalm(lmax, mmax): argument
37 ofs += lmax+1-s
43 res[0:lmax+1] = alm[0:lmax+1].real
44 res[lmax+1::2] = np.sqrt(2)*alm[lmax+1:].real
45 res[lmax+2::2] = np.sqrt(2)*alm[lmax+1:].imag
50 return ducc0.misc.vdot(compress_alm(a1, lmax), compress_alm((a2), lmax))
60 lmax, mmax = lmmax
62 nrings = lmax+1
114 nrings = lmax+1
174 mmax = lmax
[all …]
|
| H A D | test_totalconvolve.py | 62 res[0:lmax+1] = alm[0:lmax+1].real
63 res[lmax+1::2] = np.sqrt(2)*alm[lmax+1:].real
64 res[lmax+2::2] = np.sqrt(2)*alm[lmax+1:].imag
76 slm = random_alm(rng, lmax, lmax, 1)[0, :]
95 blm2 = np.zeros((nalm(lmax, lmax),))+0j
110 slm = random_alm(rng, lmax, lmax, ncomp)
123 blm2 = np.zeros((ncomp, nalm(lmax, lmax)))+0j
143 slm = random_alm(rng, lmax, lmax, ncomp).astype("c8")
157 blm2 = np.zeros((ncomp, nalm(lmax, lmax)))+0j
175 slm = random_alm(rng, lmax, lmax, 1)[0, :]
[all …]
|
| /dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/nwpw/nwpwlib/FMM/ |
| H A D | fmm.F | 592 n = (2*lmax+2)*(4*lmax+2)
647 complex*16 A(2*lmax+2,4*lmax+2)
648 complex*16 B(2*lmax+2,4*lmax+2)
686 complex*16 A(2*lmax+2,4*lmax+2)
687 complex*16 B(2*lmax+2,4*lmax+2)
688 complex*16 AB(2*lmax+2,4*lmax+2)
694 n = (2*lmax+2)*(4*lmax+2)
728 complex*16 AB(2*lmax+2,4*lmax+2)
783 n = (2*lmax+2)*(4*lmax+2)
839 complex*16 A(2*lmax+2,4*lmax+2)
[all …]
|
| /dports/science/nwchem-data/nwchem-7.0.2-release/src/nwpw/nwpwlib/FMM/ |
| H A D | fmm.F | 592 n = (2*lmax+2)*(4*lmax+2)
647 complex*16 A(2*lmax+2,4*lmax+2)
648 complex*16 B(2*lmax+2,4*lmax+2)
686 complex*16 A(2*lmax+2,4*lmax+2)
687 complex*16 B(2*lmax+2,4*lmax+2)
688 complex*16 AB(2*lmax+2,4*lmax+2)
694 n = (2*lmax+2)*(4*lmax+2)
728 complex*16 AB(2*lmax+2,4*lmax+2)
783 n = (2*lmax+2)*(4*lmax+2)
839 complex*16 A(2*lmax+2,4*lmax+2)
[all …]
|
| /dports/math/py-ducc0/ducc0-0.22.0/python/demos/ |
| H A D | sht_analysis_demo.py | 8 def nalm(lmax, mmax): argument
12 def random_alm(lmax, mmax, spin, ncomp): argument
16 res[:, 0:lmax+1].imag = 0.
20 ofs += lmax+1-s
28 nrings = lmax+1
30 nrings = lmax+2
32 nrings = 2*lmax+2
34 nrings = 2*lmax+1
35 nphi=2*lmax+2
37 alm = random_alm(lmax, lmax, spin, ncomp)
[all …]
|
| H A D | totalconvolve_demo.py | 24 def nalm(lmax, mmax): argument
36 def compress_alm(alm, lmax): argument
38 res[0:lmax+1] = alm[0:lmax+1].real
39 res[lmax+1::2] = np.sqrt(2)*alm[lmax+1:].real
40 res[lmax+2::2] = np.sqrt(2)*alm[lmax+1:].imag
44 def myalmdot(a1, a2, lmax): argument
45 return ducc0.misc.vdot(compress_alm(a1, lmax), compress_alm(a2, lmax))
47 lmax = 1024 variable
60 slm = random_alm(lmax, lmax, ncomp)
86 nth = lmax+1
[all …]
|
| /dports/databases/grass7/grass-7.8.6/imagery/i.atcorr/ |
| H A D | transform.cpp | 13 lmax = 194.3f; in EtmDN()
20 lmax = 202.4f; in EtmDN()
27 lmax = 158.6f; in EtmDN()
36 lmax = 235.0f; in EtmDN()
41 lmax = 157.5f; in EtmDN()
49 lmax = 31.76f; in EtmDN()
56 lmax = 10.932f; in EtmDN()
75 lmax = 191.6f; in EtmDN()
82 lmax = 196.5f; in EtmDN()
98 lmax = 241.1f; in EtmDN()
[all …]
|
| /dports/science/healpix/Healpix_3.50/src/healpy/healpy/ |
| H A D | sphtfunc.py | 114 lmax=lmax,
127 lmax=lmax,
139 cls = alm2cl(alms1, alms2=alms2, lmax=lmax, mmax=mmax, lmax_out=lmax, nspec=nspec)
222 lmax=lmax,
235 lmax=lmax,
402 if lmax is None or lmax < 0:
419 if lmax is None or lmax < 0:
523 if lmax is None or lmax < 0:
529 lmax=lmax,
903 lmax=lmax,
[all …]
|
| /dports/science/libint2/libint-2.7.1/include/libint2/ |
| H A D | cgshellinfo.h | 64 static void compute(int (&cartindex)[lmax+1][lmax+1][lmax+1]);
66 template <unsigned int lmax> struct CGShellOrderingGenerator<CGShellOrdering_Standard,lmax> {
67 static void compute(int (&cartindex)[lmax+1][lmax+1][lmax+1]) {
81 template <unsigned int lmax> struct CGShellOrderingGenerator<CGShellOrdering_GAMESS,lmax> {
82 static void compute(int (&cartindex)[lmax+1][lmax+1][lmax+1]) {
158 template <unsigned int lmax> struct CGShellOrderingGenerator<CGShellOrdering_ORCA,lmax> {
159 static void compute(int (&cartindex)[lmax+1][lmax+1][lmax+1]) {
378 template <unsigned int lmax> struct CGShellOrderingGenerator<CGShellOrdering_MOLDEN,lmax> {
379 static void compute(int (&cartindex)[lmax+1][lmax+1][lmax+1]) {
459 int cartindex[lmax+1][lmax+1][lmax+1];
[all …]
|
| /dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/NWints/auxil/ |
| H A D | spherical.F | 7 integer lmax variable
10 lmax = 8
11 memused = (((lmax+1)*(lmax+2))/2)*(2*lmax+1)*(lmax+1)
44 integer lmax local
45 double precision d(((lmax+1)*(lmax+2))/2, -lmax:lmax, 0:lmax)
65 call dcopy((((lmax+1)*(lmax+2))/2)*(2*lmax+1)*(lmax+1),
77 do l = 1, lmax-1
171 if (((2*lmax+1)*(lmax+1)).gt.qn_sz) stop 'qn error'
177 integer lmax local
178 double precision d(((lmax+1)*(lmax+2))/2, -lmax:lmax, 0:lmax)
[all …]
|
| /dports/science/nwchem-data/nwchem-7.0.2-release/src/NWints/auxil/ |
| H A D | spherical.F | 7 integer lmax variable
10 lmax = 8
11 memused = (((lmax+1)*(lmax+2))/2)*(2*lmax+1)*(lmax+1)
44 integer lmax local
45 double precision d(((lmax+1)*(lmax+2))/2, -lmax:lmax, 0:lmax)
65 call dcopy((((lmax+1)*(lmax+2))/2)*(2*lmax+1)*(lmax+1),
77 do l = 1, lmax-1
171 if (((2*lmax+1)*(lmax+1)).gt.qn_sz) stop 'qn error'
177 integer lmax local
178 double precision d(((lmax+1)*(lmax+2))/2, -lmax:lmax, 0:lmax)
[all …]
|
| /dports/lang/yorick/yorick-y_2_2_04/i/ |
| H A D | ylm.i | 98 y_lm = array(0., 1+lmax, 1+lmax, dimsof(z));
102 m1 = 1+lmax-l;
108 if (lmax) {
109 m = indgen(0:lmax)(,-:0:lmax) - indgen(0:lmax)(-:0:lmax,);
179 lmax = y; in ylm_lm()
183 lmax = d(2) - 1; in ylm_lm()
186 l = indgen(0:lmax)(,-:0:lmax); in ylm_lm()
286 x = transpose(ylm(lmax, x+0.*v, y, z), -1)(..,1:(1+lmax)^2,1);
287 alm = array(0., 1+lmax,1+lmax);
309 func legpol(lmax, z) in legpol() argument
[all …]
|
| /dports/science/nwchem-data/nwchem-7.0.2-release/src/NWints/hondo/ |
| H A D | hnd_init_1eelec.F | 18 integer lmax, iblk, i, mroots
26 lmax = -1
27 if (.not.bas_high_angular(basis,lmax))
28 & call errquit('int_init_1eelec: lmax error',lmax, INT_ERR)
29 lmax = lmax + 1
33 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4 * npt * max(nder*3,1)
37 mroots = (2*lmax+max(nder,0)-2)/2+1
38 maxscr = 3*lmax*lmax*mroots
43 iblk = lmax*lmax*mroots
|
| H A D | hnd_init_dip.F | 16 integer lmax, i
28 lmax = -1
29 if (.not.bas_high_angular(basis,lmax))
30 & call errquit('int_init_dso: lmax error',lmax, INT_ERR)
31 lmax = lmax + 1
35 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4
39 maxscr = 9*lmax*(lmax+1)
46 dip(i) = 1 + (i-1)*lmax*(lmax+1)
|
| H A D | hnd_init_dso.F | 17 integer lmax, iblk, i, mroots
30 lmax = -1
31 if (.not.bas_high_angular(basis,lmax))
32 & call errquit('int_init_dso: lmax error',lmax, INT_ERR)
33 lmax = lmax + 1
37 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4*npt*9
41 mroots = (2*lmax+2-2)/2+1
42 iblk = lmax*lmax*4*mroots
|
| H A D | hnd_init_pso.F | 17 integer lmax, iblk, i, mroots
29 lmax = -1
30 if (.not.bas_high_angular(basis,lmax))
31 & call errquit('int_init_pso: lmax error',lmax, INT_ERR)
32 lmax = lmax + 1
36 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4*npt*3
40 mroots = (2*lmax+2-2)/2+1
41 iblk = lmax*(lmax+1)*3*mroots
|
| /dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/NWints/hondo/ |
| H A D | hnd_init_1eelec.F | 18 integer lmax, iblk, i, mroots
26 lmax = -1
27 if (.not.bas_high_angular(basis,lmax))
28 & call errquit('int_init_1eelec: lmax error',lmax, INT_ERR)
29 lmax = lmax + 1
33 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4 * npt * max(nder*3,1)
37 mroots = (2*lmax+max(nder,0)-2)/2+1
38 maxscr = 3*lmax*lmax*mroots
43 iblk = lmax*lmax*mroots
|
| H A D | hnd_init_dip.F | 16 integer lmax, i
28 lmax = -1
29 if (.not.bas_high_angular(basis,lmax))
30 & call errquit('int_init_dso: lmax error',lmax, INT_ERR)
31 lmax = lmax + 1
35 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4
39 maxscr = 9*lmax*(lmax+1)
46 dip(i) = 1 + (i-1)*lmax*(lmax+1)
|
| H A D | hnd_init_dso.F | 17 integer lmax, iblk, i, mroots
30 lmax = -1
31 if (.not.bas_high_angular(basis,lmax))
32 & call errquit('int_init_dso: lmax error',lmax, INT_ERR)
33 lmax = lmax + 1
37 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4*npt*9
41 mroots = (2*lmax+2-2)/2+1
42 iblk = lmax*lmax*4*mroots
|
| H A D | hnd_init_pso.F | 17 integer lmax, iblk, i, mroots
29 lmax = -1
30 if (.not.bas_high_angular(basis,lmax))
31 & call errquit('int_init_pso: lmax error',lmax, INT_ERR)
32 lmax = lmax + 1
36 maxbuf = lmax*(lmax+1)*lmax*(lmax+1)/4*npt*3
40 mroots = (2*lmax+2-2)/2+1
41 iblk = lmax*(lmax+1)*3*mroots
|
| /dports/science/healpix/Healpix_3.50/src/cxx/Healpix_cxx/ |
| H A D | alm_fitsio.h | 49 void get_almsize(fitshandle &inp, int &lmax, int &mmax);
68 (fitshandle &inp, int lmax, int mmax) in read_Alm_from_fits() argument
71 read_Alm_from_fits (inp, res, lmax, mmax); in read_Alm_from_fits()
81 int lmax, int mmax, int hdunum=2);
97 int lmax, int mmax, int firsthdu=2)
110 int lmax, int mmax, PDT datatype);
114 int lmax, int mmax, PDT datatype) in write_Alm_to_fits() argument
118 write_Alm_to_fits (out, alms, lmax, mmax, datatype); in write_Alm_to_fits()
124 int lmax, int mmax, PDT datatype) in write_Alm_to_fits() argument
128 write_Alm_to_fits (out, almT, lmax, mmax, datatype); in write_Alm_to_fits()
[all …]
|
| H A D | alm_fitsio.cc | 46 inp.get_key ("MAX-LPOL",lmax); in get_almsize()
53 lmax=mmax=-1; in get_almsize()
65 if (l>lmax) lmax=l; in get_almsize()
76 get_almsize (inp, lmax, mmax); in get_almsize()
84 lmax=mmax=0; in get_almsize_pol()
89 if (tlmax>lmax) lmax=tlmax; in get_almsize_pol()
101 alms.Set(lmax, mmax); in read_Alm_from_fits()
103 int max_index = lmax*lmax + lmax + mmax + 1; in read_Alm_from_fits()
122 if ((l<=lmax) && (m<=mmax)) in read_Alm_from_fits()
141 read_Alm_from_fits(inp,alms,lmax,mmax); in read_Alm_from_fits()
[all …]
|
| /dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/contrib/marat/spectrum-improved/ |
| H A D | main.f | 34 l1=lmax(1)
35 l2=lmax(1)
37 write(*,1000) lmax(i),f(i),sigma(i)
41 tmp = lmax(i)-5*sigma(i)
43 tmp = lmax(i)+5*sigma(i)
46 call plot(n,l1,l2,f,lmax,sigma,e)
51 subroutine integrate(lmax,sigma,sum) argument
53 double precision lmax,sigma,sum local
63 x = lmax-dl
64 f = exp(-((lmax-x)/sigma)**2)/x**2
[all …]
|
| /dports/science/nwchem-data/nwchem-7.0.2-release/contrib/marat/spectrum-improved/ |
| H A D | main.f | 34 l1=lmax(1)
35 l2=lmax(1)
37 write(*,1000) lmax(i),f(i),sigma(i)
41 tmp = lmax(i)-5*sigma(i)
43 tmp = lmax(i)+5*sigma(i)
46 call plot(n,l1,l2,f,lmax,sigma,e)
51 subroutine integrate(lmax,sigma,sum)
53 double precision lmax,sigma,sum
63 x = lmax-dl
64 f = exp(-((lmax-x)/sigma)**2)/x**2
[all …]
|
| /dports/science/py-abipy/abipy-0.9.0/abipy/examples/plot/ |
| H A D | plot_efatbands.py | 29 lmax = 2 variable
32 fbnc_kpath.plot_fatbands_typeview(lmax=lmax, tight_layout=True)
35 fbnc_kpath.plot_fatbands_lview(lmax=lmax, tight_layout=True)
44 fbnc_kmesh.plot_pjdos_typeview(lmax=lmax, tight_layout=True)
47 fbnc_kmesh.plot_pjdos_lview(lmax=lmax, tight_layout=True)
52 fbnc_kpath.plot_fatbands_with_pjdos(pjdosfile=fbnc_kmesh, lmax=lmax,
56 fbnc_kpath.plot_fatbands_with_pjdos(pjdosfile=fbnc_kmesh, lmax=lmax,
|