(git:374b731)
Loading...
Searching...
No Matches
lazy.F
Go to the documentation of this file.
1!--------------------------------------------------------------------------------------------------!
2! CP2K: A general program to perform molecular dynamics simulations !
3! Copyright 2000-2024 CP2K developers group <https://cp2k.org> !
4! !
5! SPDX-License-Identifier: GPL-2.0-or-later !
6!--------------------------------------------------------------------------------------------------!
7
8! **************************************************************************************************
9!> \brief Filters for interpolating scaling functions .
10!> \author Florian Schiffmann (09.2007,fschiff)
11! **************************************************************************************************
12MODULE lazy
13
14 USE kinds, ONLY: dp
15#include "../base/base_uses.f90"
16
17 IMPLICIT NONE
18
19 PRIVATE
20
21 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'lazy'
22
23! *** Public data types ***
24
25 PUBLIC :: lazy_arrays
26
27CONTAINS
28
29! **************************************************************************************************
30!> \brief ...
31!> \param itype ...
32!> \param m ...
33!> \param ch ...
34!> \param cg ...
35!> \param cgt ...
36!> \param cht ...
37!> \par History
38!> 09.2007 created [Florian Schiffmann]
39!> \author fschiff
40! **************************************************************************************************
41 SUBROUTINE lazy_arrays(itype, m, ch, cg, cgt, cht)
42 INTEGER, INTENT(in) :: itype, m
43 REAL(kind=dp), DIMENSION(:), POINTER :: ch, cg, cgt, cht
44
45 INTEGER :: i
46
47 ALLOCATE (ch(-m:m))
48 ALLOCATE (cg(-m:m))
49 ALLOCATE (cgt(-m:m))
50 ALLOCATE (cht(-m:m))
51
52 SELECT CASE (itype)
53 CASE (8)
54
55 DO i = -m, m
56 ch(i) = 0._dp
57 cht(i) = 0._dp
58 cg(i) = 0._dp
59 cgt(i) = 0._dp
60 END DO
61
62 ! The normalization is chosen such that a constant function remains the same constant
63 ! on each level of the transform
64
65 ch(-7) = -5._dp/2048._dp
66 ch(-6) = 0._dp
67 ch(-5) = 49._dp/2048._dp
68 ch(-4) = 0._dp
69 ch(-3) = -245._dp/2048._dp
70 ch(-2) = 0._dp
71 ch(-1) = 1225._dp/2048._dp
72 ch(0) = 1._dp
73 ch(1) = 1225._dp/2048._dp
74 ch(2) = 0._dp
75 ch(3) = -245._dp/2048._dp
76 ch(4) = 0._dp
77 ch(5) = 49._dp/2048._dp
78 ch(6) = 0._dp
79 ch(7) = -5._dp/2048._dp
80 !
81 cht(0) = 1._dp
82
83 ! g coefficients from h coefficients
84 DO i = -m, m - 1
85 cg(i + 1) = cht(-i)*(-1)**(i + 1)
86 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
87 END DO
88
89 CASE (14)
90 ch = (/ &
91 0._dp, 0._dp, 0._dp, 0.0000275373458862304687_dp, 0._dp, -0.000423073768615722656_dp, 0._dp, &
92 0.00310254096984863281_dp, 0._dp, -0.0146262645721435547_dp, 0._dp, &
93 0.0511919260025024414_dp, 0._dp, -0.153575778007507324_dp, 0._dp, 0.614303112030029297_dp, &
94 1._dp, 0.614303112030029297_dp, 0._dp, -0.153575778007507324_dp, 0._dp, &
95 0.0511919260025024414_dp, 0._dp, -0.0146262645721435547_dp, 0._dp, &
96 0.00310254096984863281_dp, 0._dp, -0.000423073768615722656_dp, 0._dp, &
97 0.0000275373458862304687_dp, 0._dp, 0._dp, 0._dp &
98 /)
99
100 DO i = -m, m
101 cht(i) = 0._dp
102 cg(i) = 0._dp
103 cgt(i) = 0._dp
104 END DO
105
106 ! the normalization is chosen such that a constant function remains the same constant
107 ! on each level of the transform
108
109 cht(0) = 1._dp
110
111 ! g coefficients from h coefficients
112 DO i = -m, m - 1
113 cg(i + 1) = cht(-i)*(-1)**(i + 1)
114 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
115 END DO
116
117 CASE (16)
118 ch = (/ &
119 0._dp, 0._dp, 0._dp, -6.39259815216064453e-6_dp, 0._dp, 0.000110641121864318848_dp, 0._dp, &
120 -0.000915303826332092285_dp, 0._dp, 0.00484772026538848877_dp, 0._dp, &
121 -0.0186983495950698853_dp, 0._dp, 0.0575909167528152466_dp, 0._dp, &
122 -0.159974768757820129_dp, 0._dp, 0.617045536637306213_dp, 1._dp, 0.617045536637306213_dp, &
123 0._dp, -0.159974768757820129_dp, 0._dp, 0.0575909167528152466_dp, 0._dp, &
124 -0.0186983495950698853_dp, 0._dp, 0.00484772026538848877_dp, 0._dp, &
125 -0.000915303826332092285_dp, 0._dp, 0.000110641121864318848_dp, 0._dp, &
126 -6.39259815216064453e-6_dp, 0._dp, 0._dp, 0._dp &
127 /)
128
129 DO i = -m, m
130 cht(i) = 0._dp
131 cg(i) = 0._dp
132 cgt(i) = 0._dp
133 END DO
134
135 ! the normalization is chosen such that a constant function remains the same constant
136 ! on each level of the transform
137
138 cht(0) = 1._dp
139
140 ! g coefficients from h coefficients
141 DO i = -m, m - 1
142 cg(i + 1) = cht(-i)*(-1)**(i + 1)
143 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
144 END DO
145
146 CASE (20)
147 ch = (/ &
148 0._dp, 0._dp, 0._dp, &
149 -3.53757059201598167e-7_dp, 0._dp, 7.51213519833981991e-6_dp, 0._dp, &
150 -0.0000766237790230661631_dp, 0._dp, 0.000501001632073894143_dp, 0._dp, &
151 -0.00236837135162204504_dp, 0._dp, 0.00868402828928083181_dp, 0._dp, &
152 -0.0260520848678424954_dp, 0._dp, 0.0677354206563904881_dp, 0._dp, &
153 -0.16933855164097622_dp, 0._dp, 0.620908022683579475_dp, 1._dp, 0.620908022683579475_dp, &
154 0._dp, -0.16933855164097622_dp, 0._dp, 0.0677354206563904881_dp, 0._dp, &
155 -0.0260520848678424954_dp, 0._dp, 0.00868402828928083181_dp, 0._dp, &
156 -0.00236837135162204504_dp, 0._dp, 0.000501001632073894143_dp, 0._dp, &
157 -0.0000766237790230661631_dp, 0._dp, 7.51213519833981991e-6_dp, 0._dp, &
158 -3.53757059201598167e-7_dp, &
159 0._dp, 0._dp, 0._dp/)
160
161 !******** coefficients for wavelet transform *********************
162 DO i = -m, m
163 cht(i) = 0._dp
164 cg(i) = 0._dp
165 cgt(i) = 0._dp
166 END DO
167
168 ! the normalization is chosen such that a constant function remains the same constant
169 ! on each level of the transform
170
171 cht(0) = 1._dp
172
173 ! g coefficients from h coefficients
174 DO i = -m, m - 1
175 cg(i + 1) = cht(-i)*(-1)**(i + 1)
176 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
177 END DO
178
179 CASE (24)
180 ch = (/ &
181 0._dp, 0._dp, 0._dp, &
182 -2.00495833269087598e-8_dp, 0._dp, 5.0505855142546352e-7_dp, 0._dp, &
183 -6.14044870417274069e-6_dp, 0._dp, 0.0000480399810385279125_dp, 0._dp, &
184 -0.000272226559218324837_dp, 0._dp, 0.00119360875964957813_dp, 0._dp, &
185 -0.00423188560239395883_dp, 0._dp, 0.0125613112325027032_dp, 0._dp, &
186 -0.032300514597864094_dp, 0._dp, 0.0753678673950162192_dp, 0._dp, &
187 -0.175858357255037845_dp, 0._dp, 0.623497812086043268_dp, 1._dp, 0.623497812086043268_dp, &
188 0._dp, -0.175858357255037845_dp, 0._dp, 0.0753678673950162192_dp, 0._dp, &
189 -0.032300514597864094_dp, 0._dp, 0.0125613112325027032_dp, 0._dp, &
190 -0.00423188560239395883_dp, 0._dp, 0.00119360875964957813_dp, 0._dp, &
191 -0.000272226559218324837_dp, 0._dp, 0.0000480399810385279125_dp, 0._dp, &
192 -6.14044870417274069e-6_dp, 0._dp, 5.0505855142546352e-7_dp, 0._dp, &
193 -2.00495833269087598e-8_dp, &
194 0._dp, 0._dp, 0._dp/)
195
196 !******** coefficients for wavelet transform *********************
197 DO i = -m, m
198 cht(i) = 0._dp
199 cg(i) = 0._dp
200 cgt(i) = 0._dp
201 END DO
202
203 ! the normalization is chosen such that a constant function remains the same constant
204 ! on each level of the transform
205
206 cht(0) = 1._dp
207
208 ! g coefficients from h coefficients
209 DO i = -m, m - 1
210 cg(i + 1) = cht(-i)*(-1)**(i + 1)
211 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
212 END DO
213
214 CASE (30)
215 ch = (/ &
216 0._dp, 0._dp, 0._dp, &
217 2.78364831185484718e-10_dp, 0._dp, -8.6705490009997277e-9_dp, 0._dp, &
218 1.31098700895115883e-7_dp, 0._dp, -1.28248729136526407e-6_dp, 0._dp, &
219 9.13008809805271326e-6_dp, 0._dp, -0.000050455750015554468_dp, 0._dp, &
220 0.000225566882422478798_dp, 0._dp, -0.00083996810502084962_dp, 0._dp, &
221 0.0026652834101623113_dp, 0._dp, -0.00734972091893243418_dp, 0._dp, &
222 0.0179659844685015058_dp, 0._dp, -0.039898484988490357_dp, 0._dp, &
223 0.0837868184758297496_dp, 0._dp, -0.18261229667809048_dp, 0._dp, 0.626099302896310217_dp, &
224 1._dp, 0.626099302896310217_dp, 0._dp, -0.18261229667809048_dp, 0._dp, &
225 0.0837868184758297496_dp, 0._dp, -0.039898484988490357_dp, 0._dp, &
226 0.0179659844685015058_dp, 0._dp, -0.00734972091893243418_dp, 0._dp, &
227 0.0026652834101623113_dp, 0._dp, -0.00083996810502084962_dp, 0._dp, &
228 0.000225566882422478798_dp, 0._dp, -0.000050455750015554468_dp, 0._dp, &
229 9.13008809805271326e-6_dp, 0._dp, -1.28248729136526407e-6_dp, 0._dp, &
230 1.31098700895115883e-7_dp, 0._dp, -8.6705490009997277e-9_dp, 0._dp, &
231 2.78364831185484718e-10_dp, &
232 0._dp, 0._dp, 0._dp/)
233
234 !******** coefficients for wavelet transform *********************
235 DO i = -m, m
236 cht(i) = 0._dp
237 cg(i) = 0._dp
238 cgt(i) = 0._dp
239 END DO
240
241 ! the normalization is chosen such that a constant function remains the same constant
242 ! on each level of the transform
243
244 cht(0) = 1._dp
245
246 ! g coefficients from h coefficients
247 DO i = -m, m - 1
248 cg(i + 1) = cht(-i)*(-1)**(i + 1)
249 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
250 END DO
251
252 CASE (40)
253 ch = (/ &
254 0._dp, 0._dp, 0._dp, &
255 -2.33895335687458837e-13_dp, 0._dp, 9.61499474542229436e-12_dp, 0._dp, &
256 -1.9312403731519637e-10_dp, 0._dp, 2.52621846791090201e-9_dp, 0._dp, &
257 -2.42028027409528354e-8_dp, 0._dp, 1.81103730854716044e-7_dp, 0._dp, &
258 -1.10227332483179024e-6_dp, 0._dp, 5.61214589957214346e-6_dp, 0._dp, &
259 -0.0000244006343459658411_dp, 0._dp, 0.0000920510703104425647_dp, 0._dp, &
260 -0.000305221969976730609_dp, 0._dp, 0.000899343879450152758_dp, 0._dp, &
261 -0.00237826492565707063_dp, 0._dp, 0.00569939227746221067_dp, 0._dp, &
262 -0.0125090557778066702_dp, 0._dp, 0.0254814099177543281_dp, 0._dp, &
263 -0.0491427191270976329_dp, 0._dp, 0.0930820915230908105_dp, 0._dp, &
264 -0.189611667917407207_dp, 0._dp, 0.628712372568244948_dp, 1._dp, 0.628712372568244948_dp, &
265 0._dp, -0.189611667917407207_dp, 0._dp, 0.0930820915230908105_dp, 0._dp, &
266 -0.0491427191270976329_dp, 0._dp, 0.0254814099177543281_dp, 0._dp, &
267 -0.0125090557778066702_dp, 0._dp, 0.00569939227746221067_dp, 0._dp, &
268 -0.00237826492565707063_dp, 0._dp, 0.000899343879450152758_dp, 0._dp, &
269 -0.000305221969976730609_dp, 0._dp, 0.0000920510703104425647_dp, 0._dp, &
270 -0.0000244006343459658411_dp, 0._dp, 5.61214589957214346e-6_dp, 0._dp, &
271 -1.10227332483179024e-6_dp, 0._dp, 1.81103730854716044e-7_dp, 0._dp, &
272 -2.42028027409528354e-8_dp, 0._dp, 2.52621846791090201e-9_dp, 0._dp, &
273 -1.9312403731519637e-10_dp, 0._dp, 9.61499474542229436e-12_dp, 0._dp, &
274 -2.33895335687458837e-13_dp, &
275 0._dp, 0._dp, 0._dp/)
276
277 !******** coefficients for wavelet transform *********************
278 DO i = -m, m
279 cht(i) = 0._dp
280 cg(i) = 0._dp
281 cgt(i) = 0._dp
282 END DO
283
284 ! the normalization is chosen such that a constant function remains the same constant
285 ! on each level of the transform
286
287 cht(0) = 1._dp
288
289 ! g coefficients from h coefficients
290 DO i = -m, m - 1
291 cg(i + 1) = cht(-i)*(-1)**(i + 1)
292 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
293 END DO
294
295 CASE (50)
296 !******** coefficients for wavelet transform *********************
297
298 ch(-m:0) = (/0._dp, 0._dp, 0._dp, &
299 2.03510990661269603e-16_dp, 0._dp, -1.03963806080363471e-14_dp, 0._dp, &
300 2.60602607241444435e-13_dp, 0._dp, -4.27267065360972853e-12_dp, 0._dp, &
301 5.15325765417075794e-11_dp, 0._dp, -4.87577454971540944e-10_dp, 0._dp, &
302 3.76884194923947865e-9_dp, 0._dp, -2.44743981275102471e-8_dp, 0._dp, &
303 1.36277898664545694e-7_dp, 0._dp, -6.60874540835592558e-7_dp, 0._dp, &
304 2.82580838150391301e-6_dp, 0._dp, -0.0000107609066649189415_dp, 0._dp, &
305 0.0000368023007940227798_dp, 0._dp, -0.000113853271687194218_dp, 0._dp, &
306 0.000320647989649649023_dp, 0._dp, -0.00082693428909646327_dp, 0._dp, &
307 0.00196396893660410027_dp, 0._dp, -0.00432073166052902059_dp, 0._dp, &
308 0.008863039303649273_dp, 0._dp, -0.0170899753079935742_dp, 0._dp, &
309 0.0313316213979882195_dp, 0._dp, -0.0556296135025505121_dp, 0._dp, &
310 0.0991218567863627306_dp, 0._dp, -0.193934067625492299_dp, 0._dp, &
311 0.630285719782849972_dp, 1._dp/)
312 ch(1:m) = (/ &
313 0.630285719782849972_dp, 0._dp, -0.193934067625492299_dp, 0._dp, &
314 0.0991218567863627306_dp, 0._dp, -0.0556296135025505121_dp, 0._dp, &
315 0.0313316213979882195_dp, 0._dp, -0.0170899753079935742_dp, 0._dp, &
316 0.008863039303649273_dp, 0._dp, -0.00432073166052902059_dp, 0._dp, &
317 0.00196396893660410027_dp, 0._dp, -0.00082693428909646327_dp, 0._dp, &
318 0.000320647989649649023_dp, 0._dp, -0.000113853271687194218_dp, 0._dp, &
319 0.0000368023007940227798_dp, 0._dp, -0.0000107609066649189415_dp, 0._dp, &
320 2.82580838150391301e-6_dp, 0._dp, -6.60874540835592558e-7_dp, 0._dp, &
321 1.36277898664545694e-7_dp, 0._dp, -2.44743981275102471e-8_dp, 0._dp, &
322 3.76884194923947865e-9_dp, 0._dp, -4.87577454971540944e-10_dp, 0._dp, &
323 5.15325765417075794e-11_dp, 0._dp, -4.27267065360972853e-12_dp, 0._dp, &
324 2.60602607241444435e-13_dp, 0._dp, -1.03963806080363471e-14_dp, 0._dp, &
325 2.03510990661269603e-16_dp, &
326 0._dp, 0._dp, 0._dp/)
327 DO i = -m, m
328 cht(i) = 0._dp
329 cg(i) = 0._dp
330 cgt(i) = 0._dp
331 END DO
332
333 ! the normalization is chosen such that a constant function remains the same constant
334 ! on each level of the transform
335
336 cht(0) = 1._dp
337
338 ! g coefficients from h coefficients
339 DO i = -m, m - 1
340 cg(i + 1) = cht(-i)*(-1)**(i + 1)
341 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
342 END DO
343
344 CASE (60)
345 !******** coefficients for wavelet transform *********************
346 ch(-m:0) = (/ &
347 0._dp, 0._dp, 0._dp, &
348 -1.80960777804178984e-19_dp, 0._dp, 1.10513064480060885e-17_dp, 0._dp, &
349 -3.32141991973710259e-16_dp, 0._dp, 6.54883738891560793e-15_dp, 0._dp, &
350 -9.5279163579909433e-14_dp, 0._dp, 1.09084919935365698e-12_dp, 0._dp, &
351 -1.02354148279779304e-11_dp, 0._dp, 8.09410105920413478e-11_dp, 0._dp, &
352 -5.50587106934234749e-10_dp, 0._dp, 3.27218841600752522e-9_dp, 0._dp, &
353 -1.71999647508087864e-8_dp, 0._dp, 8.07595396530604197e-8_dp, 0._dp, &
354 -3.41497481961512632e-7_dp, 0._dp, 1.30947169656570694e-6_dp, 0._dp, &
355 -4.58013372208926991e-6_dp, 0._dp, 0.0000146880150398035208_dp, 0._dp, &
356 -0.0000433840444231233622_dp, 0._dp, 0.000118515001353520526_dp, 0._dp, &
357 -0.000300581525171972349_dp, 0._dp, 0.000710396937937769235_dp, 0._dp, &
358 -0.00157035112596770041_dp, 0._dp, 0.00325946830347917649_dp, 0._dp, &
359 -0.00638065613347741822_dp, 0._dp, 0.0118436928230099234_dp, 0._dp, &
360 -0.0209956372771539552_dp, 0._dp, 0.0359258682297967678_dp, 0._dp, &
361 -0.0604028333973506096_dp, 0._dp, 0.103355959368799932_dp, 0._dp, &
362 -0.196868494035809394_dp, 0._dp, 0.631336894666561161_dp, 1._dp/)
363 ch(1:m) = (/0.631336894666561161_dp, &
364 0._dp, -0.196868494035809394_dp, 0._dp, 0.103355959368799932_dp, 0._dp, &
365 -0.0604028333973506096_dp, 0._dp, 0.0359258682297967678_dp, 0._dp, &
366 -0.0209956372771539552_dp, 0._dp, 0.0118436928230099234_dp, 0._dp, &
367 -0.00638065613347741822_dp, 0._dp, 0.00325946830347917649_dp, 0._dp, &
368 -0.00157035112596770041_dp, 0._dp, 0.000710396937937769235_dp, 0._dp, &
369 -0.000300581525171972349_dp, 0._dp, 0.000118515001353520526_dp, 0._dp, &
370 -0.0000433840444231233622_dp, 0._dp, 0.0000146880150398035208_dp, 0._dp, &
371 -4.58013372208926991e-6_dp, 0._dp, 1.30947169656570694e-6_dp, 0._dp, &
372 -3.41497481961512632e-7_dp, 0._dp, 8.07595396530604197e-8_dp, 0._dp, &
373 -1.71999647508087864e-8_dp, 0._dp, 3.27218841600752522e-9_dp, 0._dp, &
374 -5.50587106934234749e-10_dp, 0._dp, 8.09410105920413478e-11_dp, 0._dp, &
375 -1.02354148279779304e-11_dp, 0._dp, 1.09084919935365698e-12_dp, 0._dp, &
376 -9.5279163579909433e-14_dp, 0._dp, 6.54883738891560793e-15_dp, 0._dp, &
377 -3.32141991973710259e-16_dp, 0._dp, 1.10513064480060885e-17_dp, 0._dp, &
378 -1.80960777804178984e-19_dp, &
379 0._dp, 0._dp, 0._dp/)
380
381 DO i = -m, m
382 cht(i) = 0._dp
383 cg(i) = 0._dp
384 cgt(i) = 0._dp
385 END DO
386
387 ! the normalization is chosen such that a constant function remains the same constant
388 ! on each level of the transform
389
390 cht(0) = 1._dp
391
392 ! g coefficients from h coefficients
393 DO i = -m, m - 1
394 cg(i + 1) = cht(-i)*(-1)**(i + 1)
395 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
396 END DO
397
398 CASE (100)
399 !******** coefficients for wavelet transform *********************
400 ch(-m:-m + 22) = (/ &
401 0._dp, 0._dp, 0._dp, &
402 -1.2683805626484815e-31_dp, 0._dp, 1.28158741180595538e-29_dp, 0._dp, &
403 -6.41198417717232202e-28_dp, 0._dp, 2.11779334023809847e-26_dp, 0._dp, &
404 -5.19441179715542395e-25_dp, 0._dp, 1.00911662890805933e-23_dp, 0._dp, &
405 -1.61729304931893186e-22_dp, 0._dp, 2.19924673311587861e-21_dp, 0._dp, &
406 -2.59007672484430282e-20_dp, 0._dp, 2.68351845030850743e-19_dp, 0._dp/)
407 ch(-m + 23:0) = (/ &
408 -2.47631006363911635e-18_dp, 0._dp, 2.05560047903738218e-17_dp, 0._dp, &
409 -1.54763876066192241e-16_dp, 0._dp, 1.06410357358472536e-15_dp, 0._dp, &
410 -6.72076683598079057e-15_dp, 0._dp, 3.91882394735691508e-14_dp, 0._dp, &
411 -2.1187969775075261e-13_dp, 0._dp, 1.06630135673478305e-12_dp, 0._dp, &
412 -5.01180443685748998e-12_dp, 0._dp, 2.20666419683207691e-11_dp, 0._dp, &
413 -9.1258654919835045e-11_dp, 0._dp, 3.55352205999457932e-10_dp, 0._dp, &
414 -1.30569909741619005e-9_dp, 0._dp, 4.53620646231137398e-9_dp, 0._dp, &
415 -1.49279735540769725e-8_dp, 0._dp, 4.66118357913015673e-8_dp, 0._dp, &
416 -1.38309751701521672e-7_dp, 0._dp, 3.90568525175243503e-7_dp, 0._dp, &
417 -1.05103157937523999e-6_dp, 0._dp, 2.69873794098621337e-6_dp, 0._dp, &
418 -6.61998110310293365e-6_dp, 0._dp, 0.0000155312721433740143_dp, 0._dp, &
419 -0.0000348898934935080536_dp, 0._dp, 0.0000751302114254144957_dp, 0._dp, &
420 -0.000155250095336387829_dp, 0._dp, 0.000308205854288494064_dp, 0._dp, &
421 -0.000588508297900663563_dp, 0._dp, 0.00108221904295030132_dp, 0._dp, &
422 -0.00191926718143360074_dp, 0._dp, 0.00328782888345707184_dp, 0._dp, &
423 -0.00545087420152093489_dp, 0._dp, 0.00876675750345045625_dp, 0._dp, &
424 -0.0137206712673049998_dp, 0._dp, 0.0209859998811194541_dp, 0._dp, &
425 -0.031565718829452402_dp, 0._dp, 0.0471537281279474153_dp, 0._dp, &
426 -0.0711699126403181485_dp, 0._dp, 0.112357606764076737_dp, 0._dp, &
427 -0.202867901101805219_dp, 0._dp, 0.633444670787269357_dp, 1._dp/)
428 ch(1:80) = (/0.633444670787269357_dp, &
429 0._dp, -0.202867901101805219_dp, 0._dp, 0.112357606764076737_dp, 0._dp, &
430 -0.0711699126403181485_dp, 0._dp, 0.0471537281279474153_dp, 0._dp, &
431 -0.031565718829452402_dp, 0._dp, 0.0209859998811194541_dp, 0._dp, &
432 -0.0137206712673049998_dp, 0._dp, 0.00876675750345045625_dp, 0._dp, &
433 -0.00545087420152093489_dp, 0._dp, 0.00328782888345707184_dp, 0._dp, &
434 -0.00191926718143360074_dp, 0._dp, 0.00108221904295030132_dp, 0._dp, &
435 -0.000588508297900663563_dp, 0._dp, 0.000308205854288494064_dp, 0._dp, &
436 -0.000155250095336387829_dp, 0._dp, 0.0000751302114254144957_dp, 0._dp, &
437 -0.0000348898934935080536_dp, 0._dp, 0.0000155312721433740143_dp, 0._dp, &
438 -6.61998110310293365e-6_dp, 0._dp, 2.69873794098621337e-6_dp, 0._dp, &
439 -1.05103157937523999e-6_dp, 0._dp, 3.90568525175243503e-7_dp, 0._dp, &
440 -1.38309751701521672e-7_dp, 0._dp, 4.66118357913015673e-8_dp, 0._dp, &
441 -1.49279735540769725e-8_dp, 0._dp, 4.53620646231137398e-9_dp, 0._dp, &
442 -1.30569909741619005e-9_dp, 0._dp, 3.55352205999457932e-10_dp, 0._dp, &
443 -9.1258654919835045e-11_dp, 0._dp, 2.20666419683207691e-11_dp, 0._dp, &
444 -5.01180443685748998e-12_dp, 0._dp, 1.06630135673478305e-12_dp, 0._dp, &
445 -2.1187969775075261e-13_dp, 0._dp, 3.91882394735691508e-14_dp, 0._dp, &
446 -6.72076683598079057e-15_dp, 0._dp, 1.06410357358472536e-15_dp, 0._dp, &
447 -1.54763876066192241e-16_dp, 0._dp, 2.05560047903738218e-17_dp, 0._dp, &
448 -2.47631006363911635e-18_dp, 0._dp/)
449 ch(81:m) = (/2.68351845030850743e-19_dp, 0._dp, &
450 -2.59007672484430282e-20_dp, 0._dp, 2.19924673311587861e-21_dp, 0._dp, &
451 -1.61729304931893186e-22_dp, 0._dp, 1.00911662890805933e-23_dp, 0._dp, &
452 -5.19441179715542395e-25_dp, 0._dp, 2.11779334023809847e-26_dp, 0._dp, &
453 -6.41198417717232202e-28_dp, 0._dp, 1.28158741180595538e-29_dp, 0._dp, &
454 -1.2683805626484815e-31_dp, &
455 0._dp, 0._dp, 0._dp/)
456
457 DO i = -m, m
458 cht(i) = 0._dp
459 cg(i) = 0._dp
460 cgt(i) = 0._dp
461 END DO
462
463 ! the normalization is chosen such that a constant function remains the same constant
464 ! on each level of the transform
465
466 cht(0) = 1._dp
467
468 ! g coefficients from h coefficients
469 DO i = -m, m - 1
470 cg(i + 1) = cht(-i)*(-1)**(i + 1)
471 cgt(i + 1) = ch(-i)*(-1)**(i + 1)
472 END DO
473
474 END SELECT
475
476 END SUBROUTINE lazy_arrays
477END MODULE lazy
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
lazy
Filters for interpolating scaling functions .
Definition lazy.F:12
lazy::lazy_arrays
subroutine, public lazy_arrays(itype, m, ch, cg, cgt, cht)
...
Definition lazy.F:42