bug 11079 fixed: tlist(['r','num','den','dt'] replaced by rlist
[scilab.git] / scilab / modules / overloading / macros / %p_j_s.sci
index 333a29f..12b0ab1 100644 (file)
@@ -1,17 +1,26 @@
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) INRIA
+// 
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at    
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
 function f=%p_j_s(p,s)
 // %p_j_s(p,s)  computes p.^s for p polynomial matrix in special cases
 //!
-// Copyright INRIA
+
 if s==[] then f=[],return,end
 if  or(imag(s)<>0)|or(int(s)<>s) then error('%p_j_s: integer power only'),end
 [m,n]=size(p)
 [ms,ns]=size(s)
 if ms==1&ns==1 then
   if s<0 then 
-    if or(abs(coeff(p(:)))*ones(maxi(degree(p))+1,1)==0) then
+    if or(abs(coeff(p(:)))*ones(max(degree(p))+1,1)==0) then
       error(27)
     end
-    f=tlist(['r','num','den','dt'],ones(p),p.^(-s),[])
+    f=rlist(ones(p),p.^(-s),[])
   else // this case is in fact hard coded
     f=p.^s
   end
@@ -24,7 +33,7 @@ elseif m==1&n==1 then // Element wise exponentiation p.^s with p "scalar"
   p=1/p
   num(kn)=p(2).^(-s(kn))
   den(kn)=p(3).^(-s(kn))
-  f=tlist(['r','num','den','dt'],num,den,[])
+  f=rlist(num,den,[])
 elseif ms==m&ns==n then  // Element wise exponentiation
   p=p(:);s=s(:);
   kp=find(s>=0)
@@ -32,12 +41,12 @@ elseif ms==m&ns==n then  // Element wise exponentiation
   num=p
   den=ones(s)
   num(kp)=num(kp).^s(kp)
-  if or(abs(coeff(p(kn)))*ones(maxi(degree(p(kn)))+1,1)==0) then
+  if or(abs(coeff(p(kn)))*ones(max(degree(p(kn)))+1,1)==0) then
     error(27)
   end
   num(kn)=ones(p(kn))
   den(kn)=p(kn).^(-s(kn))
-  f=tlist(['r','num','den','dt'],matrix(num,n,m),matrix(den,n,m),[])
+  f = rlist(matrix(num,n,m),matrix(den,n,m),[])
 else
   error(30)
 end