62 Commits

Author SHA1 Message Date
Daniel Bauer
35ca1b1317 version up 2021-07-21 16:47:32 +02:00
Daniel Bauer
d9a219a36a better verbose output during dataset generation and flag to ignore empty histograms 2021-07-21 16:45:58 +02:00
Daniel Bauer
2c565dee28 remove some debug statements 2021-07-19 08:32:57 +02:00
Daniel Bauer
c05becea22 version up 2021-07-17 14:44:09 +02:00
Daniel Bauer
c27dbb7ab5 test g=1 for correlation 2021-07-17 14:42:21 +02:00
Daniel Bauer
c9e18ebda9 Squashed commit of the following:
commit eaebf0dcbb259decbc0d8f5bbffa62244303f6c7
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Sat Jul 17 14:12:30 2021 +0200

    error for empty timeseries

commit 1fe5383c5079ca6c0ad102d4e2cb32d5b1227e80
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Sat Jul 17 14:09:13 2021 +0200

    refractored convdt slices calculation

commit 0b1fb7fb6a72edc50b013b59623551b2ccab6913
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Sat Jul 17 13:27:59 2021 +0200

    fix histogram building without convdt

commit f9882ca4cece57451cd2971d0993a2d3621b0cdf
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Sat Jul 17 12:53:06 2021 +0200

    fix tests not compiling

commit e6550e20bde3824f8199432b08777be41fc79fa8
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Sat Jul 17 12:43:36 2021 +0200

    refractoring

commit 18c77a9b6694d491ef23cebb3918ccda8fcc44a5
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Fri Jul 16 08:16:34 2021 +0200

    run and output for multiple datasets

commit 069f318f72207c416387007435eeff9867494e7a
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Fri Jul 16 07:58:01 2021 +0200

    cleanup io.rs

commit 10efa428a0c6bf490c9f2d7a4c1df185de402a11
Author: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
Date:   Thu Jul 15 19:27:22 2021 +0200

    parse multiple datasets with convdt
2021-07-17 14:14:25 +02:00
Daniel Bauer
0b5d0ebecd test time boundaries 2021-07-17 13:31:36 +02:00
Daniel Bauer
091b1f1382 TODO 2021-07-17 13:04:42 +02:00
Daniel Bauer
9849c00321 cargo ignore examples 2021-07-17 11:30:42 +02:00
daniel
034586d08e inlined autocorrelation calculation for better performance 2020-10-26 11:24:48 +01:00
Daniel Bauer
cc4448f1d8 Update README.md 2020-10-26 10:53:57 +01:00
Daniel Bauer
5ca708abf1 version up 2020-10-26 10:52:22 +01:00
Daniel Bauer
31ffb5958d Merge branch 'corr' of github.com:danijoo/WHAM into corr 2020-10-26 10:46:43 +01:00
Daniel Bauer
ffab55e0e1 version up 2020-10-26 10:46:23 +01:00
Daniel Bauer
69d0c7ca00 better output for uncorr 2020-10-26 10:45:22 +01:00
Daniel Bauer
8a88d7742a remove gsl 2020-10-26 10:20:41 +01:00
Daniel Bauer
18debf1a28 Update README.md 2020-10-25 22:48:38 +01:00
Daniel Bauer
a856e5d34f Update cli.yml 2020-10-25 22:41:05 +01:00
Daniel Bauer
2cc34919b0 Autocorrelation (#2)
* calculate autocorrelation stats_ineff and tau
* uncorr flag
* read timeseries into vector
* uncorrelate data
* README
* Update README.md
* Update README.md

Co-authored-by: Daniel Bauer <bauer@cbs.tu-darmstadt.de>
2020-10-25 22:40:13 +01:00
Daniel Bauer
cd57289211 Update README.md 2020-10-25 22:39:15 +01:00
Daniel Bauer
9d8423bce9 Update README.md 2020-10-25 22:36:27 +01:00
Daniel Bauer
2876c1dc3a README 2020-10-25 20:49:52 +01:00
Daniel Bauer
b9d08ae24b uncorrelate data 2020-10-25 20:36:02 +01:00
Daniel Bauer
07fc8344fe read timeseries into vector 2020-10-25 19:42:33 +01:00
Daniel Bauer
d7eea7aa03 uncorr flag 2020-10-25 19:07:33 +01:00
Daniel Bauer
3d93693eac calculate autocorrelation stats_ineff and tau 2020-10-25 18:38:25 +01:00
Daniel Bauer
38c30baa5b category 2020-10-24 20:53:07 +02:00
Daniel Bauer
b130596fe8 Merge branch 'master' of github.com:danijoo/WHAM 2020-10-24 20:44:30 +02:00
Daniel Bauer
05451a6cb2 cli 2020-10-24 20:44:27 +02:00
Daniel Bauer
bf8bf39319 edit todo 2020-10-24 20:44:18 +02:00
Daniel Bauer
6c63eda769 Update README.md 2020-10-24 20:43:54 +02:00
Daniel Bauer
1bbafeba70 lock update 2020-10-24 20:38:37 +02:00
Daniel Bauer
64f468e636 version up 2020-10-24 20:35:24 +02:00
Daniel Bauer
2f93b77879 update rand lib 2020-10-24 20:35:06 +02:00
Daniel Bauer
b148998b11 update GSL version 2020-10-24 20:21:32 +02:00
Daniel Bauer
87f3cceb6a keywords,categories 2020-10-24 20:06:28 +02:00
Daniel Bauer
7b39d603bd remove unneccessary parameter in bootstrapping function 2020-09-15 15:32:05 +02:00
Daniel Bauer
465f46b280 fix unit tests 2020-09-15 15:20:40 +02:00
Daniel Bauer
d221fc7195 lock 2020-09-15 14:49:03 +02:00
Daniel Bauer
dc218df051 calculate SE of A directly 2020-09-15 14:48:36 +02:00
Daniel Bauer
1e27cfd457 report standard error (SE) instead of standard deviations 2020-09-15 14:10:48 +02:00
Daniel Bauer
7116142c2a .gitignore 2020-03-11 09:40:31 +01:00
Daniel Bauer
fbdb451f6b print less 2020-03-11 09:40:01 +01:00
Daniel Bauer
d4119b2fae clean code (2) 2020-03-11 09:31:48 +01:00
Daniel Bauer
b9f6b24b20 clean code 2020-03-10 09:44:32 +01:00
Daniel Bauer
5390b089ac remove some code smell 2020-03-09 10:05:36 +01:00
Daniel Bauer
133234706e remove 1d test because of cargo size limit 2020-03-08 12:52:55 +01:00
Daniel Bauer
1a95f222eb cargo lock 2020-03-08 12:42:28 +01:00
Daniel Bauer
77637133a7 Merge branch 'master' of github.com:danijoo/WHAM 2020-03-08 12:41:03 +01:00
Daniel Bauer
a099a71b3b bump version 2020-03-08 12:40:09 +01:00
Daniel Bauer
a4f3b3091c test bootstrapping 2020-03-08 12:39:41 +01:00
Daniel Bauer
08e7c93327 seedable rng 2020-03-08 12:29:52 +01:00
Daniel Bauer
936cbe1c58 add pi to min/max 2020-03-08 11:42:47 +01:00
Daniel Bauer
554e250838 Update README.md 2020-03-08 00:30:16 +01:00
Daniel Bauer
cd7bea063f non-cyclic 1d test 2020-03-07 21:10:31 +01:00
Daniel Bauer
d3a70520df version up 2018-12-06 11:40:37 +01:00
Daniel Bauer
aca4265a02 fix an IOB when a timeseries value is exactly on the upper histogram boundary 2018-12-06 11:40:26 +01:00
Daniel Bauer
d758c843d8 possibility to truncate timeseries with --start/stop cli arguments 2018-12-04 16:56:43 +01:00
Daniel Bauer
79db640a7b v0.9.3 2018-11-30 13:09:05 +01:00
Daniel Bauer
99f2c9069d cleanup + version 2018-11-30 13:06:58 +01:00
Daniel Bauer
7109d27543 v0.9.2 2018-11-30 12:58:53 +01:00
Daniel Bauer
c03be1723e better error messages when reading timeseries 2018-11-30 12:58:24 +01:00
683 changed files with 11510 additions and 10621 deletions

3
.gitignore vendored
View File

@@ -1,4 +1,5 @@
target
reference
.idea
*bench.*
*bench.*
*.sublime*

402
Cargo.lock generated
View File

@@ -1,372 +1,400 @@
# This file is automatically @generated by Cargo.
# It is not intended for manual editing.
[[package]]
name = "GSL"
version = "1.0.0"
name = "addr2line"
version = "0.13.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "1b6a2d3371669ab3ca9797670853d61402b03d0b4b9ebf33d677dfa720203072"
dependencies = [
"c_vec 1.0.12 (registry+https://github.com/rust-lang/crates.io-index)",
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"pkg-config 0.3.14 (registry+https://github.com/rust-lang/crates.io-index)",
"gimli",
]
[[package]]
name = "adler"
version = "0.2.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ee2a4ec343196209d6594e19543ae87a39f96d5534d7174822a3ad825dd6ed7e"
[[package]]
name = "ansi_term"
version = "0.11.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ee49baf6cb617b853aa8d93bf420db2383fab46d314482ca2803b40d5fde979b"
dependencies = [
"winapi 0.3.6 (registry+https://github.com/rust-lang/crates.io-index)",
"winapi",
]
[[package]]
name = "arrayvec"
version = "0.4.7"
name = "assert_approx_eq"
version = "1.1.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
dependencies = [
"nodrop 0.1.13 (registry+https://github.com/rust-lang/crates.io-index)",
]
checksum = "3c07dab4369547dbe5114677b33fbbf724971019f3818172d59a97a61c774ffd"
[[package]]
name = "atty"
version = "0.2.11"
version = "0.2.14"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "d9b39be18770d11421cdb1b9947a45dd3f37e93092cbf377614828a319d5fee8"
dependencies = [
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"termion 1.5.1 (registry+https://github.com/rust-lang/crates.io-index)",
"winapi 0.3.6 (registry+https://github.com/rust-lang/crates.io-index)",
"hermit-abi",
"libc",
"winapi",
]
[[package]]
name = "autocfg"
version = "1.0.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "cdb031dd78e28731d87d56cc8ffef4a8f36ca26c38fe2de700543e627f8a464a"
[[package]]
name = "backtrace"
version = "0.3.9"
version = "0.3.53"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "707b586e0e2f247cbde68cdd2c3ce69ea7b7be43e1c5b426e37c9319c4b9838e"
dependencies = [
"backtrace-sys 0.1.24 (registry+https://github.com/rust-lang/crates.io-index)",
"cfg-if 0.1.6 (registry+https://github.com/rust-lang/crates.io-index)",
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"rustc-demangle 0.1.9 (registry+https://github.com/rust-lang/crates.io-index)",
"winapi 0.3.6 (registry+https://github.com/rust-lang/crates.io-index)",
]
[[package]]
name = "backtrace-sys"
version = "0.1.24"
source = "registry+https://github.com/rust-lang/crates.io-index"
dependencies = [
"cc 1.0.25 (registry+https://github.com/rust-lang/crates.io-index)",
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"addr2line",
"cfg-if 1.0.0",
"libc",
"miniz_oxide",
"object",
"rustc-demangle",
]
[[package]]
name = "bitflags"
version = "1.0.4"
source = "registry+https://github.com/rust-lang/crates.io-index"
[[package]]
name = "c_vec"
version = "1.0.12"
source = "registry+https://github.com/rust-lang/crates.io-index"
[[package]]
name = "cc"
version = "1.0.25"
version = "1.2.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "cf1de2fe8c75bc145a2f577add951f8134889b4795d47466a54a5c846d691693"
[[package]]
name = "cfg-if"
version = "0.1.6"
version = "0.1.10"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "4785bdd1c96b2a846b2bd7cc02e86b6b3dbf14e7e53446c4f54c92a361040822"
[[package]]
name = "cfg-if"
version = "1.0.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "baf1de4339761588bc0619e3cbc0120ee582ebb74b53b4efbf79117bd2da40fd"
[[package]]
name = "clap"
version = "2.32.0"
version = "2.33.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "37e58ac78573c40708d45522f0d80fa2f01cc4f9b4e2bf749807255454312002"
dependencies = [
"ansi_term 0.11.0 (registry+https://github.com/rust-lang/crates.io-index)",
"atty 0.2.11 (registry+https://github.com/rust-lang/crates.io-index)",
"bitflags 1.0.4 (registry+https://github.com/rust-lang/crates.io-index)",
"strsim 0.7.0 (registry+https://github.com/rust-lang/crates.io-index)",
"textwrap 0.10.0 (registry+https://github.com/rust-lang/crates.io-index)",
"unicode-width 0.1.5 (registry+https://github.com/rust-lang/crates.io-index)",
"vec_map 0.8.1 (registry+https://github.com/rust-lang/crates.io-index)",
"yaml-rust 0.3.5 (registry+https://github.com/rust-lang/crates.io-index)",
"ansi_term",
"atty",
"bitflags",
"strsim",
"textwrap",
"unicode-width",
"vec_map",
"yaml-rust",
]
[[package]]
name = "cloudabi"
version = "0.0.3"
name = "const_fn"
version = "0.4.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ce90df4c658c62f12d78f7508cf92f9173e5184a539c10bfe54a3107b3ffd0f2"
[[package]]
name = "crossbeam-channel"
version = "0.5.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "dca26ee1f8d361640700bde38b2c37d8c22b3ce2d360e1fc1c74ea4b0aa7d775"
dependencies = [
"bitflags 1.0.4 (registry+https://github.com/rust-lang/crates.io-index)",
"cfg-if 1.0.0",
"crossbeam-utils",
]
[[package]]
name = "crossbeam-deque"
version = "0.2.0"
version = "0.8.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "94af6efb46fef72616855b036a624cf27ba656ffc9be1b9a3c931cfc7749a9a9"
dependencies = [
"crossbeam-epoch 0.3.1 (registry+https://github.com/rust-lang/crates.io-index)",
"crossbeam-utils 0.2.2 (registry+https://github.com/rust-lang/crates.io-index)",
"cfg-if 1.0.0",
"crossbeam-epoch",
"crossbeam-utils",
]
[[package]]
name = "crossbeam-epoch"
version = "0.3.1"
version = "0.9.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ec0f606a85340376eef0d6d8fec399e6d4a544d648386c6645eb6d0653b27d9f"
dependencies = [
"arrayvec 0.4.7 (registry+https://github.com/rust-lang/crates.io-index)",
"cfg-if 0.1.6 (registry+https://github.com/rust-lang/crates.io-index)",
"crossbeam-utils 0.2.2 (registry+https://github.com/rust-lang/crates.io-index)",
"lazy_static 1.2.0 (registry+https://github.com/rust-lang/crates.io-index)",
"memoffset 0.2.1 (registry+https://github.com/rust-lang/crates.io-index)",
"nodrop 0.1.13 (registry+https://github.com/rust-lang/crates.io-index)",
"scopeguard 0.3.3 (registry+https://github.com/rust-lang/crates.io-index)",
"cfg-if 1.0.0",
"const_fn",
"crossbeam-utils",
"lazy_static",
"memoffset",
"scopeguard",
]
[[package]]
name = "crossbeam-utils"
version = "0.2.2"
version = "0.8.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ec91540d98355f690a86367e566ecad2e9e579f230230eb7c21398372be73ea5"
dependencies = [
"cfg-if 0.1.6 (registry+https://github.com/rust-lang/crates.io-index)",
"autocfg",
"cfg-if 1.0.0",
"const_fn",
"lazy_static",
]
[[package]]
name = "either"
version = "1.5.0"
version = "1.6.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "e78d4f1cc4ae33bbfc157ed5d5a5ef3bc29227303d595861deb238fcec4e9457"
[[package]]
name = "error-chain"
version = "0.12.0"
version = "0.12.4"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "2d2f06b9cac1506ece98fe3231e3cc9c4410ec3d5b1f24ae1c8946f0742cdefc"
dependencies = [
"backtrace 0.3.9 (registry+https://github.com/rust-lang/crates.io-index)",
"backtrace",
"version_check",
]
[[package]]
name = "fuchsia-zircon"
version = "0.3.3"
name = "getrandom"
version = "0.1.15"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "fc587bc0ec293155d5bfa6b9891ec18a1e330c234f896ea47fbada4cadbe47e6"
dependencies = [
"bitflags 1.0.4 (registry+https://github.com/rust-lang/crates.io-index)",
"fuchsia-zircon-sys 0.3.3 (registry+https://github.com/rust-lang/crates.io-index)",
"cfg-if 0.1.10",
"libc",
"wasi",
]
[[package]]
name = "fuchsia-zircon-sys"
version = "0.3.3"
name = "gimli"
version = "0.22.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "aaf91faf136cb47367fa430cd46e37a788775e7fa104f8b4bcb3861dc389b724"
[[package]]
name = "hermit-abi"
version = "0.1.17"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "5aca5565f760fb5b220e499d72710ed156fdb74e631659e99377d9ebfbd13ae8"
dependencies = [
"libc",
]
[[package]]
name = "lazy_static"
version = "1.2.0"
version = "1.4.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "e2abad23fbc42b3700f2f279844dc832adb2b2eb069b2df918f455c4e18cc646"
[[package]]
name = "libc"
version = "0.2.43"
version = "0.2.80"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "4d58d1b70b004888f764dfbf6a26a3b0342a1632d33968e4a179d8011c760614"
[[package]]
name = "memoffset"
version = "0.2.1"
version = "0.5.6"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "043175f069eda7b85febe4a74abbaeff828d9f8b448515d3151a14a3542811aa"
dependencies = [
"autocfg",
]
[[package]]
name = "nodrop"
version = "0.1.13"
name = "miniz_oxide"
version = "0.4.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "0f2d26ec3309788e423cfbf68ad1800f061638098d76a83681af979dc4eda19d"
dependencies = [
"adler",
"autocfg",
]
[[package]]
name = "num_cpus"
version = "1.8.0"
version = "1.13.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "05499f3756671c15885fee9034446956fff3f243d6077b91e5767df161f766b3"
dependencies = [
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"hermit-abi",
"libc",
]
[[package]]
name = "pkg-config"
version = "0.3.14"
name = "object"
version = "0.21.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "37fd5004feb2ce328a52b0b3d01dbf4ffff72583493900ed15f22d4111c51693"
[[package]]
name = "ppv-lite86"
version = "0.2.9"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "c36fa947111f5c62a733b652544dd0016a43ce89619538a8ef92724a6f501a20"
[[package]]
name = "rand"
version = "0.5.5"
version = "0.7.3"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "6a6b1679d49b24bbfe0c803429aa1874472f50d9b363131f0e89fc356b544d03"
dependencies = [
"cloudabi 0.0.3 (registry+https://github.com/rust-lang/crates.io-index)",
"fuchsia-zircon 0.3.3 (registry+https://github.com/rust-lang/crates.io-index)",
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"rand_core 0.2.2 (registry+https://github.com/rust-lang/crates.io-index)",
"winapi 0.3.6 (registry+https://github.com/rust-lang/crates.io-index)",
"getrandom",
"libc",
"rand_chacha",
"rand_core",
"rand_hc",
]
[[package]]
name = "rand_core"
name = "rand_chacha"
version = "0.2.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "f4c8ed856279c9737206bf725bf36935d8666ead7aa69b52be55af369d193402"
dependencies = [
"rand_core 0.3.0 (registry+https://github.com/rust-lang/crates.io-index)",
"ppv-lite86",
"rand_core",
]
[[package]]
name = "rand_core"
version = "0.3.0"
version = "0.5.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "90bde5296fc891b0cef12a6d03ddccc162ce7b2aff54160af9338f8d40df6d19"
dependencies = [
"getrandom",
]
[[package]]
name = "rand_hc"
version = "0.2.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ca3129af7b92a17112d59ad498c6f81eaf463253766b90396d39ea7a39d6613c"
dependencies = [
"rand_core",
]
[[package]]
name = "rayon"
version = "1.0.3"
version = "1.5.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "8b0d8e0819fadc20c74ea8373106ead0600e3a67ef1fe8da56e39b9ae7275674"
dependencies = [
"crossbeam-deque 0.2.0 (registry+https://github.com/rust-lang/crates.io-index)",
"either 1.5.0 (registry+https://github.com/rust-lang/crates.io-index)",
"rayon-core 1.4.1 (registry+https://github.com/rust-lang/crates.io-index)",
"autocfg",
"crossbeam-deque",
"either",
"rayon-core",
]
[[package]]
name = "rayon-core"
version = "1.4.1"
version = "1.9.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "9ab346ac5921dc62ffa9f89b7a773907511cdfa5490c572ae9be1be33e8afa4a"
dependencies = [
"crossbeam-deque 0.2.0 (registry+https://github.com/rust-lang/crates.io-index)",
"lazy_static 1.2.0 (registry+https://github.com/rust-lang/crates.io-index)",
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"num_cpus 1.8.0 (registry+https://github.com/rust-lang/crates.io-index)",
]
[[package]]
name = "redox_syscall"
version = "0.1.40"
source = "registry+https://github.com/rust-lang/crates.io-index"
[[package]]
name = "redox_termios"
version = "0.1.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
dependencies = [
"redox_syscall 0.1.40 (registry+https://github.com/rust-lang/crates.io-index)",
"crossbeam-channel",
"crossbeam-deque",
"crossbeam-utils",
"lazy_static",
"num_cpus",
]
[[package]]
name = "rustc-demangle"
version = "0.1.9"
version = "0.1.18"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "6e3bad0ee36814ca07d7968269dd4b7ec89ec2da10c4bb613928d3077083c232"
[[package]]
name = "scopeguard"
version = "0.3.3"
version = "1.1.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "d29ab0c6d3fc0ee92fe66e2d99f700eab17a8d57d1c1d3b748380fb20baa78cd"
[[package]]
name = "strsim"
version = "0.7.0"
version = "0.8.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
[[package]]
name = "termion"
version = "1.5.1"
source = "registry+https://github.com/rust-lang/crates.io-index"
dependencies = [
"libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)",
"redox_syscall 0.1.40 (registry+https://github.com/rust-lang/crates.io-index)",
"redox_termios 0.1.1 (registry+https://github.com/rust-lang/crates.io-index)",
]
checksum = "8ea5119cdb4c55b55d432abb513a0429384878c15dde60cc77b1c99de1a95a6a"
[[package]]
name = "textwrap"
version = "0.10.0"
version = "0.11.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "d326610f408c7a4eb6f51c37c330e496b08506c9457c9d34287ecc38809fb060"
dependencies = [
"unicode-width 0.1.5 (registry+https://github.com/rust-lang/crates.io-index)",
"unicode-width",
]
[[package]]
name = "unicode-width"
version = "0.1.5"
version = "0.1.8"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "9337591893a19b88d8d87f2cec1e73fad5cdfd10e5a6f349f498ad6ea2ffb1e3"
[[package]]
name = "vec_map"
version = "0.8.1"
version = "0.8.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "f1bddf1187be692e79c5ffeab891132dfb0f236ed36a43c7ed39f1165ee20191"
[[package]]
name = "version_check"
version = "0.9.2"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "b5a972e5669d67ba988ce3dc826706fb0a8b01471c088cb0b6110b805cc36aed"
[[package]]
name = "wasi"
version = "0.9.0+wasi-snapshot-preview1"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "cccddf32554fecc6acb585f82a32a72e28b48f8c4c1883ddfeeeaa96f7d8e519"
[[package]]
name = "wham"
version = "0.9.1"
version = "1.1.2"
dependencies = [
"GSL 1.0.0 (registry+https://github.com/rust-lang/crates.io-index)",
"clap 2.32.0 (registry+https://github.com/rust-lang/crates.io-index)",
"error-chain 0.12.0 (registry+https://github.com/rust-lang/crates.io-index)",
"rand 0.5.5 (registry+https://github.com/rust-lang/crates.io-index)",
"rayon 1.0.3 (registry+https://github.com/rust-lang/crates.io-index)",
"assert_approx_eq",
"clap",
"error-chain",
"rand",
"rayon",
]
[[package]]
name = "winapi"
version = "0.3.6"
version = "0.3.9"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "5c839a674fcd7a98952e593242ea400abe93992746761e38641405d28b00f419"
dependencies = [
"winapi-i686-pc-windows-gnu 0.4.0 (registry+https://github.com/rust-lang/crates.io-index)",
"winapi-x86_64-pc-windows-gnu 0.4.0 (registry+https://github.com/rust-lang/crates.io-index)",
"winapi-i686-pc-windows-gnu",
"winapi-x86_64-pc-windows-gnu",
]
[[package]]
name = "winapi-i686-pc-windows-gnu"
version = "0.4.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "ac3b87c63620426dd9b991e5ce0329eff545bccbbb34f3be09ff6fb6ab51b7b6"
[[package]]
name = "winapi-x86_64-pc-windows-gnu"
version = "0.4.0"
source = "registry+https://github.com/rust-lang/crates.io-index"
checksum = "712e227841d057c1ee1cd2fb22fa7e5a5461ae8e48fa2ca79ec42cfc1931183f"
[[package]]
name = "yaml-rust"
version = "0.3.5"
source = "registry+https://github.com/rust-lang/crates.io-index"
[metadata]
"checksum GSL 1.0.0 (registry+https://github.com/rust-lang/crates.io-index)" = "2e70167131757876cda6ab95869b57ad6b55c627ff28421f5935f55f9b97d85c"
"checksum ansi_term 0.11.0 (registry+https://github.com/rust-lang/crates.io-index)" = "ee49baf6cb617b853aa8d93bf420db2383fab46d314482ca2803b40d5fde979b"
"checksum arrayvec 0.4.7 (registry+https://github.com/rust-lang/crates.io-index)" = "a1e964f9e24d588183fcb43503abda40d288c8657dfc27311516ce2f05675aef"
"checksum atty 0.2.11 (registry+https://github.com/rust-lang/crates.io-index)" = "9a7d5b8723950951411ee34d271d99dddcc2035a16ab25310ea2c8cfd4369652"
"checksum backtrace 0.3.9 (registry+https://github.com/rust-lang/crates.io-index)" = "89a47830402e9981c5c41223151efcced65a0510c13097c769cede7efb34782a"
"checksum backtrace-sys 0.1.24 (registry+https://github.com/rust-lang/crates.io-index)" = "c66d56ac8dabd07f6aacdaf633f4b8262f5b3601a810a0dcddffd5c22c69daa0"
"checksum bitflags 1.0.4 (registry+https://github.com/rust-lang/crates.io-index)" = "228047a76f468627ca71776ecdebd732a3423081fcf5125585bcd7c49886ce12"
"checksum c_vec 1.0.12 (registry+https://github.com/rust-lang/crates.io-index)" = "aa9e1d9f7d49e289f36f19effbf3d5a5e30163ecf9c7a3c9be94d5374dec5b9a"
"checksum cc 1.0.25 (registry+https://github.com/rust-lang/crates.io-index)" = "f159dfd43363c4d08055a07703eb7a3406b0dac4d0584d96965a3262db3c9d16"
"checksum cfg-if 0.1.6 (registry+https://github.com/rust-lang/crates.io-index)" = "082bb9b28e00d3c9d39cc03e64ce4cea0f1bb9b3fde493f0cbc008472d22bdf4"
"checksum clap 2.32.0 (registry+https://github.com/rust-lang/crates.io-index)" = "b957d88f4b6a63b9d70d5f454ac8011819c6efa7727858f458ab71c756ce2d3e"
"checksum cloudabi 0.0.3 (registry+https://github.com/rust-lang/crates.io-index)" = "ddfc5b9aa5d4507acaf872de71051dfd0e309860e88966e1051e462a077aac4f"
"checksum crossbeam-deque 0.2.0 (registry+https://github.com/rust-lang/crates.io-index)" = "f739f8c5363aca78cfb059edf753d8f0d36908c348f3d8d1503f03d8b75d9cf3"
"checksum crossbeam-epoch 0.3.1 (registry+https://github.com/rust-lang/crates.io-index)" = "927121f5407de9956180ff5e936fe3cf4324279280001cd56b669d28ee7e9150"
"checksum crossbeam-utils 0.2.2 (registry+https://github.com/rust-lang/crates.io-index)" = "2760899e32a1d58d5abb31129f8fae5de75220bc2176e77ff7c627ae45c918d9"
"checksum either 1.5.0 (registry+https://github.com/rust-lang/crates.io-index)" = "3be565ca5c557d7f59e7cfcf1844f9e3033650c929c6566f511e8005f205c1d0"
"checksum error-chain 0.12.0 (registry+https://github.com/rust-lang/crates.io-index)" = "07e791d3be96241c77c43846b665ef1384606da2cd2a48730abe606a12906e02"
"checksum fuchsia-zircon 0.3.3 (registry+https://github.com/rust-lang/crates.io-index)" = "2e9763c69ebaae630ba35f74888db465e49e259ba1bc0eda7d06f4a067615d82"
"checksum fuchsia-zircon-sys 0.3.3 (registry+https://github.com/rust-lang/crates.io-index)" = "3dcaa9ae7725d12cdb85b3ad99a434db70b468c09ded17e012d86b5c1010f7a7"
"checksum lazy_static 1.2.0 (registry+https://github.com/rust-lang/crates.io-index)" = "a374c89b9db55895453a74c1e38861d9deec0b01b405a82516e9d5de4820dea1"
"checksum libc 0.2.43 (registry+https://github.com/rust-lang/crates.io-index)" = "76e3a3ef172f1a0b9a9ff0dd1491ae5e6c948b94479a3021819ba7d860c8645d"
"checksum memoffset 0.2.1 (registry+https://github.com/rust-lang/crates.io-index)" = "0f9dc261e2b62d7a622bf416ea3c5245cdd5d9a7fcc428c0d06804dfce1775b3"
"checksum nodrop 0.1.13 (registry+https://github.com/rust-lang/crates.io-index)" = "2f9667ddcc6cc8a43afc9b7917599d7216aa09c463919ea32c59ed6cac8bc945"
"checksum num_cpus 1.8.0 (registry+https://github.com/rust-lang/crates.io-index)" = "c51a3322e4bca9d212ad9a158a02abc6934d005490c054a2778df73a70aa0a30"
"checksum pkg-config 0.3.14 (registry+https://github.com/rust-lang/crates.io-index)" = "676e8eb2b1b4c9043511a9b7bea0915320d7e502b0a079fb03f9635a5252b18c"
"checksum rand 0.5.5 (registry+https://github.com/rust-lang/crates.io-index)" = "e464cd887e869cddcae8792a4ee31d23c7edd516700695608f5b98c67ee0131c"
"checksum rand_core 0.2.2 (registry+https://github.com/rust-lang/crates.io-index)" = "1961a422c4d189dfb50ffa9320bf1f2a9bd54ecb92792fb9477f99a1045f3372"
"checksum rand_core 0.3.0 (registry+https://github.com/rust-lang/crates.io-index)" = "0905b6b7079ec73b314d4c748701f6931eb79fd97c668caa3f1899b22b32c6db"
"checksum rayon 1.0.3 (registry+https://github.com/rust-lang/crates.io-index)" = "373814f27745b2686b350dd261bfd24576a6fb0e2c5919b3a2b6005f820b0473"
"checksum rayon-core 1.4.1 (registry+https://github.com/rust-lang/crates.io-index)" = "b055d1e92aba6877574d8fe604a63c8b5df60f60e5982bf7ccbb1338ea527356"
"checksum redox_syscall 0.1.40 (registry+https://github.com/rust-lang/crates.io-index)" = "c214e91d3ecf43e9a4e41e578973adeb14b474f2bee858742d127af75a0112b1"
"checksum redox_termios 0.1.1 (registry+https://github.com/rust-lang/crates.io-index)" = "7e891cfe48e9100a70a3b6eb652fef28920c117d366339687bd5576160db0f76"
"checksum rustc-demangle 0.1.9 (registry+https://github.com/rust-lang/crates.io-index)" = "bcfe5b13211b4d78e5c2cadfebd7769197d95c639c35a50057eb4c05de811395"
"checksum scopeguard 0.3.3 (registry+https://github.com/rust-lang/crates.io-index)" = "94258f53601af11e6a49f722422f6e3425c52b06245a5cf9bc09908b174f5e27"
"checksum strsim 0.7.0 (registry+https://github.com/rust-lang/crates.io-index)" = "bb4f380125926a99e52bc279241539c018323fab05ad6368b56f93d9369ff550"
"checksum termion 1.5.1 (registry+https://github.com/rust-lang/crates.io-index)" = "689a3bdfaab439fd92bc87df5c4c78417d3cbe537487274e9b0b2dce76e92096"
"checksum textwrap 0.10.0 (registry+https://github.com/rust-lang/crates.io-index)" = "307686869c93e71f94da64286f9a9524c0f308a9e1c87a583de8e9c9039ad3f6"
"checksum unicode-width 0.1.5 (registry+https://github.com/rust-lang/crates.io-index)" = "882386231c45df4700b275c7ff55b6f3698780a650026380e72dabe76fa46526"
"checksum vec_map 0.8.1 (registry+https://github.com/rust-lang/crates.io-index)" = "05c78687fb1a80548ae3250346c3db86a80a7cdd77bda190189f2d0a0987c81a"
"checksum winapi 0.3.6 (registry+https://github.com/rust-lang/crates.io-index)" = "92c1eb33641e276cfa214a0522acad57be5c56b10cb348b3c5117db75f3ac4b0"
"checksum winapi-i686-pc-windows-gnu 0.4.0 (registry+https://github.com/rust-lang/crates.io-index)" = "ac3b87c63620426dd9b991e5ce0329eff545bccbbb34f3be09ff6fb6ab51b7b6"
"checksum winapi-x86_64-pc-windows-gnu 0.4.0 (registry+https://github.com/rust-lang/crates.io-index)" = "712e227841d057c1ee1cd2fb22fa7e5a5461ae8e48fa2ca79ec42cfc1931183f"
"checksum yaml-rust 0.3.5 (registry+https://github.com/rust-lang/crates.io-index)" = "e66366e18dc58b46801afbf2ca7661a9f59cc8c5962c29892b6039b4f86fa992"
checksum = "e66366e18dc58b46801afbf2ca7661a9f59cc8c5962c29892b6039b4f86fa992"

View File

@@ -1,20 +1,26 @@
[package]
name = "wham"
version = "0.9.1"
authors = ["D. Bauer <bauer@bio.tu-darmstadt.de>"]
version = "1.1.2"
authors = ["Daniel Bauer <bauer@cbs.tu-darmstadt.de>"]
description = "An implementation of the weighted histogram analysis method"
license = "GPL-3.0"
repository = "https://github.com/danijoo/WHAM"
readme = "README.md"
categories = ["science", "command-line-utilities", "algorithms"]
keywords = ["math", "statistics", "histogram", "bioinformatics", "molecular-dynamics"]
exclude = [
"example/*"
]
[dependencies]
clap = {version="2.32.0", features=['yaml']}
error-chain = "0.12.0"
rand = "0.5.5"
GSL = "1.0.0"
rand = "0.7.*"
rayon = "1.0.3"
[profile.release]
opt-level = 2
[dev-dependencies]
assert_approx_eq = "1.1.0"
[profile.release]
opt-level = 3
[features]
default = ["GSL/v2"]

104
README.md
View File

@@ -1,4 +1,4 @@
[![Build Status](https://travis-ci.com/danijoo/WHAM.svg?branch=master)](https://travis-ci.com/danijoo/WHAM) [![crates.io](https://img.shields.io/badge/crates.io-orange.svg?longCache=true)](https://www.crates.io/crates/wham) [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.1488598.svg)](https://doi.org/10.5281/zenodo.1488598)
[![Build Status](https://travis-ci.com/danijoo/WHAM.svg?branch=master)](https://travis-ci.com/danijoo/WHAM) [![crates.io](https://img.shields.io/badge/crates.io-orange.svg?longCache=true)](https://www.crates.io/crates/wham) [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.1488597.svg)](https://doi.org/10.5281/zenodo.1488597)
@@ -12,19 +12,14 @@ from umbrella sampling simulations. For more details on the method, I suggest *R
Features
---
- Fast, especially for small systems
- Multithreaded
- Multidimensional
- Error analysis
- Multithreaded (automatically runs on all available cores)
- Multidimensional (any number of collective variables are possible)
- Autocorrelation to remove correlated samples
- Error analysis via bootstrapping
- Unit tested
Installation
---
WHAM requires the GSL library to be installed:
```bash
# on debian/ubuntu:
sudo apt-get install libgsl0-dev
```
Installation from source via cargo:
```bash
# cargo installation
@@ -38,7 +33,71 @@ Usage
wham has a convenient command line interface. You can see all options with
```wham -h```:
To run the two dimensional example (simulation of dialanine phi and psi angle):
```
wham 1.1.0
D. Bauer <bauer@cbs.tu-darmstadt.de>
wham is a fast implementation of the weighted histogram analysis method (WHAM) written in Rust. It currently supports
potential of mean force (PMF) calculations in multiple dimensions at constant temperature.
Metadata file format:
/path/to/timeseries_file1 x_1 x_2 x_N fc_1 fc_2 fc_N
/path/to/timeseries_file2 x_1 x_2 x_N fc_1 fc_2 fc_N
/path/to/timeseries_file3 x_1 x_2 x_N fc_1 fc_2 fc_N
The first column is a path to a timeseries file _relative_ to the metadata file (see below). This is followed by the
position of the umbrella potential x in N dimensions and the force constant fc in each dimension. Lines starting with a
# are treated as comments and will not be parsed.
Timeseries file format:
time x_1 x_2 x_N
time x_1 x_2 x_N
time x_1 x_2 x_N
The first column will be ignored and is followed by N reaction coordinates x.
Shipped under the GPLv3 license.
USAGE:
wham [FLAGS] [OPTIONS] --bins <BINS> --max <HIST_MAX> --file <METADATA> --min <HIST_MIN> --temperature <temperature>
FLAGS:
-c, --cyclic For periodic reaction coordinates. If this is set, the first and last coordinate bin in each
dimension are treated as neighbors for the bias calculation.
-h, --help Prints help information
-g, --uncorr Estimates statistical inefficiency of each timeseries via autocorrelation and removes correlated
samples (default is off).
-V, --version Prints version information
-v, --verbose Enables verbose output.
OPTIONS:
-b, --bins <BINS> Number of histogram bins (comma separated).
--bt <bootstrap> Number of bayesian bootstrapping runs for error analysis by assigning random
weights (defaults to 0).
--seed <bootstrap_seed> Random seed for bootstrapping runs.
--convdt <convdt> Performs WHAM for slices with the given delta in time and returns an output file
for each slice. THis is useful to check the result for convergence. Example: with
--convdt 100 and a timeseries ranging from 0-300, free energy surfaces for slices
0-100, 0-200 and 0-300 will be given returned.
--end <end> Skip rows in timeseries with an index larger than this value (defaults to 1e+20)
-i, --iterations <ITERATIONS> Stop WHAM after this many iterations without convergence (defaults to 100,000).
--max <HIST_MAX> Histogram maxima (comma separated). Also accepts "pi".
-f, --file <METADATA> Path to the metadata file.
--min <HIST_MIN> Histogram minima (comma separated for multiple dimensions). Also accepts "pi".
-o, --output <output> Free energy output file (defaults to wham.out).
--start <start> Skip rows in timeseries with an index smaller than this value (defaults to 0)
-T, --temperature <temperature> WHAM temperature in Kelvin.
-t, --tolerance <TOLERANCE> Abortion criteria for WHAM calculation. WHAM stops if abs(F_new - F_old) <
tolerance (defaults to 0.000001).
```
Examples
---
The example folder contains input and output files for two simple test systems:
- 1d_cyclic: Phi torsion angle of dialanine in vaccum
- 2d_cyclic: Phi and psi torsion angles of the same system
The command below will run the two dimensional example (simulation of dialanine phi and psi angle) and calculate the free energy based on the two collective variables
in the range of -3.14 to 3.14, with 100 bins in each dimension and periodic collective variables:
```bash
wham --max 3.14,3.14 --min -3.14,-3.14 -T 300 --bins 100,100 --cyclic -f example/2d/metadata.dat
> Supplied WHAM options: Metadata=example/2d/metadata.dat, hist_min=[-3.14, -3.14], hist_max=[3.14, 3.14], bins=[100, 100] verbose=false, tolerance=0.000001, iterations=100000, temperature=300, cyclic=true
@@ -55,7 +114,6 @@ wham --max 3.14,3.14 --min -3.14,-3.14 -T 300 --bins 100,100 --cyclic -f example
```
After convergence, final bias offsets (F) and the free energy will be dumped to stdout and the output file is written.
The output file contains the free energy and probability for each bin. Probabilities are normalized to sum to P=1.0 and
the smallest free energy is set to 0 (with other free energies based on that).
```
@@ -75,26 +133,30 @@ the smallest free energy is set to 0 (with other free energies based on that).
Error analysis
---
WHAM can perform error analysis using the bayesian bootstrapping method. Every simulation window is assumed to be an
individual set of data point. By calculating probabilities N times with randomly assigned weights for each window,
individual set of data points. By calculating probabilities N times with randomly assigned weights for each window,
one can estimate the error as standard deviation between the N bootstrapping runs. For more details see
*Van der Spoel, D. et al. (2010). g_wham—A Free Weighted Histogram Analysis Implementation Including Robust Error and
Autocorrelation Estimates, JCTC, 6(12), 3713-3720*.
To perform bayesian bootstrapping in WHAM, use the ```-bt <RUNS>``` flag to perform <RUNS> individual bootstrapping
runs. The error estimates of bin probabilities and free energy will be given as separate column (+/-) in the output file.
If no error analysis is performed, these columns are set to 0.0.
runs. The error estimates of bin probabilities and free energy will be given as standard error (SE) in a
separate column (+/-) in the output file. If no error analysis is performed, these columns are set to 0.0.
Examples
Autocorrelation analysis
---
The example folder contains input and output files for two simple test systems:
With the ```--uncorr``` flag, WHAM calculates the autocorrelation time ```tau``` for all timeseries and all collective
variables. Timeseries are then filtered based on their highest autocorrelation time to remove correlated samples from
the dataset. This reduces the number of data points but can improve the accuracy of the result.
- 1d: Phi torsion angle of dialanine in vaccum
- 2d: Phi and psi torsion angles of the same system
For filtering, the statistical inefficiency `g` is calculated: ```g = 1 + 2*tau```, and only every `g`th element of the
timeseries is used for unbiasing. A more detailed description of the method can be found in
*Chodera, J.D. et al. (2007). Use of the weighted histogram analysis method for the analysis of simulated and parallel
tempering simulations, JCTC 3(1):26-41*
TODO
---
- Autocorrelation
- Option to output histograms
- Replica exchange
License & Citing
@@ -102,7 +164,7 @@ License & Citing
WHAM is licensed under the GPL-3.0 license. Please read the LICENSE file in this
repository for more information.
There's no publication for this WHAM implementation. However, there is a citeabe DOI. If you use this software for your work, please consider citing it: *Bauer, D, WHAM - An efficient weighted histogram analysis implementation written in Rust, Zenodo. https://doi.org/10.5281/zenodo.1488597*
There's no publication for this WHAM implementation. However, there is a citeabe DOI. If you use this software for your work, please consider citing it: *Bauer, D., WHAM - An efficient weighted histogram analysis implementation written in Rust, Zenodo. https://doi.org/10.5281/zenodo.1488597*
Parts of this work, especially some perfomance optimizations and the I/O format, are inspired by the
implementation of A. Grossfield (*Grossfield, A, WHAM: the weighted histogram analysis method, http://membrane.urmc.rochester.edu/content/wham*).

View File

@@ -1,101 +0,0 @@
#coord1 Free Energy +/- Probability +/-
-3.108600 7.118166 0.000000 0.003561 0.000000
-3.045800 5.331290 0.000000 0.007290 0.000000
-2.983000 3.879241 0.000000 0.013048 0.000000
-2.920200 2.968659 0.000000 0.018796 0.000000
-2.857400 1.942494 0.000000 0.028362 0.000000
-2.794600 1.418408 0.000000 0.034994 0.000000
-2.731800 1.169398 0.000000 0.038668 0.000000
-2.669000 0.843024 0.000000 0.044073 0.000000
-2.606200 0.636021 0.000000 0.047887 0.000000
-2.543400 0.720123 0.000000 0.046299 0.000000
-2.480600 1.066434 0.000000 0.040297 0.000000
-2.417800 1.509855 0.000000 0.033734 0.000000
-2.355000 2.015549 0.000000 0.027544 0.000000
-2.292200 2.226942 0.000000 0.025306 0.000000
-2.229400 2.573176 0.000000 0.022026 0.000000
-2.166600 2.591946 0.000000 0.021861 0.000000
-2.103800 2.476835 0.000000 0.022893 0.000000
-2.041000 2.513953 0.000000 0.022555 0.000000
-1.978200 2.470032 0.000000 0.022956 0.000000
-1.915400 2.227585 0.000000 0.025299 0.000000
-1.852600 2.105310 0.000000 0.026570 0.000000
-1.789800 1.778779 0.000000 0.030286 0.000000
-1.727000 1.447127 0.000000 0.034593 0.000000
-1.664200 0.906384 0.000000 0.042968 0.000000
-1.601400 0.304263 0.000000 0.054699 0.000000
-1.538600 0.247397 0.000000 0.055960 0.000000
-1.475800 0.000000 0.000000 0.061795 0.000000
-1.413000 0.536142 0.000000 0.049843 0.000000
-1.350200 1.428046 0.000000 0.034859 0.000000
-1.287400 2.398149 0.000000 0.023627 0.000000
-1.224600 3.784647 0.000000 0.013552 0.000000
-1.161800 5.696583 0.000000 0.006297 0.000000
-1.099000 7.656158 0.000000 0.002870 0.000000
-1.036200 9.955000 0.000000 0.001142 0.000000
-0.973400 12.357407 0.000000 0.000436 0.000000
-0.910600 14.954254 0.000000 0.000154 0.000000
-0.847800 17.745327 0.000000 0.000050 0.000000
-0.785000 20.559461 0.000000 0.000016 0.000000
-0.722200 22.783560 0.000000 0.000007 0.000000
-0.659400 25.175865 0.000000 0.000003 0.000000
-0.596600 26.419980 0.000000 0.000002 0.000000
-0.533800 27.885999 0.000000 0.000001 0.000000
-0.471000 29.029334 0.000000 0.000001 0.000000
-0.408200 30.366556 0.000000 0.000000 0.000000
-0.345400 31.615252 0.000000 0.000000 0.000000
-0.282600 32.781716 0.000000 0.000000 0.000000
-0.219800 33.731378 0.000000 0.000000 0.000000
-0.157000 34.476014 0.000000 0.000000 0.000000
-0.094200 35.387300 0.000000 0.000000 0.000000
-0.031400 35.579756 0.000000 0.000000 0.000000
0.031400 35.520194 0.000000 0.000000 0.000000
0.094200 35.344763 0.000000 0.000000 0.000000
0.157000 34.886460 0.000000 0.000000 0.000000
0.219800 33.633998 0.000000 0.000000 0.000000
0.282600 32.686511 0.000000 0.000000 0.000000
0.345400 31.211777 0.000000 0.000000 0.000000
0.408200 29.680451 0.000000 0.000000 0.000000
0.471000 28.032242 0.000000 0.000001 0.000000
0.533800 26.439599 0.000000 0.000002 0.000000
0.596600 24.432467 0.000000 0.000003 0.000000
0.659400 22.309893 0.000000 0.000008 0.000000
0.722200 20.186040 0.000000 0.000019 0.000000
0.785000 18.295046 0.000000 0.000040 0.000000
0.847800 16.221232 0.000000 0.000093 0.000000
0.910600 14.251911 0.000000 0.000204 0.000000
0.973400 12.562353 0.000000 0.000402 0.000000
1.036200 11.184029 0.000000 0.000698 0.000000
1.099000 10.040677 0.000000 0.001103 0.000000
1.161800 9.369064 0.000000 0.001444 0.000000
1.224600 9.095392 0.000000 0.001612 0.000000
1.287400 9.256035 0.000000 0.001511 0.000000
1.350200 9.814068 0.000000 0.001208 0.000000
1.413000 10.945207 0.000000 0.000768 0.000000
1.475800 12.518898 0.000000 0.000409 0.000000
1.538600 14.409744 0.000000 0.000191 0.000000
1.601400 16.451965 0.000000 0.000084 0.000000
1.664200 18.589166 0.000000 0.000036 0.000000
1.727000 20.669504 0.000000 0.000016 0.000000
1.789800 22.804791 0.000000 0.000007 0.000000
1.852600 24.531690 0.000000 0.000003 0.000000
1.915400 26.182612 0.000000 0.000002 0.000000
1.978200 27.279439 0.000000 0.000001 0.000000
2.041000 28.571380 0.000000 0.000001 0.000000
2.103800 29.328969 0.000000 0.000000 0.000000
2.166600 29.897739 0.000000 0.000000 0.000000
2.229400 30.320908 0.000000 0.000000 0.000000
2.292200 30.065391 0.000000 0.000000 0.000000
2.355000 29.809137 0.000000 0.000000 0.000000
2.417800 29.340601 0.000000 0.000000 0.000000
2.480600 28.473346 0.000000 0.000001 0.000000
2.543400 27.679353 0.000000 0.000001 0.000000
2.606200 26.453160 0.000000 0.000002 0.000000
2.669000 24.424649 0.000000 0.000003 0.000000
2.731800 22.285933 0.000000 0.000008 0.000000
2.794600 19.958407 0.000000 0.000021 0.000000
2.857400 17.606291 0.000000 0.000053 0.000000
2.920200 15.420138 0.000000 0.000128 0.000000
2.983000 13.069599 0.000000 0.000328 0.000000
3.045800 11.036085 0.000000 0.000740 0.000000
3.108600 9.015657 0.000000 0.001664 0.000000

View File

@@ -0,0 +1,9 @@
Example data: Dialanine torsion angle at 300K
wham.out was produced with the folling settings:
wham --max pi --min -pi --cyclic\
--bins 100 \
-T 300 \
-f metadata.dat \
-o wham.out \
--bt 100 --seed 1234

101
example/1d_cyclic/wham.out Normal file
View File

@@ -0,0 +1,101 @@
#coord1 Free Energy +/- Probability +/-
-3.110177 7.158102 0.071438 0.003494 0.000068
-3.047345 5.365727 0.069449 0.007168 0.000135
-2.984513 3.873190 0.067495 0.013039 0.000232
-2.921681 2.953162 0.067589 0.018855 0.000343
-2.858849 1.949554 0.064640 0.028195 0.000480
-2.796017 1.391747 0.063570 0.035261 0.000584
-2.733186 1.128270 0.061710 0.039189 0.000620
-2.670354 0.839970 0.060445 0.043991 0.000667
-2.607522 0.624769 0.060762 0.047955 0.000739
-2.544690 0.663757 0.060786 0.047211 0.000731
-2.481858 1.051932 0.059774 0.040407 0.000617
-2.419026 1.463048 0.060422 0.034267 0.000527
-2.356194 1.990616 0.055242 0.027734 0.000360
-2.293363 2.190692 0.044458 0.025597 0.000210
-2.230531 2.553036 0.042144 0.022136 0.000159
-2.167699 2.572522 0.043313 0.021964 0.000165
-2.104867 2.472360 0.039799 0.022863 0.000157
-2.042035 2.517562 0.036792 0.022453 0.000203
-1.979203 2.469778 0.036324 0.022887 0.000237
-1.916372 2.223125 0.049496 0.025266 0.000565
-1.853540 2.080157 0.039205 0.026756 0.000453
-1.790708 1.793841 0.041488 0.030011 0.000558
-1.727876 1.458784 0.042362 0.034326 0.000565
-1.665044 0.914441 0.037561 0.042697 0.000682
-1.602212 0.303813 0.028701 0.054540 0.000761
-1.539380 0.268335 0.025765 0.055321 0.000771
-1.476549 0.000000 0.023371 0.061604 0.000866
-1.413717 0.537061 0.026370 0.049671 0.000715
-1.350885 1.439940 0.025620 0.034586 0.000547
-1.288053 2.391838 0.027653 0.023614 0.000393
-1.225221 3.779470 0.027093 0.013538 0.000229
-1.162389 5.685555 0.027762 0.006305 0.000111
-1.099557 7.661896 0.028996 0.002855 0.000051
-1.036726 9.946594 0.032118 0.001142 0.000021
-0.973894 12.359408 0.042801 0.000434 0.000009
-0.911062 14.954655 0.052443 0.000153 0.000004
-0.848230 17.742242 0.063261 0.000050 0.000002
-0.785398 20.555785 0.068133 0.000016 0.000001
-0.722566 22.811160 0.075518 0.000007 0.000000
-0.659734 25.178094 0.086226 0.000003 0.000000
-0.596903 26.442288 0.091783 0.000002 0.000000
-0.534071 27.897565 0.093175 0.000001 0.000000
-0.471239 29.062473 0.096027 0.000001 0.000000
-0.408407 30.384521 0.096855 0.000000 0.000000
-0.345575 31.638454 0.099956 0.000000 0.000000
-0.282743 32.817727 0.107689 0.000000 0.000000
-0.219911 33.770369 0.108946 0.000000 0.000000
-0.157080 34.505503 0.112902 0.000000 0.000000
-0.094248 35.431659 0.120856 0.000000 0.000000
-0.031416 35.615810 0.122046 0.000000 0.000000
0.031416 35.561946 0.118225 0.000000 0.000000
0.094248 35.382089 0.108652 0.000000 0.000000
0.157080 34.934827 0.118384 0.000000 0.000000
0.219911 33.673460 0.119195 0.000000 0.000000
0.282743 32.731563 0.121592 0.000000 0.000000
0.345575 31.261855 0.129698 0.000000 0.000000
0.408407 29.717377 0.139789 0.000000 0.000000
0.471239 28.076620 0.137289 0.000001 0.000000
0.534071 26.481097 0.135712 0.000002 0.000000
0.596903 24.487358 0.135156 0.000003 0.000000
0.659734 22.344251 0.131027 0.000008 0.000000
0.722566 20.241543 0.133145 0.000018 0.000001
0.785398 18.341869 0.131952 0.000039 0.000002
0.848230 16.261582 0.134676 0.000091 0.000006
0.911062 14.301801 0.134417 0.000199 0.000013
0.973894 12.603788 0.131018 0.000394 0.000024
1.036726 11.249601 0.131276 0.000678 0.000043
1.099557 10.087886 0.132498 0.001079 0.000070
1.162389 9.443303 0.132990 0.001398 0.000092
1.225221 9.152799 0.132250 0.001570 0.000104
1.288053 9.331937 0.133178 0.001462 0.000099
1.350885 9.905546 0.133357 0.001161 0.000078
1.413717 11.042050 0.133807 0.000736 0.000051
1.476549 12.598167 0.132722 0.000395 0.000027
1.539380 14.520167 0.131816 0.000183 0.000012
1.602212 16.569783 0.131773 0.000080 0.000005
1.665044 18.687390 0.132601 0.000034 0.000002
1.727876 20.775408 0.133867 0.000015 0.000001
1.790708 22.905200 0.129265 0.000006 0.000000
1.853540 24.643852 0.128894 0.000003 0.000000
1.916372 26.301740 0.130266 0.000002 0.000000
1.979203 27.372071 0.128620 0.000001 0.000000
2.042035 28.697726 0.133263 0.000001 0.000000
2.104867 29.417901 0.133513 0.000000 0.000000
2.167699 30.008351 0.130925 0.000000 0.000000
2.230531 30.406016 0.124398 0.000000 0.000000
2.293363 30.171275 0.122493 0.000000 0.000000
2.356194 29.884646 0.130278 0.000000 0.000000
2.419026 29.428153 0.130844 0.000000 0.000000
2.481858 28.546982 0.148114 0.000001 0.000000
2.544690 27.757520 0.133045 0.000001 0.000000
2.607522 26.505787 0.134364 0.000001 0.000000
2.670354 24.491866 0.115061 0.000003 0.000000
2.733186 22.320664 0.110036 0.000008 0.000000
2.796017 20.052723 0.107894 0.000020 0.000001
2.858849 17.655650 0.105964 0.000052 0.000002
2.921681 15.471590 0.107005 0.000125 0.000005
2.984513 13.138167 0.099378 0.000318 0.000012
3.047345 11.092386 0.087770 0.000722 0.000022
3.110177 9.065722 0.077092 0.001626 0.000038

View File

@@ -0,0 +1,101 @@
#coord1 Free Energy +/- Probability +/-
-3.110177 7.531315 0.000000 0.003080 0.000000
-3.047345 5.690157 0.000000 0.006443 0.000000
-2.984513 4.243063 0.000000 0.011509 0.000000
-2.921681 3.334686 0.000000 0.016564 0.000000
-2.858849 2.277349 0.000000 0.025309 0.000000
-2.796017 1.723296 0.000000 0.031604 0.000000
-2.733186 1.246264 0.000000 0.038265 0.000000
-2.670354 1.099867 0.000000 0.040578 0.000000
-2.607522 0.771910 0.000000 0.046279 0.000000
-2.544690 0.770616 0.000000 0.046303 0.000000
-2.481858 1.265507 0.000000 0.037971 0.000000
-2.419026 1.562335 0.000000 0.033711 0.000000
-2.356194 1.891577 0.000000 0.029542 0.000000
-2.293363 2.227858 0.000000 0.025816 0.000000
-2.230531 2.488355 0.000000 0.023256 0.000000
-2.167699 2.502265 0.000000 0.023127 0.000000
-2.104867 2.358037 0.000000 0.024503 0.000000
-2.042035 2.278147 0.000000 0.025301 0.000000
-1.979203 2.974067 0.000000 0.019141 0.000000
-1.916372 2.696600 0.000000 0.021393 0.000000
-1.853540 2.361827 0.000000 0.024466 0.000000
-1.790708 1.516746 0.000000 0.034332 0.000000
-1.727876 1.526829 0.000000 0.034194 0.000000
-1.665044 0.884114 0.000000 0.044244 0.000000
-1.602212 0.323912 0.000000 0.055385 0.000000
-1.539380 0.197985 0.000000 0.058253 0.000000
-1.476549 0.000000 0.000000 0.063065 0.000000
-1.413717 0.458247 0.000000 0.052481 0.000000
-1.350885 1.389410 0.000000 0.036131 0.000000
-1.288053 2.386522 0.000000 0.024225 0.000000
-1.225221 3.743253 0.000000 0.014062 0.000000
-1.162389 5.566654 0.000000 0.006770 0.000000
-1.099557 7.822800 0.000000 0.002740 0.000000
-1.036726 10.128719 0.000000 0.001087 0.000000
-0.973894 12.199246 0.000000 0.000474 0.000000
-0.911062 14.488129 0.000000 0.000189 0.000000
-0.848230 16.902310 0.000000 0.000072 0.000000
-0.785398 18.910200 0.000000 0.000032 0.000000
-0.722566 21.241681 0.000000 0.000013 0.000000
-0.659734 22.706373 0.000000 0.000007 0.000000
-0.596903 24.531129 0.000000 0.000003 0.000000
-0.534071 25.936227 0.000000 0.000002 0.000000
-0.471239 27.000262 0.000000 0.000001 0.000000
-0.408407 28.673293 0.000000 0.000001 0.000000
-0.345575 29.335203 0.000000 0.000000 0.000000
-0.282743 30.841118 0.000000 0.000000 0.000000
-0.219911 31.983859 0.000000 0.000000 0.000000
-0.157080 32.144015 0.000000 0.000000 0.000000
-0.094248 33.885395 0.000000 0.000000 0.000000
-0.031416 33.783105 0.000000 0.000000 0.000000
0.031416 34.243727 0.000000 0.000000 0.000000
0.094248 33.975567 0.000000 0.000000 0.000000
0.157080 32.994787 0.000000 0.000000 0.000000
0.219911 32.607398 0.000000 0.000000 0.000000
0.282743 31.401902 0.000000 0.000000 0.000000
0.345575 29.911670 0.000000 0.000000 0.000000
0.408407 28.603574 0.000000 0.000001 0.000000
0.471239 26.925443 0.000000 0.000001 0.000000
0.534071 25.298070 0.000000 0.000002 0.000000
0.596903 23.638560 0.000000 0.000005 0.000000
0.659734 21.156231 0.000000 0.000013 0.000000
0.722566 19.126480 0.000000 0.000029 0.000000
0.785398 17.351953 0.000000 0.000060 0.000000
0.848230 15.135525 0.000000 0.000146 0.000000
0.911062 13.188112 0.000000 0.000319 0.000000
0.973894 11.536983 0.000000 0.000618 0.000000
1.036726 10.158328 0.000000 0.001074 0.000000
1.099557 9.109036 0.000000 0.001636 0.000000
1.162389 8.282343 0.000000 0.002279 0.000000
1.225221 8.022102 0.000000 0.002530 0.000000
1.288053 8.162415 0.000000 0.002391 0.000000
1.350885 8.600135 0.000000 0.002006 0.000000
1.413717 9.837348 0.000000 0.001222 0.000000
1.476549 11.363156 0.000000 0.000663 0.000000
1.539380 13.077849 0.000000 0.000333 0.000000
1.602212 15.353594 0.000000 0.000134 0.000000
1.665044 17.565051 0.000000 0.000055 0.000000
1.727876 19.710884 0.000000 0.000023 0.000000
1.790708 21.721260 0.000000 0.000010 0.000000
1.853540 23.567649 0.000000 0.000005 0.000000
1.916372 25.008817 0.000000 0.000003 0.000000
1.979203 26.405367 0.000000 0.000002 0.000000
2.042035 28.070821 0.000000 0.000001 0.000000
2.104867 28.877213 0.000000 0.000001 0.000000
2.167699 29.378146 0.000000 0.000000 0.000000
2.230531 31.093267 0.000000 0.000000 0.000000
2.293363 30.704994 0.000000 0.000000 0.000000
2.356194 30.563093 0.000000 0.000000 0.000000
2.419026 31.215952 0.000000 0.000000 0.000000
2.481858 30.331416 0.000000 0.000000 0.000000
2.544690 29.005123 0.000000 0.000001 0.000000
2.607522 27.674618 0.000000 0.000001 0.000000
2.670354 24.911788 0.000000 0.000003 0.000000
2.733186 22.746417 0.000000 0.000007 0.000000
2.796017 20.608288 0.000000 0.000016 0.000000
2.858849 18.018280 0.000000 0.000046 0.000000
2.921681 15.949215 0.000000 0.000105 0.000000
2.984513 13.617809 0.000000 0.000268 0.000000
3.047345 11.432193 0.000000 0.000645 0.000000
3.110177 9.461494 0.000000 0.001420 0.000000

File diff suppressed because it is too large Load Diff

Some files were not shown because too many files have changed in this diff Show More