Line data Source code
1 : !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2 : !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
3 : !!!!
4 : !!!! MIT License
5 : !!!!
6 : !!!! ParaMonte: plain powerful parallel Monte Carlo library.
7 : !!!!
8 : !!!! Copyright (C) 2012-present, The Computational Data Science Lab
9 : !!!!
10 : !!!! This file is part of the ParaMonte library.
11 : !!!!
12 : !!!! Permission is hereby granted, free of charge, to any person obtaining a
13 : !!!! copy of this software and associated documentation files (the "Software"),
14 : !!!! to deal in the Software without restriction, including without limitation
15 : !!!! the rights to use, copy, modify, merge, publish, distribute, sublicense,
16 : !!!! and/or sell copies of the Software, and to permit persons to whom the
17 : !!!! Software is furnished to do so, subject to the following conditions:
18 : !!!!
19 : !!!! The above copyright notice and this permission notice shall be
20 : !!!! included in all copies or substantial portions of the Software.
21 : !!!!
22 : !!!! THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
23 : !!!! EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
24 : !!!! MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
25 : !!!! IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM,
26 : !!!! DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR
27 : !!!! OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE
28 : !!!! OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
29 : !!!!
30 : !!!! ACKNOWLEDGMENT
31 : !!!!
32 : !!!! ParaMonte is an honor-ware and its currency is acknowledgment and citations.
33 : !!!! As per the ParaMonte library license agreement terms, if you use any parts of
34 : !!!! this library for any purposes, kindly acknowledge the use of ParaMonte in your
35 : !!!! work (education/research/industry/development/...) by citing the ParaMonte
36 : !!!! library as described on this page:
37 : !!!!
38 : !!!! https://github.com/cdslaborg/paramonte/blob/main/ACKNOWLEDGMENT.md
39 : !!!!
40 : !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
41 : !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
42 :
43 : !> \brief This include file contains the body of the submodules
44 : !> [ParaDRAM_mod@Test_runSampler_smod](@ref paradram_mod@test_runsampler_smod) and
45 : !> [ParaDISE_mod@@Test_runSampler_smod](@ref paradise_mod@test_runsampler_smod).
46 : !> \author Amir Shahmoradi
47 :
48 : !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
49 :
50 : contains
51 :
52 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
53 :
54 : !> \brief
55 : !> Test the ParaDXXX sampler with no input arguments or input file.
56 4 : module function test_runSampler_1() result(assertion)
57 : implicit none
58 : logical :: assertion
59 2 : type(ParaDXXX_type) :: PD
60 2 : assertion = .true.
61 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
62 : call PD%runSampler ( ndim = 1_IK &
63 : , getLogFunc = getLogFuncMVN &
64 : , mpiFinalizeRequested = .false. &
65 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_1" &
66 2 : )
67 2 : assertion = assertion .and. .not. PD%Err%occurred
68 : #endif
69 30 : end function test_runSampler_1
70 :
71 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
72 :
73 : !> \brief
74 : !> Test the ParaDXXX sampler with an internal input file.
75 2 : module function test_runSampler_2() result(assertion)
76 2 : use String_mod, only: num2str
77 : implicit none
78 : logical :: assertion
79 2 : type(ParaDXXX_type) :: PD
80 : integer(IK) , parameter :: chainSize_ref = 100_IK
81 : integer(IK) , parameter :: userSeed_ref = 1111_IK
82 :
83 2 : assertion = .true.
84 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
85 : call PD%runSampler ( ndim = 1_IK &
86 : , getLogFunc = getLogFuncMVN &
87 : , mpiFinalizeRequested = .false. &
88 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_2" &
89 : , inputFile = ParaDXXX_NML//" randomSeed = "//num2str(userSeed_ref)//" chainSize = "//num2str(chainSize_ref)//" /" &
90 2 : )
91 2 : assertion = assertion .and. .not. PD%Err%occurred .and. PD%SpecBase%RandomSeed%userSeed==userSeed_ref .and. PD%SpecMCMC%ChainSize%val==chainSize_ref
92 : ! LCOV_EXCL_START
93 : if (Test%isDebugMode .and. .not. assertion) then
94 : write(Test%outputUnit,"(*(g0))")
95 : write(Test%outputUnit,"(*(g0))") "Image%id, PD%Err%occurred = ", PD%Image%id, PD%Err%occurred
96 : write(Test%outputUnit,"(*(g0))") "Image%id, PD%Err%msg = ", PD%Image%id, PD%Err%msg
97 : write(Test%outputUnit,"(*(g0))") "Image%id, userSeed_ref = ", PD%Image%id, userSeed_ref
98 : write(Test%outputUnit,"(*(g0))") "Image%id, PD%SpecBase%RandomSeed%userSeed = ", PD%Image%id, PD%SpecBase%RandomSeed%userSeed
99 : write(Test%outputUnit,"(*(g0))") "Image%id, chainSize_ref = ", PD%Image%id, chainSize_ref
100 : write(Test%outputUnit,"(*(g0))") "Image%id, PD%SpecMCMC%ChainSize%val = ", PD%Image%id, PD%SpecMCMC%ChainSize%val
101 : write(Test%outputUnit,"(*(g0))")
102 : end if
103 : ! LCOV_EXCL_STOP
104 : #endif
105 30 : end function test_runSampler_2
106 :
107 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
108 :
109 : !> \brief
110 : !> Test the ParaDXXX sampler with an external input file.
111 2 : module function test_runSampler_3() result(assertion)
112 : implicit none
113 : logical :: assertion
114 2 : type(ParaDXXX_type) :: PD
115 2 : assertion = .true.
116 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
117 : call PD%runSampler ( ndim = 1_IK &
118 : , getLogFunc = getLogFuncMVN &
119 : , mpiFinalizeRequested = .false. &
120 : , inputFile = Test%inDir//"/Test_ParaDXXX_mod@test_runSampler_3.in" &
121 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_3" &
122 2 : )
123 2 : assertion = assertion .and. .not. PD%Err%occurred
124 : #endif
125 32 : end function test_runSampler_3
126 :
127 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
128 :
129 : !> \brief
130 : !> Test the ParaDXXX sampler with a path to a non-existing external input file.
131 2 : module function test_runSampler_4() result(assertion)
132 : implicit none
133 : logical :: assertion
134 2 : type(ParaDXXX_type) :: PD
135 2 : assertion = .true.
136 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
137 : call PD%runSampler ( ndim = 1_IK &
138 : , getLogFunc = getLogFuncMVN &
139 : , mpiFinalizeRequested = .false. &
140 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_4" &
141 : , inputFile = " " &
142 2 : )
143 2 : assertion = assertion .and. .not. PD%Err%occurred
144 : #endif
145 32 : end function test_runSampler_4
146 :
147 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
148 :
149 : !> \brief
150 : !> Test the ParaDXXX sampler with wrong SpecBase input values.
151 2 : module function test_runSampler_5() result(assertion)
152 : implicit none
153 : logical :: assertion
154 2 : type(ParaDXXX_type) :: PD
155 2 : assertion = .true.
156 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
157 : call PD%runSampler ( ndim = 1_IK &
158 : , getLogFunc = getLogFuncMVN &
159 : , mpiFinalizeRequested = .false. &
160 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_5" &
161 : , inputFile = Test%inDir//"/Test_ParaDXXX_mod@test_runSampler_5.in" &
162 2 : )
163 2 : assertion = assertion .and. PD%Err%occurred
164 : #endif
165 4 : end function test_runSampler_5
166 :
167 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
168 :
169 : !> \brief
170 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file.
171 2 : module function test_runSampler_6() result(assertion)
172 : implicit none
173 : logical :: assertion
174 2 : type(ParaDXXX_type) :: PD1, PD2
175 2 : assertion = .true.
176 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
177 :
178 : ! Run the fresh simulation
179 :
180 : call PD1%runSampler ( ndim = 2_IK &
181 : , getLogFunc = getLogFuncMVN &
182 : , mpiFinalizeRequested = .false. &
183 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_6" &
184 : , outputRealPrecision = 16_IK &
185 : , randomSeed = 12345_IK &
186 : , sampleSize = 250_IK &
187 : , chainSize = 500_IK &
188 2 : )
189 :
190 2 : assertion = assertion .and. .not. PD1%Err%occurred
191 :
192 2 : if (.not. assertion) then
193 : ! LCOV_EXCL_START
194 : if (Test%isDebugMode) then
195 : write(Test%outputUnit,"(*(g0,:,' '))")
196 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
197 : write(Test%outputUnit,"(*(g0,:,' '))")
198 : end if
199 : return
200 : end if
201 : ! LCOV_EXCL_STOP
202 :
203 2 : call Test%Image%sync()
204 :
205 : block
206 : use System_mod, only: removeFile
207 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
208 : end block
209 :
210 2 : call Test%Image%sync()
211 :
212 : ! restart the simulation with the same configuration
213 :
214 : call PD2%runSampler ( ndim = 2_IK &
215 : , getLogFunc = getLogFuncMVN &
216 : , mpiFinalizeRequested = .false. &
217 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_6" &
218 : , outputRealPrecision = 16_IK &
219 : , randomSeed = 12345_IK &
220 : , sampleSize = 250_IK &
221 : , chainSize = 500_IK &
222 2 : )
223 :
224 2 : assertion = assertion .and. .not. PD2%Err%occurred
225 :
226 2 : if (.not. assertion) then
227 : ! LCOV_EXCL_START
228 : if (Test%isDebugMode) then
229 : write(Test%outputUnit,"(*(g0,:,' '))")
230 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
231 : write(Test%outputUnit,"(*(g0,:,' '))")
232 : end if
233 : return
234 : end if
235 : ! LCOV_EXCL_STOP
236 :
237 2 : if (PD2%Image%isLeader) then
238 : ! LCOV_EXCL_START
239 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-12_RK ) ! by default, the output precision is only 8 digits
240 : if (.not. assertion) then
241 : if (Test%isDebugMode) then
242 : write(Test%outputUnit,"(*(g0,:,' '))")
243 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
244 : write(Test%outputUnit,"(*(g0,:,' '))")
245 : end if
246 : return
247 : end if
248 : end if
249 : ! LCOV_EXCL_STOP
250 : #endif
251 60 : end function test_runSampler_6
252 :
253 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
254 :
255 : !> \brief
256 : !> Test the ParaDXXX sampler with a wrong internal input namelist group:
257 : !> + Infinity values for `domainLowerLimitVec` and `domainUpperLimitVec`.
258 2 : module function test_runSampler_7() result(assertion)
259 2 : use Constants_mod, only: IK, RK
260 : implicit none
261 : logical :: assertion
262 2 : type(ParaDXXX_type) :: PD
263 2 : assertion = .true.
264 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
265 : call PD%runSampler ( ndim = 1_IK &
266 : , getLogFunc = getLogFuncMVN &
267 : , mpiFinalizeRequested = .false. &
268 : , inputFile = "&ParaDXXX randomSeed = 1111 /" &
269 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"@SpecBase/test_runSampler_7" &
270 2 : )
271 2 : assertion = assertion .and. .not. PD%Err%occurred
272 : #endif
273 30 : end function test_runSampler_7
274 :
275 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
276 :
277 : !> \brief
278 : !> Test the ParaDXXX sampler with a wrong internal input namelist group:
279 : !> + Infinity values for `domainLowerLimitVec` and `domainUpperLimitVec`.
280 2 : module function test_runSampler_8() result(assertion)
281 2 : use Constants_mod, only: IK, RK
282 : implicit none
283 : logical :: assertion
284 2 : type(ParaDXXX_type) :: PD
285 2 : assertion = .true.
286 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
287 : call PD%runSampler ( ndim = 1_IK &
288 : , getLogFunc = getLogFuncMVN &
289 : , mpiFinalizeRequested = .false. &
290 : , inputFile = Test%inDir//"/Test_ParaDXXX_mod@test_runSampler_8.in" &
291 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"@SpecBase/test_runSampler_8" &
292 2 : )
293 2 : assertion = assertion .and. .not. PD%Err%occurred
294 : #endif
295 30 : end function test_runSampler_8
296 :
297 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
298 :
299 : !> \brief
300 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file.
301 2 : module function test_runSampler_9() result(assertion)
302 : implicit none
303 : logical :: assertion
304 2 : type(ParaDXXX_type) :: PD1, PD2
305 2 : assertion = .true.
306 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
307 :
308 : ! Run the fresh simulation as a reference run
309 :
310 : call PD1%runSampler ( ndim = 2_IK &
311 : , getLogFunc = getLogFuncMVN &
312 : , mpiFinalizeRequested = .false. &
313 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_9.ref" &
314 : , proposalModel = "uniform" &
315 : , randomSeed = 12345_IK &
316 : , sampleSize = 250_IK &
317 : , chainSize = 500_IK &
318 2 : )
319 2 : assertion = assertion .and. .not. PD1%Err%occurred
320 :
321 2 : if (.not. assertion) then
322 : ! LCOV_EXCL_START
323 : if (Test%isDebugMode) then
324 : write(Test%outputUnit,"(*(g0,:,' '))")
325 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
326 : write(Test%outputUnit,"(*(g0,:,' '))")
327 : end if
328 : return
329 : end if
330 : ! LCOV_EXCL_STOP
331 :
332 2 : call Test%Image%sync()
333 :
334 : ! Run the fresh simulation
335 :
336 : call PD1%runSampler ( ndim = 2_IK &
337 : , getLogFunc = getLogFuncMVN &
338 : , mpiFinalizeRequested = .false. &
339 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_9" &
340 : , proposalModel = "uniform" &
341 : , randomSeed = 12345_IK &
342 : , sampleSize = 250_IK &
343 : , chainSize = 500_IK &
344 2 : )
345 :
346 2 : assertion = assertion .and. .not. PD1%Err%occurred
347 :
348 2 : if (.not. assertion) then
349 : ! LCOV_EXCL_START
350 : if (Test%isDebugMode) then
351 : write(Test%outputUnit,"(*(g0,:,' '))")
352 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(2)", Test%Image%id, PD1%Err%occurred
353 : write(Test%outputUnit,"(*(g0,:,' '))")
354 : end if
355 : return
356 : end if
357 : ! LCOV_EXCL_STOP
358 :
359 2 : call Test%Image%sync()
360 :
361 : block
362 : use System_mod, only: removeFile
363 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
364 : end block
365 :
366 2 : call Test%Image%sync()
367 :
368 : ! restart the simulation with the same configuration
369 :
370 : call PD2%runSampler ( ndim = 2_IK &
371 : , getLogFunc = getLogFuncMVN &
372 : , mpiFinalizeRequested = .false. &
373 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_9" &
374 : , proposalModel = "uniform" &
375 : , randomSeed = 12345_IK &
376 : , sampleSize = 250_IK &
377 : , chainSize = 500_IK &
378 2 : )
379 :
380 2 : assertion = assertion .and. .not. PD2%Err%occurred
381 :
382 2 : if (.not. assertion) then
383 : ! LCOV_EXCL_START
384 : if (Test%isDebugMode) then
385 : write(Test%outputUnit,"(*(g0,:,' '))")
386 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred", Test%Image%id, PD2%Err%occurred
387 : write(Test%outputUnit,"(*(g0,:,' '))")
388 : end if
389 : return
390 : end if
391 : ! LCOV_EXCL_STOP
392 :
393 2 : if (PD2%Image%isLeader) then
394 : ! LCOV_EXCL_START
395 :
396 : assertion = assertion .and. all(shape(PD2%RefinedChain%LogFuncState) == shape(PD1%RefinedChain%LogFuncState))
397 :
398 : if (.not. assertion) then
399 : if (Test%isDebugMode) then
400 : write(*,"(10(g0,:,', '))")
401 : write(*,"(10(g0,:,', '))") "shape(PD1%RefinedChain%LogFuncState)", shape(PD1%RefinedChain%LogFuncState)
402 : write(*,"(10(g0,:,', '))") "shape(PD2%RefinedChain%LogFuncState)", shape(PD2%RefinedChain%LogFuncState)
403 : write(*,"(10(g0,:,', '))")
404 : end if
405 : return
406 : end if
407 :
408 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-6_RK ) ! by default, the output precision is only 8 digits
409 :
410 : if (.not. assertion) then
411 : if (Test%isDebugMode) then
412 : write(Test%outputUnit,"(*(g0,:,' '))")
413 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
414 : write(Test%outputUnit,"(*(g0,:,' '))")
415 : end if
416 : return
417 : end if
418 :
419 : end if
420 : ! LCOV_EXCL_STOP
421 : #endif
422 60 : end function test_runSampler_9
423 :
424 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
425 :
426 : !> \brief
427 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file, with an ASCII output file.
428 2 : module function test_runSampler_10() result(assertion)
429 : implicit none
430 : logical :: assertion
431 2 : type(ParaDXXX_type) :: PD1, PD2
432 2 : assertion = .true.
433 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
434 :
435 : ! Run the fresh simulation as a reference run
436 :
437 : call PD1%runSampler ( ndim = 2_IK &
438 : , getLogFunc = getLogFuncMVN &
439 : , mpiFinalizeRequested = .false. &
440 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_10.ref" &
441 : , restartFileFormat = "ascii" &
442 : , proposalModel = "uniform" &
443 : , randomSeed = 12345_IK &
444 : , sampleSize = 50_IK &
445 : , chainSize = 500_IK &
446 2 : )
447 2 : assertion = assertion .and. .not. PD1%Err%occurred
448 : if (.not. assertion) return ! LCOV_EXCL_LINE
449 :
450 2 : call Test%Image%sync()
451 :
452 : ! Run the fresh simulation
453 :
454 : call PD1%runSampler ( ndim = 2_IK &
455 : , getLogFunc = getLogFuncMVN &
456 : , mpiFinalizeRequested = .false. &
457 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_10" &
458 : , restartFileFormat = "ascii" &
459 : , proposalModel = "uniform" &
460 : , randomSeed = 12345_IK &
461 : , sampleSize = 50_IK &
462 : , chainSize = 500_IK &
463 2 : )
464 2 : assertion = assertion .and. .not. PD1%Err%occurred
465 : if (.not. assertion) return ! LCOV_EXCL_LINE
466 :
467 2 : call Test%Image%sync()
468 :
469 : block
470 : use System_mod, only: removeFile
471 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
472 : end block
473 :
474 2 : call Test%Image%sync()
475 :
476 : ! restart the simulation with the same configuration
477 :
478 : call PD2%runSampler ( ndim = 2_IK &
479 : , getLogFunc = getLogFuncMVN &
480 : , mpiFinalizeRequested = .false. &
481 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_10" &
482 : , restartFileFormat = "ascii" &
483 : , proposalModel = "uniform" &
484 : , randomSeed = 12345_IK &
485 : , sampleSize = 50_IK &
486 : , chainSize = 500_IK &
487 2 : )
488 2 : assertion = assertion .and. .not. PD2%Err%occurred
489 : if (.not. assertion) return ! LCOV_EXCL_LINE
490 :
491 : if (PD2%Image%isLeader) assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-6_RK ) ! LCOV_EXCL_LINE ! by default, the output precision is only 8 digits
492 :
493 : #endif
494 60 : end function test_runSampler_10
495 :
496 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
497 :
498 : !> \brief
499 : !> Test the ParaDXXX sampler to restart a complete simulation without an
500 : !> output sample file, with an ASCII output file but with binary chain file.
501 2 : module function test_runSampler_11() result(assertion)
502 : implicit none
503 : logical :: assertion
504 2 : type(ParaDXXX_type) :: PD1, PD2
505 2 : assertion = .true.
506 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
507 :
508 : ! Run the fresh simulation as a reference run
509 :
510 : call PD1%runSampler ( ndim = 2_IK &
511 : , getLogFunc = getLogFuncMVN &
512 : , mpiFinalizeRequested = .false. &
513 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_11.ref" &
514 : , restartFileFormat = "ascii" &
515 : , chainFileFormat = "verbose" &
516 : , proposalModel = "uniform" &
517 : , randomSeed = 12345_IK &
518 : , sampleSize = 250_IK &
519 : , chainSize = 500_IK &
520 2 : )
521 2 : assertion = assertion .and. .not. PD1%Err%occurred
522 : if (.not. assertion) return ! LCOV_EXCL_LINE
523 :
524 2 : call Test%Image%sync()
525 :
526 : ! Run the fresh simulation
527 :
528 : call PD1%runSampler ( ndim = 2_IK &
529 : , getLogFunc = getLogFuncMVN &
530 : , mpiFinalizeRequested = .false. &
531 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_11" &
532 : , restartFileFormat = "ascii" &
533 : , chainFileFormat = "verbose" &
534 : , proposalModel = "uniform" &
535 : , randomSeed = 12345_IK &
536 : , sampleSize = 250_IK &
537 : , chainSize = 500_IK &
538 2 : )
539 2 : assertion = assertion .and. .not. PD1%Err%occurred
540 : if (.not. assertion) return ! LCOV_EXCL_LINE
541 :
542 2 : call Test%Image%sync()
543 :
544 : block
545 : use System_mod, only: removeFile
546 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
547 : end block
548 :
549 2 : call Test%Image%sync()
550 :
551 : ! restart the simulation with the same configuration
552 :
553 : call PD2%runSampler ( ndim = 2_IK &
554 : , getLogFunc = getLogFuncMVN &
555 : , mpiFinalizeRequested = .false. &
556 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_11" &
557 : , restartFileFormat = "ascii" &
558 : , chainFileFormat = "verbose" &
559 : , proposalModel = "uniform" &
560 : , randomSeed = 12345_IK &
561 : , sampleSize = 250_IK &
562 : , chainSize = 500_IK &
563 2 : )
564 2 : assertion = assertion .and. .not. PD2%Err%occurred
565 : if (.not. assertion) return ! LCOV_EXCL_LINE
566 :
567 : if (PD2%Image%isLeader) assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-6_RK ) ! LCOV_EXCL_LINE ! by default, the output precision is only 8 digits
568 : #endif
569 60 : end function test_runSampler_11
570 :
571 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
572 :
573 : !> \brief
574 : !> Test the ParaDXXX sampler to restart a complete simulation without an
575 : !> output sample file, with an ASCII output file but with binary chain file.
576 2 : module function test_runSampler_12() result(assertion)
577 : implicit none
578 : logical :: assertion
579 2 : type(ParaDXXX_type) :: PD1, PD2
580 2 : assertion = .true.
581 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
582 :
583 : ! Run the fresh simulation as a reference run
584 :
585 : call PD1%runSampler ( ndim = 2_IK &
586 : , getLogFunc = getLogFuncMVN &
587 : , mpiFinalizeRequested = .false. &
588 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_12.ref" &
589 : , restartFileFormat = "ascii" &
590 : , chainFileFormat = "binary" &
591 : , proposalModel = "uniform" &
592 : , randomSeed = 12345_IK &
593 : , sampleSize = 250_IK &
594 : , chainSize = 500_IK &
595 2 : )
596 2 : assertion = assertion .and. .not. PD1%Err%occurred
597 : if (.not. assertion) return ! LCOV_EXCL_LINE
598 :
599 2 : call Test%Image%sync()
600 :
601 : ! Run the fresh simulation
602 :
603 : call PD1%runSampler ( ndim = 2_IK &
604 : , getLogFunc = getLogFuncMVN &
605 : , mpiFinalizeRequested = .false. &
606 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_12" &
607 : , restartFileFormat = "ascii" &
608 : , chainFileFormat = "binary" &
609 : , proposalModel = "uniform" &
610 : , randomSeed = 12345_IK &
611 : , sampleSize = 250_IK &
612 : , chainSize = 500_IK &
613 2 : )
614 2 : assertion = assertion .and. .not. PD1%Err%occurred
615 : if (.not. assertion) return ! LCOV_EXCL_LINE
616 :
617 2 : call Test%Image%sync()
618 :
619 : block
620 : use System_mod, only: removeFile
621 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
622 : end block
623 :
624 2 : call Test%Image%sync()
625 :
626 : ! restart the simulation with the same configuration
627 :
628 : call PD2%runSampler ( ndim = 2_IK &
629 : , getLogFunc = getLogFuncMVN &
630 : , mpiFinalizeRequested = .false. &
631 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_12" &
632 : , restartFileFormat = "ascii" &
633 : , chainFileFormat = "binary" &
634 : , proposalModel = "uniform" &
635 : , randomSeed = 12345_IK &
636 : , sampleSize = 250_IK &
637 : , chainSize = 500_IK &
638 2 : )
639 2 : assertion = assertion .and. .not. PD2%Err%occurred
640 : if (.not. assertion) return ! LCOV_EXCL_LINE
641 :
642 : if (PD2%Image%isLeader) assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-6_RK ) ! LCOV_EXCL_LINE ! by default, the output precision is only 8 digits
643 :
644 : #endif
645 60 : end function test_runSampler_12
646 :
647 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
648 :
649 : !> \brief
650 : !> Test the ParaDXXX sampler with an empty external input file.
651 : !> This should first search for the specific namelist, then for the generic namelist,
652 : !> then issues a warning and runs the simulation without setting the simulation specifications.
653 2 : module function test_runSampler_13() result(assertion)
654 : implicit none
655 : logical :: assertion
656 2 : type(ParaDXXX_type) :: PD
657 2 : assertion = .true.
658 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
659 : call PD%runSampler ( ndim = 1_IK &
660 : , getLogFunc = getLogFuncMVN &
661 : , mpiFinalizeRequested = .false. &
662 : , inputFile = Test%inDir//"/Test_ParaDXXX_mod@test_runSampler_13.in" &
663 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_13" &
664 2 : )
665 2 : assertion = assertion .and. .not. PD%Err%occurred
666 : #endif
667 32 : end function test_runSampler_13
668 :
669 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
670 :
671 : !> \brief
672 : !> Test the ParaDXXX sampler with an empty external input file.
673 : !> This should first search for the specific namelist, then for the generic namelist,
674 : !> then issues a warning and runs the simulation without setting the simulation specifications.
675 2 : module function test_runSampler_14() result(assertion)
676 : implicit none
677 : logical :: assertion
678 2 : type(ParaDXXX_type) :: PD
679 2 : assertion = .true.
680 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
681 : call PD%runSampler ( ndim = 1_IK &
682 : , getLogFunc = getLogFuncMVN &
683 : , mpiFinalizeRequested = .false. &
684 : , inputFile = "&nonsense /" &
685 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_14" &
686 2 : )
687 2 : assertion = assertion .and. .not. PD%Err%occurred
688 : #endif
689 32 : end function test_runSampler_14
690 :
691 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
692 :
693 : !> \brief
694 : !> Test whether the read method of the `ParaMCMCRefinedChain_type` class can successfully read an external sample file.
695 2 : module function test_runSampler_15() result(assertion)
696 2 : use ParaMCMCRefinedChain_mod, only: readRefinedChain, RefinedChain_type
697 : use Statistics_mod, only: flatten
698 : implicit none
699 : logical :: assertion
700 : real(RK) , parameter :: tolerance = 1.e-10_RK
701 : character(*), parameter :: DELIM = "delim"
702 : integer(IK) , parameter :: NDIM = 2_IK
703 2 : type(ParaDXXX_type) :: PD
704 2 : type(RefinedChain_type) :: RefinedChain
705 : real(RK), allocatable :: Difference(:,:)
706 : real(RK), allocatable :: FlattenedLogFuncState(:,:)
707 2 : assertion = .true.
708 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
709 :
710 : call PD%runSampler ( ndim = NDIM &
711 : , getLogFunc = getLogFuncMVN &
712 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_15" &
713 : , mpiFinalizeRequested = .false. &
714 : , outputRealPrecision = 15_IK &
715 : , outputDelimiter = DELIM &
716 : , sampleSize = 10_IK &
717 2 : )
718 2 : assertion = assertion .and. .not. PD%Err%occurred
719 : if (.not. assertion) return ! LCOV_EXCL_LINE
720 :
721 2 : if (PD%Image%isLeader) then
722 : ! LCOV_EXCL_START
723 :
724 : RefinedChain = readRefinedChain( sampleFilePath = PD%SampleFile%Path%original, delimiter = PD%SpecBase%OutputDelimiter%val, ndim = PD%nd%val )
725 :
726 : assertion = assertion .and. RefinedChain%numRefinement == 0_IK
727 :
728 : if (.not. assertion) then
729 : if (Test%isDebugMode) then
730 : write(*,"(10(g0,:,', '))")
731 : write(*,"(10(g0,:,', '))") "RefinedChain%numRefinement", RefinedChain%numRefinement
732 : write(*,"(10(g0,:,', '))")
733 : end if
734 : return
735 : end if
736 :
737 : ! NOTE: Keep in mind that `PD%RefinedChain` is a weighted chain internally, but unweighted when read from the external file.
738 : FlattenedLogFuncState = flatten(nd = PD%nd%val + 1, np = PD%RefinedChain%Count(PD%RefinedChain%numRefinement)%compact, Point = PD%RefinedChain%LogFuncState, Weight = PD%RefinedChain%Weight)
739 : assertion = assertion .and. all( shape(FlattenedLogFuncState) == shape(RefinedChain%LogFuncState) )
740 :
741 : if (.not. assertion) then
742 : if (Test%isDebugMode) then
743 : write(*,"(10(g0,:,', '))")
744 : write(*,"(10(g0,:,', '))") "shape(PD%RefinedChain%LogFuncState) ", shape(PD%RefinedChain%LogFuncState)
745 : write(*,"(10(g0,:,', '))") "shape(RefinedChain%LogFuncState) ", shape(RefinedChain%LogFuncState)
746 : write(*,"(10(g0,:,', '))")
747 : write(*,"(10(g0,:,', '))") "PD%RefinedChain%LogFuncState ", PD%RefinedChain%LogFuncState
748 : write(*,"(10(g0,:,', '))") "RefinedChain%LogFuncState ", RefinedChain%LogFuncState
749 : write(*,"(10(g0,:,', '))")
750 : end if
751 : return
752 : end if
753 :
754 : Difference = abs(FlattenedLogFuncState - RefinedChain%LogFuncState)
755 : assertion = assertion .and. all(Difference < tolerance)
756 :
757 : if (.not. assertion) then
758 : if (Test%isDebugMode) then
759 : write(*,"(10(g0,:,', '))")
760 : write(*,"(10(g0,:,', '))") "PD%RefinedChain%LogFuncState", PD%RefinedChain%LogFuncState
761 : write(*,"(10(g0,:,', '))") "RefinedChain%LogFuncState ", RefinedChain%LogFuncState
762 : write(*,"(10(g0,:,', '))") "Difference ", Difference
763 : write(*,"(10(g0,:,', '))")
764 : end if
765 : return
766 : end if
767 :
768 : end if
769 : ! LCOV_EXCL_STOP
770 :
771 : #endif
772 40 : end function test_runSampler_15
773 :
774 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
775 :
776 : !> \brief
777 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file, in multichain parallelism.
778 2 : module function test_runSampler_16() result(assertion)
779 : implicit none
780 : logical :: assertion
781 2 : type(ParaDXXX_type) :: PD1, PD2
782 2 : assertion = .true.
783 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
784 :
785 : ! Run the fresh simulation
786 :
787 : call PD1%runSampler ( ndim = 2_IK &
788 : , getLogFunc = getLogFuncMVN &
789 : , mpiFinalizeRequested = .false. &
790 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_16" &
791 : , parallelizationModel = "multi chain" &
792 : , outputRealPrecision = 16_IK &
793 : , randomSeed = 12345_IK &
794 : , sampleSize = 250_IK &
795 : , chainSize = 500_IK &
796 2 : )
797 2 : assertion = assertion .and. .not. PD1%Err%occurred
798 :
799 2 : if (.not. assertion) then
800 : ! LCOV_EXCL_START
801 : if (Test%isDebugMode) then
802 : write(Test%outputUnit,"(*(g0,:,' '))")
803 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
804 : write(Test%outputUnit,"(*(g0,:,' '))")
805 : end if
806 : return
807 : end if
808 : ! LCOV_EXCL_STOP
809 :
810 2 : call Test%Image%sync()
811 :
812 : block
813 : use System_mod, only: removeFile
814 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
815 : end block
816 :
817 2 : call Test%Image%sync()
818 :
819 : ! restart the simulation with the same configuration
820 :
821 : call PD2%runSampler ( ndim = 2_IK &
822 : , getLogFunc = getLogFuncMVN &
823 : , mpiFinalizeRequested = .false. &
824 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_16" &
825 : , parallelizationModel = "multi chain" &
826 : , outputRealPrecision = 16_IK &
827 : , randomSeed = 12345_IK &
828 : , sampleSize = 250_IK &
829 : , chainSize = 500_IK &
830 2 : )
831 :
832 2 : assertion = assertion .and. .not. PD2%Err%occurred
833 :
834 2 : if (.not. assertion) then
835 : ! LCOV_EXCL_START
836 : if (Test%isDebugMode) then
837 : write(Test%outputUnit,"(*(g0,:,' '))")
838 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
839 : write(Test%outputUnit,"(*(g0,:,' '))")
840 : end if
841 : return
842 : end if
843 : ! LCOV_EXCL_STOP
844 :
845 2 : if (PD2%Image%isLeader) then
846 : ! LCOV_EXCL_START
847 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-12_RK ) ! by default, the output precision is only 8 digits
848 : if (.not. assertion) then
849 : ! LCOV_EXCL_START
850 : if (Test%isDebugMode) then
851 : write(Test%outputUnit,"(*(g0,:,' '))")
852 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
853 : write(Test%outputUnit,"(*(g0,:,' '))")
854 : end if
855 : return
856 : end if
857 : end if
858 : ! LCOV_EXCL_STOP
859 : #endif
860 60 : end function test_runSampler_16
861 :
862 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
863 :
864 : !> \brief
865 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file, in singlechain parallelism.
866 : !> Also, ensure the delayed rejection sampling is activated by setting a low value for `adaptiveUpdatePeriod` and
867 : !> a large value for `scaleFactor`.
868 2 : module function test_runSampler_17() result(assertion)
869 2 : use Constants_mod, only: RK
870 : implicit none
871 : logical :: assertion
872 2 : type(ParaDXXX_type) :: PD1, PD2
873 : real(RK), parameter :: targetAcceptanceRate(*) = [0.2_RK, 0.23_RK]
874 2 : assertion = .true.
875 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
876 :
877 : ! Run the fresh simulation
878 :
879 : call PD1%runSampler ( ndim = 2_IK &
880 : , getLogFunc = getLogFuncMVN &
881 : , mpiFinalizeRequested = .false. &
882 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_17" &
883 : , targetAcceptanceRate = targetAcceptanceRate &
884 : , parallelizationModel = "single chain" &
885 : , greedyAdaptationCount = 30_IK & ! This must remain larger than adaptiveUpdatePeriod
886 : , delayedRejectionCount = 4_IK &
887 : , adaptiveUpdatePeriod = 1_IK &
888 : , outputRealPrecision = 16_IK &
889 : , scaleFactor = "5 * gelman" &
890 : , randomSeed = 12345_IK &
891 : , sampleSize = 250_IK &
892 : , chainSize = 500_IK &
893 2 : )
894 2 : assertion = assertion .and. .not. PD1%Err%occurred
895 :
896 2 : if (.not. assertion) then
897 : ! LCOV_EXCL_START
898 : if (Test%isDebugMode) then
899 : write(Test%outputUnit,"(*(g0,:,' '))")
900 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
901 : write(Test%outputUnit,"(*(g0,:,' '))")
902 : end if
903 : return
904 : end if
905 : ! LCOV_EXCL_STOP
906 :
907 2 : call Test%Image%sync()
908 :
909 : block
910 : use System_mod, only: removeFile
911 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
912 : end block
913 :
914 2 : call Test%Image%sync()
915 :
916 : ! restart the simulation with the same configuration
917 :
918 : call PD2%runSampler ( ndim = 2_IK &
919 : , getLogFunc = getLogFuncMVN &
920 : , mpiFinalizeRequested = .false. &
921 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_17" &
922 : , targetAcceptanceRate = targetAcceptanceRate &
923 : , parallelizationModel = "single chain" &
924 : , greedyAdaptationCount = 30_IK & ! This must remain larger than adaptiveUpdatePeriod
925 : , delayedRejectionCount = 4_IK &
926 : , adaptiveUpdatePeriod = 1_IK &
927 : , outputRealPrecision = 16_IK &
928 : , scaleFactor = "5 * gelman" &
929 : , randomSeed = 12345_IK &
930 : , sampleSize = 250_IK &
931 : , chainSize = 500_IK &
932 2 : )
933 :
934 2 : assertion = assertion .and. .not. PD2%Err%occurred
935 :
936 2 : if (.not. assertion) then
937 : ! LCOV_EXCL_START
938 : if (Test%isDebugMode) then
939 : write(Test%outputUnit,"(*(g0,:,' '))")
940 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
941 : write(Test%outputUnit,"(*(g0,:,' '))")
942 : end if
943 : return
944 : end if
945 : ! LCOV_EXCL_STOP
946 :
947 2 : if (PD2%Image%isLeader) then
948 : ! LCOV_EXCL_START
949 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-12_RK ) ! by default, the output precision is only 8 digits
950 : if (.not. assertion) then
951 : ! LCOV_EXCL_START
952 : if (Test%isDebugMode) then
953 : write(Test%outputUnit,"(*(g0,:,' '))")
954 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
955 : write(Test%outputUnit,"(*(g0,:,' '))")
956 : end if
957 : return
958 : end if
959 : end if
960 : ! LCOV_EXCL_STOP
961 : #endif
962 58 : end function test_runSampler_17
963 :
964 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
965 :
966 : !> \brief
967 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file, in singlechain parallelism.
968 : !> Also, ensure the delayed rejection sampling is activated by setting a low value for `adaptiveUpdatePeriod` and
969 : !> a large value for `scaleFactor`.
970 : !> Also, use uniform proposal, as opposed to Normal in `test_runSampler_17()`.
971 2 : module function test_runSampler_18() result(assertion)
972 2 : use Constants_mod, only: RK
973 : implicit none
974 : logical :: assertion
975 2 : type(ParaDXXX_type) :: PD1, PD2
976 : real(RK), parameter :: targetAcceptanceRate(*) = [0.2_RK, 0.23_RK]
977 2 : assertion = .true.
978 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
979 :
980 : ! Run the fresh simulation
981 :
982 : call PD1%runSampler ( ndim = 2_IK &
983 : , getLogFunc = getLogFuncMVN &
984 : , mpiFinalizeRequested = .false. &
985 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_18" &
986 : , targetAcceptanceRate = targetAcceptanceRate &
987 : , parallelizationModel = "single chain" &
988 : , greedyAdaptationCount = 30_IK & ! This must remain larger than adaptiveUpdatePeriod
989 : , delayedRejectionCount = 4_IK &
990 : , adaptiveUpdatePeriod = 1_IK &
991 : , outputRealPrecision = 16_IK &
992 : , scaleFactor = "5 * gelman" &
993 : , proposalModel = "uniform" &
994 : , randomSeed = 12345_IK &
995 : , sampleSize = 250_IK &
996 : , chainSize = 500_IK &
997 2 : )
998 2 : assertion = assertion .and. .not. PD1%Err%occurred
999 :
1000 2 : if (.not. assertion) then
1001 : ! LCOV_EXCL_START
1002 : if (Test%isDebugMode) then
1003 : write(Test%outputUnit,"(*(g0,:,' '))")
1004 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
1005 : write(Test%outputUnit,"(*(g0,:,' '))")
1006 : end if
1007 : return
1008 : end if
1009 : ! LCOV_EXCL_STOP
1010 :
1011 2 : call Test%Image%sync()
1012 :
1013 : block
1014 : use System_mod, only: removeFile
1015 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
1016 : end block
1017 :
1018 2 : call Test%Image%sync()
1019 :
1020 : ! restart the simulation with the same configuration
1021 :
1022 : call PD2%runSampler ( ndim = 2_IK &
1023 : , getLogFunc = getLogFuncMVN &
1024 : , mpiFinalizeRequested = .false. &
1025 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_18" &
1026 : , targetAcceptanceRate = targetAcceptanceRate &
1027 : , parallelizationModel = "singlechain" &
1028 : , greedyAdaptationCount = 30_IK & ! This must remain larger than adaptiveUpdatePeriod
1029 : , delayedRejectionCount = 4_IK &
1030 : , adaptiveUpdatePeriod = 1_IK &
1031 : , outputRealPrecision = 16_IK &
1032 : , scaleFactor = "5 * gelman" &
1033 : , proposalModel = "uniform" &
1034 : , randomSeed = 12345_IK &
1035 : , sampleSize = 250_IK &
1036 : , chainSize = 500_IK &
1037 2 : )
1038 :
1039 2 : assertion = assertion .and. .not. PD2%Err%occurred
1040 :
1041 2 : if (.not. assertion) then
1042 : ! LCOV_EXCL_START
1043 : if (Test%isDebugMode) then
1044 : write(Test%outputUnit,"(*(g0,:,' '))")
1045 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
1046 : write(Test%outputUnit,"(*(g0,:,' '))")
1047 : end if
1048 : return
1049 : end if
1050 : ! LCOV_EXCL_STOP
1051 :
1052 2 : if (PD2%Image%isLeader) then
1053 : ! LCOV_EXCL_START
1054 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-12_RK ) ! by default, the output precision is only 8 digits
1055 : if (.not. assertion) then
1056 : ! LCOV_EXCL_START
1057 : if (Test%isDebugMode) then
1058 : write(Test%outputUnit,"(*(g0,:,' '))")
1059 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
1060 : write(Test%outputUnit,"(*(g0,:,' '))")
1061 : end if
1062 : return
1063 : end if
1064 : end if
1065 : ! LCOV_EXCL_STOP
1066 : #endif
1067 58 : end function test_runSampler_18
1068 :
1069 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1070 :
1071 : !> \brief
1072 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file, in multichain parallelism.
1073 : !> Also, ensure the delayed rejection sampling is activated by setting a low value for `adaptiveUpdatePeriod` and
1074 : !> a large value for `scaleFactor`.
1075 : !> Also, avoid delayed rejection, as opposed to what is done in `test_runSampler_17()`.
1076 : !> Also, request the `maxCumSum` sample refinement and count method.
1077 2 : module function test_runSampler_19() result(assertion)
1078 2 : use Constants_mod, only: RK
1079 : implicit none
1080 : logical :: assertion
1081 2 : type(ParaDXXX_type) :: PD1, PD2
1082 : real(RK), parameter :: targetAcceptanceRate(*) = [0.2_RK, 0.23_RK]
1083 2 : assertion = .true.
1084 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
1085 :
1086 : ! Run the fresh simulation
1087 :
1088 : call PD1%runSampler ( ndim = 2_IK &
1089 : , getLogFunc = getLogFuncMVN &
1090 : , mpiFinalizeRequested = .false. &
1091 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_19" &
1092 : , sampleRefinementMethod = "MaxCumSumAutoCorr-median" &
1093 : , targetAcceptanceRate = targetAcceptanceRate &
1094 : , parallelizationModel = "multi chain" &
1095 : , greedyAdaptationCount = 2_IK & ! This must remain larger than adaptiveUpdatePeriod
1096 : , adaptiveUpdatePeriod = 1_IK &
1097 : , outputRealPrecision = 16_IK &
1098 : , chainFileFormat = "verbose" &
1099 : , scaleFactor = "5 * gelman" &
1100 : , proposalModel = "uniform" &
1101 : , randomSeed = 12345_IK &
1102 : , sampleSize = -50_IK &
1103 : , chainSize = 700_IK &
1104 2 : )
1105 2 : assertion = assertion .and. .not. PD1%Err%occurred
1106 :
1107 2 : if (.not. assertion) then
1108 : ! LCOV_EXCL_START
1109 : if (Test%isDebugMode) then
1110 : write(Test%outputUnit,"(*(g0,:,' '))")
1111 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
1112 : write(Test%outputUnit,"(*(g0,:,' '))")
1113 : end if
1114 : return
1115 : end if
1116 : ! LCOV_EXCL_STOP
1117 :
1118 2 : call Test%Image%sync()
1119 :
1120 : block
1121 : use System_mod, only: removeFile
1122 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
1123 : end block
1124 :
1125 2 : call Test%Image%sync()
1126 :
1127 : ! restart the simulation with the same configuration
1128 :
1129 : call PD2%runSampler ( ndim = 2_IK &
1130 : , getLogFunc = getLogFuncMVN &
1131 : , mpiFinalizeRequested = .false. &
1132 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_19" &
1133 : , sampleRefinementMethod = "MaxCumSumAutoCorr-median" &
1134 : , targetAcceptanceRate = targetAcceptanceRate &
1135 : , parallelizationModel = "multi chain" &
1136 : , greedyAdaptationCount = 2_IK & ! This must remain larger than adaptiveUpdatePeriod
1137 : , adaptiveUpdatePeriod = 1_IK &
1138 : , outputRealPrecision = 16_IK &
1139 : , chainFileFormat = "verbose" &
1140 : , scaleFactor = "5 * gelman" &
1141 : , proposalModel = "uniform" &
1142 : , randomSeed = 12345_IK &
1143 : , sampleSize = -50_IK &
1144 : , chainSize = 700_IK &
1145 2 : )
1146 :
1147 2 : assertion = assertion .and. .not. PD2%Err%occurred
1148 :
1149 2 : if (.not. assertion) then
1150 : ! LCOV_EXCL_START
1151 : if (Test%isDebugMode) then
1152 : write(Test%outputUnit,"(*(g0,:,' '))")
1153 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
1154 : write(Test%outputUnit,"(*(g0,:,' '))")
1155 : end if
1156 : return
1157 : end if
1158 : ! LCOV_EXCL_STOP
1159 :
1160 2 : if (PD2%Image%isLeader) then
1161 : ! LCOV_EXCL_START
1162 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-12_RK ) ! by default, the output precision is only 8 digits
1163 : if (.not. assertion) then
1164 : ! LCOV_EXCL_START
1165 : if (Test%isDebugMode) then
1166 : write(Test%outputUnit,"(*(g0,:,' '))")
1167 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
1168 : write(Test%outputUnit,"(*(g0,:,' '))")
1169 : end if
1170 : return
1171 : end if
1172 : end if
1173 : ! LCOV_EXCL_STOP
1174 : #endif
1175 58 : end function test_runSampler_19
1176 :
1177 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1178 :
1179 : !> \brief
1180 : !> Test the ParaDXXX sampler to restart a complete simulation without an output sample file, in multichain parallelism.
1181 : !> Also, ensure the delayed rejection sampling is activated by setting a low value for `adaptiveUpdatePeriod` and
1182 : !> a large value for `scaleFactor`.
1183 : !> Also, avoid delayed rejection, as opposed to what is done in `test_runSampler_17()`.
1184 : !> Also, request the `cutoffAutoCorr` sample refinement and count method.
1185 2 : module function test_runSampler_20() result(assertion)
1186 2 : use Constants_mod, only: RK
1187 : implicit none
1188 : logical :: assertion
1189 2 : type(ParaDXXX_type) :: PD1, PD2
1190 : real(RK), parameter :: targetAcceptanceRate(*) = [0.2_RK, 0.23_RK]
1191 2 : assertion = .true.
1192 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
1193 :
1194 : ! Run the fresh simulation
1195 :
1196 : call PD1%runSampler ( ndim = 2_IK &
1197 : , getLogFunc = getLogFuncMVN &
1198 : , mpiFinalizeRequested = .false. &
1199 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_20" &
1200 : , sampleRefinementMethod = "cutoffAutoCorr-maximum" &
1201 : , targetAcceptanceRate = targetAcceptanceRate &
1202 : , parallelizationModel = "multi chain" &
1203 : , greedyAdaptationCount = 2_IK & ! This must remain larger than adaptiveUpdatePeriod
1204 : , adaptiveUpdatePeriod = 1_IK &
1205 : , outputRealPrecision = 16_IK &
1206 : , chainFileFormat = "binary" &
1207 : , scaleFactor = "5 * gelman" &
1208 : , proposalModel = "uniform" &
1209 : , randomSeed = 12345_IK &
1210 : , sampleSize = -50_IK &
1211 : , chainSize = 700_IK &
1212 2 : )
1213 2 : assertion = assertion .and. .not. PD1%Err%occurred
1214 :
1215 2 : if (.not. assertion) then
1216 : ! LCOV_EXCL_START
1217 : if (Test%isDebugMode) then
1218 : write(Test%outputUnit,"(*(g0,:,' '))")
1219 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD1%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
1220 : write(Test%outputUnit,"(*(g0,:,' '))")
1221 : end if
1222 : return
1223 : end if
1224 : ! LCOV_EXCL_STOP
1225 :
1226 2 : call Test%Image%sync()
1227 :
1228 : block
1229 : use System_mod, only: removeFile
1230 2 : call removeFile(PD1%SampleFile%Path%original, PD1%Err) ! delete the sample file
1231 : end block
1232 :
1233 2 : call Test%Image%sync()
1234 :
1235 : ! restart the simulation with the same configuration
1236 :
1237 : call PD2%runSampler ( ndim = 2_IK &
1238 : , getLogFunc = getLogFuncMVN &
1239 : , mpiFinalizeRequested = .false. &
1240 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_20" &
1241 : , sampleRefinementMethod = "cutoffAutoCorr-maximum" &
1242 : , targetAcceptanceRate = targetAcceptanceRate &
1243 : , parallelizationModel = "multi chain" &
1244 : , greedyAdaptationCount = 2_IK & ! This must remain larger than adaptiveUpdatePeriod
1245 : , adaptiveUpdatePeriod = 1_IK &
1246 : , outputRealPrecision = 16_IK &
1247 : , chainFileFormat = "binary" &
1248 : , scaleFactor = "5 * gelman" &
1249 : , proposalModel = "uniform" &
1250 : , randomSeed = 12345_IK &
1251 : , sampleSize = -50_IK &
1252 : , chainSize = 700_IK &
1253 2 : )
1254 :
1255 2 : assertion = assertion .and. .not. PD2%Err%occurred
1256 :
1257 2 : if (.not. assertion) then
1258 : ! LCOV_EXCL_START
1259 : if (Test%isDebugMode) then
1260 : write(Test%outputUnit,"(*(g0,:,' '))")
1261 : write(Test%outputUnit,"(*(g0,:,' '))") "process, PD2%Err%occurred(1)", Test%Image%id, PD1%Err%occurred
1262 : write(Test%outputUnit,"(*(g0,:,' '))")
1263 : end if
1264 : return
1265 : end if
1266 : ! LCOV_EXCL_STOP
1267 :
1268 2 : if (PD2%Image%isLeader) then
1269 : ! LCOV_EXCL_START
1270 : assertion = assertion .and. all( abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState) < 1.e-12_RK ) ! by default, the output precision is only 8 digits
1271 : if (.not. assertion) then
1272 : ! LCOV_EXCL_START
1273 : if (Test%isDebugMode) then
1274 : write(Test%outputUnit,"(*(g0,:,' '))")
1275 : write(Test%outputUnit,"(*(g0,:,' '))") "process, Difference:", Test%Image%id, abs(PD2%RefinedChain%LogFuncState - PD1%RefinedChain%LogFuncState)
1276 : write(Test%outputUnit,"(*(g0,:,' '))")
1277 : end if
1278 : return
1279 : end if
1280 : end if
1281 : ! LCOV_EXCL_STOP
1282 : #endif
1283 58 : end function test_runSampler_20
1284 :
1285 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1286 :
1287 : !> \brief
1288 : !> Test whether the read method of the `ParaMCMCRefinedChain_type` class can successfully read an external sample file.
1289 : !> \remark
1290 : !> This is similar to test #15, except that it also verifies the functionality of the optional argument `tenabled` to `readRefinedChain()`.
1291 2 : module function test_runSampler_21() result(assertion)
1292 2 : use ParaMCMCRefinedChain_mod, only: readRefinedChain, RefinedChain_type
1293 : use Statistics_mod, only: flatten
1294 : implicit none
1295 : logical :: assertion
1296 : real(RK) , parameter :: tolerance = 1.e-10_RK
1297 : character(*), parameter :: DELIM = "delim"
1298 : integer(IK) , parameter :: NDIM = 2_IK
1299 2 : type(ParaDXXX_type) :: PD
1300 2 : type(RefinedChain_type) :: RefinedChain
1301 : real(RK), allocatable :: Difference(:,:)
1302 : real(RK), allocatable :: FlattenedLogFuncState(:,:)
1303 2 : assertion = .true.
1304 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
1305 :
1306 : call PD%runSampler ( ndim = NDIM &
1307 : , getLogFunc = getLogFuncMVN &
1308 : , mpiFinalizeRequested = .false. &
1309 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_21" &
1310 : , sampleRefineMentMethod = "batchMeans-med" &
1311 : , parallelizationModel = "multichain" &
1312 : , outputRealPrecision = 15_IK &
1313 : , outputDelimiter = DELIM &
1314 : , sampleSize = 10_IK &
1315 2 : )
1316 2 : assertion = assertion .and. .not. PD%Err%occurred
1317 : if (.not. assertion) return ! LCOV_EXCL_LINE
1318 :
1319 2 : if (PD%Image%isLeader) then
1320 :
1321 : RefinedChain = readRefinedChain ( sampleFilePath = PD%SampleFile%Path%original & ! LCOV_EXCL_LINE
1322 : , delimiter = PD%SpecBase%OutputDelimiter%val & ! LCOV_EXCL_LINE
1323 : , tenabled = .true. & ! LCOV_EXCL_LINE
1324 : , ndim = PD%nd%val & ! LCOV_EXCL_LINE
1325 2 : )
1326 :
1327 : RefinedChain = readRefinedChain ( sampleFilePath = PD%SampleFile%Path%original & ! LCOV_EXCL_LINE
1328 : , delimiter = PD%SpecBase%OutputDelimiter%val & ! LCOV_EXCL_LINE
1329 : , tenabled = .false. & ! LCOV_EXCL_LINE
1330 : , ndim = PD%nd%val & ! LCOV_EXCL_LINE
1331 10 : )
1332 :
1333 2 : assertion = assertion .and. RefinedChain%numRefinement == 0_IK
1334 :
1335 2 : if (.not. assertion) then
1336 : ! LCOV_EXCL_START
1337 : if (Test%isDebugMode) then
1338 : write(*,"(10(g0,:,', '))")
1339 : write(*,"(10(g0,:,', '))") "RefinedChain%numRefinement", RefinedChain%numRefinement
1340 : write(*,"(10(g0,:,', '))")
1341 : end if
1342 : return
1343 : end if
1344 : ! LCOV_EXCL_STOP
1345 :
1346 : ! NOTE: Keep in mind that `PD%RefinedChain` is a weighted chain internally, but unweighted when read from the external file.
1347 82 : FlattenedLogFuncState = flatten(nd = PD%nd%val + 1, np = PD%RefinedChain%Count(PD%RefinedChain%numRefinement)%compact, Point = PD%RefinedChain%LogFuncState, Weight = PD%RefinedChain%Weight)
1348 6 : assertion = assertion .and. all( shape(FlattenedLogFuncState) == shape(RefinedChain%LogFuncState) )
1349 :
1350 2 : if (.not. assertion) then
1351 : ! LCOV_EXCL_START
1352 : if (Test%isDebugMode) then
1353 : write(*,"(10(g0,:,', '))")
1354 : write(*,"(10(g0,:,', '))") "shape(PD%RefinedChain%LogFuncState) ", shape(PD%RefinedChain%LogFuncState)
1355 : write(*,"(10(g0,:,', '))") "shape(RefinedChain%LogFuncState) ", shape(RefinedChain%LogFuncState)
1356 : write(*,"(10(g0,:,', '))")
1357 : write(*,"(10(g0,:,', '))") "PD%RefinedChain%LogFuncState ", PD%RefinedChain%LogFuncState
1358 : write(*,"(10(g0,:,', '))") "RefinedChain%LogFuncState ", RefinedChain%LogFuncState
1359 : write(*,"(10(g0,:,', '))")
1360 : end if
1361 : return
1362 : end if
1363 : ! LCOV_EXCL_STOP
1364 :
1365 82 : Difference = abs(FlattenedLogFuncState - RefinedChain%LogFuncState)
1366 82 : assertion = assertion .and. all(Difference < tolerance)
1367 :
1368 2 : if (.not. assertion) then
1369 : ! LCOV_EXCL_START
1370 : if (Test%isDebugMode) then
1371 : write(*,"(10(g0,:,', '))")
1372 : write(*,"(10(g0,:,', '))") "PD%RefinedChain%LogFuncState", PD%RefinedChain%LogFuncState
1373 : write(*,"(10(g0,:,', '))") "RefinedChain%LogFuncState ", RefinedChain%LogFuncState
1374 : write(*,"(10(g0,:,', '))") "Difference ", Difference
1375 : write(*,"(10(g0,:,', '))")
1376 : end if
1377 : return
1378 : end if
1379 : ! LCOV_EXCL_STOP
1380 :
1381 : end if
1382 :
1383 : #endif
1384 40 : end function test_runSampler_21
1385 :
1386 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1387 :
1388 : !> \brief
1389 : !> Test whether an object of [ChainFileContents_type](@ref paramcmcrefinedchain_mod::chainfilecontents_type)
1390 : !> class can handle zero-count input MCMC chain for refinement.
1391 : !> \remark
1392 2 : module function test_runSampler_22() result(assertion)
1393 2 : use ParaMCMCRefinedChain_mod, only: readRefinedChain, RefinedChain_type
1394 : implicit none
1395 : logical :: assertion
1396 : real(RK) , parameter :: tolerance = 1.e-10_RK
1397 : character(*), parameter :: DELIM = "delim"
1398 : integer(IK) , parameter :: NDIM = 2_IK
1399 2 : type(ParaDXXX_type) :: PD
1400 2 : type(RefinedChain_type) :: RefinedChain
1401 2 : assertion = .true.
1402 : #if defined CODECOV_ENABLED || defined SAMPLER_TEST_ENABLED
1403 :
1404 : call PD%runSampler ( ndim = NDIM &
1405 : , getLogFunc = getLogFuncMVN &
1406 : , mpiFinalizeRequested = .false. &
1407 : , outputFileName = Test%outDir//"/"//MODULE_NAME//"/test_runSampler_22" &
1408 : , sampleRefineMentMethod = "maxCumSumAutoCorr-minimum" &
1409 : , parallelizationModel = "multi chain" &
1410 : , chainSize = NDIM + 1_IK &
1411 2 : )
1412 2 : assertion = assertion .and. .not. PD%Err%occurred
1413 : if (.not. assertion) return ! LCOV_EXCL_LINE
1414 :
1415 2 : if (PD%Image%isLeader) then
1416 : ! LCOV_EXCL_START
1417 :
1418 : ! Now, set PD%Chain%ndim and PD%Chain%Count%Compact and to zero.
1419 : ! This should lead to the use of shape(PD%Chain%State(:,:)) to infer
1420 : ! ndim and Compact chain length.
1421 : ! To test the behavior of the procedure in the presence of 1 point in the chain,
1422 : ! we will also resize the input PD%Chain%State(1:ndim,1:PD%Chain%Count%Compact) to PD%Chain%State(1:ndim,1:1)
1423 :
1424 : PD%Chain%ndim = 0
1425 : PD%Chain%Count%compact = 0
1426 : PD%Chain%State = PD%Chain%State(1:NDIM, 1:1)
1427 :
1428 : call PD%RefinedChain%get(CFC = PD%Chain, Err = PD%Err)
1429 :
1430 : assertion = assertion .and. .not. PD%Err%occurred .and. all(PD%RefinedChain%IAC==0_IK) .and. size(PD%RefinedChain%Count)==1_IK .and. PD%RefinedChain%numRefinement==0_IK
1431 : if (.not. assertion) return ! LCOV_EXCL_LINE
1432 :
1433 : if (.not. assertion) then
1434 : if (Test%isDebugMode) then
1435 : write(*,"(10(g0,:,', '))")
1436 : write(*,"(10(g0,:,', '))") "Test%Image%id, RefinedChain%numRefinement", RefinedChain%numRefinement
1437 : write(*,"(10(g0,:,', '))") "Test%Image%id, RefinedChain%IAC", RefinedChain%IAC
1438 : write(*,"(10(g0,:,', '))")
1439 : end if
1440 : return
1441 : end if
1442 :
1443 : end if
1444 : ! LCOV_EXCL_STOP
1445 :
1446 : #endif
1447 34 : end function test_runSampler_22
1448 :
1449 : !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
