new mpz_fac_ui code

[David Newman]

  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.

--27463816-1920821747-1055672314=:24428
Content-Type: TEXT/PLAIN; charset=US-ASCII

Hi,

Here's some factorial code which uses the factorization of n! and repeated
squaring. It runs about twice as fast as the original code for large enough n.
Unfortunately, 'large enough' is about 15000 on my Athlon 1.2G so a bit of
threshold handling is needed.

Could someone give me a hand with making it less machine-dependent? 
(e.g. the nextprime function and subroutines assume unsigned long is 
32-bit). 

Thanks,
Dave Newman

--27463816-1920821747-1055672314=:24428
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="fac_new.c"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.LNX.4.44.0306151118340.24428@sheba.jesus.ox.ac.uk>
Content-Description: 
Content-Disposition: attachment; filename="fac_new.c"

I2luY2x1ZGUgPHN0ZGlvLmg+DQojaW5jbHVkZSA8dGltZS5oPg0KI2luY2x1
ZGUgPGdtcC5oPg0KDQojaW5jbHVkZSAiZ21wLWltcGwuaCINCiNpbmNsdWRl
ICJsb25nbG9uZy5oIg0KI2luY2x1ZGUgIm1fc3Vicy5oIg0KDQp2b2lkIGZh
Y19uZXcgKCk7DQp1bnNpZ25lZCBsb25nIG5leHRwcmltZSAoKTsNCg0KdW5z
aWduZWQgc2hvcnQgcHJpbWVzW10gPQ0Kew0KICAzLDUsNywxMSwxMywxNywx
OSwyMywyOSwzMSwzNyw0MSw0Myw0Nyw1Myw1OSw2MSw2Nyw3MSw3Myw3OSw4
Myw4OSw5NywNCiAgMTAxLDEwMywxMDcsMTA5LDExMywxMjcsMTMxLDEzNywx
MzksMTQ5LDE1MSwxNTcsMTYzLDE2NywxNzMsMTc5LDE4MSwNCiAgMTkxLDE5
MywxOTcsMTk5LDIxMSwyMjMsMjI3LDIyOSwyMzMsMjM5LDI0MSwyNTEsMjU3
LDI2MywyNjksMjcxLDI3NywNCiAgMjgxLDI4MywyOTMsMzA3LDMxMSwzMTMs
MzE3LDMzMSwzMzcsMzQ3LDM0OSwzNTMsMzU5LDM2NywzNzMsMzc5LDM4MywN
CiAgMzg5LDM5Nyw0MDEsNDA5LDQxOSw0MjEsNDMxLDQzMyw0MzksNDQzLDQ0
OSw0NTcsNDYxLDQ2Myw0NjcsNDc5LDQ4NywNCiAgNDkxLDQ5OSw1MDMsNTA5
LDUyMSw1MjMsNTQxLDU0Nyw1NTcsNTYzLDU2OSw1NzEsNTc3LDU4Nyw1OTMs
NTk5LDYwMSwNCiAgNjA3LDYxMyw2MTcsNjE5LDYzMSw2NDEsNjQzLDY0Nyw2
NTMsNjU5LDY2MSw2NzMsNjc3LDY4Myw2OTEsNzAxLDcwOSwNCiAgNzE5LDcy
Nyw3MzMsNzM5LDc0Myw3NTEsNzU3LDc2MSw3NjksNzczLDc4Nyw3OTcsODA5
LDgxMSw4MjEsODIzLDgyNywNCiAgODI5LDgzOSw4NTMsODU3LDg1OSw4NjMs
ODc3LDg4MSw4ODMsODg3LDkwNyw5MTEsOTE5LDkyOSw5MzcsOTQxLDk0NywN
CiAgOTUzLDk2Nyw5NzEsOTc3LDk4Myw5OTEsOTk3DQp9Ow0KDQojZGVmaW5l
IE5VTUJFUl9PRl9QUklNRVMgMTY3DQoNCi8qIGZvciBzaW5nbGUgbm9uLXpl
cm8gbGltYiAqLw0KI2RlZmluZSBNUFpfU0VUXzFfTlooeixuKSAgICAgICBc
DQogIGRvIHsgICAgICAgICAgICAgICAgICAgICAgICAgIFwNCiAgICBtcHpf
cHRyICBfX3ogPSAoeik7ICAgICAgICAgXA0KICAgIEFTU0VSVCAoKG4pICE9
IDApOyAgICAgICAgICBcDQogICAgUFRSKF9feilbMF0gPSAobik7ICAgICAg
ICAgIFwNCiAgICBTSVooX196KSA9IDE7ICAgICAgICAgICAgICAgXA0KICB9
IHdoaWxlICgwKQ0KDQovKiBmb3Igc2luZ2xlIG5vbi16ZXJvIGxpbWIgKi8N
CiNkZWZpbmUgTVBaX0lOSVRfU0VUXzFfTlooeixuKSAgICAgICAgICAgICAg
ICAgIFwNCiAgZG8geyAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgIFwNCiAgICBtcHpfcHRyICBfX2l6ID0gKHopOyAgICAgICAg
ICAgICAgICAgICAgICAgIFwNCiAgICBBTExPQyhfX2l6KSA9IDE7ICAgICAg
ICAgICAgICAgICAgICAgICAgICAgIFwNCiAgICBQVFIoX19peikgPSBfX0dN
UF9BTExPQ0FURV9GVU5DX0xJTUJTICgxKTsgIFwNCiAgICBNUFpfU0VUXzFf
TlogKF9faXosIG4pOyAgICAgICAgICAgICAgICAgICAgIFwNCiAgfSB3aGls
ZSAoMCkNCg0KLyogZm9yIHNyYz4wIGFuZCBuPjAgKi8NCiNkZWZpbmUgTVBa
X01VTF8xX1BPUyhkc3Qsc3JjLG4pICAgICAgICAgICAgICAgICAgICAgICAg
XA0KICBkbyB7ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICBcDQogICAgbXB6X3B0ciAgICBfX2RzdCA9IChkc3Qp
OyAgICAgICAgICAgICAgICAgICAgICAgICAgIFwNCiAgICBtcHpfc3JjcHRy
IF9fc3JjID0gKHNyYyk7ICAgICAgICAgICAgICAgICAgICAgICAgICAgXA0K
ICAgIG1wX3NpemVfdCAgX19zaXplID0gU0laKF9fc3JjKTsgICAgICAgICAg
ICAgICAgICAgICBcDQogICAgbXBfcHRyICAgICBfX2RzdF9wOyAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgIFwNCiAgICBtcF9saW1iX3QgIF9f
YzsgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgXA0KICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICBcDQogICAgQVNTRVJUIChfX3NpemUgPiAwKTsgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgIFwNCiAgICBBU1NFUlQgKChuKSAhPSAw
KTsgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgXA0KICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICBcDQogICAgTVBaX1JFQUxMT0MgKF9fZHN0LCBfX3NpemUrMSk7ICAg
ICAgICAgICAgICAgICAgICAgIFwNCiAgICBfX2RzdF9wID0gUFRSKF9fZHN0
KTsgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgXA0KICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICBcDQogICAgX19jID0gbXBuX211bF8xIChfX2RzdF9wLCBQVFIoX19zcmMp
LCBfX3NpemUsIG4pOyAgIFwNCiAgICBfX2RzdF9wW19fc2l6ZV0gPSBfX2M7
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgXA0KICAgIFNJWihfX2Rz
dCkgPSBfX3NpemUgKyAoX19jICE9IDApOyAgICAgICAgICAgICAgICAgICBc
DQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgIFwNCiAgfSB3aGlsZSAoMCkNCg0KDQp2b2lkIGZhY19u
ZXcgKG1wel9wdHIgcmVzdWx0LCB1bnNpZ25lZCBsb25nIG4pDQp7DQogIC8q
IFdyaXRlIG4hID0gcFsxXV5yWzFdICogcFsyXV5yWzJdICogLi4uICogcFt0
XV5yW3RdLCB3aGVyZSB0aGUgcFtpXSBhcmUNCiAgICAgcHJpbWVzIChzdHJp
Y3RseSBpbmNyZWFzaW5nKSBhbmQgdGhlIHJbaV0gYXJlIHRoZWlyIGV4cG9u
ZW50cyAoZGVjcmVhc2luZykuDQoNCiAgICAgTm93IG11bHRpcGx5IHRvZ2V0
aGVyIGFsbCB0aGUgcHJpbWVzIHdpdGggdGhlIHNhbWUgZXhwb25lbnQgdXNp
bmcgdGhlDQogICAgIHN0YWNrIGFsZ29yaXRobSBvZiB0aGUgb3JpZ2luYWwg
bXB6X2ZhY191aSBmdW5jdGlvbi4gRm9yIGxvdyB2YWx1ZXMgb2YNCiAgICAg
cltpXSwgdGhlcmUgYXJlIHR5cGljYWxseSBhIGxvdCBvZiBwcmltZXMgcCB3
aXRoIGV4cG9uZW50IHJbaV0sIGZvcg0KICAgICBpbnN0YW5jZSBhbGwgcHJp
bWVzIHAgd2l0aCBuLzIgPCBwIDwgbiBoYXZlIGV4cG9uZW50IDEuDQoNCiAg
ICAgU28gd2UgaGF2ZSBuISA9IFBbMV1eUlsxXSAqIFBbMl1eUlsyXSAqIC4u
LiAqIFBbdV1eUlt1XSwNCiAgICAgd2hlcmUgdGhlIGV4cG9uZW50cyBSW2ld
IGFyZSBub3cgc3RyaWN0bHkgZGVjcmVhc2luZy4NCg0KICAgICBMZXQgRVti
XSBiZSB0aGUgcHJvZHVjdCBvZiBhbGwgdGhlIFBbaV0gd2l0aCBiaXQgYiBz
ZXQgaW4gUltpXS4gVGhlbg0KICAgICB3ZSBjYW4gY2FsY3VsYXRlIG4hIGJ5
IHJlcGVhdGVkIHNxdWFyaW5nLg0KDQogICAgIG4hID0gKCAuLi4gKCgoKHpb
MzFdXjIgKiB6WzMwXSleMiAqIHpbMjldKV4yICogLi4uICleMiAqIHpbMF0u
DQoNCiAgICAgVGhpcyBpcyBhYm91dCB0d2ljZSBhcyBmYXN0IGFzIHRoZSB2
YW5pbGxhIHN0YWNrIGFsZ29yaXRobSBmb3IgbGFyZ2Ugbi4gKi8NCg0KICAv
KiBQb3NzaWJsZSBpbXByb3ZlbWVudHM6DQoNCiAgICAgVGhlIGFsZ29yaXRo
bSBpcyBzbG93IGZvciBzbWFsbCBuLiBFaXRoZXIgbWFrZSBsb29rdXAgdGFi
bGVzLCB1c2UgdGhlDQogICAgIG9yaWdpbmFsIGFsZ29yaXRobSBiZWxvdyBh
IGNlcnRhaW4gdGhyZXNob2xkLCBvciBib3RoLg0KDQogICAgIEZpbmRpbmcg
YSBiZXR0ZXIgd2F5IHRvIGZpbmQgdGhlIHByaW1lcyA8PSBuIHdvdWxkIGdp
dmUgYSBzbWFsbCBzcGVlZHVwLg0KICAgICBJcyBzaWV2aW5nIGZlYXNpYmxl
IGZvciBsYXJnZSBuPw0KDQogICAgIE1vZGlmeSB0aGUgc3RhY2sgYWxnb3Jp
dGhtIHNvIHRoYXQgc3RhY2tbaV0gaGFzIGV4YWN0bHkgMl5pIGxpbWJzLiBU
aGlzDQogICAgIHdvdWxkIG1lYW4gdGhhdCB0aGUgc2l6ZSBvZiBzdGFja1tp
XSBpcyBjb25zdGFudC4NCiAgICAgDQogICAgIFVzZSBsb3ctbGV2ZWwgbWVt
b3J5IGFsbG9jYXRpb24gKi8NCiAgICAgDQogIC8qIEtub3duIGJ1Z3M6DQoN
CiAgICAgbiA9IDB4ZmZmZmZmPz8gKG9uIGEgMzItYml0IHN5c3RlbSkgd2ls
bCBrZWVwIGxvb2tpbmcgZm9yIHByaW1lcyBhcyB3ZQ0KICAgICBuZXZlciBo
YXZlIHAgPiBuLiBIb3dldmVyLCBzaW5jZSB0aGUgaG9zdCBzeXN0ZW0gcHJv
YmFibHkgbGFja3MgdGhlDQogICAgIHN0b3JhZ2UgdG8gZXZhbHVhdGUgKDB4
ZmZmZmZmZmYpISwgdGhpcyBpcyBhIG1pbm9yIHByb2JsZW0uDQogICovDQoN
CiAgI2RlZmluZSBNQVhfU1RBQ0tfU0laRSAzMA0KICBtcHpfdCBzdGFja1tN
QVhfU1RBQ0tfU0laRV07DQogIGludCBzdGFja19wdHIgPSAtMSwgc3RhY2tf
c2l6ZSA9IC0xOw0KICB1bnNpZ25lZCBsb25nIHRyZWVfc2l6ZSA9IDA7DQog
IA0KICB1bnNpZ25lZCBsb25nIHQsIHQwLCB0MTsNCiAgdW5zaWduZWQgbG9u
ZyBpLCBtLCBwLCBwcmltZW51bTsgDQogIHVuc2lnbmVkIGxvbmcgaW5kZXgs
IGxhc3RpbmRleDsNCiAgbXB6X3QgeSwgelszMl07DQogDQogIG1wel9zZXRf
dWkgKHJlc3VsdCwgMSk7ICANCiAgbXB6X2luaXQgKHkpOw0KICANCiAgLyog
eltpXSBhcyBpbiBwcmVhbWJsZS4gV2UgbWlnaHQgYXMgd2VsbCBpbml0aWFs
aXNlIHRoZW0gYWxsLCBidXQgMzAgd291bGQNCiAgICAgcHJvYmFibHkgYmUg
bW9yZSB0aGFuIHN1ZmZpY2llbnQuICovDQogIGZvciAoaSA9IDA7IGkgPCAz
MjsgaSsrKQ0KICAgIG1wel9pbml0X3NldF91aSAoeltpXSwgMSk7DQogIA0K
ICAvLyBzdGFydCB3aXRoIHByaW1lWzFdID0gMywgYXMgMiBpcyB0cmVhdGVk
IHNlcGFyYXRlbHkgKHNlZSBiZWxvdykuDQogIHAgPSAzOyBwcmltZW51bSA9
IDA7DQogIGxhc3RpbmRleCA9IDA7DQogIHQgPSAxOw0KICANCiAgZm9yICg7
OykNCiAgew0KICAgIC8qIHNldCBpbmRleCB0byBiZSB0aGUgZXhwb25lbnQg
b2YgcCBpbiBuIQ0KICAgICAgIGluZGV4ID0gW24vcF0gKyBbbi9wXjJdICsg
W24vcF4zXSArIC4uLiAoWy5dIGRlbm90ZXMgcm91bmRpbmcgZG93bikgKi8N
CiAgICBmb3IgKGluZGV4ID0gMCwgbSA9IG47IG07KQ0KICAgIHsNCiAgICAg
IG0gLz0gcDsNCg0KICAgICAgaW5kZXggKz0gbTsNCiAgICB9DQoNCiAgICBp
ZiAoaW5kZXggIT0gbGFzdGluZGV4IHx8IHAgPiBuKQ0KICAgIHsNCiAgICAg
IC8qIFdlIGhhdmUgYSBjb21wbGV0ZSBQW2ldIChlaXRoZXIgd2UndmUgY2hh
bmdlZCBleHBvbmVudCBvciB3ZSd2ZSBydW4NCiAgICAgICAgIG91dCBvZiBw
cmltZXMpLiBBZGQgaXQgdG8gYWxsIHRoZSBhcHByb3ByaWF0ZSB6W2pdLiAq
Lw0KICAgICAgDQogICAgICBpZiAoc3RhY2tfcHRyIDwgMCkNCiAgICAgIHsN
CiAgICAgICAgLyogc21hbGwgYW1vdW50IG9mIHNwZWVkdXAgaWYgdCBpcyBv
bmUgbGltYi4gKi8NCg0KCWZvciAobSA9IDA7IGxhc3RpbmRleDsgbSsrKQ0K
CXsNCgkgIGlmIChsYXN0aW5kZXggJiAxKSBtcHpfbXVsX3VpICh6W21dLCB6
W21dLCB0KTsNCg0KCSAgbGFzdGluZGV4ID4+PSAxOw0KCX0NCiAgICAgIH0N
CiAgICAgIGVsc2UNCiAgICAgIHsNCiAgICAgICAgLyogbXVsdGlwbHkgdG9n
ZXRoZXIgZXZlcnl0aGluZyBvbiB0aGUgc3RhY2suICovDQoJTVBaX01VTF8x
X1BPUyAoeSwgc3RhY2tbc3RhY2tfcHRyLS1dLCB0KTsNCiANCiAgICAgICAg
d2hpbGUgKHN0YWNrX3B0ciA+PSAwKQ0KICAgICAgICAgIG1wel9tdWwgKHks
IHksIHN0YWNrW3N0YWNrX3B0ci0tXSk7DQogICAgICANCiAgICAgICAgZm9y
IChtID0gMDsgbGFzdGluZGV4OyBtKyspDQogICAgICAgIHsNCiAgICAgICAg
ICBpZiAobGFzdGluZGV4ICYgMSkgbXB6X211bCAoelttXSwgelttXSwgeSk7
DQoNCiAgCSAgbGFzdGluZGV4ID4+PSAxOw0KICAgICAgICB9DQogICAgICB9
DQoNCiAgICAgIGlmIChwID4gbikgYnJlYWs7DQogICAgICANCiAgICAgIGxh
c3RpbmRleCA9IGluZGV4Ow0KICAgICAgdCA9IHA7DQogICAgICB0cmVlX3Np
emUgPSAwOw0KICAgIH0NCiAgICBlbHNlDQogICAgew0KICAgICAgdW11bF9w
cG1tICh0MSwgdDAsIHQsIHApOw0KDQogICAgICBpZiAodDEpDQogICAgICB7
DQogICAgICAgIHRyZWVfc2l6ZSsrOw0KCQ0KCWlmICghKHRyZWVfc2l6ZSAm
IDEpKQ0KCXsNCgkgIE1QWl9NVUxfMV9QT1MgKHN0YWNrW3N0YWNrX3B0cl0s
IHN0YWNrW3N0YWNrX3B0cl0sIHQpOw0KCQ0KCSAgZm9yIChpID0gNDsgISh0
cmVlX3NpemUgJiAoaSAtIDEpKTsgaSA8PD0gMSkNCgkgIHsNCgkgICAgbXB6
X211bCAoc3RhY2tbc3RhY2tfcHRyIC0gMV0sDQoJICAgICAgICAgICAgIHN0
YWNrW3N0YWNrX3B0cl0sIHN0YWNrW3N0YWNrX3B0ciAtIDFdKTsNCg0KICAg
ICAgICAgICAgc3RhY2tfcHRyLS07DQogICAgICAgICAgfQ0KCX0NCgllbHNl
DQoJew0KCSAgaWYgKCsrc3RhY2tfcHRyID4gc3RhY2tfc2l6ZSkNCgkgIHsN
CgkgICAgaWYgKHN0YWNrX3B0ciA+IE1BWF9TVEFDS19TSVpFKQ0KCSAgICAg
IGFib3J0ICgpOw0KDQogICAgICAgICAgICBNUFpfSU5JVF9TRVRfMV9OWiAo
c3RhY2tbc3RhY2tfcHRyXSwgdCk7DQogICAgICAgICAgICBzdGFja19zaXpl
Kys7DQoJICB9DQoJICBlbHNlDQoJICAgIE1QWl9TRVRfMV9OWiAoc3RhY2tb
c3RhY2tfcHRyXSwgdCk7DQoJfQ0KCSAgDQogICAgICAgIHQgPSBwOw0KICAg
ICAgfQ0KICAgICAgZWxzZQ0KICAgICAgICB0ID0gdDA7DQogICAgfQ0KICAg
ICAgDQogICAgLyogZ2V0IHRoZSBuZXh0IHByaW1lLiAqLw0KICAgIHAgPSAo
KytwcmltZW51bSA8IE5VTUJFUl9PRl9QUklNRVMpID8gcHJpbWVzW3ByaW1l
bnVtXSA6IG5leHRwcmltZShwKTsNCiAgfQ0KDQogIHdoaWxlIChzdGFja19z
aXplID49IDApDQogICAgbXB6X2NsZWFyKHN0YWNrW3N0YWNrX3NpemUtLV0p
Ow0KICAgIA0KICBtcHpfY2xlYXIgKHkpOw0KICANCiAgLyogRm9yIGxhcmdl
IG4sIGFib3V0IDgwJSBvZiB0aGUgY29tcHV0YXRpb24gaXMgaW4gdGhlIGZv
bGxvd2luZyBsb29wOw0KICAgICB0aGVyZSdzIG5vdCBtdWNoIHNjb3BlIGZv
ciBvcHRpbWl6YXRpb24uICovDQogIA0KICBmb3IgKGkgPSAwOyBpIDwgMzI7
IGkrKykNCiAgew0KICAgIG1wel9tdWwgKHJlc3VsdCwgcmVzdWx0LCByZXN1
bHQpOw0KICAgIG1wel9tdWwgKHJlc3VsdCwgcmVzdWx0LCB6WzMxLWldKTsN
CiAgDQogICAgbXB6X2NsZWFyICh6WzMxLWldKTsNCiAgfQ0KICANCiAgLyog
RmluZCB0aGUgZXhwb25lbnQgb2YgMiBpbiBuISwgbXVsdGlwbHkgcmVzdWx0
IGJ5IDJeeCBieSBzaGlmdGluZy4gVGhpcw0KICAgICBzYXZlcyBhIGxvdCBv
ZiB0aW1lISAqLw0KICANCiAgZm9yIChpbmRleCA9IDA7IG47IGluZGV4ICs9
IG4pDQogICAgbiA+Pj0gMTsNCg0KICBtcHpfbXVsXzJleHAgKHJlc3VsdCwg
cmVzdWx0LCBpbmRleCk7DQp9DQoNCnVuc2lnbmVkIGxvbmcgbmV4dHByaW1l
ICh1bnNpZ25lZCBsb25nIHApDQp7DQogIGlmIChwID09IDIpIHJldHVybiAz
Ow0KICANCiAgZm9yICg7OykNCiAgew0KICAgIHAgKz0gMjsNCg0KICAgIC8q
IFRoaXMgZ3VhcmFudGVlcyBwcmltYWxpdHkgZm9yIDMyLWJpdCB2YWx1ZXMg
b2YgcCwgc2VlDQogICAgICAgaHR0cDovL3d3dy51dG0uZWR1L3Jlc2VhcmNo
L3ByaW1lcy9wcm92ZS9wcm92ZTJfMy5odG1sICovDQoNCiAgICBpZiAoIW1f
c3BwICgyLCBwKSkgY29udGludWU7DQogICAgaWYgKCFtX3NwcCAoNywgcCkp
IGNvbnRpbnVlOw0KICAgIGlmIChtX3NwcCAoNjEsIHApKSByZXR1cm4gcDsN
CiAgfQ0KfQ0K
--27463816-1920821747-1055672314=:24428
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="m_subs.h"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.LNX.4.44.0306151118341.24428@sheba.jesus.ox.ac.uk>
Content-Description: 
Content-Disposition: attachment; filename="m_subs.h"

c3RhdGljIGlubGluZSB1bnNpZ25lZCBsb25nIG1fbXVsICh1bnNpZ25lZCBs
b25nIHgsIHVuc2lnbmVkIGxvbmcgeSwNCiAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgdW5zaWduZWQgbG9uZyBtKQ0KLy8geHkgbW9kIG0N
CnsNCiAgdW5zaWduZWQgbG9uZyB6Ow0KICBfX2FzbV9fICgibXVsbCAlJWVj
eFxuIg0KICAgICAgICAgICAiZGl2bCAlJWVieCIgOiAiPWQiICh6KSA6ICJh
IiAoeCksICJjIiAoeSksICJiIiAobSkgOiAiY2MiKTsNCiAgcmV0dXJuIHo7
DQp9DQoNCg0Kc3RhdGljIGlubGluZSB1bnNpZ25lZCBsb25nIG1fc3FyICh1
bnNpZ25lZCBsb25nIHgsIHVuc2lnbmVkIGxvbmcgbSkNCi8vIHheMiBtb2Qg
bQ0Kew0KICB1bnNpZ25lZCBsb25nIHo7DQogIF9fYXNtX18gKCJtdWxsICUl
ZWF4XG4iDQogICAgICAgICAgICJkaXZsICUlZWN4IiA6ICI9ZCIgKHopIDog
ICJhIiAoeCksICJjIiAobSkgOiAiY2MiKTsNCiAgcmV0dXJuIHo7DQp9DQoN
Cg0KdW5zaWduZWQgbG9uZyBtX3BvdyAodW5zaWduZWQgbG9uZyB4LCB1bnNp
Z25lZCBsb25nIGEsDQogICAgICAgICAgICAgICAgICAgICB1bnNpZ25lZCBs
b25nIG0pDQovLyB4XmEgbW9kIG0NCnsNCiAgdW5zaWduZWQgbG9uZyBwYXJ0
aWFsOw0KDQogIGZvciAocGFydGlhbCA9IDE7OykNCiAgew0KICAgIGlmIChh
ICYgMSkgcGFydGlhbCA9IG1fbXVsICh4LCBwYXJ0aWFsLCBtKTsNCgkJDQog
ICAgYSA+Pj0gMTsNCgkJDQogICAgaWYgKCFhKSByZXR1cm4gcGFydGlhbDsN
Cg0KICAgIHggPSBtX3NxciAoeCwgbSk7DQogIH0NCn0NCg0KDQppbnQgbV9z
cHAgKGEsIG4pDQovLyBSZXR1cm5zIDEgaWZmIG4gaXMgYSBzdHJvbmcgcHNl
dWRvcHJpbWUgdG8gYmFzZSBhDQp7DQogIHVuc2lnbmVkIGxvbmcgciwgcywg
dCwgZDsNCg0KICBmb3IgKHMgPSAwLCBkID0gbi0xOyAhKGQgJiAxKTsgcysr
LCBkID4+PSAxKTsNCg0KICB0ID0gbV9wb3cgKGEsIGQsIG4pOw0KCQ0KICBp
ZiAodCA9PSAxKSByZXR1cm4gMTsNCgkNCiAgZm9yIChyID0gMDsgciA8IHM7
IHIrKykNCiAgew0KICAgIGlmICh0ID09IG4tMSkgcmV0dXJuIDE7DQoNCiAg
ICB0ID0gbV9zcXIgKHQsIG4pOw0KICB9DQoJDQogIHJldHVybiAwOw0KfQ0K

--27463816-1920821747-1055672314=:24428--