1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38 package ffx.algorithms.groovy;
39
40 import ffx.algorithms.misc.AlgorithmsTest;
41 import org.junit.Assert;
42 import org.junit.Test;
43
44 public class BARSingleTopologyTest extends AlgorithmsTest {
45
46
47
48
49 @Test
50 public void testBAR() {
51
52 if (!ffxOpenMM) {
53 return;
54 }
55
56
57 String filepath = getResourcePath("testBar/singleTopology/dimethylphosphate.10.xyz");
58 String[] args = {"-t", "298", "--nw", "8", "--ac", "1-13", filepath};
59 binding.setVariable("args", args);
60
61 BAR bar = new BAR(binding).run();
62 algorithmsScript = bar;
63
64 System.out.println("after run");
65 double expectedFepFor = -54.3012;
66 double actualFepFor = bar.getFepFor();
67 Assert.assertEquals(expectedFepFor, actualFepFor, 0.2);
68
69 double expectedFepBack = -52.9186;
70 double actualFepBack = bar.getFepBack();
71 Assert.assertEquals(expectedFepBack, actualFepBack, 0.2);
72
73 double expectedhFor = -32.0327;
74 double actualhFor = bar.gethFor();
75 Assert.assertEquals(expectedhFor, actualhFor, 7);
76
77 double expectedhBack = 52.9186;
78 double actualhBack = bar.gethBack();
79 Assert.assertEquals(expectedhBack, actualhBack, 0.1);
80
81 double expectedhBAR = 14.6302;
82 double actualhBAR = bar.gethBAR();
83 Assert.assertEquals(expectedhBAR, actualhBAR, 1);
84
85 double expectedsFor = 0.0771;
86 double actualsFor = bar.getsFor();
87 Assert.assertEquals(expectedsFor, actualsFor, 0.05);
88
89 double expectedsBack = 0.3552;
90 double actualsBack = bar.getsBack();
91 Assert.assertEquals(expectedsBack, actualsBack, 0.001);
92
93 double expectedsBAR = 0.2280;
94 double actualsBAR = bar.getsBAR();
95 Assert.assertEquals(expectedsBAR, actualsBAR, 0.01);
96
97 double expectedBARIteration = -53.3234;
98 double actualBARIteration = bar.getBarEnergy();
99 Assert.assertEquals(expectedBARIteration, actualBARIteration, 0.1);
100
101 double expectedBARBootstrap = -53.3224;
102 double actualBARBootstrap = bar.getBarEnergyBS();
103 Assert.assertEquals(expectedBARBootstrap, actualBARBootstrap, 0.2);
104 }
105
106 @Test
107 public void testBARHelp() {
108 String[] args = {"-h"};
109 binding.setVariable("args", args);
110
111
112 BAR bar = new BAR(binding).run();
113 algorithmsScript = bar;
114 }
115
116 }