|
307 | 307 | @time simulate!(s_gpu, simulator, n_steps; n_threads=1)
|
308 | 308 | coord_diff = sum(sum(map(x -> abs.(x), s.coords .- Array(s_gpu.coords)))) / (3 * n_atoms)
|
309 | 309 | E_diff = abs(potential_energy(s) - potential_energy(s_gpu))
|
310 |
| - @info "$(rpad(name, 19)) - difference per coordinate $coord_diff - potential energy difference $E_diff" |
311 | 310 | @test coord_diff < 1e-4u"nm"
|
312 | 311 | @test E_diff < 5e-4u"kJ * mol^-1"
|
313 | 312 | end
|
|
472 | 471 | @time simulate!(s_gpu, simulator, n_steps)
|
473 | 472 | coord_diff = sum(sum(map(x -> abs.(x), s.coords .- Array(s_gpu.coords)))) / (3 * n_atoms)
|
474 | 473 | E_diff = abs(potential_energy(s) - potential_energy(s_gpu))
|
475 |
| - @info "$(rpad(inter, 19)) - difference per coordinate $coord_diff - potential energy difference $E_diff_start (start) and $E_diff" |
476 | 474 | @test coord_diff < 5e-4u"nm"
|
477 | 475 | @test E_diff < 5e-3u"kJ * mol^-1"
|
478 | 476 | end
|
|
913 | 911 | @test efficiency > 0.2 # This is a fairly arbitrary threshold, but it's a good tests for very bad cases
|
914 | 912 | @test efficiency < 1.0 # Bad acceptance rate?
|
915 | 913 | @info "Exchange Efficiency: $efficiency"
|
| 914 | + |
916 | 915 | for id in eachindex(repsys.replicas)
|
917 | 916 | mean_temp = mean(values(repsys.replicas[id].loggers.temp))
|
918 | 917 | @test (0.9 * temp) < mean_temp < (1.1 * temp)
|
@@ -1063,7 +1062,7 @@ end
|
1063 | 1062 | simulate!(deepcopy(sys), lang, 1_000; n_threads=1, rng=rng)
|
1064 | 1063 | @time simulate!(sys, lang, n_steps; n_threads=1, rng=rng)
|
1065 | 1064 |
|
1066 |
| - @test 270.0u"K" < mean(values(sys.loggers.temperature)) < 300.0u"K" |
| 1065 | + @test 260.0u"K" < mean(values(sys.loggers.temperature)) < 300.0u"K" |
1067 | 1066 | @test 50.0u"kJ * mol^-1" < mean(values(sys.loggers.total_energy )) < 120.0u"kJ * mol^-1"
|
1068 | 1067 | @test 50.0u"kJ * mol^-1" < mean(values(sys.loggers.kinetic_energy)) < 120.0u"kJ * mol^-1"
|
1069 | 1068 | @test mean(values(sys.loggers.potential_energy)) < 0.0u"kJ * mol^-1"
|
@@ -1103,7 +1102,7 @@ end
|
1103 | 1102 | simulate!(deepcopy(sys), sim, 1_000; n_threads=1, rng=rng)
|
1104 | 1103 | @time simulate!(sys, sim, n_steps; n_threads=1, rng=rng)
|
1105 | 1104 |
|
1106 |
| - @test 270.0u"K" < mean(values(sys.loggers.temperature)) < 300.0u"K" |
| 1105 | + @test 260.0u"K" < mean(values(sys.loggers.temperature)) < 300.0u"K" |
1107 | 1106 | @test 50.0u"kJ * mol^-1" < mean(values(sys.loggers.total_energy )) < 120.0u"kJ * mol^-1"
|
1108 | 1107 | @test 50.0u"kJ * mol^-1" < mean(values(sys.loggers.kinetic_energy)) < 120.0u"kJ * mol^-1"
|
1109 | 1108 | @test mean(values(sys.loggers.potential_energy)) < 0.0u"kJ * mol^-1"
|
|
0 commit comments