13 lines
273 B
Mathematica
13 lines
273 B
Mathematica
|
function [b,a] = exclusionCut(x,thesign)
|
||
|
|
|
||
|
|
zv = find(abs(x) <= 1e-5);
|
||
|
|
nz = find(abs(x) > 1e-5);
|
||
|
|
a = spalloc(1,length(x),length(zv));
|
||
|
|
b = length(x)-length(zv)-1;
|
||
|
|
if thesign == -1
|
||
|
|
a(nz) = 1;
|
||
|
|
a(zv) = -1;
|
||
|
|
elseif thesign == 1
|
||
|
|
a(zv) = 1;
|
||
|
|
a(nz) = -1;
|
||
|
|
end
|