Skip to content
Snippets Groups Projects
Commit f8cd42e2 authored by Giovanni Bussi's avatar Giovanni Bussi
Browse files

Regtest

Some test for the leptonx2 optimization.
Also some test for different functions.
parent 77340130
No related branches found
No related tags found
No related merge requests found
#! FIELDS time c1 cs1
0.000000 1.000000 1.000000
1.000000 1.000000 1.000000
2.000000 0.999987 0.999987
3.000000 0.999849 0.999847
4.000000 0.999152 0.999139
5.000000 0.996773 0.996723
6.000000 0.990427 0.990277
7.000000 0.976206 0.975834
8.000000 0.948488 0.947683
9.000000 0.900818 0.899268
10.000000 0.828379 0.825698
11.000000 0.731515 0.727320
12.000000 0.617808 0.611837
13.000000 0.500000 0.492188
14.000000 0.390634 0.381113
15.000000 0.297631 0.286656
16.000000 0.223422 0.211288
17.000000 0.166646 0.153625
18.000000 0.124277 0.110594
19.000000 0.093051 0.078880
20.000000 0.070130 0.055601
21.000000 0.053280 0.038488
22.000000 0.040834 0.025847
23.000000 0.031576 0.016444
24.000000 0.024635 0.009395
25.000000 0.019387 0.004065
26.000000 0.000000 0.000000
27.000000 0.000000 0.000000
28.000000 0.000000 0.000000
29.000000 0.000000 0.000000
#! FIELDS time parameter c1 cs1
0.000000 0 0.0000 0.0000
0.000000 1 0.0000 0.0000
0.000000 2 0.0000 0.0000
0.000000 3 0.0000 0.0000
0.000000 4 0.0000 0.0000
0.000000 5 0.0000 0.0000
0.000000 6 0.0000 0.0000
0.000000 7 0.0000 0.0000
0.000000 8 0.0000 0.0000
0.000000 9 0.0000 0.0000
0.000000 10 0.0000 0.0000
0.000000 11 0.0000 0.0000
0.000000 12 0.0000 0.0000
0.000000 13 0.0000 0.0000
0.000000 14 0.0000 0.0000
1.000000 0 0.0000 0.0000
1.000000 1 0.0000 0.0000
1.000000 2 0.0000 0.0000
1.000000 3 0.0000 0.0000
1.000000 4 0.0000 0.0000
1.000000 5 -0.0000 -0.0000
1.000000 6 0.0000 0.0000
1.000000 7 0.0000 0.0000
1.000000 8 0.0000 0.0000
1.000000 9 0.0000 0.0000
1.000000 10 0.0000 0.0000
1.000000 11 0.0000 0.0000
1.000000 12 0.0000 0.0000
1.000000 13 0.0000 0.0000
1.000000 14 0.0000 0.0000
2.000000 0 0.0000 0.0000
2.000000 1 0.0000 0.0000
2.000000 2 0.0004 0.0004
2.000000 3 0.0000 0.0000
2.000000 4 0.0000 0.0000
2.000000 5 -0.0004 -0.0004
2.000000 6 0.0000 0.0000
2.000000 7 0.0000 0.0000
2.000000 8 0.0000 0.0000
2.000000 9 0.0000 0.0000
2.000000 10 0.0000 0.0000
2.000000 11 0.0000 0.0000
2.000000 12 0.0000 0.0000
2.000000 13 0.0000 0.0000
2.000000 14 0.0001 0.0001
3.000000 0 0.0000 0.0000
3.000000 1 0.0000 0.0000
3.000000 2 0.0030 0.0031
3.000000 3 0.0000 0.0000
3.000000 4 0.0000 0.0000
3.000000 5 -0.0030 -0.0031
3.000000 6 0.0000 0.0000
3.000000 7 0.0000 0.0000
3.000000 8 0.0000 0.0000
3.000000 9 0.0000 0.0000
3.000000 10 0.0000 0.0000
3.000000 11 0.0000 0.0000
3.000000 12 0.0000 0.0000
3.000000 13 0.0000 0.0000
3.000000 14 0.0009 0.0009
4.000000 0 0.0000 0.0000
4.000000 1 0.0000 0.0000
4.000000 2 0.0127 0.0129
4.000000 3 0.0000 0.0000
4.000000 4 0.0000 0.0000
4.000000 5 -0.0127 -0.0129
4.000000 6 0.0000 0.0000
4.000000 7 0.0000 0.0000
4.000000 8 0.0000 0.0000
4.000000 9 0.0000 0.0000
4.000000 10 0.0000 0.0000
4.000000 11 0.0000 0.0000
4.000000 12 0.0000 0.0000
4.000000 13 0.0000 0.0000
4.000000 14 0.0051 0.0052
5.000000 0 0.0000 0.0000
5.000000 1 0.0000 0.0000
5.000000 2 0.0386 0.0392
5.000000 3 0.0000 0.0000
5.000000 4 0.0000 0.0000
5.000000 5 -0.0386 -0.0392
5.000000 6 0.0000 0.0000
5.000000 7 0.0000 0.0000
5.000000 8 0.0000 0.0000
5.000000 9 0.0000 0.0000
5.000000 10 0.0000 0.0000
5.000000 11 0.0000 0.0000
5.000000 12 0.0000 0.0000
5.000000 13 0.0000 0.0000
5.000000 14 0.0193 0.0196
6.000000 0 0.0000 0.0000
6.000000 1 0.0000 0.0000
6.000000 2 0.0948 0.0963
6.000000 3 0.0000 0.0000
6.000000 4 0.0000 0.0000
6.000000 5 -0.0948 -0.0963
6.000000 6 0.0000 0.0000
6.000000 7 0.0000 0.0000
6.000000 8 0.0000 0.0000
6.000000 9 0.0000 0.0000
6.000000 10 0.0000 0.0000
6.000000 11 0.0000 0.0000
6.000000 12 0.0000 0.0000
6.000000 13 0.0000 0.0000
6.000000 14 0.0569 0.0578
7.000000 0 0.0000 0.0000
7.000000 1 0.0000 0.0000
7.000000 2 0.1991 0.2022
7.000000 3 0.0000 0.0000
7.000000 4 0.0000 0.0000
7.000000 5 -0.1991 -0.2022
7.000000 6 0.0000 0.0000
7.000000 7 0.0000 0.0000
7.000000 8 0.0000 0.0000
7.000000 9 0.0000 0.0000
7.000000 10 0.0000 0.0000
7.000000 11 0.0000 0.0000
7.000000 12 0.0000 0.0000
7.000000 13 0.0000 0.0000
7.000000 14 0.1394 0.1415
8.000000 0 0.0000 0.0000
8.000000 1 0.0000 0.0000
8.000000 2 0.3664 0.3722
8.000000 3 0.0000 0.0000
8.000000 4 0.0000 0.0000
8.000000 5 -0.3664 -0.3722
8.000000 6 0.0000 0.0000
8.000000 7 0.0000 0.0000
8.000000 8 0.0000 0.0000
8.000000 9 0.0000 0.0000
8.000000 10 0.0000 0.0000
8.000000 11 0.0000 0.0000
8.000000 12 0.0000 0.0000
8.000000 13 0.0000 0.0000
8.000000 14 0.2932 0.2977
9.000000 0 0.0000 0.0000
9.000000 1 0.0000 0.0000
9.000000 2 0.5956 0.6049
9.000000 3 0.0000 0.0000
9.000000 4 0.0000 0.0000
9.000000 5 -0.5956 -0.6049
9.000000 6 0.0000 0.0000
9.000000 7 0.0000 0.0000
9.000000 8 0.0000 0.0000
9.000000 9 0.0000 0.0000
9.000000 10 0.0000 0.0000
9.000000 11 0.0000 0.0000
9.000000 12 0.0000 0.0000
9.000000 13 0.0000 0.0000
9.000000 14 0.5361 0.5444
10.000000 0 0.0000 0.0000
10.000000 1 0.0000 0.0000
10.000000 2 0.8530 0.8663
10.000000 3 0.0000 0.0000
10.000000 4 0.0000 0.0000
10.000000 5 -0.8530 -0.8663
10.000000 6 0.0000 0.0000
10.000000 7 0.0000 0.0000
10.000000 8 0.0000 0.0000
10.000000 9 0.0000 0.0000
10.000000 10 0.0000 0.0000
10.000000 11 0.0000 0.0000
10.000000 12 0.0000 0.0000
10.000000 13 0.0000 0.0000
10.000000 14 0.8530 0.8663
11.000000 0 0.0000 0.0000
11.000000 1 0.0000 0.0000
11.000000 2 1.0713 1.0880
11.000000 3 0.0000 0.0000
11.000000 4 0.0000 0.0000
11.000000 5 -1.0713 -1.0880
11.000000 6 0.0000 0.0000
11.000000 7 0.0000 0.0000
11.000000 8 0.0000 0.0000
11.000000 9 0.0000 0.0000
11.000000 10 0.0000 0.0000
11.000000 11 0.0000 0.0000
11.000000 12 0.0000 0.0000
11.000000 13 0.0000 0.0000
11.000000 14 1.1784 1.1968
12.000000 0 0.0000 0.0000
12.000000 1 0.0000 0.0000
12.000000 2 1.1806 1.1991
12.000000 3 0.0000 0.0000
12.000000 4 0.0000 0.0000
12.000000 5 -1.1806 -1.1991
12.000000 6 0.0000 0.0000
12.000000 7 0.0000 0.0000
12.000000 8 0.0000 0.0000
12.000000 9 0.0000 0.0000
12.000000 10 0.0000 0.0000
12.000000 11 0.0000 0.0000
12.000000 12 0.0000 0.0000
12.000000 13 0.0000 0.0000
12.000000 14 1.4167 1.4389
13.000000 0 0.0000 0.0000
13.000000 1 0.0000 0.0000
13.000000 2 1.1538 1.1719
13.000000 3 0.0000 0.0000
13.000000 4 0.0000 0.0000
13.000000 5 -1.1538 -1.1719
13.000000 6 0.0000 0.0000
13.000000 7 0.0000 0.0000
13.000000 8 0.0000 0.0000
13.000000 9 0.0000 0.0000
13.000000 10 0.0000 0.0000
13.000000 11 0.0000 0.0000
13.000000 12 0.0000 0.0000
13.000000 13 0.0000 0.0000
13.000000 14 1.5000 1.5234
14.000000 0 0.0000 0.0000
14.000000 1 0.0000 0.0000
14.000000 2 1.0202 1.0361
14.000000 3 0.0000 0.0000
14.000000 4 0.0000 0.0000
14.000000 5 -1.0202 -1.0361
14.000000 6 0.0000 0.0000
14.000000 7 0.0000 0.0000
14.000000 8 0.0000 0.0000
14.000000 9 0.0000 0.0000
14.000000 10 0.0000 0.0000
14.000000 11 0.0000 0.0000
14.000000 12 0.0000 0.0000
14.000000 13 0.0000 0.0000
14.000000 14 1.4282 1.4506
15.000000 0 0.0000 0.0000
15.000000 1 0.0000 0.0000
15.000000 2 0.8362 0.8493
15.000000 3 0.0000 0.0000
15.000000 4 0.0000 0.0000
15.000000 5 -0.8362 -0.8493
15.000000 6 0.0000 0.0000
15.000000 7 0.0000 0.0000
15.000000 8 0.0000 0.0000
15.000000 9 0.0000 0.0000
15.000000 10 0.0000 0.0000
15.000000 11 0.0000 0.0000
15.000000 12 0.0000 0.0000
15.000000 13 0.0000 0.0000
15.000000 14 1.2543 1.2739
16.000000 0 0.0000 0.0000
16.000000 1 0.0000 0.0000
16.000000 2 0.6506 0.6608
16.000000 3 0.0000 0.0000
16.000000 4 0.0000 0.0000
16.000000 5 -0.6506 -0.6608
16.000000 6 0.0000 0.0000
16.000000 7 0.0000 0.0000
16.000000 8 0.0000 0.0000
16.000000 9 0.0000 0.0000
16.000000 10 0.0000 0.0000
16.000000 11 0.0000 0.0000
16.000000 12 0.0000 0.0000
16.000000 13 0.0000 0.0000
16.000000 14 1.0410 1.0573
17.000000 0 0.0000 0.0000
17.000000 1 0.0000 0.0000
17.000000 2 0.4901 0.4978
17.000000 3 0.0000 0.0000
17.000000 4 0.0000 0.0000
17.000000 5 -0.4901 -0.4978
17.000000 6 0.0000 0.0000
17.000000 7 0.0000 0.0000
17.000000 8 0.0000 0.0000
17.000000 9 0.0000 0.0000
17.000000 10 0.0000 0.0000
17.000000 11 0.0000 0.0000
17.000000 12 0.0000 0.0000
17.000000 13 0.0000 0.0000
17.000000 14 0.8333 0.8463
18.000000 0 0.0000 0.0000
18.000000 1 0.0000 0.0000
18.000000 2 0.3628 0.3684
18.000000 3 0.0000 0.0000
18.000000 4 0.0000 0.0000
18.000000 5 -0.3628 -0.3684
18.000000 6 0.0000 0.0000
18.000000 7 0.0000 0.0000
18.000000 8 0.0000 0.0000
18.000000 9 0.0000 0.0000
18.000000 10 0.0000 0.0000
18.000000 11 0.0000 0.0000
18.000000 12 0.0000 0.0000
18.000000 13 0.0000 0.0000
18.000000 14 0.6530 0.6632
19.000000 0 0.0000 0.0000
19.000000 1 0.0000 0.0000
19.000000 2 0.2665 0.2707
19.000000 3 0.0000 0.0000
19.000000 4 0.0000 0.0000
19.000000 5 -0.2665 -0.2707
19.000000 6 0.0000 0.0000
19.000000 7 0.0000 0.0000
19.000000 8 0.0000 0.0000
19.000000 9 0.0000 0.0000
19.000000 10 0.0000 0.0000
19.000000 11 0.0000 0.0000
19.000000 12 0.0000 0.0000
19.000000 13 0.0000 0.0000
19.000000 14 0.5064 0.5143
20.000000 0 0.0000 0.0000
20.000000 1 0.0000 0.0000
20.000000 2 0.1956 0.1987
20.000000 3 0.0000 0.0000
20.000000 4 0.0000 0.0000
20.000000 5 -0.1956 -0.1987
20.000000 6 0.0000 0.0000
20.000000 7 0.0000 0.0000
20.000000 8 0.0000 0.0000
20.000000 9 0.0000 0.0000
20.000000 10 0.0000 0.0000
20.000000 11 0.0000 0.0000
20.000000 12 0.0000 0.0000
20.000000 13 0.0000 0.0000
20.000000 14 0.3913 0.3974
21.000000 0 0.0000 0.0000
21.000000 1 0.0000 0.0000
21.000000 2 0.1441 0.1464
21.000000 3 0.0000 0.0000
21.000000 4 0.0000 0.0000
21.000000 5 -0.1441 -0.1464
21.000000 6 0.0000 0.0000
21.000000 7 0.0000 0.0000
21.000000 8 0.0000 0.0000
21.000000 9 0.0000 0.0000
21.000000 10 0.0000 0.0000
21.000000 11 0.0000 0.0000
21.000000 12 0.0000 0.0000
21.000000 13 0.0000 0.0000
21.000000 14 0.3026 0.3074
22.000000 0 0.0000 0.0000
22.000000 1 0.0000 0.0000
22.000000 2 0.1068 0.1085
22.000000 3 0.0000 0.0000
22.000000 4 0.0000 0.0000
22.000000 5 -0.1068 -0.1085
22.000000 6 0.0000 0.0000
22.000000 7 0.0000 0.0000
22.000000 8 0.0000 0.0000
22.000000 9 0.0000 0.0000
22.000000 10 0.0000 0.0000
22.000000 11 0.0000 0.0000
22.000000 12 0.0000 0.0000
22.000000 13 0.0000 0.0000
22.000000 14 0.2350 0.2387
23.000000 0 0.0000 0.0000
23.000000 1 0.0000 0.0000
23.000000 2 0.0798 0.0810
23.000000 3 0.0000 0.0000
23.000000 4 0.0000 0.0000
23.000000 5 -0.0798 -0.0810
23.000000 6 0.0000 0.0000
23.000000 7 0.0000 0.0000
23.000000 8 0.0000 0.0000
23.000000 9 0.0000 0.0000
23.000000 10 0.0000 0.0000
23.000000 11 0.0000 0.0000
23.000000 12 0.0000 0.0000
23.000000 13 0.0000 0.0000
23.000000 14 0.1835 0.1863
24.000000 0 0.0000 0.0000
24.000000 1 0.0000 0.0000
24.000000 2 0.0601 0.0610
24.000000 3 0.0000 0.0000
24.000000 4 0.0000 0.0000
24.000000 5 -0.0601 -0.0610
24.000000 6 0.0000 0.0000
24.000000 7 0.0000 0.0000
24.000000 8 0.0000 0.0000
24.000000 9 0.0000 0.0000
24.000000 10 0.0000 0.0000
24.000000 11 0.0000 0.0000
24.000000 12 0.0000 0.0000
24.000000 13 0.0000 0.0000
24.000000 14 0.1442 0.1464
25.000000 0 0.0000 0.0000
25.000000 1 0.0000 0.0000
25.000000 2 0.0456 0.0463
25.000000 3 0.0000 0.0000
25.000000 4 0.0000 0.0000
25.000000 5 -0.0456 -0.0463
25.000000 6 0.0000 0.0000
25.000000 7 0.0000 0.0000
25.000000 8 0.0000 0.0000
25.000000 9 0.0000 0.0000
25.000000 10 0.0000 0.0000
25.000000 11 0.0000 0.0000
25.000000 12 0.0000 0.0000
25.000000 13 0.0000 0.0000
25.000000 14 0.1141 0.1159
26.000000 0 0.0000 0.0000
26.000000 1 0.0000 0.0000
26.000000 2 0.0000 0.0000
26.000000 3 0.0000 0.0000
26.000000 4 0.0000 0.0000
26.000000 5 0.0000 0.0000
26.000000 6 0.0000 0.0000
26.000000 7 0.0000 0.0000
26.000000 8 0.0000 0.0000
26.000000 9 0.0000 0.0000
26.000000 10 0.0000 0.0000
26.000000 11 0.0000 0.0000
26.000000 12 0.0000 0.0000
26.000000 13 0.0000 0.0000
26.000000 14 0.0000 0.0000
27.000000 0 0.0000 0.0000
27.000000 1 0.0000 0.0000
27.000000 2 0.0000 0.0000
27.000000 3 0.0000 0.0000
27.000000 4 0.0000 0.0000
27.000000 5 0.0000 0.0000
27.000000 6 0.0000 0.0000
27.000000 7 0.0000 0.0000
27.000000 8 0.0000 0.0000
27.000000 9 0.0000 0.0000
27.000000 10 0.0000 0.0000
27.000000 11 0.0000 0.0000
27.000000 12 0.0000 0.0000
27.000000 13 0.0000 0.0000
27.000000 14 0.0000 0.0000
28.000000 0 0.0000 0.0000
28.000000 1 0.0000 0.0000
28.000000 2 0.0000 0.0000
28.000000 3 0.0000 0.0000
28.000000 4 0.0000 0.0000
28.000000 5 0.0000 0.0000
28.000000 6 0.0000 0.0000
28.000000 7 0.0000 0.0000
28.000000 8 0.0000 0.0000
28.000000 9 0.0000 0.0000
28.000000 10 0.0000 0.0000
28.000000 11 0.0000 0.0000
28.000000 12 0.0000 0.0000
28.000000 13 0.0000 0.0000
28.000000 14 0.0000 0.0000
29.000000 0 0.0000 0.0000
29.000000 1 0.0000 0.0000
29.000000 2 0.0000 0.0000
29.000000 3 0.0000 0.0000
29.000000 4 0.0000 0.0000
29.000000 5 0.0000 0.0000
29.000000 6 0.0000 0.0000
29.000000 7 0.0000 0.0000
29.000000 8 0.0000 0.0000
29.000000 9 0.0000 0.0000
29.000000 10 0.0000 0.0000
29.000000 11 0.0000 0.0000
29.000000 12 0.0000 0.0000
29.000000 13 0.0000 0.0000
29.000000 14 0.0000 0.0000
......@@ -8,5 +8,16 @@ PRINT ...
DUMPDERIVATIVES ARG=c,cs FILE=deriv FMT=%8.4f
# this is just to make sure that stretching is correct with the optimized x2 implementation
c1: COORDINATION GROUPA=1 GROUPB=2 SWITCH={MATHEVAL R_0=1.3 D_MAX=2.6 FUNC=1/(1+x2^3) NOSTRETCH}
cs1: COORDINATION GROUPA=1 GROUPB=2 SWITCH={CUSTOM R_0=1.3 D_MAX=2.6 FUNC=1/(1+x2^3)}
PRINT ...
ARG=c1,cs1
FILE=COLVAR1
... PRINT
DUMPDERIVATIVES ARG=c1,cs1 FILE=deriv1 FMT=%8.4f
ENDPLUMED
#! FIELDS time c d dfast e f g h
0.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
1.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
2.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
3.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
4.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
5.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
6.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
7.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
8.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
9.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
10.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
11.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
12.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
13.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
14.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
15.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
16.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
17.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
18.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
19.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
20.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
21.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
22.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
23.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
24.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
25.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
26.000000 171.1815 171.2953 171.2953 182.5448 182.6585 857.1425 857.1425
27.000000 11555.3194 11555.3194 11555.3194 11555.9968 11555.9968 9813.4170 9813.4170
28.000000 11555.9993 11555.9993 11555.9993 11556.0000 11556.0000 10952.2272 10952.2272
29.000000 11555.7324 11555.7324 11555.7324 11555.9995 11555.9995 10083.8540 10083.8540
30.000000 11555.9999 11555.9999 11555.9999 11556.0000 11556.0000 11126.5294 11126.5294
31.000000 11551.3675 11551.3676 11551.3676 11555.8899 11555.8899 9262.6433 9262.6433
32.000000 11554.8705 11554.8705 11554.8705 11555.9900 11555.9900 9704.2327 9704.2327
33.000000 11555.8052 11555.8052 11555.8052 11555.9998 11555.9998 10036.1770 10036.1770
34.000000 11553.8594 11553.8594 11553.8594 11555.9738 11555.9738 9431.1577 9431.1577
35.000000 11555.9485 11555.9485 11555.9485 11556.0000 11556.0000 10418.1634 10418.1634
36.000000 11553.9347 11553.9347 11553.9347 11555.9756 11555.9756 9466.8778 9466.8778
37.000000 11554.7972 11554.7972 11554.7972 11555.9914 11555.9914 9641.7889 9641.7889
38.000000 11555.8020 11555.8020 11555.8020 11555.9998 11555.9998 10035.9562 10035.9562
39.000000 11555.4367 11555.4367 11555.4367 11555.9983 11555.9983 9790.0490 9790.0490
40.000000 11555.9104 11555.9104 11555.9104 11556.0000 11556.0000 10230.8314 10230.8314
41.000000 11555.6978 11555.6978 11555.6978 11555.9995 11555.9995 9969.0934 9969.0934
42.000000 11548.9710 11548.9711 11548.9711 11555.8145 11555.8145 8973.7830 8973.7830
43.000000 11555.5133 11555.5133 11555.5133 11555.9984 11555.9984 9896.2166 9896.2166
44.000000 11555.4781 11555.4781 11555.4781 11555.9985 11555.9985 9788.6446 9788.6446
45.000000 11555.4807 11555.4807 11555.4807 11555.9980 11555.9980 9886.2067 9886.2067
46.000000 11555.1715 11555.1716 11555.1716 11555.9951 11555.9951 9746.5077 9746.5077
47.000000 11555.0278 11555.0279 11555.0279 11555.9938 11555.9938 9663.8156 9663.8156
48.000000 11555.8818 11555.8818 11555.8818 11555.9999 11555.9999 10206.8420 10206.8420
49.000000 11550.7984 11550.7984 11550.7984 11555.8733 11555.8733 9243.9690 9243.9690
50.000000 11552.3158 11552.3158 11552.3158 11555.9456 11555.9456 9153.2704 9153.2704
51.000000 11555.9024 11555.9024 11555.9024 11556.0000 11556.0000 10294.5076 10294.5076
52.000000 11553.8360 11553.8360 11553.8360 11555.9785 11555.9785 9341.5602 9341.5602
53.000000 11550.5747 11550.5748 11550.5748 11555.8742 11555.8742 9101.5537 9101.5537
54.000000 11555.8688 11555.8688 11555.8688 11555.9999 11555.9999 10148.0927 10148.0927
55.000000 11555.9098 11555.9098 11555.9098 11556.0000 11556.0000 10260.4494 10260.4494
56.000000 11555.9968 11555.9968 11555.9968 11556.0000 11556.0000 10772.3540 10772.3540
57.000000 11555.9475 11555.9475 11555.9475 11556.0000 11556.0000 10328.1273 10328.1273
58.000000 11555.2189 11555.2190 11555.2190 11555.9956 11555.9956 9797.4731 9797.4731
59.000000 11554.3484 11554.3485 11554.3485 11555.9833 11555.9833 9545.7559 9545.7559
60.000000 11555.9989 11555.9989 11555.9989 11556.0000 11556.0000 10908.7613 10908.7613
61.000000 11554.1871 11554.1871 11554.1871 11555.9787 11555.9787 9479.4950 9479.4950
62.000000 11553.9961 11553.9961 11553.9961 11555.9772 11555.9772 9441.4599 9441.4599
63.000000 11553.1389 11553.1389 11553.1389 11555.9552 11555.9552 9341.9431 9341.9431
64.000000 11554.0967 11554.0968 11554.0968 11555.9775 11555.9775 9529.9313 9529.9313
65.000000 11555.9954 11555.9954 11555.9954 11556.0000 11556.0000 10722.0256 10722.0256
66.000000 11555.6342 11555.6342 11555.6342 11555.9992 11555.9992 9924.2156 9924.2156
67.000000 11553.6744 11553.6744 11553.6744 11555.9706 11555.9706 9387.7469 9387.7469
68.000000 11555.7428 11555.7428 11555.7428 11555.9997 11555.9997 9968.6967 9968.6967
69.000000 11552.3911 11552.3911 11552.3911 11555.9355 11555.9355 9240.7544 9240.7544
70.000000 11555.9997 11555.9997 11555.9997 11556.0000 11556.0000 11018.4173 11018.4173
71.000000 11555.0612 11555.0612 11555.0612 11555.9936 11555.9936 9740.3220 9740.3220
72.000000 11555.9224 11555.9224 11555.9224 11556.0000 11556.0000 10349.6091 10349.6091
73.000000 11555.9231 11555.9231 11555.9231 11556.0000 11556.0000 10299.3941 10299.3941
74.000000 11555.7558 11555.7558 11555.7558 11555.9997 11555.9997 9982.0367 9982.0367
75.000000 11551.0097 11551.0098 11551.0098 11555.8819 11555.8819 9173.7783 9173.7783
76.000000 11554.1278 11554.1279 11554.1279 11555.9794 11555.9794 9509.6803 9509.6803
77.000000 11555.7544 11555.7544 11555.7544 11555.9997 11555.9997 10021.2115 10021.2115
78.000000 11553.2694 11553.2694 11553.2694 11555.9644 11555.9644 9388.2415 9388.2415
79.000000 11555.9087 11555.9087 11555.9087 11556.0000 11556.0000 10253.3549 10253.3549
80.000000 11555.2733 11555.2733 11555.2733 11555.9964 11555.9964 9793.0432 9793.0432
81.000000 11554.4183 11554.4184 11554.4184 11555.9865 11555.9865 9492.8793 9492.8793
82.000000 11554.0377 11554.0377 11554.0377 11555.9786 11555.9786 9505.6322 9505.6322
83.000000 11555.8667 11555.8667 11555.8667 11555.9999 11555.9999 10172.6161 10172.6161
84.000000 11555.9784 11555.9784 11555.9784 11556.0000 11556.0000 10476.2075 10476.2075
85.000000 11551.7929 11551.7930 11551.7930 11555.9225 11555.9225 9164.1025 9164.1025
86.000000 11555.3497 11555.3497 11555.3497 11555.9976 11555.9976 9734.5934 9734.5934
87.000000 11555.5810 11555.5810 11555.5810 11555.9991 11555.9991 9854.4028 9854.4028
88.000000 11554.7942 11554.7942 11554.7942 11555.9882 11555.9882 9788.3750 9788.3750
89.000000 11555.2354 11555.2354 11555.2354 11555.9969 11555.9969 9676.7742 9676.7742
90.000000 11554.4651 11554.4651 11554.4651 11555.9841 11555.9841 9605.0962 9605.0962
91.000000 11555.2909 11555.2909 11555.2909 11555.9971 11555.9971 9710.0986 9710.0986
92.000000 11555.8838 11555.8838 11555.8838 11556.0000 11556.0000 10164.7871 10164.7871
93.000000 11554.7613 11554.7613 11554.7613 11555.9908 11555.9908 9563.8592 9563.8592
94.000000 11555.8958 11555.8958 11555.8958 11555.9999 11555.9999 10340.5430 10340.5430
95.000000 11555.8504 11555.8504 11555.8504 11555.9999 11555.9999 10188.6328 10188.6328
96.000000 11553.4806 11553.4806 11553.4806 11555.9709 11555.9709 9306.7080 9306.7080
97.000000 11554.3927 11554.3927 11554.3927 11555.9860 11555.9860 9488.4639 9488.4639
98.000000 11553.9367 11553.9367 11553.9367 11555.9757 11555.9757 9486.2964 9486.2964
99.000000 11555.8015 11555.8015 11555.8015 11555.9998 11555.9998 10061.1075 10061.1075
100.000000 11555.4155 11555.4155 11555.4155 11555.9978 11555.9978 9858.3232 9858.3232
101.000000 11555.2815 11555.2815 11555.2815 11555.9968 11555.9968 9727.2081 9727.2081
102.000000 11549.7314 11549.7315 11549.7315 11555.8361 11555.8361 9060.2674 9060.2674
103.000000 11555.9984 11555.9984 11555.9984 11556.0000 11556.0000 10858.0931 10858.0931
104.000000 11555.3982 11555.3982 11555.3982 11555.9979 11555.9979 9756.6344 9756.6344
105.000000 11555.5026 11555.5027 11555.5027 11555.9983 11555.9983 9887.7851 9887.7851
106.000000 11551.1266 11551.1266 11551.1266 11555.9018 11555.9018 9097.0913 9097.0913
107.000000 11555.7818 11555.7818 11555.7818 11555.9997 11555.9997 10098.5292 10098.5292
......@@ -5,3 +5,21 @@ PRINT ARG=c,d FILE=colvar FMT=%8.4f STRIDE=1
FLUSH STRIDE=1
s: COMBINE ARG=c,d PERIODIC=NO
RESTRAINT ARG=s KAPPA=0 SLOPE=1 AT=0
dfast: COORDINATION GROUPA=1-108 GROUPB=1-108 SWITCH={CUSTOM FUNC=1/(1+x2^3) R_0=1}
e: COORDINATION GROUPA=1-108 GROUPB=1-108 R_0=1 NN=7 D_0=0.1
f: COORDINATION GROUPA=1-108 GROUPB=1-108 SWITCH={CUSTOM FUNC=1/(1+x^7) R_0=1 D_0=0.1}
g: COORDINATION GROUPA=1-108 GROUPB=1-108 SWITCH={EXP R_0=1}
h: COORDINATION GROUPA=1-108 GROUPB=1-108 SWITCH={CUSTOM FUNC=exp(-x) R_0=1}
PRINT ARG=c,d,dfast,e,f,g,h FILE=check1 FMT=%8.4f
# this would check derivatives:
# DUMPDERIVATIVES ARG=c,d,dfast FILE=der0 FMT=%8.4f
# DUMPDERIVATIVES ARG=e,f FILE=der1 FMT=%8.4f
# DUMPDERIVATIVES ARG=g,h FILE=der2 FMT=%8.4f
# to detect errors we just add and subtract biases
# this should add zero forces
all: COMBINE ARG=c,d,dfast,e,f,g,h COEFFICIENTS=1,-0.5,-0.5,1,-1,1,-1 PERIODIC=NO
RESTRAINT ARG=all KAPPA=0 SLOPE=0.001 AT=0
# SLOPE is very small in order to avoid numerical errors to appear
# it should be sufficient to detect errors
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment