[Rd] Bug in stars function (PR#738)

tdye@lava.net tdye@lava.net
Sun, 19 Nov 2000 08:57:14 +0100 (MET)


  This message is in MIME format.  The first part should be readable text,
  while the remaining parts are likely unreadable without MIME-aware tools.
  Send mail to mime@docserver.cac.washington.edu for more info.

--1465796871-645734172-974620647=:18687
Content-Type: TEXT/PLAIN; charset=US-ASCII

Hi all,

I've discovered and (hopefully) fixed a bug in the stars function I wrote
some time ago.  A corrected function is attached.

Please let me know if you have any questions or comments.

Tom

Thomas S. Dye, Ph.D.                                    http://www.tsdye.com
Home: 812A 19th Avenue, Honolulu, Hawaii 96816.  (808) 739-1367 or 387-9352.
Work: International Archaeological Research Institute, Inc., 2081 Young St.,
      Honolulu, Hawaii 96826. Voice (808) 946-2548; Fax 943-0716.

--1465796871-645734172-974620647=:18687
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="new-stars.r"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.BSI.4.10.10011182157270.18687@malasada.lava.net>
Content-Description: bug fix
Content-Disposition: attachment; filename="new-stars.r"

ZnVuY3Rpb24gKHgsIGZ1bGwgPSBUUlVFLCBzY2FsZSA9IFRSVUUsIHJhZGl1
cyA9IFRSVUUsIGxhYmVscyA9IGRpbW5hbWVzKHgpW1sxXV0sIA0KICAgICAg
ICAgICAgbG9jYXRpb25zID0gTlVMTCwgeGxpbWl0ID0gTlVMTCwgeWxpbWl0
ID0gTlVMTCwgbGVuID0gMSwgDQogICAgICAgICAgICBjb2xvcnMgPSBOVUxM
LCBrZXkubG9jID0gTlVMTCwga2V5LmxhYmVscyA9IE5VTEwsIGRyYXcuc2Vn
bWVudHMgPSBGQUxTRSwgDQogICAgICAgICAgICBkcmF3LmF4ZXMgPSBGQUxT
RSwgLi4uKSANCnsNCiAgaWYgKGlzLmRhdGEuZnJhbWUoeCkpIA0KICAgIHgg
PC0gYXMubWF0cml4KHgpDQogIGVsc2UgaWYgKCFpcy5tYXRyaXgoeCkpIA0K
ICAgIHN0b3AoInggbXVzdCBiZSBhIG1hdHJpeCBvciBhIGRhdGEgZnJhbWUi
KQ0KICBpZiAoIWlzLm51bWVyaWMoeCkpIA0KICAgIHN0b3AoImRhdGEgaW4g
eCBtdXN0IGJlIG51bWVyaWMiKQ0KICBuLmxvYyA8LSBucm93KHgpDQogIG4u
c2VnIDwtIG5jb2woeCkNCiAgZGVnIDwtIHBpLzE4MA0KICBzZWcuY29sb3Jz
IDwtIGlmICghaXMubnVsbChjb2xvcnMpKSANCiAgICBjb2xvcnMNCiAgZWxz
ZSAxOm4uc2VnDQogIGlmIChpcy5udWxsKGxvY2F0aW9ucykpIHsNCiAgICBt
YXQuZGltIDwtIGNlaWxpbmcoc3FydChuLmxvYykpDQogICAgdGVtcC5sb2Mu
MSA8LSByZXAoeCA9IDIuMSAqIDE6bWF0LmRpbSwgdGltZXMgPSBtYXQuZGlt
LCANCiAgICAgICAgICAgICAgICAgICAgICBsZW5ndGggPSBuLmxvYykNCiAg
ICB0ZW1wLmxvYy4yIDwtIHJlcCh4ID0gMi4xICogbWF0LmRpbToxLCByZXAo
eCA9IG1hdC5kaW0sIA0KICAgICAgICAgICAgICAgICAgICAgICAgdGltZXMg
PSBtYXQuZGltKSwgbGVuZ3RoID0gbi5sb2MpDQogICAgbG9jIDwtIG1hdHJp
eChkYXRhID0gYyh0ZW1wLmxvYy4xLCB0ZW1wLmxvYy4yKSwgbmNvbCA9IDIp
DQogIH0NCiAgZWxzZSB7DQogICAgaWYgKCFpcy5tYXRyaXgobG9jYXRpb25z
KSB8fCBuY29sKGxvY2F0aW9ucykgIT0gMikgDQogICAgICBzdG9wKCJsb2Nh
dGlvbnMgbXVzdCBiZSBhIDItY29sdW1uIG1hdHJpeC4iKQ0KICAgIGxvYyA8
LSAuQWxpYXMobG9jYXRpb25zKQ0KICB9DQogIGlmIChuLmxvYyAhPSBucm93
KGxvYykpIA0KICAgIHN0b3AoIm51bWJlciBvZiByb3dzIG9mIGxvY2F0aW9u
cyBhbmQgeCBtdXN0IGJlIGVxdWFsLiIpDQogIGFuZ2xlcyA8LSBpZiAoZnVs
bCkgDQogICAgc2VxKDAsIDIgKiBwaSwgbGVuZ3RoID0gbi5zZWcgKyAxKVst
KG4uc2VnICsgMSldDQogIGVsc2UgaWYgKGRyYXcuc2VnbWVudHMpIA0KICAg
IHNlcSgwLCBwaSwgbGVuZ3RoID0gbi5zZWcgKyAxKVstKG4uc2VnICsgMSld
DQogIGVsc2Ugc2VxKDAsIHBpLCBsZW5ndGggPSBuLnNlZykNCiAgaWYgKGxl
bmd0aChhbmdsZXMpICE9IG4uc2VnKSANCiAgICBzdG9wKCJsZW5ndGgoYW5n
bGVzKSBtdXN0IGJlIHRoZSBzYW1lIGFzIG5jb2woeCkiKQ0KICB4W2lzLm5h
KHgpXSA8LSAwDQogIGlmIChzY2FsZSkgDQogICAgeCA8LSBzd2VlcCh4LCAy
LCBhcHBseSh4LCAyLCBtYXgpLCBGVU4gPSAiLyIpDQojIENvbHVtbnMgb2Yg
MHMgd2lsbCBwdXQgTkFzIGluIHgsIG5leHQgbGluZSBnZXRzIHJpZCBvZiB0
aGVtDQogIHhbaXMubmEoeCldIDwtIDANCiAgeCA8LSB4ICogbGVuDQogIHRl
bXAueGxpbSA8LSBpZiAoaXMubnVsbCh4bGltaXQpKSANCiAgICByYW5nZShs
b2NbLCAxXSArIG1heCh4KSwgbG9jWywgMV0gLSBtYXgoeCkpDQogIGVsc2Ug
eGxpbWl0DQogIHRlbXAueWxpbSA8LSBpZiAoaXMubnVsbCh5bGltaXQpKSAN
CiAgICByYW5nZShsb2NbLCAyXSArIG1heCh4KSwgbG9jWywgMl0gLSBtYXgo
eCkpDQogIGVsc2UgeWxpbWl0DQogIG9wYXIgPC0gcGFyKG5vLnJlYWRvbmx5
ID0gVFJVRSkNCiAgb24uZXhpdChwYXIob3BhcikpDQogIHBsb3QoMCwgdHlw
ZSA9ICJuIiwgLi4uLCB4bGltID0gdGVtcC54bGltLCB5bGltID0gdGVtcC55
bGltLCANCiAgICAgICB4bGFiID0gIiIsIHlsYWIgPSAiIiwgYXNwID0gMSwg
YXhlcyA9IGRyYXcuYXhlcykNCiAgaWYgKGRyYXcuc2VnbWVudHMpIHsNCiAg
ICBmb3IgKGkgaW4gMTpuLmxvYykgew0KICAgICAgcG9seS54IDwtIE5BDQog
ICAgICBwb2x5LnkgPC0gTkENCiAgICAgIHN0YXJ0LnguY29vcmQgPC0geFtp
LCBdICogY29zKGFuZ2xlcykgKyBsb2NbaSwgMV0NCiAgICAgIHN0YXJ0Lnku
Y29vcmQgPC0geFtpLCBdICogc2luKGFuZ2xlcykgKyBsb2NbaSwgMl0NCiAg
ICAgIGZvciAoaiBpbiAxOm4uc2VnKSB7DQogICAgICAgIHBvbHkueCA8LSBj
KHBvbHkueCwgbG9jW2ksIDFdLCBzdGFydC54LmNvb3JkW2pdKQ0KICAgICAg
ICBwb2x5LnkgPC0gYyhwb2x5LnksIGxvY1tpLCAyXSwgc3RhcnQueS5jb29y
ZFtqXSkNCiAgICAgICAgbmV4dC5hbmdsZSA8LSBpZiAoaiA8IG4uc2VnKSAN
CiAgICAgICAgICBhbmdsZXNbaiArIDFdDQogICAgICAgIGVsc2UgKGlmIChm
dWxsKSANCiAgICAgICAgICAgICAgMzYwDQogICAgICAgIGVsc2UgMTgwKSAq
IGRlZw0KICAgICAgICBrIDwtIHNlcShmcm9tID0gYW5nbGVzW2pdLCB0byA9
IG5leHQuYW5nbGUsIGJ5ID0gZGVnKQ0KICAgICAgICBwb2x5LnggPC0gYyhw
b2x5LngsIHhbaSwgal0gKiBjb3MoaykgKyBsb2NbaSwgDQogICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAxXSwg
TkEpDQogICAgICAgIHBvbHkueSA8LSBjKHBvbHkueSwgeFtpLCBqXSAqIHNp
bihrKSArIGxvY1tpLCANCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgIDJdLCBOQSkNCiAgICAgIH0NCiAgICAg
IHBhcihsd2QgPSAwLjI1KQ0KICAgICAgcG9seWdvbihwb2x5LngsIHBvbHku
eSwgY29sID0gc2VnLmNvbG9ycykNCiAgICAgIHBhcihsd2QgPSAxKQ0KICAg
ICAgaWYgKCFpcy5udWxsKGxhYmVscykpIA0KICAgICAgICB0ZXh0KGxvY1tp
LCAxXSwgbG9jW2ksIDJdIC0gaWYgKGZ1bGwpIA0KICAgICAgICAgICAgIG1h
eCh4KQ0KICAgICAgICBlbHNlIDAuMSAqIG1heCh4KSwgbGFiZWxzW2ldLCBj
ZXggPSAwLjUsIGFkaiA9IGMoMC41LCANCiAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIDEpKQ0KICAgIH0NCiAg
fQ0KICBlbHNlIHsNCiAgICBmb3IgKGkgaW4gMTpuLmxvYykgew0KICAgICAg
dGVtcC54LmNvb3JkIDwtIHhbaSwgXSAqIGNvcyhhbmdsZXMpICsgbG9jW2ks
IDFdDQogICAgICB0ZW1wLnkuY29vcmQgPC0geFtpLCBdICogc2luKGFuZ2xl
cykgKyBsb2NbaSwgMl0NCiAgICAgIGlmIChyYWRpdXMpIHsNCiAgICAgICAg
cGFyKGx3ZCA9IDAuMjUpDQogICAgICAgIHNlZ21lbnRzKHJlcChsb2NbaSwg
MV0sIG4uc2VnKSwgcmVwKGxvY1tpLCAyXSwgDQogICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgIG4uc2VnKSwgdGVtcC54LmNv
b3JkLCB0ZW1wLnkuY29vcmQpDQogICAgICAgIHBhcihsd2QgPSAxKQ0KICAg
ICAgfQ0KICAgICAgbGluZXMoYyh0ZW1wLnguY29vcmQsIHRlbXAueC5jb29y
ZFsxXSksIGModGVtcC55LmNvb3JkLCANCiAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgIHRlbXAueS5jb29yZFsxXSks
IGx3ZCA9IDAuMjUpDQogICAgICBpZiAoIWlzLm51bGwobGFiZWxzKSkgDQog
ICAgICAgIHRleHQobG9jW2ksIDFdLCBsb2NbaSwgMl0gLSBpZiAoZnVsbCkg
DQogICAgICAgICAgICAgbWF4KHgpDQogICAgICAgIGVsc2UgMC4xICogbWF4
KHgpLCBsYWJlbHNbaV0sIGNleCA9IDAuNSwgYWRqID0gYygwLjUsIA0KICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgMSkpDQogICAgfQ0KICB9DQogIGlmICghaXMubnVsbChrZXkubG9jKSkg
ew0KICAgIGlmIChkcmF3LnNlZ21lbnRzKSB7DQogICAgICBrZXkueCA8LSBO
QQ0KICAgICAga2V5LnkgPC0gTkENCiAgICAgIGtleS54LmNvb3JkIDwtIGNv
cyhhbmdsZXMpICogbGVuICsga2V5LmxvY1sxXQ0KICAgICAga2V5LnkuY29v
cmQgPC0gc2luKGFuZ2xlcykgKiBsZW4gKyBrZXkubG9jWzJdDQogICAgICBm
b3IgKGogaW4gMTpuLnNlZykgew0KICAgICAgICBrZXkueCA8LSBjKGtleS54
LCBrZXkubG9jWzFdLCBrZXkueC5jb29yZFtqXSkNCiAgICAgICAga2V5Lnkg
PC0gYyhrZXkueSwga2V5LmxvY1syXSwga2V5LnkuY29vcmRbal0pDQogICAg
ICAgIGsgPC0gYW5nbGVzW2pdICsgZGVnDQogICAgICAgIG5leHQuYW5nbGUg
PC0gaWYgKGogPCBuLnNlZykgDQogICAgICAgICAgYW5nbGVzW2ogKyAxXQ0K
ICAgICAgICBlbHNlIChpZiAoZnVsbCkgDQogICAgICAgICAgICAgIDM2MA0K
ICAgICAgICBlbHNlIDE4MCkgKiBkZWcNCiAgICAgICAgd2hpbGUgKGsgPCBu
ZXh0LmFuZ2xlKSB7DQogICAgICAgICAga2V5LnggPC0gYyhrZXkueCwgbGVu
ICogY29zKGspICsga2V5LmxvY1sxXSkNCiAgICAgICAgICBrZXkueSA8LSBj
KGtleS55LCBsZW4gKiBzaW4oaykgKyBrZXkubG9jWzJdKQ0KICAgICAgICAg
IGsgPC0gayArIGRlZw0KICAgICAgICB9DQogICAgICAgIGtleS54IDwtIGMo
a2V5LngsIGxlbiAqIGNvcyhuZXh0LmFuZ2xlKSArIGtleS5sb2NbMV0sIA0K
ICAgICAgICAgICAgICAgICAgIE5BKQ0KICAgICAgICBrZXkueSA8LSBjKGtl
eS55LCBsZW4gKiBzaW4obmV4dC5hbmdsZSkgKyBrZXkubG9jWzJdLCANCiAg
ICAgICAgICAgICAgICAgICBOQSkNCiAgICAgIH0NCiAgICAgIHBhcihsd2Qg
PSAwLjI1KQ0KICAgICAgcG9seWdvbihrZXkueCwga2V5LnksIGNvbCA9IHNl
Zy5jb2xvcnMpDQogICAgICBwYXIobHdkID0gMSkNCiAgICB9DQogICAgZWxz
ZSB7DQogICAgICB0ZW1wLnguY29vcmQgPC0gY29zKGFuZ2xlcykgKiBsZW4g
KyBrZXkubG9jWzFdDQogICAgICB0ZW1wLnkuY29vcmQgPC0gc2luKGFuZ2xl
cykgKiBsZW4gKyBrZXkubG9jWzJdDQogICAgICBwYXIobHdkID0gMC4yNSkN
CiAgICAgIGlmIChyYWRpdXMpIA0KICAgICAgICBzZWdtZW50cyhyZXAoa2V5
LmxvY1sxXSwgbi5zZWcpLCByZXAoa2V5LmxvY1syXSwgDQogICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBuLnNlZyksIHRl
bXAueC5jb29yZCwgdGVtcC55LmNvb3JkKQ0KICAgICAgbGluZXMoYyh0ZW1w
LnguY29vcmQsIHRlbXAueC5jb29yZFsxXSksIGModGVtcC55LmNvb3JkLCAN
CiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgIHRlbXAueS5jb29yZFsxXSkpDQogICAgICBwYXIobHdkID0gMSkNCiAg
ICB9DQogICAgaWYgKGlzLm51bGwoa2V5LmxhYmVscykpIA0KICAgICAga2V5
LmxhYmVscyA8LSBkaW1uYW1lcyh4KVtbMl1dDQogICAgbGFiLmFuZ2wgPC0g
YW5nbGVzICsgaWYgKGRyYXcuc2VnbWVudHMpIA0KICAgICAgKGFuZ2xlc1sy
XSAtIGFuZ2xlc1sxXSkvMg0KICAgIGVsc2UgMA0KICAgIGxhYmVsLnguY29v
cmQgPC0gY29zKGxhYi5hbmdsKSAqIGxlbiAqIDEuMSArIGtleS5sb2NbMV0N
CiAgICBsYWJlbC55LmNvb3JkIDwtIHNpbihsYWIuYW5nbCkgKiBsZW4gKiAx
LjEgKyBrZXkubG9jWzJdDQogICAgZm9yIChrIGluIDE6bi5zZWcpIHsNCiAg
ICAgIHRleHQuYWRqIDwtIGlmIChsYWIuYW5nbFtrXSA8ICg5MCAqIGRlZykg
fHwgbGFiLmFuZ2xba10gPiANCiAgICAgICAgICAgICAgICAgICAgICAoMjcw
ICogZGVnKSkgDQogICAgICAgIDANCiAgICAgIGVsc2UgaWYgKGxhYi5hbmds
W2tdID4gKDkwICogZGVnKSAmJiBsYWIuYW5nbFtrXSA8IA0KICAgICAgICAg
ICAgICAgKDI3MCAqIGRlZykpIA0KICAgICAgICAxDQogICAgICBlbHNlIDAu
NQ0KICAgICAgaWYgKGxhYi5hbmdsW2tdIDw9ICg5MCAqIGRlZykpIA0KICAg
ICAgICB0ZXh0LmFkaiA8LSBjKHRleHQuYWRqLCAwLjUgKiAoMSAtIGxhYi5h
bmdsW2tdLyg5MCAqIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICBkZWcpKSkNCiAgICAgIGVsc2Ug
aWYgKGxhYi5hbmdsW2tdID4gKDkwICogZGVnKSAmIGxhYi5hbmdsW2tdIDw9
IA0KICAgICAgICAgICAgICAgKDI3MCAqIGRlZykpIA0KICAgICAgICB0ZXh0
LmFkaiA8LSBjKHRleHQuYWRqLCAobGFiLmFuZ2xba10gLSAoOTAgKiANCiAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
IGRlZykpLygxODAgKiBkZWcpKQ0KICAgICAgZWxzZSBpZiAobGFiLmFuZ2xb
a10gPiAoMjcwICogZGVnKSkgDQogICAgICAgIHRleHQuYWRqIDwtIGModGV4
dC5hZGosIDEgLSAoMC41ICogKGxhYi5hbmdsW2tdIC0gDQogICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICgyNzAgKiBkZWcp
KS8oOTAgKiBkZWcpKSkNCiAgICAgIHRleHQuZGVmYXVsdCh4ID0gbGFiZWwu
eC5jb29yZFtrXSwgeSA9IGxhYmVsLnkuY29vcmRba10sIA0KICAgICAgICAg
ICAgICAgICAgIGxhYmVscyA9IGtleS5sYWJlbHNba10sIGNleCA9IDAuNSwg
YWRqID0gdGV4dC5hZGopDQogICAgfQ0KICB9DQogIGludmlzaWJsZSgpDQp9
DQo=
--1465796871-645734172-974620647=:18687--

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._