From octave-sources-request at bevo dot che dot wisc dot edu Mon Jan 15 17:17:31 2001 Subject: Upper Triangle/Lower Triangle matrix by CHOL From: "Cai Jianming" To: Date: Tue, 16 Jan 2001 07:15:35 +0800 This is a multi-part message in MIME format. ------=_NextPart_000_0021_01C07F8C.1E998360 Content-Type: multipart/alternative; boundary="----=_NextPart_001_0022_01C07F8C.1E998360" ------=_NextPart_001_0022_01C07F8C.1E998360 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hi, The attached file will allow CHOL to return a upper triangle/lower = triangle matrix, so that ldiv can be done more efficiently. Any other = operations will get passed to octave_matrix. I.e. typical code could be: c =3D chol(a) x =3D c\ (c' \ b) To compile: mkoctfile *.cc -o chol.oct and copy chol.oct into your libexec/octave/site/oct directory. Regards, Jianming ------=_NextPart_001_0022_01C07F8C.1E998360 Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable

Hi,

The attached file will = allow CHOL to=20 return a upper triangle/lower triangle matrix, so that ldiv can be done = more=20 efficiently. Any other operations will get passed to octave_matrix. I.e. = typical=20 code could be:

c =3D chol(a)

x =3D c\ (c' \ b)

To compile:

mkoctfile *.cc -o=20 chol.oct

and copy chol.oct into your=20 libexec/octave/site/oct directory.

Regards,

Jianming

------=_NextPart_001_0022_01C07F8C.1E998360-- ------=_NextPart_000_0021_01C07F8C.1E998360 Content-Type: application/x-compressed; name="ov-re-tri.tgz" Content-Transfer-Encoding: base64 Content-Disposition: attachment; filename="ov-re-tri.tgz" H4sIAHLtYjoAA+1Z4VLbSBLOX/QUHd8WJxPb2CaBCsbZaA3J+ooFCpskVG2dapDHWBdZ8o5kg4/l nv26Z0b2yBawt5ts7q7UVWBpuqenp6f7m56RN4qCmuc9+5pUb9Truy9fPqsj7e2+yvwivdp5ufes vtfcqb96tbfTbKB8o7m3+wzqX9UqTdM4YQLg2WefhY+54Sn+/yhtb1lWJ5rMhX89SsDulKHx+vVu hf7vwd+iUQgfa3DEkii0rP7Ij2HoBxzwd8JEAtEQTr2EzXjNstQDsYaCc4ijYXLDBG/BPJqCx0IQ fODHifCvpgmKJcDCwXYkYBwN/OEcG6xpOOACkhGHhItxTNrp5f3JBbznIRcsgLPpVeB7cOx7PIw5 MLSDWuIRH8DVnMStdzR6T48O7yLUyhI/ClvAfeQLmHER4zs0K4DD2ywhCwVEE5Iqo1lzK2DJUjAz t+UUBuCH0r5RNEGLR6gHJ3XjBwFccZjGfDgNKoCS8LHb//H0om85J5fw0Tk/d076ly2UTEYRcvmM Kz3+eBL4qBbtFixM5jT/n47OOz+ivPND97jbv0SDrXfd/slRrwfvTs/BgTPnvN/tXBw753B2cX52 2juqAfQ4f8Jz1lB6Hh004AnzgxgneYnrFKNJwQBGNFvBPe7P0CAGHkbI08thsSAKr+XEdFi0INam yLDpnJ5ddk/eo4XdIYRRUoEb4WMwJJGUoZWzclauAq9eQ5+jezicBczjUIXelDru7NQr8EMUJyT0 kwNQbzYajWpjp75XgYueg7Pa2rasv/ihF0wHHErRrCp4FVewNiqZ7R1Uftv58fSY2pfNg6uAp61m M65tWB0Eg6w0FyIS2aZr4U94nG2LvKQaXf0j2zhNaBFoGH6LsRfCLPIHgCGG8BQELlrsRpPYptZy y7IOj95dnLiHx4dge4jgFWDiOq5AiD8YUxULoFTdqmIe3frhMAJ8/jn82XqLhifzCR+GcHccsQHD 6cG7aeiRl+/h7h5IGdhvZ0zcsfuy7OP5mJa30EEOjz/PYci8JBL+P+XSkEQnGk8opWkJV6QqIFWJ +0oaPfF8POY4Gw8mUewnGF8YgUM/xNUkXWOGvFvQBmB8YMJKxlt/iJORT/r3Ozj/ex/OoY1Z8F1N cng4gFSKnpd9/CH5gR4V85ZRNNGz+CtsgUAtbKnDYGs9aW/5uvBiqWzdoasjGevujAVT7iIYJZg6 Cb7hOtEKJnJZECzacpVqAQ+vMUPssuIPwdb8521owK+/pqsIb6BZRgmAO/kfYCJQmzuN2TUHu0Rr VUIliodDTkW4GJma7q0V42h8bYVd1z1TE4Vi1ER0E0vbdLun270omI7DeGE28rDZ9WMXExPBqg3y 16UhtG0YjxgBobecaKbHAdTV9NZMXxN9syL6kwoTO6O5huKCM0wWXBzkZX2nu4z1dFSgabdoPUSk 6znITHYx9RJetjburA01YYqCFr0QJsgQB3tckc1l2Z4uJq1euw1NWk2ZgPiCM9hQBm2ooUg95sFV hClHJrsBpiTCbRs1BzFvpXLKJEOiDHeKA9bGhl5fnM3+vuDXGH1cpB5AFRurCFJe6MXw2FjEDjrC bpRxbGWV67E4ORhEiO/8DdjLGRridRIP+Q0sTbDJJzVafHecLlEFTBMvJhMuclThgsyvuDueks8X Rt7TD0dnqFe5LDq69kGPgPvIOpZQXmzc6yxQKsww8SLK79v8SOko5krApF3+YOAoLV8kfv7AquUs 01fweMarcid0bwRWCNLtWaDAl7IBWiuIgG1Yp7ZaLTiOPCw7PjDh094V78PWlmJgFcn3ofPixaLl KBwoNlUA37rM/q+lZUH09Q6BT5z/dpsv6/r819zdbaBco9nY3SvOf38GPVQaWyZoGwDvRbifpJvp JoIX7SsSRpPy/rLcUEhhj8tSwE6oULp/SClorcuWTejnaOtrbf0a/s/R+C/DTlkoZ0QUYm+ZHTw8 r/CF6AJ2Vva0rQTPvWWJQnrLzmq0yAXhdMwFIXAU6mOjO9SFdXZ2stMmMDVix+n13c7pyQfXOX+/ 7gfYlKiYa1qK3LNMLYObmbRUQiwZ4w69MOOlRywF5QuQZpgeeWKGrXU3Zzb9RMxd3NawtPQR/ZcK DNdneyvgx40qM9dHFC0LQb1ptNv1stp6FtWtUmIWuEaJq5lY5WqmZMvaTJ7424mY8kXFTEPZ1No+ jm6wnkk3OTxV26RStOstEAehwP8vXuDGnTK8tnjRaIF3EHr4X7KkLh3houKVn7frNbJd7r9qbLMe vMIS97N8Xm62iwJH2aRqrFWbPLLJHHhprEc2meZ+MZu0DbKMsEsnWDgwCiYWXk8DhF2FI6W0isJB n0vt5eWhRsVBTtyrX7mOD9QM8ghtxiGLY/86hJyElMe2TfAHtxUwIW4TxCjORGgaikqXLXuQTGth wno2ZHOBhTFWT/yhXFtBHwKf6ptsiteWSspLBMbyk+KVYqCsUWDNA/IECXaETuZsvAlRXFFBPhEu k8engRvPw4TdpscQ08IVF0hlNqlY6y3dEcVwcABkD9WvMirheyjJhxLsQ0kmT6lMUiXoL8KiRJ3R DYGP6Bxp6JU3HqdnpE7PXc7cgFFi58PoprE6MFt6T2a60v1DF3vbwcCfSbUVDQiGeilB+ns5A6RR k/GQGjazY35Ct84aK4i9JnVJUs01KYlw9icJUjbm5KXEMkR8leo6yWgOUK2Ckkul0osC7QQzPG2d Q+ouwBzg0+MDyINA/MuUbuvGi0x+YhRrcRR37FR/BS71mEscRy9RYOu+FEBmHJuwuwJw4GhdGudy kdnRIyvAQ2BThyaIp+P2pYS81hIgferjHyA6+hodSQz/qp9Q1C9vObavejiyK7G2iSUWtzP31spB 9PHpmQj+H04vnVmVIP0N9QFRrT49Rbkv+YZjvshcM4GzBv9m0JgY7qRZ6fR63fcy6xXc/pbcNKun Sl5CoutN8F5LNBMsViM4xdTce1lpT/ek13eOj10FVtHENfBqaRpNRDfL4dZ6jbgY4wGbhQ/3MrrJ +ct+CsHMLhkPVIAEMkMufYzdWRy6/JdHFSjfKV/QGnVPjtzTTt/5cOSiutOO0z89N1dBX1dLsQ9u //LsqHvoIpZ3+t3Tk57rnBy6h07fMbtUoIT/ZVnwrc9H/+9knPq+2hhPnP+bL3eb+vzf2Gs0dun8 X9+pF+f/P4MW5/8DlX7bWIgM/eva6I21xgqiqv5MlcMc31avWMxzeVhtVmfcy+UZ39JyuPqTWg4n /bKWb+eY0Q6az158f8vjTR7oM3t4dsij8ptuWx/gP9Qt9hhugw9xMS8Rc3O5E3bN8/tNRDWaJpNp fr94Pk7YVS5rlt6pEtfyAoR5/REXkX/S0i1rZ6ZWKosHCu7esOAzF9i2vQ2HdDGsvgrqXvGECSx/ OiSvehkXE3gemKivyZn9BvdU1UzXQTycjhfHHVVayEqpXa/IF1kVthtUTagT5W+8uJKVy9rtUSur Ie+W6rCcGUffXshrpfQuRQlkrmPgt17BtNZ1w+N3KWo0Pwyk6xczXFSahnK4M253WvdyLFnb/N5T co4rjJOWOfJyrN99HkUd2HfGEr4v/WtMVWo/POocO+frlUmG90g5Yt23AIM4E/UDn10XFUlBBRVU UEEFFVRQQQUVVFBBBRVUUEEFFVRQQQUVVFBBBRVUUEEFFVRQQQUVVFBBBRVUUEHfkv4NvE0FcQBQ AAAA ------=_NextPart_000_0021_01C07F8C.1E998360-- _________________________________________________________ Do You Yahoo!? Get your free at yahoo dot com address at http://mail.yahoo.com ------------------------------------------------------------- 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 -------------------------------------------------------------