annotate lin_solve_dgr_2.m @ 0:a4a7e3405062
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 |
|
| rev |
line source |
|
dan@0
|
1 % solve multiple 2nd degree linear systems at once - no pivoting!
|
|
dan@0
|
2 function dg2 = lin_solve_dgr_2(A,b,R)
|
|
dan@0
|
3 nrm = (-A(2,1,:) .* A(1,2,:) + A(1,1,:) .* A(2,2,:));
|
|
dan@0
|
4 dg2 = zeros(2,R);
|
|
dan@0
|
5 dg2(1,:) = -(b(2,1,:) .* A(1,2,:) - b(1,1,:) .* A(2,2,:)) ./ nrm;
|
|
dan@0
|
6 dg2(2,:) = -(b(1,1,:) .* A(2,1,:) - b(2,1,:) .* A(1,1,:)) ./ nrm;
|
|
dan@0
|
7 end |
