From help-request at octave dot org Sun Jul 17 23:53:44 2005 Subject: re: bicubic spline interpolation in octave From: Hoxide Ma To: help at octave dot org Date: Mon, 18 Jul 2005 12:51:44 +0800 (CST) --0-26952336-1121662304=:11106 Content-Type: text/plain; charset=gb2312 Content-Transfer-Encoding: 8bit Content-Id: Content-Disposition: inline new version of bicubic.m : now you can try it, by run as follow: c:\> octave bincubic_t.m :) A student of Math in Soochow University in China. Intrested in Math, Python, Twisted, Prolog and .NET/MONO. ___________________________________________________________ ÑÅ»¢Ãâ·ÑGÓÊÏ䣭No.1µÄ·À¶¾·ÀÀ¬»ø³¬´óÓÊÏä http://cn.mail.yahoo.com/?id=77072 --0-26952336-1121662304=:11106 Content-Type: application/octet-stream; name="bicubic.m" Content-Transfer-Encoding: base64 Content-Description: 1097144256-bicubic.m Content-Disposition: attachment; filename="bicubic.m" IyMgQ29weXJpZ2h0IChDKSAyMDA1IEhveGlkZQ0KIyMNCiMjIFRoaXMgcHJv Z3JhbSBpcyBmcmVlIHNvZnR3YXJlOyB5b3UgY2FuIHJlZGlzdHJpYnV0ZSBp dCBhbmQvb3IgbW9kaWZ5DQojIyBpdCB1bmRlciB0aGUgdGVybXMgb2YgdGhl IEdOVSBHZW5lcmFsIFB1YmxpYyBMaWNlbnNlIGFzIHB1Ymxpc2hlZCBieQ0K IyMgdGhlIEZyZWUgU29mdHdhcmUgRm91bmRhdGlvbjsgZWl0aGVyIHZlcnNp b24gMiBvZiB0aGUgTGljZW5zZSwgb3INCiMjIChhdCB5b3VyIG9wdGlvbikg YW55IGxhdGVyIHZlcnNpb24uDQojIw0KIyMgVGhpcyBwcm9ncmFtIGlzIGRp c3RyaWJ1dGVkIGluIHRoZSBob3BlIHRoYXQgaXQgd2lsbCBiZSB1c2VmdWws DQojIyBidXQgV0lUSE9VVCBBTlkgV0FSUkFOVFk7IHdpdGhvdXQgZXZlbiB0 aGUgaW1wbGllZCB3YXJyYW50eSBvZg0KIyMgTUVSQ0hBTlRBQklMSVRZIG9y IEZJVE5FU1MgRk9SIEEgUEFSVElDVUxBUiBQVVJQT1NFLiAgU2VlIHRoZQ0K IyMgR05VIEdlbmVyYWwgUHVibGljIExpY2Vuc2UgZm9yIG1vcmUgZGV0YWls cy4NCiMjDQojIyBZb3Ugc2hvdWxkIGhhdmUgcmVjZWl2ZWQgYSBjb3B5IG9m IHRoZSBHTlUgR2VuZXJhbCBQdWJsaWMgTGljZW5zZQ0KIyMgYWxvbmcgd2l0 aCB0aGlzIHByb2dyYW07IGlmIG5vdCwgd3JpdGUgdG8gdGhlIEZyZWUgU29m dHdhcmUNCiMjIEZvdW5kYXRpb24sIEluYy4sIDU5IFRlbXBsZSBQbGFjZSwg U3VpdGUgMzMwLCBCb3N0b24sIE1BICAwMjExMS0xMzA3ICBVU0ENCg0KIyMg VGhpcyBwcm9ncmFtIGJhc2VkIG9uIGludGVycDIubSB3aGF0IGlzIGluIHRo ZSBNYXRsYWINCiMjIENVQklDIDItRCBiaWN1YmljIGRhdGEgaW50ZXJwb2xh dGlvbi4NCiMjICAgQ1VCSUMoLi4uKSBpcyB0aGUgc2FtZSBhcyBMSU5FQVIo Li4uLikgZXhjZXB0IHRoYXQgaXQgdXNlcw0KIyMgICBiaWN1YmljIGludGVy cG9sYXRpb24uDQojIyAgIA0KIyMgICBUaGlzIGZ1bmN0aW9uIG5lZWRzIGFi b3V0IDctOCB0aW1lcyBTSVpFKFhJKSBtZW1vcnkgdG8gYmUgYXZhaWxhYmxl Lg0KIyMNCiMjICAgU2VlIGFsc28gTElORUFSLg0KIyMgICBDbGF5IE0uIFRo b21wc29uIDQtMjYtOTEsIHJldmlzZWQgNy0zLTkxLCAzLTIyLTkzIGJ5IENN VC4NCg0KDQojIyAgIEJhc2VkIG9uICJDdWJpYyBDb252b2x1dGlvbiBJbnRl cnBvbGF0aW9uIGZvciBEaWdpdGFsIEltYWdlDQojIyAgIFByb2Nlc3Npbmci LCBSb2JlcnQgRy4gS2V5cywgSUVFRSBUcmFucy4gb24gQWNvdXN0aWNzLCBT cGVlY2gsIGFuZA0KIyMgICBTaWduYWwgUHJvY2Vzc2luZywgVm9sLiAyOSwg Tm8uIDYsIERlYy4gMTk4MSwgcHAuIDExNTMtMTE2MC4NCiMjICAgDQojIyAg IFlvdSBjYW4gYXNsbyBzZWUsICJFRkZJQ0lFTlQgSU1BR0UgTUFHTklGSUNB VElPTiBCWSBCSUNVQklDDQojIyAgIFNQTElORSBJTlRFUlBPTEFUSU9OIiwg TWljaGFlbCBKLiBBcmFtaW5pIA0KIyMgICBUaGUgdXJsIGlzIDogDQojIyAg IGh0dHA6Ly9tZW1iZXJzLmJlbGxhdGxhbnRpYy5uZXQvJTdFdnplMnZydmEv ZGVzaWduLmh0bWwgICANCg0KDQpmdW5jdGlvbiBGID0gYmljdWJpYyhhcmcx LGFyZzIsYXJnMyxhcmc0LGFyZzUpDQoNCmlmIG5hcmdpbj09MSwgIyBjdWJp Yyh6KSwgRXhwYW5kIFoNCiAgW25yb3dzLG5jb2xzXSA9IHNpemUoYXJnMSk7 DQogIHMgPSAxOi41Om5jb2xzOyBzaXpzID0gc2l6ZShzKTsNCiAgdCA9ICgx Oi41Om5yb3dzKSc7IHNpenQgPSBzaXplKHQpOw0KDQplbHNlaWYgbmFyZ2lu PT0yLCAjIGN1YmljKHosbiksIEV4cGFuZCBaIG4gdGltZXMNCiAgW25yb3dz LG5jb2xzXSA9IHNpemUoYXJnMSk7DQogIG50aW1lcyA9IGZsb29yKGFyZzIp Ow0KICBzID0gMToxLygyXm50aW1lcyk6bmNvbHM7IHNpenMgPSBzaXplKHMp Ow0KICB0ID0gKDE6MS8oMl5udGltZXMpOm5yb3dzKSc7IHNpenQgPSBzaXpl KHQpOw0KDQplbHNlaWYgbmFyZ2luPT0zLCAjIGN1YmljKHoscyx0KSwgTm8g WCBvciBZIHNwZWNpZmllZC4NCiAgW25yb3dzLG5jb2xzXSA9IHNpemUoYXJn MSk7DQogIHMgPSBhcmcyOyB0ID0gYXJnMzsNCg0KZWxzZWlmIG5hcmdpbj09 NCwNCiAgZXJyb3IoJ1dyb25nIG51bWJlciBvZiBpbnB1dCBhcmd1bWVudHMu Jyk7DQoNCmVsc2VpZiBuYXJnaW49PTUsICMgY3ViaWMoeCx5LHoscyx0KSwg WCBhbmQgWSBzcGVjaWZpZWQuDQogIFtucm93cyxuY29sc10gPSBzaXplKGFy ZzMpOw0KICBteCA9IHByb2Qoc2l6ZShhcmcxKSk7IG15ID0gcHJvZChzaXpl KGFyZzIpKTsNCiAgaWYgYW55KFtteCBteV0gfj0gW25jb2xzIG5yb3dzXSkg JiAuLi4NCiAgICAgfmlzZXF1YWwoc2l6ZShhcmcxKSxzaXplKGFyZzIpLHNp emUoYXJnMykpDQogICAgZXJyb3IoJ1RoZSBsZW5ndGhzIG9mIHRoZSBYIGFu ZCBZIHZlY3RvcnMgbXVzdCBtYXRjaCBaLicpOw0KICBlbmQNCiAgaWYgYW55 KFtucm93cyBuY29sc108WzMgM10pLCBlcnJvcignWiBtdXN0IGJlIGF0IGxl YXN0IDMtYnktMy4nKTsgZW5kDQogIHMgPSAxICsgKGFyZzQtYXJnMSgxKSkv KGFyZzEobXgpLWFyZzEoMSkpKihuY29scy0xKTsNCiAgdCA9IDEgKyAoYXJn NS1hcmcyKDEpKS8oYXJnMihteSktYXJnMigxKSkqKG5yb3dzLTEpOw0KICAN CmVuZA0KDQppZiBhbnkoW25yb3dzIG5jb2xzXTxbMyAzXSksIGVycm9yKCda IG11c3QgYmUgYXQgbGVhc3QgMy1ieS0zLicpOyBlbmQNCiMjIGlmIH5pc2Vx dWFsKHNpemUocyksc2l6ZSh0KSksIA0KIyMgICBlcnJvcignWEkgYW5kIFlJ IG11c3QgYmUgdGhlIHNhbWUgc2l6ZS4nKTsNCiMjIGVuZA0KDQojIENoZWNr IGZvciBvdXQgb2YgcmFuZ2UgdmFsdWVzIG9mIHMgYW5kIHNldCB0byAxDQpz b3V0ID0gZmluZCgoczwxKXwocz5uY29scykpOw0KaWYgbGVuZ3RoKHNvdXQp PjAsIHMoc291dCkgPSBvbmVzKHNpemUoc291dCkpOyBlbmQNCg0KIyBDaGVj ayBmb3Igb3V0IG9mIHJhbmdlIHZhbHVlcyBvZiB0IGFuZCBzZXQgdG8gMQ0K dG91dCA9IGZpbmQoKHQ8MSl8KHQ+bnJvd3MpKTsNCmlmIGxlbmd0aCh0b3V0 KT4wLCB0KHRvdXQpID0gb25lcyhzaXplKHRvdXQpKTsgZW5kDQoNCiMgTWF0 cml4IGVsZW1lbnQgaW5kZXhpbmcNCiMjICBuZHggPSBmbG9vcih0KStmbG9v cihzLTEpKihucm93cysyKTsNCg0KIyBDb21wdXRlIGludGVwb2xhdGlvbiBw YXJhbWV0ZXJzLCBjaGVjayBmb3IgYm91bmRhcnkgdmFsdWUuDQoNCm5kcyA9 IGZsb29yKHMpOw0KaWYgaXNlbXB0eShzKSwgZCA9IHM7IGVsc2UgZCA9IGZp bmQocz09bmNvbHMpOyBlbmQNCnMoOikgPSAocyAtIGZsb29yKHMpKTsNCmlm IGxlbmd0aChkKT4wLCBzKGQpID0gcyhkKSsxOyBuZHMoZCkgPSBuZHMoZCkt MTsgZW5kDQoNCm5kdCA9IGZsb29yKHQpOw0KIyBDb21wdXRlIGludGVwb2xh dGlvbiBwYXJhbWV0ZXJzLCBjaGVjayBmb3IgYm91bmRhcnkgdmFsdWUuDQpp ZiBpc2VtcHR5KHQpLCBkID0gdDsgZWxzZSBkID0gZmluZCh0PT1ucm93cyk7 IGVuZA0KdCg6KSA9ICh0IC0gZmxvb3IodCkpOw0KaWYgbGVuZ3RoKGQpPjAs IHQoZCkgPSB0KGQpKzE7IG5kdChkKSA9IG5kdChkKS0xOyBlbmQNCmQgPSBb XTsNCg0KDQppZiBuYXJnaW49PTUsDQogICMgRXhwYW5kIHogc28gaW50ZXJw b2xhdGlvbiBpcyB2YWxpZCBhdCB0aGUgYm91bmRhcmllcy4NCiAgenogPSB6 ZXJvcyhzaXplKGFyZzMpKzIpOw0KICB6eigxLDI6bmNvbHMrMSkgPSAzKmFy ZzMoMSw6KS0zKmFyZzMoMiw6KSthcmczKDMsOik7DQogIHp6KDI6bnJvd3Mr MSwyOm5jb2xzKzEpID0gYXJnMzsNCiAgenoobnJvd3MrMiwyOm5jb2xzKzEp ID0gMyphcmczKG5yb3dzLDopLTMqYXJnMyhucm93cy0xLDopK2FyZzMobnJv d3MtMiw6KTsNCiAgenooOiwxKSA9IDMqenooOiwyKS0zKnp6KDosMykrenoo Oiw0KTsNCiAgenooOixuY29scysyKSA9IDMqenooOixuY29scysxKS0zKnp6 KDosbmNvbHMpK3p6KDosbmNvbHMtMSk7DQogIG5yb3dzID0gbnJvd3MrMjsg bmNvbHMgPSBuY29scysyOw0KZWxzZQ0KICAjIEV4cGFuZCB6IHNvIGludGVy cG9sYXRpb24gaXMgdmFsaWQgYXQgdGhlIGJvdW5kYXJpZXMuDQogIHp6ID0g emVyb3Moc2l6ZShhcmcxKSsyKTsNCiAgenooMSwyOm5jb2xzKzEpID0gMyph cmcxKDEsOiktMyphcmcxKDIsOikrYXJnMSgzLDopOw0KICB6eigyOm5yb3dz KzEsMjpuY29scysxKSA9IGFyZzE7DQogIHp6KG5yb3dzKzIsMjpuY29scysx KSA9IDMqYXJnMShucm93cyw6KS0zKmFyZzEobnJvd3MtMSw6KSthcmcxKG5y b3dzLTIsOik7DQogIHp6KDosMSkgPSAzKnp6KDosMiktMyp6eig6LDMpK3p6 KDosNCk7DQogIHp6KDosbmNvbHMrMikgPSAzKnp6KDosbmNvbHMrMSktMyp6 eig6LG5jb2xzKSt6eig6LG5jb2xzLTEpOw0KICBucm93cyA9IG5yb3dzKzI7 IG5jb2xzID0gbmNvbHMrMjsNCmVuZA0KDQojIE5vdyBpbnRlcnBvbGF0ZSB1 c2luZyBjb21wdXRhdGlvbmFsbHkgZWZmaWNpZW50IGFsZ29yaXRobS4NCnQw ID0gKCgyLXQpLip0LTEpLip0Ow0KdDEgPSAoMyp0LTUpLip0Lip0KzI7DQp0 MiA9ICgoNC0zKnQpLip0KzEpLip0Ow0KdCg6KSA9ICh0LTEpLip0Lip0Ow0K DQoNCmZvciBpID0gMTpsZW5ndGgobmRzKQ0KICBzcyA9IHMoaSkgOw0KICBk cyA9IG5kcyhpKTsNCiAgRig6LCBpKSA9ICggenoobmR0LCBkcykuKnQwICsg enoobmR0KzEsIGRzKS4qdDEgKyB6eihuZHQrMiwgZHMpLip0MiArIHp6KG5k dCszLCBkcykuKnQgKSAuLi4NCiAgICAgICogKCgoMi1zcykgKnNzLTEpICpz cykgLi4uDQogICAgICArICggenoobmR0LCBkcysxKS4qdDAgKyB6eihuZHQr MSwgZHMrMSkuKnQxICsgenoobmR0KzIsIGRzKzEpLip0MiArIHp6KG5kdCsz LCBkcysxKS4qdCApIC4uLg0KICAgICAgKiAoKDMqc3MtNSkqc3MuKnNzKzIp IC4uLg0KICAgICAgKyAoIHp6KG5kdCwgZHMrMikuKnQwICsgenoobmR0KzEs IGRzKzIpLip0MSArIHp6KG5kdCsyLCBkcysyKS4qdDIgKyB6eihuZHQrMywg ZHMrMikuKnQgKSAuLi4NCiAgICAgICogKCgoNC0zKnNzKSpzcysxKSpzcykg Li4uDQogICAgICArICggenoobmR0LCBkcyszKS4qdDAgKyB6eihuZHQrMSwg ZHMrMykuKnQxICsgenoobmR0KzIsIGRzKzMpLip0MiArIHp6KG5kdCszLCBk cyszKS4qdCApIC4uLg0KICAgICAgKiAoKHNzLTEpKnNzKnNzKTsNCmVuZGZv cg0KDQpGKDopID0gRi80Ow0KIyMgIE5vdyBzZXQgb3V0IG9mIHJhbmdlIHZh bHVlcyB0byBOYU4uDQppZiBsZW5ndGgoc291dCk+MCwgRihzb3V0KSA9IE5h TjsgZW5kDQppZiBsZW5ndGgodG91dCk+MCwgRih0b3V0KSA9IE5hTjsgZW5k DQoNCmVuZGZ1bmN0aW9uDQoNCg== --0-26952336-1121662304=:11106 Content-Type: application/octet-stream; name="bicubic_t.m" Content-Transfer-Encoding: base64 Content-Description: 4221707798-bicubic_t.m Content-Disposition: attachment; filename="bicubic_t.m" eiA9IFsxICAyICAzICA0DQogICAgIDQgIDMgIDIgIDENCiAgICAgMyAgMiAg MSAgMg0KICAgICBdOw0KW256LCBjel0gPSBzaXplKHopOw0KDQp4ID0gbGlu c3BhY2UoMCwgMSwgY3opOw0KeSA9IGxpbnNwYWNlKDAsIDEsIG56KTsNCnp6 ID0gYmljdWJpYyh6LCAyKTsNCg0KW256eiwgY3p6XSA9IHNpemUoenopOw0K eHggPSBsaW5zcGFjZSgwLCAxLCBjenopOw0KeXkgPSBsaW5zcGFjZSgwLCAx LCBuenopOw0KDQptZXNoKHgsIHksIHopOw0KaW5wdXQoIkVudGVyLi4uIik7 DQpncmF3KCJob2xkIG9uIik7DQptZXNoKHh4LCB5eSwgenopOw0KaW5wdXQo IkVudGVyLi4uIik7DQo= --0-26952336-1121662304=:11106-- ------------------------------------------------------------- Octave is freely available under the terms of the GNU GPL. Octave's home on the web: http://www.octave.org How to fund new projects: http://www.octave.org/funding.html Subscription information: http://www.octave.org/archive.html -------------------------------------------------------------