Mercurial > repos > public > sbplib
comparison minors.m @ 886:8894e9c49e40 feature/timesteppers
Merge with default for latest changes
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 15 Nov 2018 16:36:21 -0800 |
| parents | eaf557023fbe |
| children |
comparison
equal
deleted
inserted
replaced
| 816:b5e5b195da1e | 886:8894e9c49e40 |
|---|---|
| 1 function [minor, sub] = minors(A) | 1 function [minor, sub] = minors(A, verbose) |
| 2 default_arg('verbose', true); | |
| 2 [n, m] = size(A); | 3 [n, m] = size(A); |
| 3 | 4 |
| 4 if n ~= m | 5 if n ~= m |
| 5 error('A must be square'); | 6 error('A must be square'); |
| 6 end | 7 end |
| 16 sub{end + 1} = A(C(i,:),C(i,:)); | 17 sub{end + 1} = A(C(i,:),C(i,:)); |
| 17 end | 18 end |
| 18 end | 19 end |
| 19 | 20 |
| 20 for i = 1:length(sub) | 21 for i = 1:length(sub) |
| 21 fprintf('%d:\n', ks{i}); | 22 if verbose |
| 22 disp(sub{i}) | 23 fprintf('%d:\n', ks{i}); |
| 24 disp(sub{i}) | |
| 25 end | |
| 23 | 26 |
| 24 minor(i) = det(sub{i}); | 27 minor(i) = det(sub{i}); |
| 25 end | 28 end |
| 26 end | 29 end |
| 27 | 30 |