sosfilt.cc
"Eric Chassande-Mottin" <[email protected]> Mon, 25 Feb 2008 12:13:04 +0100
| Newsgroups | gmane.comp.gnu.octave.sources |
|---|---|
| Message-ID | <[email protected]> |
------=_Part_8466_26758615.1203937984075 Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: 7bit Content-Disposition: inline here is a possible implementation for sosfilt(). I have checked that it gives consistent results as compared to matlab's version. eric ------=_Part_8466_26758615.1203937984075 Content-Type: text/x-c++src; name=sosfilt.cc Content-Transfer-Encoding: base64 X-Attachment-Id: f_fd2xju5z1 Content-Disposition: attachment; filename=sosfilt.cc Ci8vIHNvc2ZpbHQuY2MgLS0gU2Vjb25kIG9yZGVyIElJUiBmaWx0ZXJpbmcKCi8vIFRoaXMgZmls ZSBtYXkgYmUgdXNlZCB1bmRlciB0aGUgdGVybXMgZGVmaW5lZCBieSB0aGUgR05VIAovLyBHZW5l cmFsIFB1YmxpYyBMaWNlbnNlIFt3cml0ZSB0byB0aGUgRnJlZSBTb2Z0d2FyZSBGb3VuZGF0aW9u LCAKLy8gNTkgVGVtcGxlIFBsYWNlIC0gU3VpdGUgMzMwLCBCb3N0b24sIE1BICAwMjExMS0xMzA3 LCBVU0FdCi8vCi8vIFdJVEhPVVQgQU5ZIFdBUlJBTlRZOyB3aXRob3V0IGV2ZW4gdGhlIGltcGxp ZWQgd2FycmFudHkgCi8vIG9mIE1FUkNIQU5UQUJJTElUWSBvciBGSVRORVNTIEZPUiBBIFBBUlRJ Q1VMQVIgUFVSUE9TRS4gIAovLwovLyBDb3B5cmlnaHQgKEMpIDIwMDggRXJpYyBDaGFzc2FuZGUt TW90dGluLCBDTlJTIChGcmFuY2UpCgovLwoKI2luY2x1ZGUgPG9jdGF2ZS9jb25maWcuaD4KCiNp bmNsdWRlIDxpb3N0cmVhbT4KI2luY2x1ZGUgPHNzdHJlYW0+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7 CiNpbmNsdWRlIDxzdHJpbmc+CgojaW5jbHVkZSA8b2N0YXZlL2RlZnVuLWRsZC5oPgojaW5jbHVk ZSA8b2N0YXZlL2Vycm9yLmg+CiNpbmNsdWRlIDxvY3RhdmUvb2N0LW9iai5oPgojaW5jbHVkZSA8 b2N0YXZlL3BhZ2VyLmg+CiNpbmNsdWRlIDxvY3RhdmUvc3ltdGFiLmg+CiNpbmNsdWRlIDxvY3Rh dmUvdmFyaWFibGVzLmg+CgojaW5jbHVkZSA8dmVjdG9yPgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoK CkRFRlVOX0RMRCAoc29zZmlsdCwgYXJncywsCiAgIlNlY29uZCBvcmRlciBJSVIgZmlsdGVyaW5n XG4iKQp7CiAgb2N0YXZlX3ZhbHVlX2xpc3QgcmV0dmFsOwogIAogIGludCBuYXJnaW4gPSBhcmdz Lmxlbmd0aCAoKTsKICAKICBpZiAobmFyZ2luIDwgMikKICAgIHsKICAgICAgZXJyb3IoIlVzYWdl OiAgeSA9IHNvc2ZpbHQoc29zLHgpIik7CiAgICAgIHJldHVybiByZXR2YWw7CiAgICB9CgogIE1h dHJpeCBzb3MoIGFyZ3MoMCkubWF0cml4X3ZhbHVlKCkgKTsKCiAgQ29sdW1uVmVjdG9yIHgoIGFy Z3MoMSkudmVjdG9yX3ZhbHVlKCkgKTsKICBpbnQgbj14Lmxlbmd0aCgpOwoKICBpZiAoc29zLmNv bHVtbnMoKSE9NikKICAgIHsKICAgICAgZXJyb3IoIlNlY29uZC1vcmRlciBtYXRyaXggbXVzdCBi ZSBhIG5vbi1lbXB0eSBMeDYgbWF0cml4Iik7CiAgICAgIHJldHVybiByZXR2YWw7CiAgICB9Cgog IENvbHVtblZlY3RvciB5KG4sMC4wKTsKICAKICBmb3IgKGludCBqPTA7IGo8c29zLnJvd3MoKTsg aisrKQogICAgewoKICAgICAgZG91YmxlIHYwPTAuMCwgdjE9MC4wLCB2Mj0wLjA7CgogICAgICBk b3VibGUgYjA9c29zKGosMCk7CiAgICAgIGRvdWJsZSBiMT1zb3MoaiwxKTsKICAgICAgZG91Ymxl IGIyPXNvcyhqLDIpOwogICAgICBkb3VibGUgYTE9c29zKGosNCk7CiAgICAgIGRvdWJsZSBhMj1z b3Moaiw1KTsKCiAgICAgIGZvciAoaW50IGk9MDsgaTxuOyBpKyspCgl7CgkgIHYwPXgoaSktYTEq djEtYTIqdjI7CgkgIHkoaSk9YjAqdjArYjEqdjErYjIqdjI7CgkgIHYyPXYxOwoJICB2MT12MDsK CX0KCiAgICAgIHg9eTsKCiAgICB9CgogIHJldHZhbCgwKT15OwogIAogIHJldHVybiByZXR2YWw7 Cn0KCg== ------=_Part_8466_26758615.1203937984075 Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit Content-Disposition: inline _______________________________________________ Octave-sources mailing list [email protected] https://www.cae.wisc.edu/mailman/listinfo/octave-sources ------=_Part_8466_26758615.1203937984075--