From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2048/DEPENDENCIES/xdync.sci | 37 +++++++++++++++++++++++++++++++++++++
 1 file changed, 37 insertions(+)
 create mode 100755 2048/DEPENDENCIES/xdync.sci

(limited to '2048/DEPENDENCIES/xdync.sci')

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
-- 
cgit