diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2048/DEPENDENCIES/xdync.sci | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '2048/DEPENDENCIES/xdync.sci')
| -rwxr-xr-x | 2048/DEPENDENCIES/xdync.sci | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/2048/DEPENDENCIES/xdync.sci b/2048/DEPENDENCIES/xdync.sci new file mode 100755 index 000000000..e082e587f --- /dev/null +++ b/2048/DEPENDENCIES/xdync.sci @@ -0,0 +1,37 @@ +// function [Y,degY,X,degX,B,degB,A,degA] = xdync(N,degN,D,degD,C,degC,gap)
+// given coefficient matrix in T1, primary redundant row information sel,
+// solves XD + YN = C
+
+// calling order changed on 16 April 2005. Old order:
+// function [B,degB,A,degA,Y,degY,X,degX] = xdync(N,degN,D,degD,C,degC,gap)
+
+function [Y,degY,X,degX,B,degB,A,degA] = xdync(N,degN,D,degD,C,degC,gap)
+if argn(2) == 6
+ gap = 1.0e+8;
+end
+
+[F,degF] = rowjoin(D,degD,N,degN);
+
+[Frows,Fbcols] = polsize(F,degF); //Fbcols = block columns
+
+[B,degB,A,degA,S,sel,degT1,Fbcols] = left_prm(N,degN,D,degD,3,gap);
+//if issoln(D,degD,C,degC,B,degB,A,degA)
+ [Crows,Ccols] = size(C);
+ [Srows,Scols] = size(S);
+ S = clean(S);
+ S = S(mtlb_logical(sel),:);
+ T2 =[];
+
+ for i = 1:Crows,
+ Saug = seshft(S,C(i,:),0);
+ b = cindep(Saug);
+ b = move_sci(b,find(sel),Srows);
+ T2 =[T2; b];
+ end
+
+ [X,degX,Y,degY] = colsplit(T2,degT1,Fbcols,Frows-Fbcols);
+
+ [X,degX] = clcoef(X,degX);
+ [Y,degY] = clcoef(Y,degY);
+ Y = clean(Y); X = clean(X);
+endfunction
|