Path: utzoo!utgpu!cs.utexas.edu!wuarchive!zaphod.mps.ohio-state.edu!rpi!bu.edu!ncsuvx!news From: arg@ccvr1.ncsu.edu (Ron Gallant) Newsgroups: alt.sources Subject: nlmdl Part 4 of 6 Message-ID: <1990Dec22.230111.9649@ncsuvx.ncsu.edu> Date: 22 Dec 90 23:01:11 GMT Reply-To: arg@ccvr1.ncsu.edu (Ron Gallant) Organization: NCSU Computing Center Lines: 1152 #! /bin/sh # This is a shell archive. Remove anything before this line, then unpack # it by saving it into a file and typing "sh file". To overwrite existing # files, type "sh file -c". You can also feed this as standard input via # unshar, or by typing "sh 'ch5eg1/electa.dat' <<'END_OF_FILE' X 1 1 0.056731 0.280382 0.662888 0.46931 X 2 1 0.103444 0.252128 0.644427 0.79539 X 3 1 0.158353 0.270089 0.571558 0.45756 X 4 1 0.108075 0.305072 0.586853 0.94713 X 5 1 0.083921 0.211656 0.704423 1.22054 X 6 1 0.112165 0.290532 0.597302 0.93181 X 7 1 0.071274 0.240518 0.688208 1.79152 X 8 1 0.076510 0.210503 0.712987 0.51442 X 9 1 0.066173 0.202999 0.730828 0.78407 X 10 1 0.094836 0.270281 0.634883 1.01354 X 11 1 0.078501 0.293953 0.627546 0.83854 X 12 1 0.059530 0.228752 0.711718 1.53957 X 13 1 0.208982 0.328053 0.462965 1.06694 X 14 1 0.083702 0.297272 0.619027 0.82437 X 15 1 0.138705 0.358329 0.502966 0.80712 X 16 1 0.111378 0.322564 0.566058 0.53169 X 17 1 0.092919 0.259633 0.647448 0.85439 X 18 1 0.039353 0.158205 0.802442 1.93326 X 19 1 0.066577 0.247454 0.685970 1.37160 X 20 2 0.102844 0.244335 0.652821 0.92766 X 21 2 0.125485 0.230305 0.644210 1.80934 X 22 2 0.154316 0.235135 0.610549 2.41501 X 23 2 0.165714 0.276980 0.557305 0.84658 X 24 2 0.145370 0.173112 0.681518 1.60788 X 25 2 0.184467 0.268865 0.546668 0.73838 X 26 2 0.162269 0.280939 0.556792 0.81116 X 27 2 0.112016 0.220850 0.667133 2.01503 X 28 2 0.226863 0.257833 0.515304 2.32035 X 29 2 0.118028 0.219830 0.662142 2.40172 X 30 2 0.137761 0.345117 0.517122 0.57141 X 31 2 0.079115 0.257319 0.663566 0.94474 X 32 2 0.185022 0.265051 0.549928 1.63778 X 33 2 0.144524 0.276133 0.579343 0.75816 X 34 2 0.201734 0.241966 0.556300 1.00136 X 35 2 0.094890 0.227651 0.677459 1.11384 X 36 2 0.102843 0.264515 0.632642 1.07185 X 37 2 0.107760 0.214232 0.678009 1.53659 X 38 2 0.156552 0.236422 0.607026 0.24099 X 39 2 0.088431 0.222746 0.688822 0.58066 X 40 2 0.146236 0.301884 0.551880 2.52983 X 41 3 0.080802 0.199005 0.720192 1.14741 X 42 3 0.100711 0.387758 0.511531 0.97934 X 43 3 0.073483 0.335280 0.591237 1.09361 X 44 3 0.059455 0.259823 0.680722 2.19468 X 45 3 0.076195 0.378371 0.545434 1.98221 X 46 3 0.076926 0.325032 0.598042 1.78194 X 47 3 0.086052 0.339653 0.574295 3.24274 X 48 3 0.069359 0.278369 0.652272 0.47593 X 49 3 0.071265 0.273866 0.654869 1.38369 X 50 3 0.100562 0.306247 0.593191 1.57831 X 51 3 0.050203 0.294285 0.655513 2.16900 X 52 3 0.059627 0.311932 0.628442 2.11575 X 53 3 0.081433 0.328604 0.589962 0.35681 X 54 3 0.075762 0.285972 0.638265 1.55275 X 55 3 0.042910 0.372337 0.584754 1.06305 X 56 3 0.086846 0.340184 0.572970 4.02013 X 57 3 0.102537 0.335535 0.561928 0.60712 X 58 3 0.068766 0.310782 0.620452 1.15334 X 59 3 0.058405 0.307111 0.634485 2.43797 X 60 4 0.055227 0.300839 0.643934 0.10082 X 61 4 0.107435 0.273937 0.618628 0.69302 X 62 4 0.105958 0.291205 0.602837 1.12592 X 63 4 0.132278 0.279429 0.588293 1.84425 X 64 4 0.094195 0.328866 0.576940 1.57972 X 65 4 0.115259 0.401079 0.483663 1.27034 X 66 4 0.150229 0.317866 0.531905 0.56330 X 67 4 0.168780 0.307669 0.523551 3.43139 X 68 4 0.118222 0.318080 0.563698 1.00979 X 69 4 0.103394 0.307671 0.588936 2.08458 X 70 4 0.124007 0.362115 0.513879 1.30410 X 71 4 0.197987 0.280130 0.521884 3.48146 X 72 4 0.108083 0.337004 0.554913 0.53206 X 73 5 0.088798 0.232568 0.678634 3.28987 X 74 5 0.100508 0.272139 0.627353 0.32678 X 75 5 0.127303 0.298519 0.574178 0.52452 X 76 5 0.109718 0.228172 0.662109 0.36622 X 77 5 0.130080 0.231037 0.638883 0.63788 X 78 5 0.148562 0.323579 0.527859 1.42239 X 79 5 0.106306 0.252137 0.641556 0.93535 X 80 5 0.080877 0.214172 0.704951 1.26243 X 81 5 0.081810 0.135665 0.782525 1.51472 X 82 5 0.131749 0.278338 0.589913 2.07858 X 83 5 0.059180 0.254533 0.686287 1.60681 X 84 5 0.078620 0.267252 0.654128 1.54706 X 85 5 0.090220 0.293831 0.615949 2.61162 X 86 5 0.086916 0.193967 0.719117 2.96418 X 87 5 0.132383 0.230489 0.637127 0.26912 X 88 5 0.085560 0.252321 0.662120 0.42554 X 89 5 0.071368 0.276238 0.652393 1.01926 X 90 5 0.061196 0.245025 0.693780 1.53807 X 91 5 0.086608 0.233981 0.679411 0.75711 X 92 5 0.105628 0.305471 0.588901 0.83647 X 93 5 0.078158 0.202536 0.719307 1.92096 X 94 5 0.048632 0.216807 0.734560 1.57795 X 95 5 0.094527 0.224344 0.681128 0.83216 X 96 5 0.092809 0.209154 0.698037 1.39364 X 97 5 0.035751 0.166231 0.798018 1.72697 X 98 5 0.065205 0.205058 0.729736 2.04120 X 99 5 0.092561 0.193848 0.713591 2.04708 X 100 5 0.063119 0.234114 0.702767 3.43969 X 101 5 0.091186 0.224488 0.684326 2.66918 X 102 5 0.047291 0.262623 0.690086 2.71072 X 103 5 0.081575 0.206400 0.712025 3.36803 X 104 5 0.108165 0.243650 0.648185 0.65682 X 105 5 0.079534 0.320450 0.600017 0.95523 X 106 5 0.084828 0.247189 0.667984 0.61441 X 107 5 0.063747 0.210343 0.725910 1.85034 X 108 5 0.081108 0.249960 0.668932 2.11274 X 109 5 0.089942 0.206601 0.703457 1.54120 X 110 5 0.046717 0.224784 0.728499 3.54351 X 111 5 0.114925 0.272279 0.612796 2.61769 X 112 5 0.115055 0.264415 0.620530 3.00236 X 113 5 0.081511 0.223870 0.694618 1.74166 X 114 5 0.109658 0.343593 0.546750 1.17640 X 115 5 0.114263 0.304761 0.580976 0.74566 X 116 5 0.115089 0.226412 0.658499 1.30392 X 117 5 0.040622 0.198986 0.760392 2.13339 X 118 5 0.073245 0.238522 0.688234 2.83039 X 119 5 0.087954 0.287450 0.624596 1.62179 X 120 5 0.091967 0.206131 0.701902 2.18534 X 121 5 0.142746 0.302939 0.554315 0.26503 X 122 5 0.117972 0.253811 0.628217 0.05082 X 123 5 0.071573 0.248324 0.680103 0.42740 X 124 5 0.073628 0.290586 0.635786 0.47979 X 125 5 0.121075 0.350781 0.528145 0.59551 X 126 5 0.077335 0.339358 0.583307 0.47506 X 127 5 0.074766 0.167202 0.758032 2.11867 X 128 5 0.208580 0.331363 0.460058 1.13621 X 129 5 0.080195 0.210619 0.709185 2.61204 X 130 5 0.066156 0.204118 0.729726 1.45227 X 131 5 0.112282 0.252638 0.635080 0.79071 X 132 5 0.041310 0.093106 0.865584 1.30697 X 133 5 0.102675 0.297009 0.600316 0.93691 X 134 5 0.102902 0.270832 0.626266 0.98718 X 135 5 0.118932 0.250104 0.630964 1.40085 X 136 5 0.139760 0.322394 0.537846 1.78710 X 137 5 0.121616 0.214626 0.663758 8.46237 X 138 5 0.065701 0.263818 0.670481 1.58663 X 139 5 0.034029 0.175181 0.790790 2.62535 X 140 5 0.074476 0.194744 0.730780 4.29430 X 141 5 0.059568 0.229705 0.710727 0.65404 X 142 5 0.088128 0.295546 0.616326 0.41292 X 143 5 0.075522 0.213622 0.710856 2.02370 X 144 5 0.057089 0.195720 0.747190 1.76998 X 145 5 0.096331 0.301692 0.601977 0.99891 X 146 5 0.120824 0.250280 0.628896 0.27942 X 147 6 0.034529 0.193456 0.772015 0.91673 X 148 6 0.026971 0.180848 0.792181 1.15617 X 149 6 0.045271 0.141894 0.812835 1.57107 X 150 6 0.067708 0.219302 0.712990 1.24515 X 151 6 0.079335 0.230693 0.689972 1.70748 X 152 6 0.022703 0.178896 0.798401 1.79959 X 153 6 0.043053 0.157142 0.799805 4.61665 X 154 6 0.057157 0.245931 0.696912 0.59504 X 155 6 0.063229 0.136192 0.800579 1.42499 X 156 6 0.076873 0.214209 0.708918 1.34371 X 157 6 0.027353 0.124894 0.847753 2.74908 X 158 6 0.067823 0.146994 0.785183 1.84628 X 159 6 0.056388 0.189185 0.754428 3.82472 X 160 6 0.036841 0.194994 0.768165 1.18199 X 161 6 0.059160 0.138681 0.802158 2.07338 X 162 6 0.051980 0.215700 0.732320 0.80376 X 163 6 0.027300 0.145072 0.827628 1.52316 X 164 6 0.014790 0.179619 0.805591 3.17526 X 165 6 0.047865 0.167561 0.784574 3.30794 X 166 6 0.115629 0.231381 0.652990 0.72456 X 167 7 0.104970 0.147525 0.747505 0.50274 X 168 7 0.119254 0.187409 0.693337 1.22571 X 169 7 0.042564 0.112839 0.844596 2.13534 X 170 7 0.096756 0.150178 0.753066 5.56011 X 171 7 0.063013 0.168422 0.768565 3.11725 X 172 7 0.080060 0.143934 0.776006 0.99796 X 173 7 0.097493 0.173391 0.729116 0.67859 X 174 7 0.102526 0.220954 0.676520 0.79027 X 175 7 0.085538 0.195686 0.718776 2.24498 X 176 7 0.068733 0.166248 0.765019 2.01993 X 177 7 0.094915 0.140119 0.764966 4.07330 X 178 7 0.076163 0.132046 0.791792 3.66432 X 179 7 0.099943 0.176885 0.723172 0.40768 X 180 7 0.081494 0.175082 0.743425 1.09065 X 181 7 0.196026 0.299348 0.504626 1.35008 X 182 7 0.093173 0.235816 0.671011 1.06138 X 183 7 0.172293 0.173032 0.654675 0.99219 X 184 7 0.067736 0.159600 0.772663 3.69199 X 185 7 0.102033 0.171697 0.726271 2.36676 X 186 7 0.067977 0.151109 0.780914 1.84563 X 187 8 0.071073 0.238985 0.689942 0.18316 X 188 8 0.049453 0.286788 0.663759 2.23986 X 189 8 0.062748 0.255129 0.682123 3.48084 X 190 8 0.032376 0.154905 0.812719 7.26135 X 191 8 0.055055 0.225296 0.719648 1.68814 X 192 8 0.037829 0.179051 0.783120 1.13804 X 193 8 0.020102 0.172396 0.807502 1.40894 X 194 8 0.021917 0.149092 0.828992 3.47472 X 195 8 0.047590 0.174735 0.777675 3.37689 X 196 8 0.063446 0.235823 0.700731 3.14810 X 197 8 0.034719 0.159398 0.805883 3.21710 X 198 8 0.055428 0.200488 0.744084 1.13941 X 199 8 0.058074 0.254823 0.687103 2.55414 X 200 8 0.060719 0.209763 0.729518 0.29071 X 201 8 0.045681 0.206177 0.748142 1.21336 X 202 8 0.040151 0.263161 0.696688 1.02370 X 203 8 0.072230 0.281460 0.646310 1.40580 X 204 8 0.064366 0.269816 0.665819 0.97704 X 205 8 0.035993 0.191422 0.772585 2.09909 X 206 9 0.091638 0.215290 0.693073 1.03679 X 207 9 0.072171 0.236658 0.691171 2.36788 X 208 9 0.056187 0.195345 0.748468 3.45908 X 209 9 0.095888 0.229586 0.674526 3.63796 X 210 9 0.069809 0.219558 0.710633 2.56887 X 211 9 0.142920 0.223801 0.633279 2.00319 X 212 9 0.087323 0.196401 0.716276 2.40644 X 213 9 0.064517 0.218711 0.716772 2.58552 X 214 9 0.086882 0.194778 0.718341 8.94023 X 215 9 0.067463 0.219228 0.713309 3.75275 X 216 9 0.105610 0.230661 0.663730 0.34082 X 217 9 0.138992 0.283123 0.577885 1.62649 X 218 9 0.081364 0.186967 0.731670 2.31678 X 219 9 0.114535 0.221751 0.663714 1.77709 X 220 9 0.069940 0.280622 0.649438 1.38765 X 221 9 0.073137 0.143219 0.783643 3.46442 X 222 9 0.096326 0.243241 0.660434 1.74696 X 223 9 0.083284 0.202951 0.713765 1.28613 X 224 9 0.179133 0.299403 0.521465 1.15897 END_OF_FILE if test 16800 -ne `wc -c <'ch5eg1/electa.dat'`; then echo shar: \"'ch5eg1/electa.dat'\" unpacked with wrong size! fi # end of 'ch5eg1/electa.dat' fi if test -f 'ch5eg1/electc2.dat' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'ch5eg1/electc2.dat'\" else echo shar: Extracting \"'ch5eg1/electc2.dat'\" \(16800 characters\) sed "s/^X//" >'ch5eg1/electc2.dat' <<'END_OF_FILE' X 1 0 0 1 1 0 0.700 X 2 1 0 0 1 1.320 0.700 X 3 1 0 0 1 1.320 0.700 X 4 1 0 0 1 1.320 2.495 X 5 1 0 0 0 1.320 3.590 X 6 1 0 0 1 0 1.795 X 7 1 0 0 1 0 1.795 X 8 1 0 0 0 1.320 0.700 X 9 1 0 0 1 1.320 0.700 X 10 1 0 0 0 1.985 1.795 X 11 1 0 0 1 2.640 0.700 X 12 0 0 1 1 1.985 1.795 X 13 1 0 0 1 1.320 1.795 X 14 1 0 0 1 1.320 1.795 X 15 1 0 0 1 1.320 1.795 X 16 0 0 1 1 0 1.795 X 17 1 0 0 1 1.985 1.795 X 18 0 0 1 1 0 0.700 X 19 1 0 0 1 1.320 1.795 X 20 1 0 0 1 0 1.795 X 21 1 0 0 1 1.320 0.700 X 22 1 0 0 1 1.320 0.700 X 23 1 0 0 1 3.305 1.795 X 24 0 1 0 1 0 0.700 X 25 0 0 1 1 1.320 0.700 X 26 1 0 0 1 1.985 0.700 X 27 1 0 0 1 1.320 1.795 X 28 1 0 0 1 1.320 1.795 X 29 1 0 0 1 0 1.795 X 30 0 0 1 1 0 0.700 X 31 1 0 0 1 1.320 0.700 X 32 1 0 0 1 0 1.795 X 33 1 0 0 1 1.320 3.590 X 34 1 0 0 1 0 1.795 X 35 1 0 0 1 1.320 1.795 X 36 1 0 0 1 1.320 0.700 X 37 1 0 0 1 1.320 1.795 X 38 1 0 0 1 0 0.700 X 39 1 0 0 1 0 1.795 X 40 1 0 0 1 1.320 1.400 X 41 1 0 0 1 1.985 1.795 X 42 1 0 0 1 2.640 0.700 X 43 1 0 0 1 1.320 1.795 X 44 1 0 0 1 3.970 1.795 X 45 1 0 0 1 1.985 1.795 X 46 1 0 0 0 0 1.795 X 47 1 0 0 1 0 1.795 X 48 1 0 0 1 0 0.700 X 49 1 0 0 1 0 2.495 X 50 1 0 0 1 1.985 1.795 X 51 0 0 1 1 1.985 0.700 X 52 1 0 0 0 1.320 1.795 X 53 0 0 1 1 0 0.700 X 54 1 0 0 1 1.320 1.795 X 55 1 0 0 1 1.320 1.795 X 56 1 0 0 1 1.320 1.795 X 57 1 0 0 1 1.320 1.795 X 58 1 0 0 1 1.320 1.795 X 59 1 0 0 1 1.320 1.795 X 60 1 0 0 0 0 0.700 X 61 1 0 0 1 1.320 0.700 X 62 1 0 0 1 1.320 0.700 X 63 1 0 0 1 1.320 0.700 X 64 0 1 0 1 0 1.795 X 65 0 0 1 1 0 0.700 X 66 0 0 1 1 0 1.795 X 67 1 0 0 1 1.320 1.795 X 68 1 0 0 1 1.320 0.700 X 69 1 0 0 1 0 1.400 X 70 1 0 0 1 1.320 1.795 X 71 1 0 0 1 1.985 1.795 X 72 1 0 0 0 1.320 1.795 X 73 1 0 0 0 0 1.795 X 74 1 0 0 0 0 0.700 X 75 1 0 0 0 1.320 0.700 X 76 1 0 0 1 0 0.700 X 77 0 1 0 0 1.320 0.700 X 78 1 0 0 1 0 1.795 X 79 1 0 0 0 1.320 1.795 X 80 1 0 0 1 1.320 1.795 X 81 0 0 1 1 1.320 2.495 X 82 1 0 0 1 1.985 1.795 X 83 1 0 0 1 1.320 0.700 X 84 1 0 0 1 1.985 1.795 X 85 1 0 0 1 0 1.795 X 86 1 0 0 1 1.320 2.495 X 87 1 0 0 0 0 1.795 X 88 1 0 0 0 3.305 0.700 X 89 1 0 0 1 1.985 0.700 X 90 1 0 0 1 1.985 1.795 X 91 1 0 0 1 0 1.795 X 92 1 0 0 1 1.985 1.795 X 93 1 0 0 1 1.320 1.795 X 94 1 0 0 1 0 1.795 X 95 0 0 1 1 0 1.795 X 96 1 0 0 1 1.985 0.700 X 97 1 0 0 0 1.320 0.700 X 98 1 0 0 1 1.320 0.700 X 99 1 0 0 1 2.640 1.795 X 100 1 0 0 1 1.320 1.795 X 101 1 0 0 1 1.320 1.795 X 102 1 0 0 1 0 2.495 X 103 1 0 0 1 1.320 1.795 X 104 1 0 0 1 1.320 0.700 X 105 1 0 0 1 1.320 0.700 X 106 1 0 0 1 0 1.795 X 107 1 0 0 1 1.320 0.700 X 108 1 0 0 1 1.320 0.700 X 109 1 0 0 1 0 1.795 X 110 1 0 0 1 1.320 0.700 X 111 1 0 0 1 1.320 1.795 X 112 1 0 0 1 3.970 1.795 X 113 0 0 1 1 0 1.795 X 114 1 0 0 1 1.320 1.795 X 115 0 0 1 1 0 0.700 X 116 1 0 0 1 0 1.795 X 117 1 0 0 1 1.320 1.795 X 118 1 0 0 1 1.985 0.700 X 119 1 0 0 1 0 0.700 X 120 0 0 1 1 1.320 1.795 X 121 0 0 1 1 0 1.795 X 122 1 0 0 1 0 0.700 X 123 1 0 0 1 0 0.700 X 124 1 0 0 1 0 1.795 X 125 1 0 0 0 1.320 1.795 X 126 1 0 0 1 0 0.700 X 127 1 0 0 0 1.320 1.795 X 128 1 0 0 1 0 2.495 X 129 0 0 1 1 1.320 1.795 X 130 1 0 0 1 0 1.795 X 131 1 0 0 1 0 0.700 X 132 1 0 0 0 1.320 1.795 X 133 1 0 0 1 1.320 0.700 X 134 1 0 0 1 1.320 0.700 X 135 0 0 1 1 0 0.700 X 136 1 0 0 1 1.985 1.795 X 137 1 0 0 1 1.985 1.400 X 138 0 1 0 1 0 1.795 X 139 0 0 1 1 0 0.700 X 140 0 1 0 1 7.265 1.795 X 141 1 0 0 1 1.320 0.700 X 142 0 1 0 1 0 1.795 X 143 0 1 0 1 0 0.700 X 144 0 0 1 1 0 0.700 X 145 0 0 1 1 0 1.795 X 146 1 0 0 0 0 1.795 X 147 1 0 0 0 0 0.700 X 148 1 0 0 0 0 1.795 X 149 1 0 0 0 1.320 1.795 X 150 0 1 0 1 0 1.795 X 151 1 0 0 1 2.640 1.795 X 152 1 0 0 1 0 1.795 X 153 1 0 0 0 1.320 1.795 X 154 1 0 0 1 1.320 0.700 X 155 0 0 1 1 0 0.700 X 156 1 0 0 1 1.320 0.700 X 157 0 0 1 1 1.320 1.795 X 158 1 0 0 1 3.970 2.495 X 159 1 0 0 1 0 0.700 X 160 1 0 0 1 1.320 1.795 X 161 1 0 0 1 1.320 1.795 X 162 1 0 0 1 1.320 0.700 X 163 1 0 0 1 0 1.795 X 164 1 0 0 1 0 1.795 X 165 1 0 0 1 2.640 3.590 X 166 1 0 0 1 1.320 0.700 X 167 0 1 0 1 0 0.700 X 168 1 0 0 1 2.640 0.700 X 169 0 0 1 1 0 1.795 X 170 -9999 -9999 -9999 -9999 -9999 -9999 X 171 1 0 0 1 1.320 1.795 X 172 -9999 -9999 -9999 -9999 -9999 -9999 X 173 0 0 1 1 1.320 0.700 X 174 1 0 0 1 1.320 0.700 X 175 1 0 0 1 1.985 1.795 X 176 0 0 1 1 0 0.700 X 177 1 0 0 1 3.970 1.795 X 178 1 0 0 1 1.985 1.795 X 179 0 0 1 0 0 0.700 X 180 1 0 0 1 1.320 1.795 X 181 1 0 0 1 1.985 0.700 X 182 1 0 0 1 1.985 1.795 X 183 1 0 0 1 1.320 1.795 X 184 1 0 0 1 1.320 1.795 X 185 1 0 0 1 0 2.495 X 186 1 0 0 1 0 1.795 X 187 1 0 0 0 1.320 0 X 188 1 0 0 1 1.320 1.795 X 189 1 0 0 1 0 1.795 X 190 1 0 0 0 1.985 2.495 X 191 1 0 0 1 0 1.795 X 192 1 0 0 1 0 1.795 X 193 0 0 1 0 0 0 X 194 0 0 1 1 1.320 1.795 X 195 1 0 0 1 1.320 0.700 X 196 1 0 0 1 1.985 1.795 X 197 1 0 0 1 1.985 1.795 X 198 0 0 1 1 0 1.795 X 199 1 0 0 1 1.320 1.795 X 200 1 0 0 0 0 0.700 X 201 0 1 0 1 1.985 0 X 202 0 1 0 1 0 1.795 X 203 1 0 0 1 1.320 0.700 X 204 1 0 0 1 0 0.700 X 205 0 0 1 1 1.320 0.700 X 206 1 0 0 0 1.320 0.700 X 207 0 0 1 1 0 0.700 X 208 1 0 0 1 1.320 2.495 X 209 1 0 0 1 1.320 1.795 X 210 1 0 0 1 0 0.700 X 211 0 1 0 1 0 0.700 X 212 1 0 0 1 0 1.795 X 213 1 0 0 1 1.320 1.795 X 214 1 0 0 1 1.985 1.795 X 215 1 0 0 1 1.320 1.795 X 216 1 0 0 0 0 0.700 X 217 1 0 0 1 1.320 1.795 X 218 1 0 0 1 1.985 1.795 X 219 1 0 0 0 0 0.700 X 220 1 0 0 1 1.320 1.795 X 221 1 0 0 1 1.320 1.795 X 222 0 1 0 1 0 1.795 X 223 1 0 0 1 1.320 1.795 X 224 1 0 0 1 0 1.795 END_OF_FILE if test 16800 -ne `wc -c <'ch5eg1/electc2.dat'`; then echo shar: \"'ch5eg1/electc2.dat'\" unpacked with wrong size! fi # end of 'ch5eg1/electc2.dat' fi if test -f 'status.cc' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'status.cc'\" else echo shar: Extracting \"'status.cc'\" \(16929 characters\) sed "s/^X//" >'status.cc' <<'END_OF_FILE' X/* ---------------------------------------------------------------------------- X Xnlmdl: status.cc X Xnlmdl is a C++ implementation of the statistical methods in A. Ronald XGallant, "Nonlinear Statistical Models," New York: John Wiley and Sons, X1987, ISBN 0-471-80260-3, using a matrix class realmat that is distributed Xwith it. The header files nlmdl.h and realmat.h describe the use of the Xprogram and matrix class, respectively. X XCopyright (C) 1990. X XA. Ronald Gallant XP.O. Box 5513 XRaleigh NC 27650-5513 XUSA X XPermission to use, copy, modify, and distribute this software and its Xdocumentation for any purpose and without fee is hereby granted, provided Xthat the above copyright notice appear in all copies and that both that Xcopyright notice and this permission notice appear in supporting Xdocumentation. X XThis software is provided "as is" without any expressed or implied warranty. X X---------------------------------------------------------------------------- */ X/* status is a class used by nlmdl */ X X#include "status.h" X Xstatus::status() X{ X switches[0]='\0'; X method[0]='\0'; X n=0; X M=0; X K=0; X p=0; X itheta=0; X ivar=0; X vartype[0]='\0'; X MA=0; X weights[0]='\0'; X tol=0; X eps=0; X detail[0]='\0'; X rank=0; X df[0]='\0'; X obj=0; X starting=STARTING_FILENAME; X ending=ENDING_FILENAME; X} X Xstatus::~status() { } X X Xint status::from(char* filename) X{ X X#ifdef GNU_GPP_COMPILER X X#ifdef USE_ATT_STYLE_IO_WITH_GNU X filebuf ib; X if( ib.open(filename, input) == 0 ){ X cerr << "Error, status::from, Cannot open " << filename << "\n"; X exit(1); X } X istream ifrom(&ib); X#endif X X#ifdef USE_GNU_STYLE_IO_WITH_GNU X istream ifrom(filename,io_readonly,a_useonly); X if ( !ifrom ) { X cerr << "Error, status::from, Cannot open " << filename << "\n"; X exit(1); X } X#endif X X#endif X X#ifdef TURBO_CPP_COMPILER X ifstream ifrom(filename); X if ( !ifrom ) { X cerr << "Error, status::from, Cannot open " << filename << "\n"; X exit(1); X } X#endif X X char junk[MAX_STATUS_LINE]; X char newline; X INTEGER i,len; X double x = 0; X X if (!ifrom.get(switches,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (switches) in " << filename << "\n"; X exit(1); X } X X if (!(ifrom>>method)) return 0; X if ( strcmp(method,"SUR") != 0 X && strcmp(method,"TSLS") != 0 X && strcmp(method,"GMM") != 0 ) { X cerr << "Error, status::from, Bad data (method) in " << filename << "\n"; X exit(1); X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>n)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>M)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>K)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>p)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>itheta)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>ivar)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>vartype)) return 0; X if ( strcmp(vartype,"homoskedastic") != 0 X && strcmp(vartype,"heteroskedastic") != 0 ) { X cerr << "Error, status::from, Bad data (vartype) in " << filename << "\n"; X exit(1); X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>MA)) return 0; X if ( MA<0 || MA>=n ) { X cerr << "Error, status::from, Bad data (MA) in " << filename << "\n"; X exit(1); X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>weights)) return 0; X if ( strcmp(weights,"none") != 0 && strcmp(weights,"Parzen") != 0 ) { X cerr << "Error, status::from, Bad data (weights) in " << filename << "\n"; X exit(1); X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>x)) X return 0; X else X tol = (REAL)x; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>x)) X return 0; X else X eps = (REAL)x; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!(ifrom>>detail)) return 0; X if ( strcmp(detail,"none") != 0 X && strcmp(detail,"minimal") != 0 X && strcmp(detail,"full") != 0 ) { X cerr << "Error, status::from, Bad data (detail) in " << filename << "\n"; X exit(1); X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X X if (!ifrom.get(newline)) return 0; X X for (i=1; i<=5; i++) { X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X if (!ifrom.get(newline)) return 0; X } X X#ifdef GNU_GPP_COMPILER X if (!(ifrom>>x)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (theta) in " << filename << "\n"; X exit(1); X } X#endif X X#ifdef TURBO_CPP_COMPILER X float y; //This nonsense is because Turbo can't read a null line into X x=0; //a float. Precision is lost. I should find a better way. X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X if (sscanf(junk, "%e", &x) <= 0) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (theta) in " << filename << "\n"; X exit(1); X } X x = (double)y; X#endif X X theta.resize(p,(INTEGER)1,0); X theta[1] = (REAL)x; X X if (p==1) return 0; X for (i=2; i<=p; i++) { X if(!(ifrom>>x)) { X cerr << "Error, status::from, Bad data (theta) in " << filename << "\n"; X exit(1); X } X else { X theta[i] = (REAL)x; X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (theta) in " << filename << "\n"; X exit(1); X } X } X X if (!ifrom.get(newline)) return 0; X X#ifdef GNU_GPP_COMPILER X if (!(ifrom>>x)) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X#endif X X#ifdef TURBO_CPP_COMPILER X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) return 0; X if (sscanf(junk, "%e", &x) <= 0) return 0; X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X x = (double)y; X#endif X X if ( strcmp(method,"SUR") == 0 || strcmp(method,"TSLS") == 0) { X var.resize(M,M); X var[1] = (REAL)x; X if (M==1) return 0; X for (i=2; i<=M*M; i++) { X if(!(ifrom>>x)) { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X else { X var[i] = (REAL)x; X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X } X } X else if ( strcmp(method,"GMM") == 0 ) { X len = M*K; X var.resize(len,len); X var[1] = (REAL)x; X if (len==1) return 0; X for (i=2; i<=len*len; i++) { X if(!(ifrom>>x)) { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X else { X var[i] = (REAL)x; X } X if (!ifrom.get(junk,MAX_STATUS_LINE,'\n')) { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X } X } X else { X cerr << "Error, status::from, Bad data (var) in " << filename << "\n"; X exit(1); X } X X ifrom.close(); X return 0; X X} X Xint status::to(char* filename) X{ X X#ifdef GNU_GPP_COMPILER X X#ifdef USE_ATT_STYLE_IO_WITH_GNU X filebuf ob; X if( ob.open(filename, output) == 0 ) { X cerr << "Error, status::to, Cannot open " << filename << "\n"; X exit(1); X } X ostream oto(&ob); X#endif X X#ifdef USE_GNU_STYLE_IO_WITH_GNU X ostream oto(filename,io_writeonly,a_create); X if ( !oto ) { X cerr << "Error, status::to, Cannot open " << filename << "\n"; X exit(1); X } X#endif X X#endif X X#ifdef TURBO_CPP_COMPILER X ofstream oto(filename); X if ( !oto ) { X cerr << "Error, status::to, Cannot open " << filename << "\n"; X exit(1); X } X#endif X X INTEGER i,j,len; X char temp[MAX_STATUS_LINE]; X X if (strlen(switches) == 0) { X if (!(oto << "\tThis line is passed to class model as a string.\n")) { X cerr << "Error, status::to, Error writing to " << filename << "\n"; X exit(1); X } X } X else { X if (!(oto << switches << "\n")) { X cerr << "Error, status::to, Error writing to " << filename << "\n"; X exit(1); X } X } X X sprintf(temp,"%-19s", method); X if (!(oto<0 is unwise.\n")) { X cerr << "Error, status::to, Error writing to " << filename << "\n"; X exit(1); X } X X sprintf(temp,"%-19.6e",tol); X if (!(oto<0 is unwise.\n"; X X sprintf(temp,"%-19.6e",tol); X cout< 0 && cols(var) > 0 ){ X if ( strcmp(method,"SUR") == 0 || strcmp(method,"TSLS") == 0){ X for (j=1; j<=M; j++){ X for (i=1; i<=M; i++){ X sprintf(temp,"%25.16e % 20.8f ",var.check2(i,j),var.elem(i,j)); X cout< 0 && cols(theta) > 0 ){ X for (i=1; i<=p; i++){ X sprintf(temp,"%25.16e % 20.8f ",theta[i],theta[i]); X cout< 0 && cols(D) > 0 ){ X for (i=1; i<=p; i++){ X sprintf(temp,"%25.16e % 20.8f ",D[i],D[i]); X cout<