Stéphane Mottelet [Mon, 21 Jan 2019 09:47:37 +0000 (10:47 +0100)]
comparing rlist(1,0) with %inf or poly(%inf,'x','coeff')

Change-Id: Ice3915bb191763ce6ee2ef05acc40dde089c161a

index 42604b4..a5f13b9 100644 (file)
@@ -14,7 +14,7 @@ function [r]=%p_n_r(l1,l2)
//r%p_n_r(l1,l2) <=>r= (l1<>l2 ) l1 polynomial l2 rational
//!

-    r=degree(l2("den"))==0
+    r=degree(l2("den"))<=0
if r then r=l2("num")./coeff(l2("den"))==l1,end
r=~r
endfunction
index 6b0fba3..03ed6e4 100644 (file)
@@ -14,6 +14,6 @@ function [r]=%p_o_r(l1,l2)
//r=%p_o_r(l1,l2) <=> r=(l1==l2) l1 polynomial  l2 rational
//!

-    r=degree(l2("den"))==0
+    r=degree(l2("den"))<=0
if r then r=l2("num")./coeff(l2("den"))==l1,end
endfunction
index e171343..40b2989 100644 (file)
@@ -13,7 +13,7 @@
function [r]=%r_n_p(l1,l2)
// l1<>l2

-    r=degree(l1("den"))==0
+    r=degree(l1("den"))<=0
if r then r=l1("num")./coeff(l1("den"))==l2,end
r=~r
endfunction
index 6582ff9..45e3fe2 100644 (file)
@@ -13,7 +13,7 @@
function [r]=%r_n_s(l1,l2)
// l1<>l2     rational<>constant

-    r=degree(l1("num"))<=0&degree(l1("den"))==0
+    r=degree(l1("num"))<=0&degree(l1("den"))<=0
if r then r=coeff(l1("num"))./coeff(l1("den"))==l2,end
r=~r
endfunction
index fa5b787..6027289 100644 (file)
@@ -14,6 +14,6 @@ function r=%r_o_p(l1,l2)
// l1==l2  with l1 rational and  l2 polynomial
//!

-    r=degree(l1("den"))==0
+    r=degree(l1("den"))<=0
if r then r=l1("num")./coeff(l1("den"))==l2,end
endfunction
index 1ec9240..b166a1b 100644 (file)
@@ -13,6 +13,6 @@
function r=%r_o_s(l1,l2)
// l1==l2 with l1 rational and l2 constant
//!
-    r=degree(l1("num"))<=0&degree(l1("den"))==0
+    r=degree(l1("num"))<=0&degree(l1("den"))<=0
if r then r=coeff(l1("num"))./coeff(l1("den"))==l2,end
endfunction
index 7ebcc63..7d82a39 100644 (file)
@@ -15,7 +15,7 @@ function [r]=%s_n_r(l1,l2)
// et l2 une fraction rationnelle
//!

-    r=degree(l2("num"))<=0&degree(l2("den"))==0
+    r=degree(l2("num"))<=0&degree(l2("den"))<=0
if r then r=coeff(l2("num"))./coeff(l2("den"))==l1,end
r=~r
endfunction
index 5aee133..72a2929 100644 (file)
@@ -13,6 +13,6 @@
function [r]=%s_o_r(l1,l2)
//%s_o_r(l1,l2) constant==rational
//!
-    r=degree(l2("num"))<=0&degree(l2("den"))==0
+    r=degree(l2("num"))<=0&degree(l2("den"))<=0
if r then r=coeff(l2("num"))./coeff(l2("den"))==l1,end
endfunction
deleted file mode 100644 (file)
index 20212d3..0000000
+++ /dev/null
@@ -1,341 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2009 - DIGITEO
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-// <-- CLI SHELL MODE -->
-s=poly(0,"s");
-z = s^0;
-h=1/(2*s);
-assert_checkequal(size(h), [1 1]);
-assert_checkequal(h.num, 1*z);
-assert_checkequal(h.den, 2*s);
-//basic operations with 2D matrices
-H=[h h];
-assert_checkequal(H.num, [1 1].*z);
-assert_checkequal(H.den, [2*s 2*s]);
-H=[h;h];
-assert_checkequal(H.num, [1;1].*z);
-assert_checkequal(H.den, [2*s;2*s]);
-H=[h 3];
-assert_checkequal(H.num, [1*z 3]);
-assert_checkequal(H.den, [2*s 1]);
-H=[h;5];
-assert_checkequal(H.num, [1*z;5]);
-assert_checkequal(H.den, [2*s;1]);
-H=[3 h];
-assert_checkequal(H.num, [3 1*z]);
-assert_checkequal(H.den, [1 2*s]);
-H=[5;h];
-assert_checkequal(H.num, [5;1*z]);
-assert_checkequal(H.den, [1;2*s]);
-H=[3 h;s 2];
-assert_checkequal(H(1,1), rlist(3,1));
-assert_checkequal(H(1,2), h);
-assert_checkequal(H(1, [2 1]), [h rlist(3,1)]);
-assert_checkequal(H(1, [2 2]), [h h]);
-assert_checkequal(H(1,:), [rlist(3,1) h]);
-assert_checkequal(H(:,1), [rlist(3,1); s]);
-assert_checkequal(H([2 1], :), [s rlist(2,1);rlist(3,1) h]);
-assert_checkequal(H([1 1], :), [rlist(3,1) h;rlist(3,1) h]);
-assert_checkequal(matrix(H, -1, 1), [3*z;s;h;2*z]);
-H=h;H(1,3)=1/s;
-assert_checkequal(H.num, [1 0 1].*z);
-assert_checkequal(H.den, [2*s 1 s]);
-H=h;H(3,1)=1/s;
-assert_checkequal(H.num, [1 0 1]'.*z);
-assert_checkequal(H.den, [2*s 1 s]');
-H=h;H(1,3)=1.5;
-assert_checkequal(H.num, [1 0 1.5].*z);
-assert_checkequal(H.den, [2*s 1 1]);
-H=h;H(3,1)=1.5;
-assert_checkequal(H.num, [1 0 1.5]'.*z);
-assert_checkequal(H.den, [2*s 1 1]');
-H=1.5;H(1,3)=1/s;
-assert_checkequal(H.num, [1.5 0 1].*z);
-assert_checkequal(H.den, [1 1 s]);
-H=1.5;H(3,1)=1/s;
-assert_checkequal(H.num, [1.5 0 1]'.*z);
-assert_checkequal(H.den, [1 1 s]');
-H=[h s;1 h];H(:,1)=[];
-assert_checkequal(H, [s;h]);
-H=[h s;1 h];H(2,:)=[];
-assert_checkequal(H, [h s]);
-H=h+h;c=coeff(H.den,1);
-assert_checkequal(H.num/c, 1*z);
-assert_checkequal(H.den/c, s);
-H=h+1;
-assert_checkequal(H.num, 1+2*s);
-assert_checkequal(H.den, 2*s);
-H=1+h;
-assert_checkequal(H.num, 1+2*s);
-assert_checkequal(H.den, 2*s);
-H=h+[];
-assert_checkequal(H, []);
-H=[]+h;
-assert_checkequal(H, []);
-H=h+s;
-assert_checkequal(H.num, 1+2*s^2);
-assert_checkequal(H.den, 2*s);
-H=s+h;
-assert_checkequal(H.num, 1+2*s^2);
-assert_checkequal(H.den, 2*s);
-H=h-h;
-assert_checkequal(H.num, 0*z);
-assert_checkequal(H.den, 1*s^0);
-H=h-1;
-assert_checkequal(H.num, 1-2*s);
-assert_checkequal(H.den, 2*s);
-H=1-h;
-assert_checkequal(H.num, -1+2*s);
-assert_checkequal(H.den, 2*s);
-H=[h h+1]-1;
-assert_checkequal(H(1,1), h-1);
-assert_checkequal(H(1,2), h);
-H=[h h+1]-2*h;
-assert_checkequal(H.num, [-2 -2+4*s]);
-assert_checkequal(H.den, [4*s 4*s]);
-H=-2*h+[h h+1];
-assert_checkequal(H.num, [-2 -2+4*s]);
-assert_checkequal(H.den, [4*s 4*s]);
-// *
-H=h*h;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 4*s^2);
-H=h*2;
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2*s);
-H=2*h;
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2*s);
-H=h*s;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*z);
-H=s*h;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*z);
-H=[h h+1]*h;
-assert_checkequal(H, [h*h (h+1)*h]);
-H=h*[h h+1];
-assert_checkequal(H, [h*h, (h+1)*h]);
-H=[h h+1]*2;
-assert_checkequal(H, [h*2, (h+1)*2]);
-H=2*[h h+1];
-assert_checkequal(H, [h*2, (h+1)*2]);
-H=[h h+1]*s;
-assert_checkequal(H, [h*s, (h+1)*s]);
-H=s*[h h+1];
-assert_checkequal(H, [h*s, (h+1)*s]);
-H=[h 1;s 3]*[1;h];
-assert_checkequal(H, [h+h;s+3*h]);
-H=[h 1]*[h 1;s 3];
-assert_checkequal(H, [h*h+s, h+3]);
-H=[h 1;s 3]*[1;2];
-assert_checkequal(H, [h+2;s+6]);
-H=[2 1]*[h 1;s 3];
-assert_checkequal(H, [2*h+s, 2+3]);
-// .*
-H=h.*h;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 4*s^2);
-H=h.*2;
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2*s);
-H=2 .*h;
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2*s);
-H=h.*s;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*z);
-H=s.*h;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*z);
-H=[h h+1].*h;
-assert_checkequal(H, [h*h (h+1)*h]);
-H=h.*[h h+1];
-assert_checkequal(H, [h*h (h+1)*h]);
-H=[h h+1].*2;
-assert_checkequal(H, [h*2, (h+1)*2]);
-H=2 .*[h h+1];
-assert_checkequal(H, [h*2, (h+1)*2]);
-H=[h h+1].*s;
-assert_checkequal(H, [h*s, (h+1)*s]);
-H=s.*[h h+1];
-assert_checkequal(H, [h*s, (h+1)*s]);
-H=[3 h;s 2].*[3 h;s 2];
-assert_checkequal(H, [9 h*h; s*s 4]);
-// /
-H=h/2;
-assert_checkequal(H.num, 0.5*z);
-assert_checkequal(H.den, 2*s);
-H=h/s;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*s^2);
-H=h/(h+1);
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2+4*s);
-H=[h h-1]/2;
-assert_checkequal(H, [h/2 (h-1)/2]);
-H=[h h-1]/s;
-assert_checkequal(H, [h/s (h-1)/s]);
-H=1/h;
-assert_checkequal(H.num, 2*s);
-assert_checkequal(H.den, 1*z);
-H=[1 2]/h;
-assert_checkequal(H, [1/h 2/h]);
-H=[s+1 s-2]/h;
-assert_checkequal(H, [(s+1)/h (s-2)/h]);
-H=[h+1 (h+1)*(h-1)]/h;
-assert_checkequal(H, [(h+1)/h ((h+1)*(h-1))/h]);
-H=(eye(2,2)/[3 h;s 2])*[3 h;s 2];
-assert_checkequal(coeff(H.num)./(coeff(H.den, 0) + coeff(H.den, 2)), eye(2,2));
-// ./
-H=h./2;
-assert_checkequal(H.num, 0.5*z);
-assert_checkequal(H.den, 2*s);
-H=h./s;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*s^2);
-H=h./(h+1);
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2+4*s);
-H=[h h-1]./2;
-assert_checkequal(H, [h/2 (h-1)/2]);
-H=[h h-1]./s;
-assert_checkequal(H, [h/s (h-1)/s]);
-H=1 ./h;
-assert_checkequal(H.num, 2*s);
-assert_checkequal(H.den, 1*z);
-H=[1 2]./h;
-assert_checkequal(H, [1/h 2/h]);
-H=[s+1 s-2]./h;
-assert_checkequal(H, [(s+1)/h (s-2)/h]);
-H=[h+1 (h+1)*(h-1)]./h;
-assert_checkequal(H, [(h+1)/h ((h+1)*(h-1))/h]);
-H=[3 h;s 2]./[3 h;s 2];
-assert_checkequal(coeff(H.num)./coeff(H.den), ones(2,2));
-// .\
-H=2 .\h;
-assert_checkequal(H.num, 0.5*z);
-assert_checkequal(H.den, 2*s);
-H=s.\h;
-assert_checkequal(H.num, 1*z);
-assert_checkequal(H.den, 2*s^2);
-H=(h+1).\h;
-assert_checkequal(H.num, 2*z);
-assert_checkequal(H.den, 2+4*s);
-H=2 .\[h h-1];
-assert_checkequal(H, [h/2 (h-1)/2]);
-H=s.\[h h-1];
-assert_checkequal(H, [h/s (h-1)/s]);
-H=1 ./h;
-assert_checkequal(H.num, 2*s);
-assert_checkequal(H.den, 1*z);
-H=h.\[1 2];
-assert_checkequal(H, [1/h 2/h]);
-H=h.\[s+1 s-2];
-assert_checkequal(H, [(s+1)/h (s-2)/h]);
-H=h.\[h+1 (h+1)*(h-1)];
-assert_checkequal(H, [(h+1)/h ((h+1)*(h-1))/h]);
-H=[3 h;s 2]./[3 h;s 2];
-assert_checkequal(coeff(H.num)./coeff(H.den), ones(2,2));
-// hypermatrices of rationnals
-clear H;H(1,1,2)=h;
-assert_checkequal(H(1,1,1), rlist(0, 1, []));
-assert_checkequal(H(1,1,2), h);
-H(2,1,2)=h+1;
-assert_checkequal(size(H), [2 1 2]);
-assert_checkequal(H(:,1,1), rlist([0; 0], [1; 1]));
-assert_checkequal(H([2 2], 1, 2), [h+1;h+1]);
-clear H;H(1,1,2)=h;H(2,1,1)=1;
-assert_checkequal(H(:,1,1), rlist([0; 1], [1; 1]));
-assert_checkequal(H(:,1,2), [h;0]);
-clear H;H(1,1,2)=h;H(1,1,:)=3;
-assert_checkequal(H(1,1,1), rlist(3,1));
-assert_checkequal(H(1,1,2), rlist(3,1));
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;
-assert_checkequal(H(1:2,1,1), rlist([%s;%s], [1;1]));
-clear H;H(2,2,2)=s;H(2,1,1)=h;
-assert_checkequal(H(:,:,1), [0 0;h 0]);
-assert_checkequal(H(:,:,2), rlist([0 0;0 s], [1 1;1 1]));
-clear H;H(2,2,2)=8;H(2,1,1)=h;
-assert_checkequal(H(:,:,1), [0 0;h 0]);
-assert_checkequal(H(:,:,2), rlist([0 0;0 8], [1 1;1 1]));
-clear H;H(2,2,2)=8;H(2,1,1)=h;
-H(:,:,1)=[];
-assert_checkequal(H, rlist([0 0;0 8], [1 1;1 1]));
-clear H;H(2,2,2)=8;H(2,1,1)=h;
-H2=H(2,:,:);H(1,:,:)=[];
-assert_checkequal(H, H2);
-clear H;H(2,2,2)=8;H(2,1,1)=h;
-H2=H(:,2,:);H(:,1,:)=[];
-assert_checkequal(H, H2);
-clear H;H(2,2,2)=h;H=H+1;
-assert_checkequal(H(:,:,1), rlist([1 1;1 1], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [1 1;1 h+1]);
-clear H;H(2,2,2)=h;H=1+H;
-assert_checkequal(H(:,:,1), rlist([1 1;1 1], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [1 1;1 h+1]);
-clear H;H(2,2,2)=h;H=H+s;
-assert_checkequal(H(:,:,1), rlist([s s;s s], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [s s;s h+s]);
-clear H;H(2,2,2)=h;H=s+H;
-assert_checkequal(H(:,:,1), rlist([s s; s s], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [s s;s h+s]);
-clear H;H(2,2,2)=h;H=H-1;
-assert_checkequal(H(:,:,1), -1*rlist([1 1;1 1], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [-1 -1;-1 h-1]);
-clear H;H(2,2,2)=h;H=1-H;
-assert_checkequal(H(:,:,1), rlist([1 1;1 1], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [1 1;1 1-h]);
-clear H;H(2,2,2)=h;H=H-s;
-assert_checkequal(H(:,:,1), rlist([-s -s;-s -s], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [-s -s;-s h-s]);
-clear H;H(2,2,2)=h;H=s-H;
-assert_checkequal(H(:,:,1), rlist([s s;s s], [1 1;1 1]));
-assert_checkequal(H(:,:,2), [s s;s -h+s]);
-clear H;H(2,2,2)=s;H(2,1,1)=h;H=H+H;
-assert_checkequal(H(:,:,1), [0 0;h+h 0]);
-assert_checkequal(H(:,:,2), rlist([0 0;0 s+s], [1 1;1 1]));
-clear H;H(2,2,2)=s;H(2,1,1)=h;H=H-H;
-assert_checkequal(H(:,:,1), [0 0;h-h 0]);
-assert_checkequal(H(:,:,2), rlist([0 0;0 0], [1 1;1 1]));
-// *
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H*h;
-assert_checkequal(H(:,:,1), [s*h;s*h]);
-assert_checkequal(H(:,:,2), [h*h; 0/1]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H*2;
-assert_checkequal(H(:,:,1), rlist([s*2;s*2], [1;1]));
-assert_checkequal(H(:,:,2), [h*2;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=2*H;
-assert_checkequal(H(:,:,1), rlist([s*2;s*2], [1;1]));
-assert_checkequal(H(:,:,2), [h*2;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H*s;
-assert_checkequal(H(:,:,1), rlist([s*s;s*s], [1;1]));
-assert_checkequal(H(:,:,2), [h*s;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=s*H;
-assert_checkequal(H(:,:,1), rlist([s*s;s*s], [1;1]));
-assert_checkequal(H(:,:,2), [h*s;0]);
-// .*
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H.*h;
-assert_checkequal(H(:,:,1), [s*h;s*h]);
-assert_checkequal(H(:,:,2), [h*h;0/1]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=h.*H;
-assert_checkequal(H(:,:,1), [s*h;s*h]);
-assert_checkequal(H(:,:,2), [h*h;0/1]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H.*2;
-assert_checkequal(H(:,:,1), rlist([s*2;s*2], [1;1]));
-assert_checkequal(H(:,:,2), [h*2;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=2 .*H;
-assert_checkequal(H(:,:,1), rlist([s*2;s*2], [1;1]));
-assert_checkequal(H(:,:,2), [h*2;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H.*s;
-assert_checkequal(H(:,:,1), rlist([s*s;s*s], [1;1]));
-assert_checkequal(H(:,:,2), [h*s;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=s.*H;
-assert_checkequal(H(:,:,1), rlist([s*s;s*s], [1;1]));
-assert_checkequal(H(:,:,2), [h*s;0]);
-clear H;H(1,1,2)=h;H(1:2,1,1)=%s;H=H.*H;
-assert_checkequal(H(:,:,1), rlist([s*s;s*s], [1;1]));
-assert_checkequal(H(:,:,2), [h*h;0]);
index 3b438a5..b963e1f 100644 (file)
@@ -5,7 +5,8 @@
//  This file is distributed under the same license as the Scilab package.
// =============================================================================
// <-- CLI SHELL MODE -->
-
+// <-- NO CHECK REF -->
+//

s=poly(0,"s");
z = s^0;
@@ -334,6 +335,18 @@ assert_checkequal(H, [(h+1)/h ((h+1)*(h-1))/h]);
H=[3 h;s 2]./[3 h;s 2];
assert_checkequal(coeff(H.num)./coeff(H.den), ones(2,2));

+//comparison
+//---------------
+assert_checktrue(rlist(1,0)==%inf)
+assert_checktrue(%inf==rlist(1,0))
+assert_checktrue(poly(%inf,"x","coeff")==rlist(1,0))
+assert_checktrue(rlist(1,0)==poly(%inf,"x","coeff"))
+assert_checkfalse(rlist(1,0)<>%inf)
+assert_checkfalse(%inf<>rlist(1,0))
+assert_checkfalse(poly(%inf,"x","coeff")<>rlist(1,0))
+assert_checkfalse(rlist(1,0)<>poly(%inf,"x","coeff"))
+
+
// hypermatrices of rationnals

clear H;H(1,1,2)=h;