# HG changeset patch # User chemteam # Date 1599861357 0 # Node ID 6f2029791c946a16246c6d2438174b3ef608d21b # Parent d9f8cc3258f93524a5ca29565ea0902ee33bbf60 "planemo upload for repository https://github.com/galaxycomputationalchemistry/galaxy-tools-compchem/ commit 79589d149a8ff2791d4f71d28b155011672db827" diff -r d9f8cc3258f9 -r 6f2029791c94 NEQGamma.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/NEQGamma.py Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,210 @@ +#!/usr/bin/env python + +# coding: utf-8 +# This script is a modified version of a script written +# by Steffen Wolf under the GPL v3.0. +# The original version can be accessed at +# https://github.com/moldyn/dcTMD/blob/master/NEQGamma.py + +import argparse +import json +import sys + +import numpy as np + +import pandas as pd + +import scipy +import scipy.integrate +from scipy.ndimage.filters import gaussian_filter + + +def get_file_names(list_file): + with open(list_file) as f: + return [line for line in f.read().split('\n') if line] + + +def NEQGamma(file_names, output, output_frict, vel, T, av, sigma): + N = len(file_names) + RT = 0.0083144598 * T + + sys.stdout.write("reading data...\n") + + # read in initial data to get length of necessary array + test_file = pd.read_csv(file_names[0], delim_whitespace=True, + header=None, skiprows=17, dtype=float) + length_data = len(test_file[0].values) + full_force_set = np.zeros((N, length_data)) + x = np.zeros(length_data) + t = np.zeros(length_data) + t = test_file[0].values + x = test_file[0].values * vel + + # read in data + for i in range(0, N): + current_file_name = file_names[i] + sys.stdout.write("reading file {}\n".format(current_file_name)) + input_file_data = pd.read_csv(current_file_name, delim_whitespace=True, + header=None, skiprows=17, dtype=float) + full_force_set[i, :] = input_file_data[1].values + + # preprocessing + # * force average: calculate $\left< f_c (t)\right>_N$. + # **Important:** this is an ensemble average over the trajectory ensemble + # $N$, not the time average over $t$ + av_force = np.zeros(length_data) + av_forceintegral = np.zeros(length_data) + for i in range(length_data): + av_force[i] = np.mean(full_force_set[:, i]) + av_forceintegral[1:] = scipy.integrate.cumtrapz(av_force, x) + + # calculate $\delta f_c(t) = f_c(t) - \left< f_c (t) \right>_N$ for all $t$ + sys.stdout.write("calculating fluctuations...\n") + delta_force_set = np.zeros((N, length_data)) + for i in range(length_data): + delta_force_set[:, i] = full_force_set[:, i] - av_force[i] + + # evaluation + # * optimized algorithm for numerical evaluation: + # * integrate: $\int_0^t dt' \delta f_c(t')$ for all $t'$ + # * multiply by $\delta f_c(t)$ to yield + # $\int_0^t dt'\delta f_c(t) \delta f_c(t')$ for $t$ + # with all $t' \leq t$ each + # * then calculate the ensemble average + # $\left< \int_0^t dt' \delta f_c(t) \delta f_c(t') \right>$ + int_delta_force_set = np.zeros((N, length_data)) + for n in range(N): + int_delta_force_set[n, 1:] = scipy.integrate.cumtrapz( + delta_force_set[n, :], t) + + sys.stdout.write("averaging and integrating...\n") + intcorr = np.zeros((N, length_data)) + + for n in range(N): + for i in range(length_data): + intcorr[n, i] = delta_force_set[n, i] * int_delta_force_set[n, i] + if i % 1000 == 0: + sys.stdout.write("Trajectory {:2d} {:3.1f} % done\r".format( + n + 1, (i / length_data) * 100)) + + # shape of $\int_0^t dt' \delta f_c(t) \delta f_c(t')$: + sys.stdout.write("final average...\n") + av_intcorr = np.zeros(length_data) + for i in range(length_data): + av_intcorr[i] = np.mean(intcorr[:, i]) / RT + + # autocorrelation function evaluation: + # * calculate $\left< \delta f_c(t) \delta f_c(t') \right>$ + # for the last $t$ + + corr_set = np.zeros((N, length_data)) + autocorr_set = np.zeros(length_data) + + sys.stdout.write("calculating and processing ACF...\n") + for n in range(N): + for i in range(length_data): + corr_set[n, i] = delta_force_set[ + n, i] * delta_force_set[n, length_data - 1] + + for i in range(length_data): + autocorr_set[i] = np.mean(corr_set[:, i]) + + # * Gauss filter: + sys.stdout.write("applying Gauss filter...\n") + blurr = sigma + blurred = np.zeros(length_data) + blurred = gaussian_filter(av_intcorr, sigma=blurr) + + # * sliding window average: + sys.stdout.write("applying sliding window average...\n") + window = av + runn_av = np.zeros(length_data) + runn_av = np.convolve(av_intcorr, np.ones((window,)) / window, mode='same') + + # * $W_{diss}$ from integration: + wdiss = np.zeros(length_data) + wdiss[1:] = scipy.integrate.cumtrapz(av_intcorr, x) * vel + + sys.stdout.write("writing output...\n") + dist = open(output, "w") + frict = open(output_frict, "w") + + dist.write( + "#x force_integral frict_coeff wdiss corrected_force_integral\n") + for i in range(length_data): + dist.write("{:15.8f} {:20.8f} {:20.8f} {:20.8f} {:20.8f}\n".format( + x[i], av_forceintegral[i], av_intcorr[i], wdiss[i], + av_forceintegral[i] - wdiss[i])) + + frict.write("""#x ACF frict_coeff """ + """gauss_filtered_frict_coeff av_window_frict_coeff\n""") + for i in range(length_data): + frict.write("{:15.8f} {:20.8f} {:20.8f} {:20.8f} {:20.8f}\n".format( + x[i], autocorr_set[i], av_intcorr[i], blurred[i], runn_av[i])) + + dist.close() + frict.close() + + sys.stdout.write("Done!\n") + + +def main(): + parser = argparse.ArgumentParser(description="""dcTMD friciton correction + (please cite: Wolf, S., Stock, G. Targeted Molecular Dynamics + Calculations of Free Energy Profiles Using a Nonequilibrium + Friction Correction. J. Chem. Theory Comput. 2018, 14(12), 6175-6182, + DOI: 10.1021/acs.jctc.8b00835). Integrates a constraint force file via + trapezoid rule, calculates the NEQ memory friction kernel and friction + factors, and performs a friction correction. First column: reaction + coordinate in nm calculated via t * vel. Second column: force integral, + i.e. the work profile. Third column: friction factors. Fourth column: + trapezoid integral (final value) of friction work along reaction + coordinate. Fourth column: friction corrected work profile. ATTENTION: + Use with python3 or higher!""") + parser.add_argument('-i', metavar='', type=str, + help="""List of xvg constraint force files prefix + as given by Gromacs mdrun -pf option before running + number.""") + parser.add_argument('-o', metavar='', type=str, + help="""file to write x, dG(x), friction coefficeint by + integration (time resolved), and the friction-corrected + dG(x).""") + parser.add_argument('-ofrict', metavar='', + type=str, + help="""file to write x, ACF, friction coefficeint by + integration (time resolved), gauss filtered friction + coefficient, and slide window averaged friction.""") + parser.add_argument('-vel', metavar='', type=float, + help="""pull velocity in nm/ps for converting simulation + time t into distance x""") + parser.add_argument('-T', metavar='', type=float, + help='temperature in K') + parser.add_argument('-av', metavar='', type=int, + help="""size of averaging window for displaying + Gamma(x) (recommended: 4 to 20 per 100 data points)""") + parser.add_argument('-sigma', metavar='', type=int, + help="""sigma value for Gauss filter for displaying + Gamma(x) (recommended: 4 per 100 data points)""") + parser.add_argument('-json', metavar='', type=str, + help='JSON file defining cluster membership') + + args = parser.parse_args() + + file_names = get_file_names(args.i) + if args.json: + with open(args.json) as f: + j = json.load(f) + file_names_dct = {n: [file_names[int(m)] for m in j[n]] for n in j} + + for cluster in file_names_dct: + NEQGamma(file_names_dct[cluster], + 'cluster{}_{}'.format(cluster, args.o), + 'cluster{}_{}'.format(cluster, args.ofrict), + args.vel, args.T, args.av, args.sigma) + else: + NEQGamma(file_names, args.o, args.ofrict, args.vel, + args.T, args.av, args.sigma) + + +if __name__ == "__main__": + main() diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_clusters.json --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_clusters.json Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,1 @@ +{"0": ["0", "1"], "1": ["0", "2"]} diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_ofrict.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_ofrict.txt Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,102 @@ +#x ACF frict_coeff gauss_filtered_frict_coeff av_window_frict_coeff + 0.00000000 439653.22756622 0.00000000 39.85704130 -164.30975524 + 0.00000200 41861.51118060 -309.80490686 39.90128147 -42.50284457 + 0.00000400 -4760.58383990 -1433.31122481 39.98937507 132.75776205 + 0.00000600 -52063.97568219 -1302.42266399 40.10796653 273.24484877 + 0.00000800 -94137.43307581 -696.00283593 40.26787861 322.90057440 + 0.00001000 -109240.44374298 59.67135417 40.47438205 310.49706801 + 0.00001200 -91296.29474097 200.26534006 40.74119359 300.71687201 + 0.00001400 -42279.73264591 -211.00314767 41.08344474 307.77939795 + 0.00001600 29083.13222726 -299.86397145 41.51045417 289.59363364 + 0.00001800 100957.95082594 706.27695160 42.02391473 213.62523715 + 0.00002000 150492.32652407 2436.13821354 42.62185633 115.80351000 + 0.00002200 163610.52308989 3505.21213240 43.30185229 72.52474633 + 0.00002400 137798.21546909 2809.74173435 44.06395945 119.32042200 + 0.00002600 86552.97296922 993.11451262 44.91354630 219.23877027 + 0.00002800 30065.89230936 -248.07012791 45.85833909 214.16642537 + 0.00003000 -15271.83960767 -195.60392001 46.90428227 107.42336841 + 0.00003200 -41143.37576909 141.25051892 48.05595140 -10.23715919 + 0.00003400 -45603.75654834 -363.71528632 49.31891708 -41.13071129 + 0.00003600 -31933.14574422 -1519.36792976 50.67717643 11.13789839 + 0.00003800 -11580.24564270 -1956.43454302 52.14035913 47.86749694 + 0.00004000 1795.96943076 -865.57527338 53.69638462 -13.42971260 + 0.00004200 -6069.40567724 626.10860652 55.34660090 -151.70147237 + 0.00004400 -36605.44103603 565.05574074 57.08967355 -284.06117371 + 0.00004600 -88409.03954732 -1403.86956218 58.91185508 -352.46208669 + 0.00004800 -119187.64656330 -2830.86397506 60.79652934 -358.13736161 + 0.00005000 -126958.63383213 -2293.53919791 62.72537859 -329.02117947 + 0.00005200 -105144.23378691 -417.60570188 64.67974189 -286.64057459 + 0.00005400 -57578.18477311 834.36904591 66.64059843 -231.16675138 + 0.00005600 7642.76668942 434.72799961 68.58976595 -160.41908396 + 0.00005800 73314.58537446 -519.66723923 70.51136345 -97.11322986 + 0.00006000 117239.42596973 -329.29698187 72.38916035 -71.95507125 + 0.00006200 122335.84589389 858.01810556 74.20116994 -72.69735215 + 0.00006400 90700.59252948 1441.72347492 75.92201577 -50.06936146 + 0.00006600 40928.26316446 879.60901414 77.52654008 55.63484244 + 0.00006800 -3181.66487546 334.25351485 78.99192675 197.17614864 + 0.00007000 -31227.72623817 652.00817770 80.29696328 378.84255238 + 0.00007200 -43642.08135556 1250.72698300 81.42019653 507.21229626 + 0.00007400 -37578.91493866 1051.23806218 82.34208812 533.01944777 + 0.00007600 -9313.59450148 -253.25084775 83.04041098 497.57963010 + 0.00007800 35606.42432047 -1453.27137086 83.53076487 484.49556872 + 0.00008000 74825.06321773 -880.42089139 83.79365831 512.38603041 + 0.00008200 87450.95751560 1078.66842032 83.81028733 508.57269634 + 0.00008400 84964.54130532 2679.13981866 83.55962529 417.18682766 + 0.00008600 33723.38650302 1426.95656188 83.03964822 275.01069883 + 0.00008800 26888.00526112 802.46409967 82.26457780 159.07974622 + 0.00009000 15764.22367820 273.85567979 81.25039304 86.68872404 + 0.00009200 -1018.04050367 98.53732822 80.00383608 19.95207715 + 0.00009400 -18778.21955924 125.57269254 78.52437737 -55.26075096 + 0.00009600 -25696.85381088 173.04677203 76.81790266 -95.82294253 + 0.00009800 -19020.27963144 38.14199457 74.90534804 -40.08283716 + 0.00010000 -13853.36236626 -405.56366316 72.81359766 112.25327202 + 0.00010200 -22360.23253615 -969.69926808 70.56500733 263.50051564 + 0.00010400 -39525.81005506 -1401.79910164 68.18385377 320.37515633 + 0.00010600 -47887.23125912 -1439.01003806 65.70767268 323.43291333 + 0.00010800 -38703.50620624 -1113.56692890 63.18337712 284.23466280 + 0.00011000 -15756.47051989 -682.72476001 60.65377002 277.77044438 + 0.00011200 10148.75573217 -253.52957929 58.15278984 250.27515323 + 0.00011400 32489.76783313 239.99423089 55.71193773 178.79248872 + 0.00011600 48846.14416151 861.55125965 53.36798832 70.38368037 + 0.00011800 60211.15804332 1593.45081266 51.19169051 -55.24610237 + 0.00012000 64597.66151097 2144.52398093 49.22113114 -170.30467384 + 0.00012200 58795.32774769 2216.16123424 47.49918918 -248.93156307 + 0.00012400 63226.50291119 2740.29495857 46.05677747 -270.24081648 + 0.00012600 7112.67333129 642.99155129 44.93086250 -206.96775261 + 0.00012800 4737.67676382 673.17973142 44.16608448 -50.37244372 + 0.00013000 -10629.83066809 -276.05014319 43.80477020 154.04681427 + 0.00013200 -24445.08596782 -1331.11596207 43.87860095 320.32723893 + 0.00013400 -38089.44576760 -2042.60347448 44.40990189 378.66585808 + 0.00013600 -53888.48385453 -2339.54888278 45.42258092 320.31441637 + 0.00013800 -65525.93331796 -2263.02943488 46.94717499 192.76961767 + 0.00014000 -66605.80950981 -1978.10144773 49.01173919 46.49351043 + 0.00014200 -59461.87405093 -1395.88433631 51.62964197 -100.43225450 + 0.00014400 -50082.83777468 -136.33782409 54.79891493 -258.31822307 + 0.00014600 -36265.73051644 1692.89613965 58.50910094 -447.27992440 + 0.00014800 -11065.11218743 2974.81823090 62.74931610 -601.98154580 + 0.00015000 24574.77396976 2642.88373319 67.50830566 -716.60080679 + 0.00015200 53704.30727291 913.24280375 72.76914377 -754.65637852 + 0.00015400 62641.21173742 -927.03460327 78.50486320 -713.36174226 + 0.00015600 50475.22167026 -1689.34471443 84.67941051 -602.77004154 + 0.00015800 27415.48879178 -1332.07133205 91.25069147 -429.82867913 + 0.00016000 659.52027909 -793.99131765 98.16114225 -214.90982911 + 0.00016200 -27358.17056213 -941.55813718 105.35683648 -1.93703392 + 0.00016400 -36359.81639946 -1038.93906805 112.78930315 141.63166343 + 0.00016600 -103478.83628142 -2451.04087672 120.41222030 155.07119820 + 0.00016800 -98122.08801466 -1619.20548854 128.15466148 39.58301004 + 0.00017000 -91470.31369949 -1037.16157762 135.94071127 -128.72347744 + 0.00017200 -70167.76251890 -505.22323695 143.69516807 -240.17685720 + 0.00017400 -43298.45234300 169.23053990 151.35171103 -223.11938998 + 0.00017600 -13147.13005351 1119.27836540 158.84649112 -87.45006406 + 0.00017800 13812.88135432 2035.34756564 166.10612536 85.27603300 + 0.00018000 29405.90836664 2281.35445610 173.04701036 209.91332866 + 0.00018200 27785.78282904 1475.48961051 179.59291666 253.65425754 + 0.00018400 16400.24253711 132.45287134 185.68841711 300.73216440 + 0.00018600 5032.51740518 -616.86762340 191.28914198 352.67911780 + 0.00018800 -684.42414340 -391.31151887 196.34017758 475.23116164 + 0.00019000 -2606.74806622 413.81613798 200.77390679 556.19143606 + 0.00019200 -4228.75142469 1254.39214820 204.53041970 608.04951495 + 0.00019400 -2944.53596867 1786.35191510 207.57544804 633.31067679 + 0.00019600 4064.05842526 1765.17722675 209.89395091 624.84914980 + 0.00019800 11878.81314410 1160.67458116 211.46967661 568.88523153 + 0.00020000 12188.47862387 80.82726004 212.27097983 467.11785325 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_ofrict1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_ofrict1.txt Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,102 @@ +#x ACF frict_coeff gauss_filtered_frict_coeff av_window_frict_coeff + 0.00000000 -131993.41016925 0.00000000 -31.41159270 -71.62375840 + 0.00000200 16345.04469682 -1345.26646099 -31.23092179 16.85099383 + 0.00000400 29113.39120702 -1454.28761686 -30.88206666 43.12225366 + 0.00000600 27656.43517472 -264.14035591 -30.38447804 -24.61744170 + 0.00000800 22081.05883970 570.71230720 -29.75116213 -124.82229209 + 0.00001000 9926.50084650 482.67878751 -28.98362546 -148.88674802 + 0.00001200 -5957.22964270 -306.50010800 -28.07154561 -44.25804572 + 0.00001400 -21071.08253355 -678.79838871 -27.00311927 120.09671342 + 0.00001600 -30206.93805045 129.24578885 -25.77749754 210.63856614 + 0.00001800 -30292.41328536 1433.88087902 -24.40797512 147.07777675 + 0.00002000 -21089.95150080 1769.49504443 -22.91494801 -14.98005248 + 0.00002200 -5123.22129025 525.42519661 -21.31752461 -139.14697408 + 0.00002400 14086.19434495 -1354.79390711 -19.63003766 -71.97480155 + 0.00002600 31413.15161600 -2004.09700782 -17.86212814 78.50640859 + 0.00002800 42209.76106550 -481.28911857 -16.02210724 139.74809987 + 0.00003000 42689.48379010 2092.57404588 -14.11799868 71.98635054 + 0.00003200 32140.35686690 3287.09518295 -12.15586303 -24.51897016 + 0.00003400 13442.16919797 1810.83705441 -10.14236807 -39.10004376 + 0.00003600 -9231.10226605 -1271.21578790 -8.09221902 20.27669074 + 0.00003800 -29408.62646120 -3241.15658461 -6.02052658 45.55495337 + 0.00004000 -41627.31520850 -2483.33843193 -3.94447769 -27.79664184 + 0.00004200 -42137.13334250 -1.82301038 -1.86595819 -140.52649008 + 0.00004400 -30226.42411600 1555.33658588 0.21546106 -185.57828610 + 0.00004600 -11806.47895017 960.69346965 2.28847531 -122.35443556 + 0.00004800 9446.43397092 -784.52267942 4.34049778 -10.01090329 + 0.00005000 24611.29240972 -1447.42762641 6.35945779 66.73347194 + 0.00005200 30168.44061035 -598.12158006 8.33652167 81.49338936 + 0.00005400 25606.37528575 508.73630130 10.26729402 79.79530240 + 0.00005600 13340.52480087 634.81104146 12.15255270 108.71674708 + 0.00005800 -585.05665750 -33.15102509 13.99740074 159.07844440 + 0.00006000 -9852.88813967 -485.10192035 15.80514484 172.56136499 + 0.00006200 -10897.77796887 -375.61072392 17.57305023 116.90827164 + 0.00006400 -3757.41661050 -90.31689620 19.29877158 9.99710432 + 0.00006600 8790.19739798 242.77363760 20.98538739 -93.21873071 + 0.00006800 21442.96733022 1053.59838602 22.64082363 -139.65830590 + 0.00007000 28256.12418230 2387.77239428 24.27376738 -102.30945883 + 0.00007200 26367.06998670 3253.13344370 25.89205555 -17.58643540 + 0.00007400 15591.68871872 2389.26594789 27.50466561 70.86828133 + 0.00007600 -1617.68759163 -263.98184134 29.12475256 149.70059805 + 0.00007800 -20116.12291430 -2971.49817282 30.77024110 225.25149523 + 0.00008000 -32644.46446817 -3596.40029897 32.44887403 287.44322395 + 0.00008200 -34179.74570220 -2140.04635679 34.15733687 311.78904082 + 0.00008400 -24006.60727937 -508.98011481 35.89593384 298.54506507 + 0.00008600 -10185.53972520 31.90196581 37.68950778 269.11231530 + 0.00008800 6633.72365407 -37.54573797 39.58337067 228.28577067 + 0.00009000 19302.20585842 247.03284223 41.62183096 147.64655413 + 0.00009200 25999.47877012 1170.97275460 43.82919976 -6.36815521 + 0.00009400 25523.09074980 2085.38263558 46.21203223 -211.35490584 + 0.00009600 18932.09640609 2145.82898503 48.77454558 -376.68186517 + 0.00009800 9083.65137417 1210.68354946 51.53129119 -404.49112334 + 0.00010000 12.98279885 1.81441699 54.50576457 -282.00895288 + 0.00010200 -4694.86843894 -640.49023892 57.72139479 -110.55048155 + 0.00010400 -5249.24937729 -678.97189161 61.19985042 6.48124027 + 0.00010600 -4692.49854412 -573.75725497 64.96661400 73.78939041 + 0.00010800 -4842.03807640 -559.18594489 69.05073720 130.52160886 + 0.00011000 -6444.87970700 -692.52179240 73.47622377 195.90087001 + 0.00011200 -8761.15817605 -846.60156885 78.25904164 204.62808903 + 0.00011400 -11120.94207525 -917.27323881 83.41281019 105.85238279 + 0.00011600 -12487.28078237 -820.16700462 88.95902510 -71.43589299 + 0.00011800 -10595.44252406 -521.85476379 94.92547715 -240.74696485 + 0.00012000 -4328.87693103 -167.23087237 101.32849735 -329.27887177 + 0.00012200 5117.95222015 200.58807967 108.17078638 -331.96556974 + 0.00012400 15534.72579872 837.18288805 115.44812571 -273.77032314 + 0.00012600 15372.19160032 1166.54633476 123.16577183 -143.46280361 + 0.00012800 13199.26981327 1270.03948501 131.34284212 84.53268162 + 0.00013000 3889.72485947 421.57722256 139.99916669 366.97625240 + 0.00013200 -7609.01072867 -804.54137008 149.14109773 588.88717287 + 0.00013400 -16485.08328873 -1460.38288014 158.75662442 642.10586107 + 0.00013600 -19751.70263360 -1240.39245202 168.82466596 513.59518033 + 0.00013800 -14600.56951255 -559.95458905 179.32804534 301.55264267 + 0.00014000 -1954.06550150 -51.91954230 190.25821438 121.81791909 + 0.00014200 14688.98071130 523.41469312 201.60978325 25.53662766 + 0.00014400 28909.92796387 1927.17849885 213.36848057 -16.98627890 + 0.00014600 35444.70227900 3986.15244966 225.50204946 -59.14096672 + 0.00014800 31753.99162127 5089.68547076 237.96585087 -114.19100870 + 0.00015000 19066.65231695 3745.69661691 250.71602572 -172.95897513 + 0.00015200 1033.28363425 217.77219521 263.71514409 -202.99390816 + 0.00015400 -17523.12332680 -3487.48685363 276.92589903 -176.40811711 + 0.00015600 -30732.51592565 -5061.01775777 290.30153120 -85.63379745 + 0.00015800 -34893.34868250 -4116.54923544 303.77601724 59.26312152 + 0.00016000 -28812.43829975 -2092.85670093 317.27736371 231.49793081 + 0.00016200 -15930.13034772 -649.87005150 330.75219166 385.25655146 + 0.00016400 -201.63686952 -5.91086837 344.15499751 446.82662611 + 0.00016600 2135.67041083 65.54549512 357.42478232 347.74472500 + 0.00016800 2769.88369483 94.68015637 370.48840402 101.23584678 + 0.00017000 -5563.88881282 -179.12143808 383.27051138 -180.13175856 + 0.00017200 -16655.34783006 -272.82554896 395.69587138 -336.42665297 + 0.00017400 -25743.73899387 355.10351308 407.69294890 -255.80822779 + 0.00017600 -30711.43964525 1657.54592740 419.18969315 47.13427846 + 0.00017800 -29741.31364855 2884.74159668 430.11177201 421.77851602 + 0.00018000 -22528.97726730 3023.25287064 440.38933481 698.92434815 + 0.00018200 -11097.63474715 1754.81618622 449.96430373 797.37749633 + 0.00018400 328.48617200 -54.45952336 458.78677146 829.87099890 + 0.00018600 5849.23744020 -944.02511479 466.80248629 830.16654232 + 0.00018800 3474.25156172 -537.66663608 473.94489081 826.88926756 + 0.00019000 -3995.39107607 619.79872886 480.14686318 822.15525974 + 0.00019200 -11036.49454635 1830.14069868 485.36230147 831.11133165 + 0.00019400 -14001.78241787 2571.36327134 489.57390139 844.75260910 + 0.00019600 -12028.77674547 2431.86699351 492.77536795 826.99743344 + 0.00019800 -6620.62218248 1426.36740716 494.94559802 744.12013707 + 0.00020000 563.32886370 -123.79373740 496.04678396 599.88305724 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_ofrict2.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_ofrict2.txt Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,102 @@ +#x ACF frict_coeff gauss_filtered_frict_coeff av_window_frict_coeff + 0.00000000 820049.57805000 0.00000000 -88.35169831 -280.48699846 + 0.00000200 -4358.85795000 -133.28849699 -88.22972067 -88.19151898 + 0.00000400 -87816.09622050 -2381.86180009 -87.98019601 206.78227553 + 0.00000600 -131399.91872250 -2484.15267240 -87.60563052 420.75577467 + 0.00000800 -157962.78872850 -1272.80096612 -87.09850727 415.54857404 + 0.00001000 -138547.98040200 423.68679663 -86.46203739 220.05544591 + 0.00001200 -75149.24467650 831.83995256 -85.69381797 -19.39421284 + 0.00001400 17412.26728800 -230.42745195 -84.78750605 -184.44489621 + 0.00001600 115480.21126350 -952.91141368 -83.74024527 -254.60297579 + 0.00001800 188259.22998000 590.17608290 -82.55443295 -258.51588475 + 0.00002000 212009.57662500 3845.90958964 -81.23547895 -201.73752699 + 0.00002200 179721.96976500 5899.47589006 -79.79027794 -63.40787886 + 0.00002400 98900.16009300 4279.46998290 -78.22207384 124.26696357 + 0.00002600 -2220.79417500 -104.14401258 -76.52588032 336.60454536 + 0.00002800 -90011.60591700 -3909.86256270 -74.69455570 356.26981781 + 0.00003000 -137133.80762700 -4788.99317502 -72.72531275 198.84495402 + 0.00003200 -133164.08762535 -3301.01366740 -70.62158367 8.31697208 + 0.00003400 -84350.83000350 -1403.16159163 -68.38850295 -56.74371545 + 0.00003600 -5905.70960400 -78.25817913 -66.04439500 18.20705271 + 0.00003800 72166.91716950 1135.56715521 -63.58639145 89.69932255 + 0.00004000 120505.51429200 2766.59296252 -61.02332288 7.12530570 + 0.00004200 114021.51919200 3620.20835160 -58.35161851 -226.91946783 + 0.00004400 49243.56428550 1864.88983572 -55.57114618 -478.78785643 + 0.00004600 -55557.34158750 -2090.84722338 -52.69410000 -615.02734999 + 0.00004800 -147217.76437650 -4421.29824187 -49.73405580 -593.86618618 + 0.00005000 -198300.11521500 -3386.87284209 -46.70469036 -455.78564162 + 0.00005200 -192130.49563500 -469.37379811 -43.61960647 -279.26540999 + 0.00005400 -131119.41096000 1268.58791117 -40.49219486 -129.46297215 + 0.00005600 -30325.18489485 476.93398323 -37.33374606 -38.70994611 + 0.00005800 75656.48544150 -1061.30425422 -34.15245877 -32.81832328 + 0.00006000 146410.57420500 -834.98588086 -30.95661661 -112.68532949 + 0.00006200 154533.14829000 862.10811813 -27.76062817 -209.98160347 + 0.00006400 102247.69060950 1554.68011169 -24.58254792 -222.45823938 + 0.00006600 16240.64970450 319.07926355 -21.43985136 -77.52748335 + 0.00006800 -64351.78993350 -1148.25167151 -18.34979444 150.13100453 + 0.00007000 -112082.10661935 -1258.58854242 -15.32954316 388.33926875 + 0.00007200 -119224.23239265 -304.96491060 -12.39596315 478.34792093 + 0.00007400 -82381.22600550 411.89892915 -9.56277546 365.52349202 + 0.00007600 -4786.62582450 39.57427761 -6.84009826 164.45536646 + 0.00007800 93291.24579900 -461.77296900 -4.19942429 46.40314537 + 0.00008000 168592.66270500 820.66748274 -1.64143933 56.87529179 + 0.00008200 185712.16779000 3370.67563342 0.82140500 86.93796554 + 0.00008400 154195.97020350 4763.50495633 3.16818819 36.57899653 + 0.00008600 63075.57527115 2462.32253430 5.39209560 -59.25580198 + 0.00008800 8220.18044505 342.86704241 7.50012190 -96.67053627 + 0.00009000 -39127.79252655 -1586.69979832 9.50320322 -51.91590947 + 0.00009200 -75160.38225660 -2725.86237646 11.40764252 12.14852125 + 0.00009400 -91725.42820080 -2752.77459999 13.21381024 46.03722897 + 0.00009600 -79914.56742300 -1884.11043862 14.92948464 71.63113822 + 0.00009800 -45095.51373600 -851.86132583 16.57879150 154.22708975 + 0.00010000 -14017.68385650 -233.73240573 18.18817508 297.93359459 + 0.00010200 -9179.10851145 -145.07126219 19.77382639 384.18705955 + 0.00010400 -24921.70569600 -362.01585846 21.35197643 299.36149265 + 0.00010600 -34947.44092650 -429.21542217 22.95315194 100.04281508 + 0.00010800 -25252.93753500 -253.15913551 24.61507174 -94.88386254 + 0.00011000 2475.19015500 22.70007201 26.36654762 -172.11459139 + 0.00011200 35222.93315850 372.80924371 28.22156283 -130.44753477 + 0.00011400 64497.75738300 923.77711421 30.19195988 0.84946559 + 0.00011600 84900.31497900 1691.49330817 32.30327803 152.96661078 + 0.00011800 90975.51849000 2412.35712776 34.61832038 250.48097149 + 0.00012000 77534.31040200 2545.73678188 37.17714167 252.21706474 + 0.00012200 44742.66858000 1674.16429547 40.01343067 171.04221469 + 0.00012400 19512.94726350 777.13140500 43.14629587 56.04393382 + 0.00012600 -36653.09358330 -1436.21101812 46.59846989 -26.86155891 + 0.00012800 -32854.20668700 -1201.74753470 50.39600331 -40.13377747 + 0.00013000 -21818.38497900 -753.35866591 54.55818681 -7.50287868 + 0.00013200 -2973.95105400 -99.92236918 59.09130865 4.44251935 + 0.00013400 8564.88182730 289.56830372 63.99700240 -34.78999920 + 0.00013600 1934.85723000 66.17677573 69.28874967 -91.56918241 + 0.00013800 -24497.09191800 -817.13946083 74.99494252 -123.72241529 + 0.00014000 -61646.19561450 -1857.22940676 81.14055026 -132.89858804 + 0.00014200 -101891.36043540 -2445.03687958 87.72554000 -155.00683647 + 0.00014400 -132974.69188650 -2020.12571316 94.72632160 -212.73380806 + 0.00014600 -137663.38560000 -694.65979328 102.11343515 -303.03427244 + 0.00014800 -101707.77133650 399.45884033 109.86826252 -358.25769139 + 0.00015000 -29564.63915700 261.60803247 117.98170165 -377.49069744 + 0.00015200 51251.58621300 -411.84112715 126.43906153 -370.25980842 + 0.00015400 113180.97832800 -211.80655002 135.20886135 -370.18806573 + 0.00015600 138567.47375250 1048.42865059 144.24815584 -378.03850570 + 0.00015800 127159.72719000 2228.83367262 153.50681670 -350.28522288 + 0.00016000 82819.23176700 2103.57181338 162.93055802 -236.12084910 + 0.00016200 17812.41872085 519.62486357 172.45987082 -39.62300894 + 0.00016400 -36119.00680350 -1028.87788255 182.04578703 165.67530387 + 0.00016600 -110519.23112100 -2540.67939713 191.65166966 290.61170437 + 0.00016800 -106921.60660530 -1586.40765558 201.20240370 307.48058498 + 0.00017000 -76446.40564200 -608.74088555 210.62157891 271.62818673 + 0.00017200 -23322.73389000 -98.48751533 219.84153297 271.64664462 + 0.00017400 29707.45251000 132.55950431 228.81209036 334.73930429 + 0.00017600 74300.37906600 621.24243204 237.48578959 407.50929217 + 0.00017800 98741.98990950 1466.14801470 245.79589327 425.12644069 + 0.00018000 93913.58523300 2072.72739658 253.65590701 373.98926679 + 0.00018200 59684.08905900 1660.92937662 260.99073544 290.53799094 + 0.00018400 15614.32887000 478.60229674 267.76055483 264.55674776 + 0.00018600 -11599.31914500 -357.28218104 273.94457599 316.00064189 + 0.00018800 -10596.93693600 -317.58912475 279.50299875 443.03461174 + 0.00019000 8761.45959900 261.97719041 284.36758829 522.35499452 + 0.00019200 27201.03162900 850.01206613 288.46962428 552.79203880 + 0.00019400 36952.05022500 1243.59320758 291.77115090 557.71641457 + 0.00019600 38399.36686650 1400.77162096 294.26533754 551.08843935 + 0.00019800 30865.52299050 1206.09019462 295.95010795 520.02631775 + 0.00020000 10694.24856900 434.54629645 296.80389411 446.71891701 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_outp.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_outp.txt Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,102 @@ +#x force_integral frict_coeff wdiss corrected_force_integral + 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 + 0.00000200 -0.00392321 -309.80490686 -0.00000031 -0.00392290 + 0.00000400 -0.00600385 -1433.31122481 -0.00000205 -0.00600180 + 0.00000600 -0.00840497 -1302.42266399 -0.00000479 -0.00840018 + 0.00000800 -0.01081599 -696.00283593 -0.00000679 -0.01080921 + 0.00001000 -0.01276543 59.67135417 -0.00000742 -0.01275801 + 0.00001200 -0.01395301 200.26534006 -0.00000716 -0.01394585 + 0.00001400 -0.01428272 -211.00314767 -0.00000717 -0.01427554 + 0.00001600 -0.01390386 -299.86397145 -0.00000769 -0.01389617 + 0.00001800 -0.01311337 706.27695160 -0.00000728 -0.01310609 + 0.00002000 -0.01220919 2436.13821354 -0.00000414 -0.01220505 + 0.00002200 -0.01144642 3505.21213240 0.00000181 -0.01144822 + 0.00002400 -0.01105948 2809.74173435 0.00000812 -0.01106760 + 0.00002600 -0.01123787 993.11451262 0.00001192 -0.01124979 + 0.00002800 -0.01205259 -248.07012791 0.00001267 -0.01206525 + 0.00003000 -0.01340318 -195.60392001 0.00001222 -0.01341541 + 0.00003200 -0.01506203 141.25051892 0.00001217 -0.01507420 + 0.00003400 -0.01675665 -363.71528632 0.00001195 -0.01676860 + 0.00003600 -0.01826220 -1519.36792976 0.00001006 -0.01827227 + 0.00003800 -0.01947399 -1956.43454302 0.00000659 -0.01948058 + 0.00004000 -0.02036612 -865.57527338 0.00000377 -0.02036989 + 0.00004200 -0.02092478 626.10860652 0.00000353 -0.02092831 + 0.00004400 -0.02108110 565.05574074 0.00000472 -0.02108581 + 0.00004600 -0.02089789 -1403.86956218 0.00000388 -0.02090177 + 0.00004800 -0.02047142 -2830.86397506 -0.00000036 -0.02047107 + 0.00005000 -0.01987024 -2293.53919791 -0.00000548 -0.01986476 + 0.00005200 -0.01919238 -417.60570188 -0.00000819 -0.01918419 + 0.00005400 -0.01839938 834.36904591 -0.00000777 -0.01839160 + 0.00005600 -0.01743683 434.72799961 -0.00000650 -0.01743032 + 0.00005800 -0.01635230 -519.66723923 -0.00000659 -0.01634571 + 0.00006000 -0.01530848 -329.29698187 -0.00000744 -0.01530104 + 0.00006200 -0.01449254 858.01810556 -0.00000691 -0.01448563 + 0.00006400 -0.01404057 1441.72347492 -0.00000461 -0.01403596 + 0.00006600 -0.01404177 879.60901414 -0.00000229 -0.01403949 + 0.00006800 -0.01456203 334.25351485 -0.00000108 -0.01456096 + 0.00007000 -0.01559948 652.00817770 -0.00000009 -0.01559939 + 0.00007200 -0.01700283 1250.72698300 0.00000181 -0.01700464 + 0.00007400 -0.01850442 1051.23806218 0.00000412 -0.01850854 + 0.00007600 -0.01981895 -253.25084775 0.00000491 -0.01982387 + 0.00007800 -0.02071374 -1453.27137086 0.00000321 -0.02071694 + 0.00008000 -0.02104305 -880.42089139 0.00000087 -0.02104392 + 0.00008200 -0.02077116 1078.66842032 0.00000107 -0.02077223 + 0.00008400 -0.01997438 2679.13981866 0.00000483 -0.01997921 + 0.00008600 -0.01898980 1426.95656188 0.00000894 -0.01899874 + 0.00008800 -0.01814123 802.46409967 0.00001117 -0.01815240 + 0.00009000 -0.01764277 273.85567979 0.00001224 -0.01765501 + 0.00009200 -0.01769785 98.53732822 0.00001261 -0.01771046 + 0.00009400 -0.01829684 125.57269254 0.00001284 -0.01830967 + 0.00009600 -0.01933400 173.04677203 0.00001314 -0.01934713 + 0.00009800 -0.02061627 38.14199457 0.00001335 -0.02062962 + 0.00010000 -0.02175551 -405.56366316 0.00001298 -0.02176849 + 0.00010200 -0.02224141 -969.69926808 0.00001161 -0.02225302 + 0.00010400 -0.02170288 -1401.79910164 0.00000923 -0.02171211 + 0.00010600 -0.02011306 -1439.01003806 0.00000639 -0.02011945 + 0.00010800 -0.01776066 -1113.56692890 0.00000384 -0.01776450 + 0.00011000 -0.01509030 -682.72476001 0.00000204 -0.01509235 + 0.00011200 -0.01258934 -253.52957929 0.00000111 -0.01259045 + 0.00011400 -0.01068666 239.99423089 0.00000109 -0.01068776 + 0.00011600 -0.00962881 861.55125965 0.00000220 -0.00963100 + 0.00011800 -0.00937474 1593.45081266 0.00000465 -0.00937939 + 0.00012000 -0.00963019 2144.52398093 0.00000839 -0.00963857 + 0.00012200 -0.01000354 2216.16123424 0.00001275 -0.01001629 + 0.00012400 -0.01012789 2740.29495857 0.00001771 -0.01014560 + 0.00012600 -0.01002803 642.99155129 0.00002109 -0.01004912 + 0.00012800 -0.00980763 673.17973142 0.00002241 -0.00983004 + 0.00013000 -0.00958833 -276.05014319 0.00002280 -0.00961113 + 0.00013200 -0.00971501 -1331.11596207 0.00002120 -0.00973621 + 0.00013400 -0.01040330 -2042.60347448 0.00001782 -0.01042112 + 0.00013600 -0.01175235 -2339.54888278 0.00001344 -0.01176579 + 0.00013800 -0.01367748 -2263.02943488 0.00000884 -0.01368632 + 0.00014000 -0.01586775 -1978.10144773 0.00000460 -0.01587234 + 0.00014200 -0.01781525 -1395.88433631 0.00000122 -0.01781647 + 0.00014400 -0.01894530 -136.33782409 -0.00000031 -0.01894499 + 0.00014600 -0.01884929 1692.89613965 0.00000125 -0.01885054 + 0.00014800 -0.01750485 2974.81823090 0.00000591 -0.01751076 + 0.00015000 -0.01532342 2642.88373319 0.00001153 -0.01533495 + 0.00015200 -0.01295535 913.24280375 0.00001509 -0.01297044 + 0.00015400 -0.01100340 -927.03460327 0.00001507 -0.01101847 + 0.00015600 -0.00985920 -1689.34471443 0.00001246 -0.00987166 + 0.00015800 -0.00966688 -1332.07133205 0.00000944 -0.00967631 + 0.00016000 -0.01029419 -793.99131765 0.00000731 -0.01030150 + 0.00016200 -0.01132573 -941.55813718 0.00000557 -0.01133130 + 0.00016400 -0.01208447 -1038.93906805 0.00000359 -0.01208806 + 0.00016600 -0.01209668 -2451.04087672 0.00000010 -0.01209679 + 0.00016800 -0.01109858 -1619.20548854 -0.00000397 -0.01109462 + 0.00017000 -0.00915750 -1037.16157762 -0.00000662 -0.00915087 + 0.00017200 -0.00686287 -505.22323695 -0.00000817 -0.00685471 + 0.00017400 -0.00489099 169.23053990 -0.00000850 -0.00488249 + 0.00017600 -0.00385466 1119.27836540 -0.00000721 -0.00384744 + 0.00017800 -0.00413820 2035.34756564 -0.00000406 -0.00413414 + 0.00018000 -0.00577393 2281.35445610 0.00000026 -0.00577419 + 0.00018200 -0.00843098 1475.48961051 0.00000402 -0.00843500 + 0.00018400 -0.01147212 132.45287134 0.00000562 -0.01147774 + 0.00018600 -0.01409490 -616.86762340 0.00000514 -0.01410004 + 0.00018800 -0.01561112 -391.31151887 0.00000413 -0.01561525 + 0.00019000 -0.01567912 413.81613798 0.00000415 -0.01568328 + 0.00019200 -0.01440515 1254.39214820 0.00000582 -0.01441097 + 0.00019400 -0.01231202 1786.35191510 0.00000886 -0.01232088 + 0.00019600 -0.01014927 1765.17722675 0.00001241 -0.01016169 + 0.00019800 -0.00865789 1160.67458116 0.00001534 -0.00867323 + 0.00020000 -0.00832999 80.82726004 0.00001658 -0.00834657 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_outp1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_outp1.txt Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,102 @@ +#x force_integral frict_coeff wdiss corrected_force_integral + 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 + 0.00000200 -0.00028893 -1345.26646099 -0.00000135 -0.00028759 + 0.00000400 -0.00232537 -1454.28761686 -0.00000414 -0.00232122 + 0.00000600 -0.00534240 -264.14035591 -0.00000586 -0.00533654 + 0.00000800 -0.00892032 570.71230720 -0.00000556 -0.00891476 + 0.00001000 -0.01233174 482.67878751 -0.00000450 -0.01232723 + 0.00001200 -0.01484120 -306.50010800 -0.00000433 -0.01483688 + 0.00001400 -0.01592271 -678.79838871 -0.00000531 -0.01591740 + 0.00001600 -0.01540730 129.24578885 -0.00000586 -0.01540143 + 0.00001800 -0.01350838 1433.88087902 -0.00000430 -0.01350408 + 0.00002000 -0.01074543 1769.49504443 -0.00000110 -0.01074434 + 0.00002200 -0.00782545 525.42519661 0.00000120 -0.00782665 + 0.00002400 -0.00551646 -1354.79390711 0.00000037 -0.00551683 + 0.00002600 -0.00443259 -2004.09700782 -0.00000299 -0.00442960 + 0.00002800 -0.00480792 -481.28911857 -0.00000547 -0.00480244 + 0.00003000 -0.00643049 2092.57404588 -0.00000386 -0.00642663 + 0.00003200 -0.00878093 3287.09518295 0.00000152 -0.00878245 + 0.00003400 -0.01123762 1810.83705441 0.00000661 -0.01124423 + 0.00003600 -0.01326588 -1271.21578790 0.00000715 -0.01327304 + 0.00003800 -0.01459317 -3241.15658461 0.00000264 -0.01459581 + 0.00004000 -0.01524085 -2483.33843193 -0.00000308 -0.01523777 + 0.00004200 -0.01546401 -1.82301038 -0.00000557 -0.01545845 + 0.00004400 -0.01558410 1555.33658588 -0.00000401 -0.01558009 + 0.00004600 -0.01603192 960.69346965 -0.00000150 -0.01603042 + 0.00004800 -0.01694582 -784.52267942 -0.00000132 -0.01694450 + 0.00005000 -0.01809375 -1447.42762641 -0.00000355 -0.01809019 + 0.00005200 -0.01916352 -598.12158006 -0.00000560 -0.01915792 + 0.00005400 -0.01967108 508.73630130 -0.00000569 -0.01966539 + 0.00005600 -0.01920233 634.81104146 -0.00000455 -0.01919778 + 0.00005800 -0.01764643 -33.15102509 -0.00000394 -0.01764248 + 0.00006000 -0.01531761 -485.10192035 -0.00000446 -0.01531314 + 0.00006200 -0.01285302 -375.61072392 -0.00000532 -0.01284769 + 0.00006400 -0.01095150 -90.31689620 -0.00000579 -0.01094571 + 0.00006600 -0.01011817 242.77363760 -0.00000564 -0.01011253 + 0.00006800 -0.01052398 1053.59838602 -0.00000434 -0.01051964 + 0.00007000 -0.01200084 2387.77239428 -0.00000090 -0.01199995 + 0.00007200 -0.01412820 3253.13344370 0.00000474 -0.01413295 + 0.00007400 -0.01634020 2389.26594789 0.00001038 -0.01635058 + 0.00007600 -0.01802042 -263.98184134 0.00001251 -0.01803293 + 0.00007800 -0.01864988 -2971.49817282 0.00000927 -0.01865915 + 0.00008000 -0.01803190 -3596.40029897 0.00000271 -0.01803460 + 0.00008200 -0.01641443 -2140.04635679 -0.00000303 -0.01641140 + 0.00008400 -0.01424356 -508.98011481 -0.00000568 -0.01423789 + 0.00008600 -0.01233851 31.90196581 -0.00000616 -0.01233236 + 0.00008800 -0.01108021 -37.54573797 -0.00000616 -0.01107405 + 0.00009000 -0.01041675 247.03284223 -0.00000595 -0.01041080 + 0.00009200 -0.01057237 1170.97275460 -0.00000453 -0.01056784 + 0.00009400 -0.01152362 2085.38263558 -0.00000128 -0.01152234 + 0.00009600 -0.01304294 2145.82898503 0.00000295 -0.01304590 + 0.00009800 -0.01473765 1210.68354946 0.00000631 -0.01474396 + 0.00010000 -0.01613022 1.81441699 0.00000752 -0.01613774 + 0.00010200 -0.01683142 -640.49023892 0.00000688 -0.01683830 + 0.00010400 -0.01665238 -678.97189161 0.00000556 -0.01665794 + 0.00010600 -0.01558814 -573.75725497 0.00000431 -0.01559245 + 0.00010800 -0.01375773 -559.18594489 0.00000318 -0.01376091 + 0.00011000 -0.01139273 -692.52179240 0.00000193 -0.01139466 + 0.00011200 -0.00886230 -846.60156885 0.00000039 -0.00886269 + 0.00011400 -0.00659599 -917.27323881 -0.00000138 -0.00659461 + 0.00011600 -0.00490657 -820.16700462 -0.00000311 -0.00490346 + 0.00011800 -0.00384287 -521.85476379 -0.00000446 -0.00383842 + 0.00012000 -0.00322160 -167.23087237 -0.00000514 -0.00321646 + 0.00012200 -0.00279559 200.58807967 -0.00000511 -0.00279048 + 0.00012400 -0.00221567 837.18288805 -0.00000407 -0.00221159 + 0.00012600 -0.00179223 1166.54633476 -0.00000207 -0.00179016 + 0.00012800 -0.00161866 1270.03948501 0.00000037 -0.00161903 + 0.00013000 -0.00151181 421.57722256 0.00000206 -0.00151387 + 0.00013200 -0.00185055 -804.54137008 0.00000168 -0.00185223 + 0.00013400 -0.00284118 -1460.38288014 -0.00000059 -0.00284059 + 0.00013600 -0.00463146 -1240.39245202 -0.00000329 -0.00462817 + 0.00013800 -0.00718450 -559.95458905 -0.00000509 -0.00717941 + 0.00014000 -0.01016260 -51.91954230 -0.00000570 -0.01015689 + 0.00014200 -0.01298551 523.41469312 -0.00000523 -0.01298028 + 0.00014400 -0.01501735 1927.17849885 -0.00000278 -0.01501457 + 0.00014600 -0.01576224 3986.15244966 0.00000313 -0.01576537 + 0.00014800 -0.01501718 5089.68547076 0.00001221 -0.01502939 + 0.00015000 -0.01296828 3745.69661691 0.00002104 -0.01298932 + 0.00015200 -0.01017811 217.77219521 0.00002501 -0.01020312 + 0.00015400 -0.00739771 -3487.48685363 0.00002174 -0.00741945 + 0.00015600 -0.00530829 -5061.01775777 0.00001319 -0.00532148 + 0.00015800 -0.00432459 -4116.54923544 0.00000401 -0.00432860 + 0.00016000 -0.00449294 -2092.85670093 -0.00000220 -0.00449074 + 0.00016200 -0.00550412 -649.87005150 -0.00000494 -0.00549918 + 0.00016400 -0.00660743 -5.91086837 -0.00000560 -0.00660183 + 0.00016600 -0.00753781 65.54549512 -0.00000554 -0.00753227 + 0.00016800 -0.00787258 94.68015637 -0.00000538 -0.00786720 + 0.00017000 -0.00715284 -179.12143808 -0.00000546 -0.00714738 + 0.00017200 -0.00581344 -272.82554896 -0.00000591 -0.00580753 + 0.00017400 -0.00439586 355.10351308 -0.00000583 -0.00439003 + 0.00017600 -0.00348190 1657.54592740 -0.00000382 -0.00347808 + 0.00017800 -0.00349891 2884.74159668 0.00000072 -0.00349964 + 0.00018000 -0.00462675 3023.25287064 0.00000663 -0.00463339 + 0.00018200 -0.00676588 1754.81618622 0.00001141 -0.00677729 + 0.00018400 -0.00947283 -54.45952336 0.00001311 -0.00948594 + 0.00018600 -0.01198297 -944.02511479 0.00001211 -0.01199508 + 0.00018800 -0.01351134 -537.66663608 0.00001063 -0.01352197 + 0.00019000 -0.01359850 619.79872886 0.00001071 -0.01360921 + 0.00019200 -0.01230380 1830.14069868 0.00001316 -0.01231696 + 0.00019400 -0.01014874 2571.36327134 0.00001756 -0.01016630 + 0.00019600 -0.00786581 2431.86699351 0.00002257 -0.00788837 + 0.00019800 -0.00618981 1426.36740716 0.00002643 -0.00621624 + 0.00020000 -0.00567906 -123.79373740 0.00002773 -0.00570679 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/NEQGamma_outp2.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/NEQGamma_outp2.txt Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,102 @@ +#x force_integral frict_coeff wdiss corrected_force_integral + 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 + 0.00000200 -0.00330406 -133.28849699 -0.00000013 -0.00330393 + 0.00000400 -0.00636445 -2381.86180009 -0.00000265 -0.00636180 + 0.00000600 -0.00965354 -2484.15267240 -0.00000751 -0.00964603 + 0.00000800 -0.01252890 -1272.80096612 -0.00001127 -0.01251763 + 0.00001000 -0.01442164 423.68679663 -0.00001212 -0.01440951 + 0.00001200 -0.01503189 831.83995256 -0.00001086 -0.01502103 + 0.00001400 -0.01441631 -230.42745195 -0.00001026 -0.01440604 + 0.00001600 -0.01302548 -952.91141368 -0.00001145 -0.01301403 + 0.00001800 -0.01151471 590.17608290 -0.00001181 -0.01150290 + 0.00002000 -0.01045748 3845.90958964 -0.00000737 -0.01045010 + 0.00002200 -0.01022110 5899.47589006 0.00000237 -0.01022347 + 0.00002400 -0.01098400 4279.46998290 0.00001255 -0.01099655 + 0.00002600 -0.01275202 -104.14401258 0.00001673 -0.01276875 + 0.00002800 -0.01533740 -3909.86256270 0.00001271 -0.01535011 + 0.00003000 -0.01834052 -4788.99317502 0.00000401 -0.01834454 + 0.00003200 -0.02122997 -3301.01366740 -0.00000408 -0.02122589 + 0.00003400 -0.02350381 -1403.16159163 -0.00000878 -0.02349503 + 0.00003600 -0.02483672 -78.25817913 -0.00001026 -0.02482646 + 0.00003800 -0.02517675 1135.56715521 -0.00000921 -0.02516755 + 0.00004000 -0.02469464 2766.59296252 -0.00000530 -0.02468934 + 0.00004200 -0.02365645 3620.20835160 0.00000108 -0.02365753 + 0.00004400 -0.02230644 1864.88983572 0.00000657 -0.02231301 + 0.00004600 -0.02092224 -2090.84722338 0.00000634 -0.02092858 + 0.00004800 -0.01977588 -4421.29824187 -0.00000017 -0.01977571 + 0.00005000 -0.01901760 -3386.87284209 -0.00000798 -0.01900962 + 0.00005200 -0.01861994 -469.37379811 -0.00001183 -0.01860811 + 0.00005400 -0.01835163 1268.58791117 -0.00001103 -0.01834059 + 0.00005600 -0.01796265 476.93398323 -0.00000929 -0.01795336 + 0.00005800 -0.01738252 -1061.30425422 -0.00000987 -0.01737265 + 0.00006000 -0.01676130 -834.98588086 -0.00001177 -0.01674953 + 0.00006200 -0.01633255 862.10811813 -0.00001174 -0.01632081 + 0.00006400 -0.01629663 1554.68011169 -0.00000933 -0.01628730 + 0.00006600 -0.01682112 319.07926355 -0.00000745 -0.01681367 + 0.00006800 -0.01803550 -1148.25167151 -0.00000828 -0.01802722 + 0.00007000 -0.01990021 -1258.58854242 -0.00001069 -0.01988953 + 0.00007200 -0.02209227 -304.96491060 -0.00001225 -0.02208002 + 0.00007400 -0.02412258 411.89892915 -0.00001214 -0.02411044 + 0.00007600 -0.02554865 39.57427761 -0.00001169 -0.02553696 + 0.00007800 -0.02611824 -461.77296900 -0.00001212 -0.02610613 + 0.00008000 -0.02580973 820.66748274 -0.00001176 -0.02579797 + 0.00008200 -0.02480289 3370.67563342 -0.00000756 -0.02479533 + 0.00008400 -0.02346738 4763.50495633 0.00000057 -0.02346795 + 0.00008600 -0.02222273 2462.32253430 0.00000780 -0.02223053 + 0.00008800 -0.02150421 342.86704241 0.00001060 -0.02151481 + 0.00009000 -0.02163462 -1586.69979832 0.00000936 -0.02164398 + 0.00009200 -0.02259376 -2725.86237646 0.00000504 -0.02259881 + 0.00009400 -0.02410201 -2752.77459999 -0.00000043 -0.02410158 + 0.00009600 -0.02583460 -1884.11043862 -0.00000507 -0.02582953 + 0.00009800 -0.02750085 -851.86132583 -0.00000781 -0.02749304 + 0.00010000 -0.02870505 -233.73240573 -0.00000889 -0.02869616 + 0.00010200 -0.02898468 -145.07126219 -0.00000927 -0.02897541 + 0.00010400 -0.02805692 -362.01585846 -0.00000978 -0.02804714 + 0.00010600 -0.02599487 -429.21542217 -0.00001057 -0.02598430 + 0.00010800 -0.02318061 -253.15913551 -0.00001125 -0.02316936 + 0.00011000 -0.02011980 22.70007201 -0.00001148 -0.02010832 + 0.00011200 -0.01731324 372.80924371 -0.00001109 -0.01730216 + 0.00011400 -0.01517354 923.77711421 -0.00000979 -0.01516375 + 0.00011600 -0.01393412 1691.49330817 -0.00000718 -0.01392695 + 0.00011800 -0.01359860 2412.35712776 -0.00000307 -0.01359553 + 0.00012000 -0.01397801 2545.73678188 0.00000189 -0.01397989 + 0.00012200 -0.01476766 1674.16429547 0.00000611 -0.01477377 + 0.00012400 -0.01567924 777.13140500 0.00000856 -0.01568779 + 0.00012600 -0.01639225 -1436.21101812 0.00000790 -0.01640015 + 0.00012800 -0.01675034 -1201.74753470 0.00000526 -0.01675560 + 0.00013000 -0.01683481 -753.35866591 0.00000331 -0.01683811 + 0.00013200 -0.01677711 -99.92236918 0.00000245 -0.01677956 + 0.00013400 -0.01680665 289.56830372 0.00000264 -0.01680929 + 0.00013600 -0.01717171 66.17677573 0.00000300 -0.01717471 + 0.00013800 -0.01805922 -817.13946083 0.00000225 -0.01806146 + 0.00014000 -0.01950682 -1857.22940676 -0.00000043 -0.01950640 + 0.00014200 -0.02128489 -2445.03687958 -0.00000473 -0.02128016 + 0.00014400 -0.02288252 -2020.12571316 -0.00000920 -0.02287333 + 0.00014600 -0.02372178 -694.65979328 -0.00001191 -0.02370987 + 0.00014800 -0.02349327 399.45884033 -0.00001221 -0.02348107 + 0.00015000 -0.02231618 261.60803247 -0.00001154 -0.02230464 + 0.00015200 -0.02058259 -411.84112715 -0.00001169 -0.02057090 + 0.00015400 -0.01869749 -211.80655002 -0.00001232 -0.01868517 + 0.00015600 -0.01700933 1048.42865059 -0.00001148 -0.01699785 + 0.00015800 -0.01583020 2228.83367262 -0.00000820 -0.01582200 + 0.00016000 -0.01534495 2103.57181338 -0.00000387 -0.01534108 + 0.00016200 -0.01544411 519.62486357 -0.00000125 -0.01544286 + 0.00016400 -0.01569072 -1028.87788255 -0.00000176 -0.01568896 + 0.00016600 -0.01528460 -2540.67939713 -0.00000533 -0.01527927 + 0.00016800 -0.01372341 -1586.40765558 -0.00000945 -0.01371395 + 0.00017000 -0.01111279 -608.74088555 -0.00001165 -0.01110114 + 0.00017200 -0.00787247 -98.48751533 -0.00001236 -0.00786012 + 0.00017400 -0.00473024 132.55950431 -0.00001232 -0.00471792 + 0.00017600 -0.00244342 621.24243204 -0.00001157 -0.00243186 + 0.00017800 -0.00158672 1466.14801470 -0.00000948 -0.00157723 + 0.00018000 -0.00237524 2072.72739658 -0.00000594 -0.00236930 + 0.00018200 -0.00458287 1660.92937662 -0.00000221 -0.00458066 + 0.00018400 -0.00756423 478.60229674 -0.00000007 -0.00756416 + 0.00018600 -0.01037347 -357.28218104 0.00000005 -0.01037352 + 0.00018800 -0.01208003 -317.58912475 -0.00000062 -0.01207941 + 0.00019000 -0.01212748 261.97719041 -0.00000068 -0.01212680 + 0.00019200 -0.01054721 850.01206613 0.00000043 -0.01054764 + 0.00019400 -0.00795757 1243.59320758 0.00000253 -0.00796010 + 0.00019600 -0.00530655 1400.77162096 0.00000517 -0.00531172 + 0.00019800 -0.00351460 1206.09019462 0.00000778 -0.00352238 + 0.00020000 -0.00315052 434.54629645 0.00000942 -0.00315994 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/pull1.xvg --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/pull1.xvg Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,118 @@ +# This file was created Sat Jan 11 15:16:55 2020 +# Created by: +# :-) GROMACS - gmx mdrun, 2019.1 (-: +# +# Executable: /usr/local/bin/gmx +# Data prefix: /usr/local +# Working dir: /data/share/staging/6673623/working +# Command line: +# gmx mdrun -nt 16 -deffnm outp +# gmx mdrun is part of G R O M A C S: +# +# God Rules Over Mankind, Animals, Cosmos and Such +# +@ title "Pull force" +@ xaxis label "Time (ps)" +@ yaxis label "Force (kJ/mol/nm)" +@TYPE xy +0.0000 6024.09 +0.0020 -1440.45 +0.0040 -2511.27 +0.0060 -2897.63 +0.0080 -2775.86 +0.0100 -1984.12 +0.0120 -692.582 +0.0140 749.852 +0.0160 1926.04 +0.0180 2521.88 +0.0200 2405.94 +0.0220 1618.47 +0.0240 312.892 +0.0260 -1146.03 +0.0280 -2331.23 +0.0300 -2868.38 +0.0320 -2634.84 +0.0340 -1742.36 +0.0360 -463.328 +0.0380 764.038 +0.0400 1581.21 +0.0420 1724.85 +0.0440 1203.93 +0.0460 119.215 +0.0480 -933.681 +0.0500 -1649.19 +0.0520 -1728.6 +0.0540 -1128.9 +0.0560 -43.2839 +0.0580 1061.76 +0.0600 1706.84 +0.0620 1632.03 +0.0640 886.953 +0.0660 -265.667 +0.0680 -1413.95 +0.0700 -2156.87 +0.0720 -2271.91 +0.0740 -1707.92 +0.0760 -561.06 +0.0780 847.304 +0.0800 1993.62 +0.0820 2439.33 +0.0840 2183.08 +0.0860 1162.58 +0.0880 245.372 +0.0900 -674.664 +0.0920 -1389.64 +0.0940 -1732.39 +0.0960 -1659.95 +0.0980 -1215.14 +0.1000 -560.687 +0.1020 56.7457 +0.1040 541.266 +0.1060 941.849 +0.1080 1290.27 +0.1100 1550.28 +0.1120 1620.82 +0.1140 1483.18 +0.1160 1200.91 +0.1180 835.327 +0.1200 414.748 +0.1220 -21.986 +0.1240 -268.239 +0.1260 -610.519 +0.1280 -419.701 +0.1300 -193.452 +0.1320 11.4089 +0.1340 13.1112 +0.1360 -276.636 +0.1380 -829.048 +0.1400 -1451.56 +0.1420 -1907.91 +0.1440 -1960.87 +0.1460 -1495.45 +0.1480 -590.751 +0.1500 498.442 +0.1520 1444.86 +0.1540 2030.3 +0.1560 2092.26 +0.1580 1656.44 +0.1600 859.3 +0.1620 14.6479 +0.1640 -438.284 +0.1660 -573.584 +0.1680 32.1338 +0.1700 805.322 +0.1720 1470.23 +0.1740 1733.74 +0.1760 1558.83 +0.1780 971.192 +0.1800 103.252 +0.1820 -825.594 +0.1840 -1427.63 +0.1860 -1342.79 +0.1880 -578.404 +0.1900 513.202 +0.1920 1414.83 +0.1940 1795.16 +0.1960 1584.51 +0.1980 877.231 +0.2000 -111.269 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/pull2.xvg --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/pull2.xvg Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,118 @@ +# This file was created Sat Jan 11 15:14:02 2020 +# Created by: +# :-) GROMACS - gmx mdrun, 2019.1 (-: +# +# Executable: /usr/local/bin/gmx +# Data prefix: /usr/local +# Working dir: /data/share/staging/6673627/working +# Command line: +# gmx mdrun -nt 16 -deffnm outp +# gmx mdrun is part of G R O M A C S: +# +# GROup of MAchos and Cynical Suckers +# +@ title "Pull force" +@ xaxis label "Time (ps)" +@ yaxis label "Force (kJ/mol/nm)" +@TYPE xy +0.0000 -5098.38 +0.0020 -63.1292 +0.0040 -58.0201 +0.0060 -567.151 +0.0080 -915.192 +0.0100 -1147.66 +0.0120 -1194.57 +0.0140 -1025.71 +0.0160 -619.358 +0.0180 -30.7206 +0.0200 628.788 +0.0220 1186.76 +0.0240 1499.87 +0.0260 1501.01 +0.0280 1225.59 +0.0300 728.864 +0.0320 73.4782 +0.0340 -609.651 +0.0360 -1241.19 +0.0380 -1714.09 +0.0400 -1926.53 +0.0420 -1825.85 +0.0440 -1343.11 +0.0460 -875.662 +0.0480 -137.674 +0.0500 424.689 +0.0520 813.554 +0.0540 1028.83 +0.0560 1080.86 +0.0580 1012.46 +0.0600 876.583 +0.0620 713.725 +0.0640 570.333 +0.0660 475.042 +0.0680 392.949 +0.0700 224.142 +0.0720 -50.0798 +0.0740 -394.081 +0.0760 -697.375 +0.0780 -847.788 +0.0800 -757.177 +0.0820 -440.838 +0.0840 160.155 +0.0860 304.292 +0.0880 804.365 +0.0900 951.843 +0.0920 801.215 +0.0940 418.322 +0.0960 -64.6304 +0.0980 -449.703 +0.1000 -559.593 +0.1020 -338.869 +0.1040 98.9362 +0.1060 546.434 +0.1080 882.254 +0.1100 1007.2 +0.1120 882.558 +0.1140 546.07 +0.1160 148.665 +0.1180 -57.5016 +0.1200 49.9737 +0.1220 409.28 +0.1240 1040.8 +0.1260 684.824 +0.1280 692.54 +0.1300 134.317 +0.1320 -629.767 +0.1340 -1376.01 +0.1360 -1941.02 +0.1380 -2059.37 +0.1400 -1616.22 +0.1420 -670.138 +0.1440 475.235 +0.1460 1491.31 +0.1480 2085.01 +0.1500 2105.1 +0.1520 1531.93 +0.1540 553.708 +0.1560 -497.426 +0.1580 -1283.86 +0.1600 -1568.59 +0.1620 -1327.71 +0.1640 -455.275 +0.1660 -393.621 +0.1680 265.539 +0.1700 336.479 +0.1720 66.7614 +0.1740 -435.565 +0.1760 -1029.08 +0.1780 -1534.97 +0.1800 -1795.16 +0.1820 -1760.74 +0.1840 -1399.95 +0.1860 -849.902 +0.1880 -285.645 +0.1900 176.529 +0.1920 484.836 +0.1940 615.295 +0.1960 570.901 +0.1980 319.342 +0.2000 -63.7999 diff -r d9f8cc3258f9 -r 6f2029791c94 test-data/pull3.xvg --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/pull3.xvg Fri Sep 11 21:55:57 2020 +0000 @@ -0,0 +1,118 @@ +# This file was created Sat Jan 11 11:17:19 2020 +# Created by: +# :-) GROMACS - gmx mdrun, 2019.1 (-: +# +# Executable: /usr/local/bin/gmx +# Data prefix: /usr/local +# Working dir: /data/share/staging/6671882/working +# Command line: +# gmx mdrun -nt 16 -deffnm outp +# gmx mdrun is part of G R O M A C S: +# +# Gyas ROwers Mature At Cryogenic Speed +# +@ title "Pull force" +@ xaxis label "Time (ps)" +@ yaxis label "Force (kJ/mol/nm)" +@TYPE xy +0.0000 -9835.61 +0.0020 -1356.15 +0.0040 -812.913 +0.0060 -356.365 +0.0080 279.129 +0.0100 695.388 +0.0120 760.799 +0.0140 413.1 +0.0160 -307.339 +0.0180 -1119.04 +0.0200 -1694.31 +0.0220 -1857.34 +0.0240 -1599.83 +0.0260 -1103.08 +0.0280 -590.412 +0.0300 -216.222 +0.0320 -59.4561 +0.0340 -111.021 +0.0360 -349.112 +0.0380 -631.665 +0.0400 -749.358 +0.0420 -480.318 +0.0440 251.563 +0.0460 1193.69 +0.0480 1913.5 +0.0500 2185.92 +0.0520 1987.19 +0.0540 1406.94 +0.0560 543.203 +0.0580 -401.431 +0.0600 -1124.73 +0.0620 -1356.63 +0.0640 -1090.51 +0.0660 -579.76 +0.0680 -169.391 +0.0700 10.7899 +0.0720 33.8781 +0.0740 -114.673 +0.0760 -468.487 +0.0780 -956.942 +0.0800 -1266.95 +0.0820 -1152.33 +0.0840 -799.059 +0.0860 -57.2971 +0.0880 86.3943 +0.0900 82.0647 +0.0920 63.9564 +0.0940 41.5732 +0.0960 -114.408 +0.0980 -342.996 +0.1000 -289.586 +0.1020 234.269 +0.1040 1023.25 +0.1060 1617.73 +0.1080 1778.66 +0.1100 1502.41 +0.1120 939.611 +0.1140 235.798 +0.1160 -441.056 +0.1180 -924.133 +0.1200 -1084.76 +0.1220 -887.306 +0.1240 -645.618 +0.1260 98.3492 +0.1280 215.697 +0.1300 228.514 +0.1320 68.9249 +0.1340 -152.533 +0.1360 -314.056 +0.1380 -355.276 +0.1400 -259.327 +0.1420 62.6616 +0.1440 610.851 +0.1460 1166.95 +0.1480 1376.27 +0.1500 1070.22 +0.1520 453.658 +0.1540 -158.612 +0.1560 -587.625 +0.1580 -802.82 +0.1600 -742.418 +0.1620 -329.843 +0.1640 260.255 +0.1660 1563.85 +0.1680 2099.99 +0.1700 2283.79 +0.1720 1921.29 +0.1740 1159.2 +0.1760 121.866 +0.1780 -938.471 +0.1800 -1713.03 +0.1820 -1979.88 +0.1840 -1729.61 +0.1860 -1118.46 +0.1880 -373.46 +0.1900 343.756 +0.1920 888.764 +0.1940 1080.51 +0.1960 841.869 +0.1980 280.294 +0.2000 -318.095