libsequence  1.9.5
Grantham.cc
1 /*
2 
3  Copyright (C) 2003-2009 Kevin Thornton, krthornt[]@[]uci.edu
4 
5  Remove the brackets to email me.
6 
7  This file is part of libsequence.
8 
9  libsequence is free software: you can redistribute it and/or modify
10  it under the terms of the GNU General Public License as published by
11  the Free Software Foundation, either version 3 of the License, or
12  (at your option) any later version.
13 
14  libsequence is distributed in the hope that it will be useful,
15  but WITHOUT ANY WARRANTY; without even the implied warranty of
16  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17  GNU General Public License for more details.
18 
19  You should have received a copy of the GNU General Public License
20  long with libsequence. If not, see <http://www.gnu.org/licenses/>.
21 
22 */
23 
24 #include <Sequence/Grantham.hpp>
25 #include <cctype>
26 #include <limits>
27 namespace Sequence
28 {
29  Grantham::Grantham (const double stop) :
30  stopweight(stop)
31  {
32  D['T'-38]['M'-38]=81.03969279;
33  D['T'-38]['N'-38]=64.77428784;
34  D['T'-38]['P'-38]=37.56362476;
35  D['T'-38]['Q'-38]=41.10476252;
36  D['T'-38]['R'-38]=70.97102635;
37  D['T'-38]['S'-38]=57.74902766;
38  D['T'-38]['T'-38]=0.;
39  D['T'-38]['V'-38]=69.49621703;
40  D['T'-38]['W'-38]=150.752701;
41  D['T'-38]['Y'-38]=92.24764162;
42  D['E'-38]['*'-38]=stopweight;
43  D['E'-38]['A'-38]=106.7125514;
44  D['E'-38]['C'-38]=169.4374593;
45  D['E'-38]['D'-38]=44.60463645;
46  D['E'-38]['E'-38]=0.;
47  D['E'-38]['F'-38]=140.2120632;
48  D['E'-38]['G'-38]=97.85190048;
49  D['E'-38]['H'-38]=40.79415368;
50  D['E'-38]['I'-38]=134.1621451;
51  D['E'-38]['K'-38]=56.86795522;
52  D['E'-38]['L'-38]=138.3430673;
53  D['E'-38]['M'-38]=126.0849231;
54  D['E'-38]['N'-38]=40.85895231;
55  D['E'-38]['P'-38]=93.73085383;
56  D['E'-38]['Q'-38]=29.27370624;
57  D['E'-38]['R'-38]=54.01884479;
58  D['E'-38]['S'-38]=79.78756921;
59  D['E'-38]['T'-38]=65.50131254;
60  D['E'-38]['V'-38]=121.3285063;
61  D['E'-38]['W'-38]=160.4382217;
62  D['E'-38]['Y'-38]=122.7768651;
63  D['P'-38]['*'-38]=stopweight;
64  D['K'-38]['*'-38]=stopweight;
65  D['P'-38]['A'-38]=26.87433834;
66  D['P'-38]['C'-38]=168.5906208;
67  D['P'-38]['D'-38]=107.9096634;
68  D['P'-38]['E'-38]=93.73085383;
69  D['P'-38]['F'-38]=113.7271541;
70  D['P'-38]['G'-38]=41.62910033;
71  D['P'-38]['H'-38]=76.27712277;
72  D['P'-38]['I'-38]=95.3760884;
73  D['F'-38]['*'-38]=stopweight;
74  D['P'-38]['K'-38]=102.7139942;
75  D['P'-38]['L'-38]=97.77619789;
76  D['P'-38]['M'-38]=86.59479895;
77  D['P'-38]['N'-38]=90.15739727;
78  D['P'-38]['P'-38]=0.;
79  D['K'-38]['A'-38]=105.5709011;
80  D['P'-38]['Q'-38]=75.1358782;
81  D['P'-38]['R'-38]=102.7148709;
82  D['K'-38]['C'-38]=201.5793361;
83  D['P'-38]['S'-38]=73.35256887;
84  D['K'-38]['D'-38]=101.4568985;
85  D['P'-38]['T'-38]=37.56362476;
86  D['K'-38]['E'-38]=56.86795522;
87  D['K'-38]['F'-38]=102.141346;
88  D['P'-38]['V'-38]=67.78659846;
89  D['K'-38]['G'-38]=126.4579761;
90  D['P'-38]['W'-38]=182.7653305;
91  D['K'-38]['H'-38]=32.40288701;
92  D['K'-38]['I'-38]=101.6123327;
93  D['P'-38]['Y'-38]=109.6155723;
94  D['A'-38]['*'-38]=stopweight;
95  D['K'-38]['K'-38]=0.;
96  D['K'-38]['L'-38]=106.3355126;
97  D['K'-38]['M'-38]=94.4900972;
98  D['K'-38]['N'-38]=93.88268262;
99  D['K'-38]['P'-38]=102.7139942;
100  D['K'-38]['Q'-38]=53.22834015;
101  D['F'-38]['A'-38]=112.5814896;
102  D['K'-38]['R'-38]=26.00394626;
103  D['K'-38]['S'-38]=120.5325805;
104  D['F'-38]['C'-38]=204.3884874;
105  D['K'-38]['T'-38]=77.74105567;
106  D['F'-38]['D'-38]=176.5261931;
107  D['F'-38]['E'-38]=140.2120632;
108  D['K'-38]['V'-38]=96.9973994;
109  D['F'-38]['F'-38]=0.;
110  D['K'-38]['W'-38]=157.923034;
111  D['F'-38]['G'-38]=153.125519;
112  D['F'-38]['H'-38]=99.99505513;
113  D['K'-38]['Y'-38]=84.78644482;
114  D['F'-38]['I'-38]=21.27701376;
115  D['F'-38]['K'-38]=102.141346;
116  D['F'-38]['L'-38]=21.82392074;
117  D['F'-38]['M'-38]=28.52784356;
118  D['F'-38]['N'-38]=158.1121886;
119  D['F'-38]['P'-38]=113.7271541;
120  D['A'-38]['A'-38]=0.;
121  D['F'-38]['Q'-38]=115.5869365;
122  D['F'-38]['R'-38]=97.03272116;
123  D['A'-38]['C'-38]=195.0037439;
124  D['F'-38]['S'-38]=154.8079023;
125  D['A'-38]['D'-38]=125.7485902;
126  D['F'-38]['T'-38]=102.8587191;
127  D['A'-38]['E'-38]=106.7125514;
128  D['A'-38]['F'-38]=112.5814896;
129  D['F'-38]['V'-38]=49.93518986;
130  D['A'-38]['G'-38]=59.99535652;
131  D['F'-38]['W'-38]=142.6771181;
132  D['A'-38]['H'-38]=85.49373612;
133  D['A'-38]['I'-38]=93.66243338;
134  D['V'-38]['*'-38]=stopweight;
135  D['F'-38]['Y'-38]=21.60966489;
136  D['A'-38]['K'-38]=105.5709011;
137  D['A'-38]['L'-38]=96.18707223;
138  D['A'-38]['M'-38]=84.43958222;
139  D['A'-38]['N'-38]=110.418026;
140  D['A'-38]['P'-38]=26.87433834;
141  D['A'-38]['Q'-38]=90.76116997;
142  D['A'-38]['R'-38]=111.2645916;
143  D['A'-38]['S'-38]=99.13249708;
144  D['A'-38]['T'-38]=58.0233276;
145  D['A'-38]['V'-38]=64.43025059;
146  D['A'-38]['W'-38]=201.5127248;
147  D['A'-38]['Y'-38]=111.5882137;
148  D['Q'-38]['*'-38]=stopweight;
149  D['V'-38]['A'-38]=64.43025059;
150  D['V'-38]['C'-38]=191.232501;
151  D['V'-38]['D'-38]=152.013397;
152  D['V'-38]['E'-38]=121.3285063;
153  D['V'-38]['F'-38]=49.93518986;
154  D['V'-38]['G'-38]=108.7873884;
155  D['V'-38]['H'-38]=83.892947;
156  D['V'-38]['I'-38]=29.60907063;
157  D['L'-38]['*'-38]=stopweight;
158  D['V'-38]['K'-38]=96.9973994;
159  D['V'-38]['L'-38]=31.78478998;
160  D['V'-38]['M'-38]=21.52179938;
161  D['V'-38]['N'-38]=132.8777227;
162  D['V'-38]['P'-38]=67.78659846;
163  D['V'-38]['Q'-38]=96.32571878;
164  D['Q'-38]['A'-38]=90.76116997;
165  D['V'-38]['R'-38]=95.79707649;
166  D['V'-38]['S'-38]=123.0339581;
167  D['Q'-38]['C'-38]=154.2308431;
168  D['V'-38]['T'-38]=69.49621703;
169  D['Q'-38]['D'-38]=61.28447445;
170  D['Q'-38]['E'-38]=29.27370624;
171  D['V'-38]['V'-38]=0.;
172  D['V'-38]['W'-38]=162.8552264;
173  D['Q'-38]['F'-38]=115.5869365;
174  D['Q'-38]['G'-38]=87.1663871;
175  D['V'-38]['Y'-38]=54.66278181;
176  D['Q'-38]['H'-38]=24.08398867;
177  D['Q'-38]['I'-38]=108.566243;
178  D['G'-38]['*'-38]=stopweight;
179  D['Q'-38]['K'-38]=53.22834015;
180  D['Q'-38]['L'-38]=112.4414772;
181  D['Q'-38]['M'-38]=100.8991435;
182  D['Q'-38]['N'-38]=45.75219661;
183  D['Q'-38]['P'-38]=75.1358782;
184  D['Q'-38]['Q'-38]=0.;
185  D['L'-38]['A'-38]=96.18707223;
186  D['Q'-38]['R'-38]=42.81393193;
187  D['Q'-38]['S'-38]=68.19787424;
188  D['L'-38]['C'-38]=197.4289714;
189  D['Q'-38]['T'-38]=41.10476252;
190  D['L'-38]['D'-38]=171.7572022;
191  D['L'-38]['E'-38]=138.3430673;
192  D['Q'-38]['V'-38]=96.32571878;
193  D['L'-38]['F'-38]=21.82392074;
194  D['Q'-38]['W'-38]=141.5619982;
195  D['L'-38]['G'-38]=137.6916434;
196  D['L'-38]['H'-38]=98.69301471;
197  D['Q'-38]['Y'-38]=98.78325151;
198  D['L'-38]['I'-38]=4.855121198;
199  D['L'-38]['K'-38]=106.3355126;
200  D['L'-38]['L'-38]=0.;
201  D['L'-38]['M'-38]=14.30316651;
202  D['L'-38]['N'-38]=152.3309797;
203  D['L'-38]['P'-38]=97.77619789;
204  D['G'-38]['A'-38]=59.99535652;
205  D['L'-38]['Q'-38]=112.4414772;
206  D['L'-38]['R'-38]=101.8802928;
207  D['G'-38]['C'-38]=158.2317729;
208  D['L'-38]['S'-38]=144.0792346;
209  D['G'-38]['D'-38]=93.76755083;
210  D['L'-38]['T'-38]=92.35427468;
211  D['G'-38]['E'-38]=97.85190048;
212  D['G'-38]['F'-38]=153.125519;
213  D['L'-38]['V'-38]=31.78478998;
214  D['G'-38]['G'-38]=0.;
215  D['L'-38]['W'-38]=150.0095229;
216  D['G'-38]['H'-38]=97.53340215;
217  D['G'-38]['I'-38]=135.4188097;
218  D['L'-38]['Y'-38]=35.6773121;
219  D['G'-38]['K'-38]=126.4579761;
220  D['G'-38]['L'-38]=137.6916434;
221  D['G'-38]['M'-38]=126.9448727;
222  D['G'-38]['N'-38]=79.34585612;
223  D['G'-38]['P'-38]=41.62910033;
224  D['G'-38]['Q'-38]=87.1663871;
225  D['G'-38]['R'-38]=125.1291926;
226  D['G'-38]['S'-38]=55.26736276;
227  D['G'-38]['T'-38]=59.15645198;
228  D['G'-38]['V'-38]=108.7873884;
229  D['G'-38]['W'-38]=198.7740455;
230  D['G'-38]['Y'-38]=146.9262927;
231  D['W'-38]['*'-38]=stopweight;
232  D['R'-38]['*'-38]=stopweight;
233  D['W'-38]['A'-38]=201.5127248;
234  D['W'-38]['C'-38]=127.4031665;
235  D['W'-38]['D'-38]=175.3689846;
236  D['W'-38]['E'-38]=160.4382217;
237  D['W'-38]['F'-38]=142.6771181;
238  D['W'-38]['G'-38]=198.7740455;
239  D['W'-38]['H'-38]=147.2364297;
240  D['W'-38]['I'-38]=149.826083;
241  D['W'-38]['K'-38]=157.923034;
242  D['M'-38]['*'-38]=stopweight;
243  D['W'-38]['L'-38]=150.0095229;
244  D['W'-38]['M'-38]=152.3965508;
245  D['W'-38]['N'-38]=159.7688325;
246  D['W'-38]['P'-38]=182.7653305;
247  D['W'-38]['Q'-38]=141.5619982;
248  D['R'-38]['A'-38]=111.2645916;
249  D['W'-38]['R'-38]=132.5874819;
250  D['W'-38]['S'-38]=157.855071;
251  D['R'-38]['C'-38]=179.5348377;
252  D['W'-38]['T'-38]=150.752701;
253  D['R'-38]['D'-38]=95.81352768;
254  D['W'-38]['V'-38]=162.8552264;
255  D['R'-38]['E'-38]=54.01884479;
256  D['W'-38]['W'-38]=0.;
257  D['R'-38]['F'-38]=97.03272116;
258  D['R'-38]['G'-38]=125.1291926;
259  D['W'-38]['Y'-38]=128.9732994;
260  D['R'-38]['H'-38]=28.81922224;
261  D['R'-38]['I'-38]=97.58656354;
262  D['H'-38]['*'-38]=stopweight;
263  D['R'-38]['K'-38]=26.00394626;
264  D['R'-38]['L'-38]=101.8802928;
265  D['R'-38]['M'-38]=91.63827667;
266  D['R'-38]['N'-38]=85.11393761;
267  D['R'-38]['P'-38]=102.7148709;
268  D['M'-38]['A'-38]=84.43958222;
269  D['R'-38]['Q'-38]=42.81393193;
270  D['R'-38]['R'-38]=0.;
271  D['R'-38]['S'-38]=109.213153;
272  D['M'-38]['C'-38]=195.5543318;
273  D['M'-38]['D'-38]=160.0267338;
274  D['R'-38]['T'-38]=70.97102635;
275  D['M'-38]['E'-38]=126.0849231;
276  D['R'-38]['V'-38]=95.79707649;
277  D['M'-38]['F'-38]=28.52784356;
278  D['M'-38]['G'-38]=126.9448727;
279  D['R'-38]['W'-38]=132.5874819;
280  D['M'-38]['H'-38]=86.34392904;
281  D['M'-38]['I'-38]=10.12098632;
282  D['R'-38]['Y'-38]=77.10764416;
283  D['C'-38]['*'-38]=stopweight;
284  D['M'-38]['K'-38]=94.4900972;
285  D['M'-38]['L'-38]=14.30316651;
286  D['M'-38]['M'-38]=0.;
287  D['M'-38]['N'-38]=141.1526466;
288  D['M'-38]['P'-38]=86.59479895;
289  D['M'-38]['Q'-38]=100.8991435;
290  D['H'-38]['A'-38]=85.49373612;
291  D['M'-38]['R'-38]=91.63827667;
292  D['H'-38]['C'-38]=173.8426819;
293  D['M'-38]['S'-38]=134.863884;
294  D['M'-38]['T'-38]=81.03969279;
295  D['H'-38]['D'-38]=81.2379796;
296  D['H'-38]['E'-38]=40.79415368;
297  D['M'-38]['V'-38]=21.52179938;
298  D['H'-38]['F'-38]=99.99505513;
299  D['H'-38]['G'-38]=97.53340215;
300  D['M'-38]['W'-38]=152.3965508;
301  D['H'-38]['H'-38]=0.;
302  D['H'-38]['I'-38]=94.33752851;
303  D['M'-38]['Y'-38]=35.22269546;
304  D['H'-38]['K'-38]=32.40288701;
305  D['H'-38]['L'-38]=98.69301471;
306  D['H'-38]['M'-38]=86.34392904;
307  D['H'-38]['N'-38]=68.35485568;
308  D['H'-38]['P'-38]=76.27712277;
309  D['C'-38]['A'-38]=195.0037439;
310  D['H'-38]['Q'-38]=24.08398867;
311  D['H'-38]['R'-38]=28.81922224;
312  D['C'-38]['C'-38]=0.;
313  D['H'-38]['S'-38]=88.93552755;
314  D['C'-38]['D'-38]=153.5743572;
315  D['H'-38]['T'-38]=46.75286971;
316  D['C'-38]['E'-38]=169.4374593;
317  D['C'-38]['F'-38]=204.3884874;
318  D['H'-38]['V'-38]=83.892947;
319  D['C'-38]['G'-38]=158.2317729;
320  D['H'-38]['W'-38]=147.2364297;
321  D['C'-38]['H'-38]=173.8426819;
322  D['C'-38]['I'-38]=197.2497963;
323  D['H'-38]['Y'-38]=83.32844466;
324  D['C'-38]['K'-38]=201.5793361;
325  D['C'-38]['L'-38]=197.4289714;
326  D['C'-38]['M'-38]=195.5543318;
327  D['C'-38]['N'-38]=138.7664683;
328  D['C'-38]['P'-38]=168.5906208;
329  D['C'-38]['Q'-38]=154.2308431;
330  D['C'-38]['R'-38]=179.5348377;
331  D['C'-38]['S'-38]=111.6724957;
332  D['C'-38]['T'-38]=148.9294552;
333  D['C'-38]['V'-38]=191.232501;
334  D['C'-38]['W'-38]=127.4031665;
335  D['S'-38]['*'-38]=stopweight;
336  D['C'-38]['Y'-38]=193.7247005;
337  D['N'-38]['*'-38]=stopweight;
338  D['S'-38]['A'-38]=99.13249708;
339  D['S'-38]['C'-38]=111.6724957;
340  D['S'-38]['D'-38]=65.47081516;
341  D['S'-38]['E'-38]=79.78756921;
342  D['S'-38]['F'-38]=154.8079023;
343  D['S'-38]['G'-38]=55.26736276;
344  D['S'-38]['H'-38]=88.93552755;
345  D['S'-38]['I'-38]=141.7979605;
346  D['I'-38]['*'-38]=stopweight;
347  D['S'-38]['K'-38]=120.5325805;
348  D['*'-38]['*'-38]=0.;
349  D['S'-38]['L'-38]=144.0792346;
350  D['S'-38]['M'-38]=134.863884;
351  D['S'-38]['N'-38]=46.23977508;
352  D['S'-38]['P'-38]=73.35256887;
353  D['N'-38]['A'-38]=110.418026;
354  D['S'-38]['Q'-38]=68.19787424;
355  D['S'-38]['R'-38]=109.213153;
356  D['N'-38]['C'-38]=138.7664683;
357  D['S'-38]['S'-38]=0.;
358  D['N'-38]['D'-38]=23.005355;
359  D['S'-38]['T'-38]=57.74902766;
360  D['N'-38]['E'-38]=40.85895231;
361  D['N'-38]['F'-38]=158.1121886;
362  D['S'-38]['V'-38]=123.0339581;
363  D['N'-38]['G'-38]=79.34585612;
364  D['S'-38]['W'-38]=157.855071;
365  D['N'-38]['H'-38]=68.35485568;
366  D['S'-38]['Y'-38]=143.1073844;
367  D['N'-38]['I'-38]=148.9140957;
368  D['D'-38]['*'-38]=stopweight;
369  D['N'-38]['K'-38]=93.88268262;
370  D['N'-38]['L'-38]=152.3309797;
371  D['N'-38]['M'-38]=141.1526466;
372  D['N'-38]['N'-38]=0.;
373  D['N'-38]['P'-38]=90.15739727;
374  D['N'-38]['Q'-38]=45.75219661;
375  D['I'-38]['A'-38]=93.66243338;
376  D['N'-38]['R'-38]=85.11393761;
377  D['N'-38]['S'-38]=46.23977508;
378  D['*'-38]['A'-38]=stopweight;
379  D['I'-38]['C'-38]=197.2497963;
380  D['N'-38]['T'-38]=64.77428784;
381  D['I'-38]['D'-38]=168.0809621;
382  D['*'-38]['C'-38]=stopweight;
383  D['I'-38]['E'-38]=134.1621451;
384  D['N'-38]['V'-38]=132.8777227;
385  D['*'-38]['D'-38]=stopweight;
386  D['I'-38]['F'-38]=21.27701376;
387  D['N'-38]['W'-38]=159.7688325;
388  D['*'-38]['E'-38]=stopweight;
389  D['I'-38]['G'-38]=135.4188097;
390  D['*'-38]['F'-38]=stopweight;
391  D['I'-38]['H'-38]=94.33752851;
392  D['N'-38]['Y'-38]=142.2293696;
393  D['*'-38]['G'-38]=stopweight;
394  D['I'-38]['I'-38]=0.;
395  D['*'-38]['H'-38]=stopweight;
396  D['*'-38]['I'-38]=stopweight;
397  D['I'-38]['K'-38]=101.6123327;
398  D['I'-38]['L'-38]=4.855121198;
399  D['*'-38]['K'-38]=stopweight;
400  D['I'-38]['M'-38]=10.12098632;
401  D['*'-38]['L'-38]=stopweight;
402  D['I'-38]['N'-38]=148.9140957;
403  D['*'-38]['M'-38]=stopweight;
404  D['*'-38]['N'-38]=stopweight;
405  D['I'-38]['P'-38]=95.3760884;
406  D['D'-38]['A'-38]=125.7485902;
407  D['I'-38]['Q'-38]=108.566243;
408  D['*'-38]['P'-38]=stopweight;
409  D['I'-38]['R'-38]=97.58656354;
410  D['D'-38]['C'-38]=153.5743572;
411  D['*'-38]['Q'-38]=stopweight;
412  D['I'-38]['S'-38]=141.7979605;
413  D['D'-38]['D'-38]=0.;
414  D['*'-38]['R'-38]=stopweight;
415  D['I'-38]['T'-38]=89.2829021;
416  D['D'-38]['E'-38]=44.60463645;
417  D['*'-38]['S'-38]=stopweight;
418  D['D'-38]['F'-38]=176.5261931;
419  D['*'-38]['T'-38]=stopweight;
420  D['I'-38]['V'-38]=29.60907063;
421  D['D'-38]['G'-38]=93.76755083;
422  D['I'-38]['W'-38]=149.826083;
423  D['D'-38]['H'-38]=81.2379796;
424  D['*'-38]['V'-38]=stopweight;
425  D['D'-38]['I'-38]=168.0809621;
426  D['*'-38]['W'-38]=stopweight;
427  D['Y'-38]['*'-38]=stopweight;
428  D['I'-38]['Y'-38]=33.04769862;
429  D['D'-38]['K'-38]=101.4568985;
430  D['*'-38]['Y'-38]=stopweight;
431  D['D'-38]['L'-38]=171.7572022;
432  D['D'-38]['M'-38]=160.0267338;
433  D['D'-38]['N'-38]=23.005355;
434  D['D'-38]['P'-38]=107.9096634;
435  D['D'-38]['Q'-38]=61.28447445;
436  D['D'-38]['R'-38]=95.81352768;
437  D['D'-38]['S'-38]=65.47081516;
438  D['D'-38]['T'-38]=85.07613668;
439  D['D'-38]['V'-38]=152.013397;
440  D['D'-38]['W'-38]=175.3689846;
441  D['D'-38]['Y'-38]=159.9374066;
442  D['T'-38]['*'-38]=stopweight;
443  D['Y'-38]['A'-38]=111.5882137;
444  D['Y'-38]['C'-38]=193.7247005;
445  D['Y'-38]['D'-38]=159.9374066;
446  D['Y'-38]['E'-38]=122.7768651;
447  D['Y'-38]['F'-38]=21.60966489;
448  D['Y'-38]['G'-38]=146.9262927;
449  D['Y'-38]['H'-38]=83.32844466;
450  D['Y'-38]['I'-38]=33.04769862;
451  D['Y'-38]['K'-38]=84.78644482;
452  D['Y'-38]['L'-38]=35.6773121;
453  D['Y'-38]['M'-38]=35.22269546;
454  D['Y'-38]['N'-38]=142.2293696;
455  D['Y'-38]['P'-38]=109.6155723;
456  D['T'-38]['A'-38]=58.0233276;
457  D['Y'-38]['Q'-38]=98.78325151;
458  D['Y'-38]['R'-38]=77.10764416;
459  D['Y'-38]['S'-38]=143.1073844;
460  D['T'-38]['C'-38]=148.9294552;
461  D['Y'-38]['T'-38]=92.24764162;
462  D['T'-38]['D'-38]=85.07613668;
463  D['T'-38]['E'-38]=65.50131254;
464  D['Y'-38]['V'-38]=54.66278181;
465  D['T'-38]['F'-38]=102.8587191;
466  D['Y'-38]['W'-38]=128.9732994;
467  D['T'-38]['G'-38]=59.15645198;
468  D['T'-38]['H'-38]=46.75286971;
469  D['Y'-38]['Y'-38]=0.;
470  D['T'-38]['I'-38]=89.2829021;
471  D['T'-38]['K'-38]=77.74105567;
472  D['T'-38]['L'-38]=92.35427468;
473  }
474 
475  double
476  Grantham::operator() (char aa1, char aa2) const
483  {
484  return D[std::toupper(aa1) - 38][std::toupper(aa2) - 38];
485  }
486 }
Sequence
The namespace in which this library resides.
Grantham.hpp
Grantham&#39;s distances (Sequence::Grantham)
Sequence::Grantham::operator()
double operator()(char aa1, char aa2) const
Definition: Grantham.cc:476