2 September 2021 9:26:32.904 AM bdmlib_test(): FORTRAN90 version Test bdmlib(). MIXTURE_READ - Note: Number of records read was 38 Number of goofy records was 0 Numeric key for amino acid abbreviations. I Amino Acid Symbol 1 A Alanine 2 C Cysteine 3 D Aspartic acid 4 E Glutamic acid 5 F Phenylalanine 6 G Glycine 7 H Histidine 8 I Isoleucine 9 K Lysine 10 L Leucine 11 M Methionine 12 N Asparagine 13 P Proline 14 Q Glutamine 15 R Arginine 16 S Serine 17 T Threonine 18 V Valine 19 W Tryptophan 20 Y Tyrosine Number of components = 9 Compon: 1 2 3 4 5 6 7 8 9 Weight: 0.1781 0.0566 0.0960 0.0781 0.0835 0.0904 0.1145 0.0682 0.2346 1 A 0.2707 0.0215 0.5615 0.0701 0.0411 0.1156 0.0935 0.4522 0.0052 2 C 0.0398 0.0103 0.0454 0.0111 0.0148 0.0374 0.0047 0.1146 0.0040 3 D 0.0176 0.0117 0.4384 0.0195 0.0056 0.0124 0.3873 0.0625 0.0067 4 E 0.0164 0.0109 0.7642 0.0947 0.0102 0.0182 0.3478 0.1157 0.0061 5 F 0.0143 0.3857 0.0874 0.0132 0.1536 0.0518 0.0108 0.2842 0.0035 6 G 0.1319 0.0164 0.2591 0.0480 0.0078 0.0173 0.1059 0.1402 0.0169 7 H 0.0124 0.0762 0.2149 0.0770 0.0072 0.0049 0.0498 0.1004 0.0036 8 I 0.0226 0.0353 0.1459 0.0329 0.2996 0.7969 0.0150 0.5502 0.0022 9 K 0.0204 0.0139 0.7622 0.5766 0.0108 0.0171 0.0943 0.1440 0.0050 10 L 0.0307 0.0935 0.2473 0.0723 0.9994 0.2859 0.0278 0.7006 0.0060 11 M 0.0153 0.0220 0.1187 0.0282 0.2102 0.0758 0.0100 0.2766 0.0015 12 N 0.0483 0.0286 0.4416 0.0804 0.0061 0.0145 0.1879 0.1186 0.0042 13 P 0.0538 0.0131 0.1748 0.0377 0.0130 0.0151 0.0500 0.0975 0.0091 14 Q 0.0207 0.0230 0.5308 0.1850 0.0198 0.0114 0.1100 0.1267 0.0036 15 R 0.0236 0.0189 0.4655 0.5068 0.0145 0.0127 0.0387 0.1436 0.0066 16 S 0.2161 0.0292 0.5834 0.0737 0.0120 0.0275 0.1195 0.2790 0.0032 17 T 0.1472 0.0182 0.4456 0.0716 0.0358 0.0883 0.0658 0.3585 0.0037 18 V 0.0654 0.0361 0.2271 0.0425 0.1801 0.9443 0.0254 0.6617 0.0030 19 W 0.0038 0.0718 0.0295 0.0113 0.0127 0.0044 0.0032 0.0615 0.0028 20 Y 0.0096 0.4196 0.1211 0.0287 0.0265 0.0167 0.0187 0.1994 0.0027 Sum 1.1806 1.3558 6.6644 2.0814 2.0810 2.5682 1.7661 4.9877 0.0995 test01(): Generate a (nonrandom) sequence of 10 isoleucine results in a row. Exact component weights: 1 0.178091 2 0.565910E-01 3 0.960191E-01 4 0.781233E-01 5 0.834977E-01 6 0.904123E-01 7 0.114468 8 0.682132E-01 9 0.234585 Initial ALPHA: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 Initial estimated component weights: 1 0.111111 2 0.111111 3 0.111111 4 0.111111 5 0.111111 6 0.111111 7 0.111111 8 0.111111 9 0.111111 Current ALPHA: 1 1.02823 2 1.03844 3 1.03230 4 1.02334 5 1.21239 6 1.45770 7 1.01250 8 1.16273 9 1.03238 Estimated component weights: 1 0.102823 2 0.103844 3 0.103230 4 0.102334 5 0.121239 6 0.145770 7 0.101250 8 0.116273 9 0.103238 Current ALPHA: 1 1.05081 2 1.06947 3 1.05822 4 1.04192 5 1.41260 6 1.97645 7 1.02234 8 1.30984 9 1.05837 Estimated component weights: 1 0.955279E-01 2 0.972245E-01 3 0.962021E-01 4 0.947197E-01 5 0.128418 6 0.179677 7 0.929396E-01 8 0.119076 9 0.962151E-01 Current ALPHA: 1 1.06942 2 1.09526 3 1.07967 4 1.05717 5 1.60080 6 2.54392 7 1.03035 8 1.44355 9 1.07986 Estimated component weights: 1 0.891182E-01 2 0.912713E-01 3 0.899721E-01 4 0.880979E-01 5 0.133400 6 0.211994 7 0.858626E-01 8 0.120296 9 0.899886E-01 Current ALPHA: 1 1.08515 2 1.11719 3 1.09784 4 1.07003 5 1.77798 6 3.15070 7 1.03706 8 1.56596 9 1.09808 Estimated component weights: 1 0.834734E-01 2 0.859380E-01 3 0.844491E-01 4 0.823104E-01 5 0.136768 6 0.242361 7 0.797740E-01 8 0.120459 9 0.844679E-01 Current ALPHA: 1 1.09873 2 1.13622 3 1.11355 4 1.08110 5 1.94529 6 3.78961 7 1.04280 8 1.67886 9 1.11383 Estimated component weights: 1 0.784806E-01 2 0.811585E-01 3 0.795392E-01 4 0.772215E-01 5 0.138949 6 0.270687 7 0.744860E-01 8 0.119919 9 0.795596E-01 Current ALPHA: 1 1.11063 2 1.15298 3 1.12735 4 1.09079 5 2.10384 6 4.45523 7 1.04781 8 1.78370 9 1.12767 Estimated component weights: 1 0.740422E-01 2 0.768652E-01 3 0.751567E-01 4 0.727191E-01 5 0.140256 6 0.297016 7 0.698537E-01 8 0.118913 9 0.751783E-01 Current ALPHA: 1 1.12122 2 1.16793 3 1.13964 4 1.09938 5 2.25462 6 5.14335 7 1.05222 8 1.88165 9 1.13999 Estimated component weights: 1 0.700760E-01 2 0.729958E-01 3 0.712274E-01 4 0.687112E-01 5 0.140914 6 0.321459 7 0.657640E-01 8 0.117603 9 0.712497E-01 Current ALPHA: 1 1.13073 2 1.18142 3 1.15070 4 1.10709 5 2.39850 6 5.85066 7 1.05618 8 1.97365 9 1.15108 Estimated component weights: 1 0.665134E-01 2 0.694953E-01 3 0.676881E-01 4 0.651229E-01 5 0.141088 6 0.344157 7 0.621280E-01 8 0.116097 9 0.677109E-01 Current ALPHA: 1 1.13936 2 1.19370 3 1.16075 4 1.11408 5 2.53620 6 6.57455 7 1.05974 8 2.06047 9 1.16116 Estimated component weights: 1 0.632976E-01 2 0.663164E-01 3 0.644858E-01 4 0.618931E-01 5 0.140900 6 0.365253 7 0.588746E-01 8 0.114470 9 0.645089E-01 Current ALPHA: 1 1.14725 2 1.20495 3 1.16994 4 1.12046 5 2.66837 6 7.31291 7 1.06299 8 2.14274 9 1.17038 Estimated component weights: 1 0.603816E-01 2 0.634186E-01 3 0.615760E-01 4 0.589714E-01 5 0.140441 6 0.384890 7 0.559470E-01 8 0.112776 9 0.615992E-01 test02(): Generate many random events. We should be able to approximate the exact component weights. Exact component weights: 1 0.178091 2 0.565910E-01 3 0.960191E-01 4 0.781233E-01 5 0.834977E-01 6 0.904123E-01 7 0.114468 8 0.682132E-01 9 0.234585 Initial ALPHA: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 Initial estimated component weights: 1 0.111111 2 0.111111 3 0.111111 4 0.111111 5 0.111111 6 0.111111 7 0.111111 8 0.111111 9 0.111111 Event 1 Current ALPHA: 1 1.03099 2 1.00784 3 1.01637 4 1.69625 5 1.00026 6 1.00040 7 1.16669 8 1.00125 9 1.07995 Estimated component weights: 1 0.103099 2 0.100784 3 0.101637 4 0.169625 5 0.100026 6 0.100040 7 0.116669 8 0.100125 9 0.107995 Event 2 Current ALPHA: 1 1.86007 2 1.00804 3 1.04569 4 1.69927 5 1.00059 6 1.01398 7 1.17162 8 1.12078 9 1.07997 Estimated component weights: 1 0.169097 2 0.916397E-01 3 0.950626E-01 4 0.154479 5 0.909626E-01 6 0.921802E-01 7 0.106511 8 0.101889 9 0.981789E-01 Event 3 Current ALPHA: 1 1.86031 2 1.00829 3 1.09503 4 1.69951 5 1.00390 6 1.02138 7 1.17167 8 2.05994 9 1.07997 Estimated component weights: 1 0.155026 2 0.840244E-01 3 0.912525E-01 4 0.141626 5 0.836579E-01 6 0.851150E-01 7 0.976391E-01 8 0.171662 9 0.899973E-01 Event 4 Current ALPHA: 1 1.89937 2 1.01098 3 1.09634 4 1.70285 5 1.29509 6 1.37897 7 1.17250 8 2.36346 9 1.08045 Estimated component weights: 1 0.146105 2 0.777677E-01 3 0.843337E-01 4 0.130989 5 0.996225E-01 6 0.106074 7 0.901920E-01 8 0.181804 9 0.831116E-01 Event 5 Current ALPHA: 1 2.02511 2 1.01485 3 1.09639 4 1.70705 5 1.29556 6 1.37942 7 1.18678 8 2.36374 9 1.93110 Estimated component weights: 1 0.144651 2 0.724891E-01 3 0.783139E-01 4 0.121932 5 0.925397E-01 6 0.985300E-01 7 0.847703E-01 8 0.168839 9 0.137936 Event 6 Current ALPHA: 1 2.13594 2 1.02752 3 1.16753 4 1.74016 5 1.35644 6 1.38795 7 1.21319 8 3.03988 9 1.93138 Estimated component weights: 1 0.142396 2 0.685013E-01 3 0.778357E-01 4 0.116011 5 0.904295E-01 6 0.925301E-01 7 0.808793E-01 8 0.202659 9 0.128759 Event 7 Current ALPHA: 1 2.19052 2 1.03115 3 1.16757 4 1.74398 5 1.35740 6 1.38841 7 1.22020 8 3.04015 9 2.86062 Estimated component weights: 1 0.136908 2 0.644466E-01 3 0.729733E-01 4 0.108999 5 0.848377E-01 6 0.867759E-01 7 0.762625E-01 8 0.190009 9 0.178789 Event 8 Current ALPHA: 1 2.20803 2 1.05008 3 1.16761 4 1.74869 5 1.77209 6 1.39884 7 1.22235 8 3.04482 9 3.38750 Estimated component weights: 1 0.129884 2 0.617697E-01 3 0.686830E-01 4 0.102864 5 0.104240 6 0.822848E-01 7 0.719028E-01 8 0.179107 9 0.199265 Event 9 Current ALPHA: 1 2.24116 2 1.05226 3 1.16763 4 1.75094 5 1.77282 6 1.39911 7 1.22650 8 3.04497 9 4.34460 Estimated component weights: 1 0.124509 2 0.584590E-01 3 0.648685E-01 4 0.972745E-01 5 0.984902E-01 6 0.777286E-01 7 0.681386E-01 8 0.169165 9 0.241367 Event 10 Current ALPHA: 1 2.28269 2 1.05349 3 1.20746 4 1.76736 5 1.77289 6 1.39918 7 2.11444 8 3.04825 9 4.35423 Estimated component weights: 1 0.120142 2 0.554469E-01 3 0.635507E-01 4 0.930189E-01 5 0.933102E-01 6 0.736411E-01 7 0.111286 8 0.160434 9 0.229170 Event 100 Current ALPHA: 1 21.3749 2 8.95503 3 4.63682 4 12.7703 5 7.86991 6 8.61096 7 11.7578 8 11.3270 9 21.6973 Estimated component weights: 1 0.196100 2 0.821563E-01 3 0.425396E-01 4 0.117159 5 0.722011E-01 6 0.789996E-01 7 0.107869 8 0.103918 9 0.199058 Event 200 Current ALPHA: 1 34.4465 2 15.6518 3 11.3350 4 29.0368 5 12.7849 6 17.5014 7 25.5737 8 18.1236 9 44.5463 Estimated component weights: 1 0.164816 2 0.748889E-01 3 0.542345E-01 4 0.138932 5 0.611719E-01 6 0.837387E-01 7 0.122362 8 0.867157E-01 9 0.213140 Event 300 Current ALPHA: 1 48.3599 2 22.6987 3 19.2354 4 36.2869 5 18.7754 6 25.4487 7 39.3333 8 25.2848 9 73.5768 Estimated component weights: 1 0.156504 2 0.734586E-01 3 0.622503E-01 4 0.117433 5 0.607619E-01 6 0.823584E-01 7 0.127292 8 0.818278E-01 9 0.238113 Event 400 Current ALPHA: 1 64.4545 2 28.2282 3 27.0089 4 45.9565 5 26.6684 6 33.7430 7 51.5825 8 29.8647 9 101.493 Estimated component weights: 1 0.157590 2 0.690175E-01 3 0.660363E-01 4 0.112363 5 0.652039E-01 6 0.825012E-01 7 0.126119 8 0.730189E-01 9 0.248150 Event 500 Current ALPHA: 1 79.0510 2 34.2127 3 34.0886 4 59.0129 5 31.7912 6 40.7571 7 66.2022 8 36.3440 9 127.540 Estimated component weights: 1 0.155307 2 0.672154E-01 3 0.669717E-01 4 0.115939 5 0.624581E-01 6 0.800728E-01 7 0.130063 8 0.714028E-01 9 0.250570 Event 600 Current ALPHA: 1 98.7020 2 38.8284 3 42.5945 4 67.2570 5 38.2531 6 50.1103 7 79.3531 8 43.8279 9 150.074 Estimated component weights: 1 0.162072 2 0.637577E-01 3 0.699418E-01 4 0.110438 5 0.628129E-01 6 0.822830E-01 7 0.130301 8 0.719670E-01 9 0.246426 Event 700 Current ALPHA: 1 112.888 2 43.3981 3 50.4021 4 78.6663 5 45.6683 6 60.8942 7 91.3692 8 54.4425 9 171.271 Estimated component weights: 1 0.159222 2 0.612103E-01 3 0.710891E-01 4 0.110954 5 0.644122E-01 6 0.858875E-01 7 0.128871 8 0.767878E-01 9 0.241567 Event 800 Current ALPHA: 1 126.701 2 51.2223 3 56.8434 4 88.4260 5 50.9693 6 68.9147 7 101.681 8 61.9700 9 202.272 Estimated component weights: 1 0.156614 2 0.633156E-01 3 0.702638E-01 4 0.109303 5 0.630029E-01 6 0.851850E-01 7 0.125688 8 0.766008E-01 9 0.250027 Event 900 Current ALPHA: 1 144.597 2 54.7744 3 67.2465 4 96.0825 5 60.7902 6 79.1509 7 116.146 8 70.0087 9 220.204 Estimated component weights: 1 0.159073 2 0.602579E-01 3 0.739786E-01 4 0.105701 5 0.668759E-01 6 0.870748E-01 7 0.127773 8 0.770173E-01 9 0.242248 Event 1000 Current ALPHA: 1 159.400 2 59.6180 3 75.2611 4 105.247 5 67.1469 6 90.3037 7 128.539 8 78.2125 9 245.272 Estimated component weights: 1 0.157978 2 0.590862E-01 3 0.745898E-01 4 0.104308 5 0.665480E-01 6 0.894982E-01 7 0.127393 8 0.775149E-01 9 0.243084 bdmlib_test(): Normal end of execution. 2 September 2021 9:26:32.937 AM