Mercurial > hg > ddm

diff lin_solve_dgr_2.m @ 0:a4a7e3405062

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import DDM code by Sašo Muševič
author Dan Stowell <dan.stowell@elec.qmul.ac.uk>
date Thu, 09 May 2013 20:04:15 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/lin_solve_dgr_2.m	Thu May 09 20:04:15 2013 +0100
@@ -0,0 +1,7 @@
+% solve multiple 2nd degree linear systems at once - no pivoting!
+function dg2 = lin_solve_dgr_2(A,b,R)
+nrm =  (-A(2,1,:) .* A(1,2,:) + A(1,1,:) .* A(2,2,:));
+  dg2 = zeros(2,R);
+  dg2(1,:) = -(b(2,1,:) .* A(1,2,:) - b(1,1,:) .* A(2,2,:)) ./ nrm;
+  dg2(2,:) = -(b(1,1,:) .* A(2,1,:) - b(2,1,:) .* A(1,1,:)) ./ nrm;
+end
\ No newline at end of file