/
LinearPartition-E.patch
84 lines (78 loc) · 14.8 KB
/
LinearPartition-E.patch
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
diff --git a/Makefile b/Makefile
index 49a7aca..a4c6e20 100755
--- a/Makefile
+++ b/Makefile
@@ -6,16 +6,17 @@
################################
CC=g++
-DEPS=src/bpp.cpp src/LinearPartition.h src/Utils/energy_parameter.h src/Utils/feature_weight.h src/Utils/intl11.h src/Utils/intl21.h src/Utils/intl22.h src/Utils/utility_v.h src/Utils/utility.h
+DEPS=src/bpp.cpp src/LinearPartition.h src/Utils/energy_parameter.h src/Utils/feature_weight.h src/Utils/feature_weight_e.h src/Utils/intl11.h src/Utils/intl21.h src/Utils/intl22.h src/Utils/utility_v.h src/Utils/utility.h
CFLAGS=-std=c++11 -O3
.PHONY : clean linearpartition
-objects=bin/linearpartition_v bin/linearpartition_c
+objects=bin/linearpartition_v bin/linearpartition_c bin/linearpartition_e
linearpartition: src/LinearPartition.cpp $(DEPS)
chmod +x linearpartition draw_bpp_plot draw_heatmap
mkdir -p bin
$(CC) src/LinearPartition.cpp $(CFLAGS) -Dlpv -o bin/linearpartition_v
$(CC) src/LinearPartition.cpp $(CFLAGS) -o bin/linearpartition_c
+ $(CC) src/LinearPartition.cpp $(CFLAGS) -Dlpe -o bin/linearpartition_e
clean:
-rm $(objects)
\ No newline at end of file
diff --git a/src/Utils/feature_weight_e.h b/src/Utils/feature_weight_e.h
new file mode 100644
index 0000000..5bbb480
--- /dev/null
+++ b/src/Utils/feature_weight_e.h
@@ -0,0 +1,36 @@
+/*
+ *feature_weight_e.h*
+ Feature weights corresponding to EternaFold parameters.
+
+ author: Hannah Wayment-Steele
+ last updated: 08/2020
+*/
+
+
+#ifndef FASTCKY_W
+#define FASTCKY_W
+double multi_base = -0.1678166237;
+double multi_unpaired = -0.1409575697;
+double multi_paired = -1.310997242;
+double external_unpaired = 0.03755355437;
+double external_paired = -1.85767603;
+double base_pair[25] = {0.0,0.0,0.0,0.5327179945,0.0,0.0,0.0,2.122319587,0.0,0.0,0.0,2.122319587,0.0,-0.1412363884,0.0,0.5327179945,0.0,-0.1412363884,0.0,0.0,0.0,0.0,0.0,0.0,0.0};
+double internal_1x1_nucleotides[25] = {0.1837156984,0.2226675311,-0.07277379208,-0.09832713256,0.0,0.2226675311,-0.1797487299,0.03108743734,-0.1597067764,0.0,-0.07277379208,0.03108743734,0.09609737908,-0.1602012398,0.0,-0.09832713256,-0.1597067764,-0.1602012398,0.2459086952,0.0,0.0,0.0,0.0,0.0,0.0};
+double helix_stacking[625] = {0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.02764044118,0.0,0.0,0.0,0.4696145329,0.0,0.0,0.0,0.776993008,0.0,-0.1104463524,0.0,0.2864048526,0.0,0.4540812927,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.4468032603,0.0,0.0,0.0,1.029274484,0.0,0.0,0.0,0.3929919526,0.0,-0.3459365498,0.0,0.4696145329,0.0,0.6287656792,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.6836104681,0.0,0.0,0.0,0.2732504364,0.0,0.0,0.0,1.029274484,0.0,0.3765137104,0.0,0.4696145329,0.0,0.4999765137,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.1829584868,0.0,0.0,0.0,0.4999765137,0.0,0.0,0.0,0.6287656792,0.0,0.3173690216,0.0,0.4540812927,0.0,-0.0446909317,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.5401517598,0.0,0.0,0.0,0.6836104681,0.0,0.0,0.0,0.4468032603,0.0,-0.135286073,0.0,-0.02764044118,0.0,-0.1829584868,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.135286073,0.0,0.0,0.0,0.3765137104,0.0,0.0,0.0,-0.3459365498,0.0,0.451152827,0.0,-0.1104463524,0.0,0.3173690216,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0};
+double terminal_mismatch[625] = {0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.1586543853,-0.2902029646,-0.5468962138,-0.3528072539,0.0,0.269564011,-0.0418594662,0.004971123862,-0.1827469251,0.0,0.2995183939,-0.2012943035,-0.315652769,-0.3808878513,0.0,0.1919565022,0.05728586921,-0.02952347661,0.4718429954,0.0, 0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.09659272509,-0.1200138647,-0.2002381523,0.06715107298,0.0,0.2017879572,-0.04392645962,0.04095549347,0.232043656,0.0,0.6152857958,-0.8105483681,0.002807443234,-0.5036549226,0.0,-0.1424978609,0.08430415936,-0.3245449555,0.1919151072,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.4091277793,-0.4170328747,-0.06698906591,-0.1061375196,0.0,-0.008837262959,-0.2164860416,-0.3015065085,0.1299463362,0.0,0.8245824224,-0.1512801456,0.06736172344,-0.2483687541,0.0,-0.3357792151,0.1235612671,-0.2161815311,0.3472537064,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.255998158,0.03250369155,-0.01800555649,-0.1354303386,0.0,0.05490470564,-0.05390857182,-0.06432237365,-0.06201131973,0.0,0.08776852488,-0.05139215521,-0.106355729,0.00269328788,0.0,-0.05591162221,-0.04981094068,-0.2077227846,-0.1310127278,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.05141211647,0.1621047955,-0.05925449632,-0.3058060186,0.0,0.1510784363,0.03393908063,0.04203397301,-0.006430315313,0.0,0.1309705558,-0.005036020835,0.2390136618,-0.2910511859,0.0,-0.3399086304,-0.1941963663,-0.3600733328,0.1602979785,0.0, 0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.1404687483,0.049444791,-0.186536165,0.1357420138,0.0,-0.1621545287,-0.1294710799,0.02763302429,0.1625389898,0.0,0.4013932363,-0.009318323141,-0.1705770853,-0.01749332984,0.0,-0.1102881768,-0.2890183689,-0.07852966316,-0.142694209,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0};
+double bulge_0x1_nucleotides[5] = {-0.1798668952,0.2158414683,-0.1360090628,0.268592261,0.0};
+double helix_closing[25] = {0.0,0.0,0.0,-1.602693751,0.0,0.0,0.0,-1.199564201,0.0,0.0,0.0,-1.483170995,0.0,-1.769088545,0.0,-1.332114253,0.0,-1.145755799,0.0,0.0,0.0,0.0,0.0,0.0,0.0};
+double dangle_left[125] = {0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.3993589451,-0.2408605898,-0.06731640015,-0.008245797416,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.4525488126,-0.4430884346,-0.01683782679,-0.1094610743,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.4888081949,-0.2396554394,0.05406231407,-0.5504680093,0.0,0.0,0.0,0.0,0.0,0.0,-0.1657104844,-0.157383192,-0.1468725266,-0.3085738355,0.0,0.0,0.0,0.0,0.0,0.0,-0.4749381413,-0.2370847341,-0.4723950953,-0.252040466,0.0,0.0,0.0,0.0,0.0,0.0,-0.2673277995,-0.08994780847,0.04311136881,-0.1593516373,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0};
+double dangle_right[125] = {0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.1092729217,-0.3493057962,-0.3860230963,0.01334883966,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.2849809668,0.06753686781,-0.1296601497,-0.2581790752,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,-0.09267992263,-0.1006649989,-0.2346025594,-0.4124780161,0.0,0.0,0.0,0.0,0.0,0.0,-0.181273482,-0.2376225672,-0.09744253428,-0.2669369855,0.0,0.0,0.0,0.0,0.0,0.0,-0.1720675441,-0.4276012742,-0.2796275066,-0.3404954982,0.0,0.0,0.0,0.0,0.0,0.0,-0.09730378801,-0.1402949362,-0.2132405912,-0.07119569275,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0};
+double internal_explicit[21] = {0.0,0.0,0.0,0.0,0.0,0.1087190703,-0.1094777733,-0.2558696183,-0.08235905722,0.0,0.08550718036,0.04422915709,0.0574689997,0.0,0.0,0.03963042785,0.2419582591,0.0,0.0,0.0,-0.06591705055};
+double hairpin_length[31] = {-2.833040024,-4.197257197,-4.049986717,-3.023332787,-1.39930522,-1.054570104,-1.086421309,-0.982678683,-1.41138261,-1.639503618,-1.851258516,-2.136310756,-2.34848385,-2.639172048,-2.892747884,-3.013902282,-3.092678922,-3.116771955,-3.161669086,-3.216435677,-3.261627968,-3.277618192,-3.276476008,-3.264809733,-3.242474502,-3.211948442,-3.177487653,-3.144491766,-3.110972356,-3.076751282,-3.045836654};
+double bulge_length[31] = {0.0,-0.920316622,-2.009191015,-2.896090024,-3.50500818,-3.875673028,-4.095815732,-4.208515353,-4.252726461,-4.478718995,-4.740450924,-4.926752683,-5.061734723,-5.127826263,-5.214171072,-5.264829236,-5.313353555,-5.356344107,-5.392593013,-5.419682264,-5.453101725,-5.470635696,-5.482155819,-5.482497612,-5.483202129,-5.484268521,-5.479068025,-5.47652413,-5.47395731,-5.471109108,-5.470155293};
+double internal_length[31] = {0.0,0.0,-0.345020502,-0.798760074,-1.143021873,-1.316921234,-1.452690695,-1.535487563,-1.683273183,-1.780441353,-1.986254543,-2.335789199,-2.527779016,-2.681245657,-2.75290067,-2.824705992,-2.883486824,-2.958477596,-3.016059015,-3.035815717,-3.057453413,-3.076034873,-3.107772289,-3.130278855,-3.141304061,-3.14202303,-3.144856457,-3.156706136,-3.164174509,-3.16948557,-3.172132284};
+double internal_symmetric_length[16] = {0.0,0.089130266,0.069541462,-0.035554523,-0.180280936,-0.259090297,-0.425844035,-0.524984188,-0.572013829,-0.636356186,-0.665844794,-0.684341631,-0.697466116,-0.704242001,-0.705848654,-0.706663366};
+double internal_asymmetry[29] = {0.0,-0.434150768,-0.988647308,-1.501240421,-1.844922335,-1.986256761,-2.041066577,-2.171950668,-2.280710602,-2.377744933,-2.44651694,-2.534021026,-2.592057002,-2.63684331,-2.673547186,-2.693239044,-2.704551064,-2.731772066,-2.754116508,-2.769911389,-2.785011722,-2.787995996,-2.782538517,-2.781756278,-2.782490795,-2.783816609,-2.784628174,-2.785411015,-2.785471292};
+double hairpin_length_at_least[31] = {-2.833040024,-1.364217173,0.1472704798,1.02665393,1.624027567,0.3447351159,-0.03185120443,0.1037426259,-0.4287039276,-0.228121008,-0.2117548975,-0.2850522401,-0.2121730939,-0.2906881976,-0.2535758364,-0.1211543976,-0.07877664033,-0.02409303342,-0.04489713102,-0.05476659036,-0.04519229163,-0.0159902242,0.001142184893,0.01166627441,0.02233523082,0.03052606024,0.03446078876,0.03299588703,0.03351941082,0.03422107388,0.03091462776};
+double bulge_length_at_least[31] = {0,-0.9203166218,-1.088874393,-0.8868990087,-0.6089181569,-0.3706648473,-0.2201427043,-0.1126996205,-0.04421110857,-0.2259925335, -0.2617319299, -0.1863017584, -0.1349820405, -0.06609153928, -0.08634480938, -0.05065816443, -0.04852431806, -0.0429905529, -0.0362489059, -0.02708925061, -0.03341946057, -0.01753397155, -0.01152012338, -0.0003417922264, -0.0007045177054, -0.001066391869, 0.005200496185, 0.002543895308, 0.002566819644, 0.00284820206, 0.0009538152694};
+double internal_length_at_least[31] = {0.0,0.0,-0.345020502,-0.4537395723,-0.344261799,-0.1738993611,-0.1357694609,-0.08279686733,-0.1477856205,-0.09716816966, -0.20581319, -0.3495346562, -0.1919898171, -0.1534666412, -0.07165501224, -0.07180532248, -0.05878083178, -0.07499077249, -0.05758141914, -0.01975670154, -0.02163769616, -0.01858146005, -0.03173741584, -0.02250656646, -0.01102520591, -0.0007189686725, -0.002833426606, -0.01184967974, -0.007468372851, -0.00531106075, -0.002646713724};
+double internal_symmetric_length_at_least[16] = {0.0,0.08913026599,-0.01958880431,-0.1050959847,-0.1447264125,-0.07880936196,-0.1667537379,-0.09914015307,-0.04702964029,-0.06434235756, -0.02948860748, -0.01849683765, -0.01312448468, -0.006775884947, -0.00160665253, -0.0008147124093};
+double internal_asymmetry_at_least[29] = {0.0,-0.434150768,-0.5544965399,-0.5125931132,-0.343681914,-0.1413344261,-0.05480981619,-0.1308840911,-0.1087599338,-0.09703433049, -0.06877200709, -0.08750408608, -0.05803597568, -0.04478630849, -0.03670387583, -0.01969185841, -0.01131202, -0.02722100209, -0.02234444178, -0.01579488106, -0.01510033223, -0.002984274291, 0.005457479108, 0.0007822391233, -0.0007345172669, -0.001325813753, -0.0008115654138, -0.0007828410789, -6.027669541e-05};
+#endif
\ No newline at end of file
diff --git a/src/Utils/utility.h b/src/Utils/utility.h
index 0601a4d..a5866b1 100755
--- a/src/Utils/utility.h
+++ b/src/Utils/utility.h
@@ -13,7 +13,11 @@
#include <cstring>
#include <assert.h>
+#ifdef lpe
+#include "feature_weight_e.h"
+#else
#include "feature_weight.h"
+#endif
#define INF 1000000007