1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - 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 [R,Q]=pdiv(P1,P2)
// Element wise euclidan division of a polynomial matrix
// by a polynomial
// This is just a loop for the primitive pppdiv
//!
[lhs,rhs]=argn(0);
[n,m]=size(P1);
[n1,m1]=size(P2);
// Special case for constant matrices
if type(P1)==1&type(P2)==1 then
Q=P1./P2;R=0*P1;
if lhs==1 then R=Q; end
return;
end
R=[],Q=[]
if n1==1 & m1==1 then
for l=1:n,
for k=1:m,
[rlk,qlk]=pppdiv(P1(l,k),P2),R(l,k)=rlk;Q(l,k)=qlk;
end;
end
if lhs==1 then R=Q;end
return;
end
for l=1:n,
for k=1:m,
[rlk,qlk]=pppdiv(P1(l,k),P2(l,k)),R(l,k)=rlk;Q(l,k)=qlk;
end;
end
if lhs==1 then R=Q; end
endfunction
|
