Ticket #494: commands_file_of_Singular_with_PID_38385

File commands_file_of_Singular_with_PID_38385, 66.0 KB (added by barakat, 10 years ago)
Line 
1option(noredefine);option(redSB);LIB "matrix.lib";LIB "involut.lib";LIB "nctools.lib";LIB "poly.lib";LIB "random.lib";LIB "finvar.lib";
2int homalg_variable_i; int homalg_variable_j; int homalg_variable_k; list homalg_variable_l;
3
4
5proc IsDiagonalMatrix (matrix m)
6{
7  int min=nrows(m);
8  if (min>ncols(m))
9  {
10    min=ncols(m);
11  }
12  matrix z[nrows(m)][ncols(m)];
13  matrix c = m;
14  for (int i=1; i<=min; i++)
15  {
16    c[i,i]=0;
17  }
18  return(c==z);
19}
20
21
22proc ZeroRows (module m)
23{
24  list l;
25  int s = 1;
26  for (int i=1;i<=ncols(m);i++)
27  {
28    if (m[i]==0)
29    {
30      l[s]=i; s++;
31    }
32  }
33  if (size(l)==0)
34  {
35    return("[]"));
36  }
37  return(string(l));
38}
39
40
41proc ZeroColumns (matrix n)
42{
43  matrix m=module(transpose(n));
44  list l;
45  int s = 1;
46  for (int i=1;i<=ncols(m);i++)
47  {
48    if (m[i]==0)
49    {
50      l[s]=i; s++;
51    }
52  }
53  if (size(l)==0)
54  {
55    return("[]"));
56  }
57  return(string(l));
58}
59
60
61proc PositionOfFirstNonZeroEntryPerRow (matrix M)
62{
63  int b = 1;
64  intmat m[1][ncols(M)];
65  for (int i=1; i<=ncols(M); i++)
66  {
67    for (int j=1; j<=nrows(M); j++)
68    {
69      if ( M[j,i] <> 0 ) { m[1,i] = j; break; }
70    }
71    if ( b && i > 1 ) { if ( m[1,i] <> m[1,i-1] ) { b = 0; } } // Singular is strange
72  }
73  if ( b ) { return(m[1,1]); } else { return(m); }
74}
75
76
77proc PositionOfFirstNonZeroEntryPerColumn (matrix M)
78{
79  int b = 1;
80  intmat m[1][nrows(M)];
81  for (int j=1; j<=nrows(M); j++)
82  {
83    for (int i=1; i<=ncols(M); i++)
84    {
85      if ( M[j,i] <> 0 ) { m[1,j] = i; break; }
86    }
87    if ( b && j > 1 ) { if ( m[1,j] <> m[1,j-1] ) { b = 0; } } // Singular is strange
88  }
89  if ( b ) { return(m[1,1]); } else { return(m); }
90}
91
92
93proc IndicatorMatrixOfNonZeroEntries(matrix M)
94{
95  intmat m[ncols(M)][nrows(M)];
96  for (int i=1; i<=ncols(M); i++)
97  {
98    for (int j=1; j<=nrows(M); j++)
99    {
100      m[i,j] = ( M[j,i] <> 0 );
101    }
102  }
103  return(m);
104}
105
106
107proc MatrixOfSymbols (matrix m)
108{
109  int i; int j; poly e;
110  int r=nrows(m);
111  int c=ncols(m);
112  matrix n[r][c]=0;
113  for(i=1;i<=r;i++)
114  {
115    for(j=1;j<=c;j++)
116    {
117      e=m[i,j];
118      if(e!=0)
119      { n[i,j]=e-jet(e,deg(e)-1); }
120    }
121  }
122  return(n);
123}
124
125
126proc BasisOfRowModule (matrix M)
127{
128  return(std(M));
129}
130
131
132proc BasisOfColumnModule (matrix M)
133{
134  return(Involution(BasisOfRowModule(Involution(M))));
135}
136
137
138proc SyzygiesGeneratorsOfRows (matrix M)
139{
140  return(SyzForHomalg(M));
141}
142
143
144proc SyzygiesGeneratorsOfColumns (matrix M)
145{
146  return(Involution(SyzForHomalg(Involution(M))));
147}
148
149
150proc ReducedSyzygiesGeneratorsOfRows (matrix M)
151{
152  return(ReducedSyzForHomalg(M));
153}
154
155
156proc ReducedSyzygiesGeneratorsOfColumns (matrix M)
157{
158  return(Involution(ReducedSyzForHomalg(Involution(M))));
159}
160
161
162proc BasisOfRowsCoeff (matrix M)
163{
164  matrix B = BasisOfRowModule(M);
165  matrix T = lift(M,B);
166  list l = B,T;
167  return(l)
168}
169
170
171proc BasisOfColumnsCoeff (matrix M)
172{
173  list l = BasisOfRowsCoeff(Involution(M));
174  matrix B = l[1];
175  matrix T = l[2];
176  l = Involution(B),Involution(T);
177  return(l);
178}
179
180
181proc DecideZeroRows (matrix A, module B)
182{
183  attrib(B,"isSB",1);
184  return(reduce(A,B));
185}
186
187
188proc DecideZeroColumns (matrix A, matrix B)
189{
190  return(Involution(DecideZeroRows(Involution(A),Involution(B))));
191}
192
193
194proc DecideZeroRowsEffectively (matrix A, module B)
195{
196  attrib(B,"isSB",1);
197  matrix M = reduce(A,B);
198  matrix T = lift(B,M-A);
199  list l = M,T;
200  return(l);
201}
202
203
204proc DecideZeroColumnsEffectively (matrix A, matrix B)
205{
206  list l = DecideZeroRowsEffectively(Involution(A),Involution(B));
207  matrix B = l[1];
208  matrix T = l[2];
209  l = Involution(B),Involution(T);
210  return(l);
211}
212
213
214proc RelativeSyzygiesGeneratorsOfRows (matrix M1, matrix M2)
215{
216  return(BasisOfRowModule(modulo(M1, M2)));
217}
218
219
220proc RelativeSyzygiesGeneratorsOfColumns (matrix M1, matrix M2)
221{
222  return(Involution(RelativeSyzygiesGeneratorsOfRows(Involution(M1),Involution(M2))));
223}
224
225
226proc PartiallyReducedBasisOfRowModule (matrix M)
227{
228  return(mstd(M)[2]);
229}
230
231
232proc PartiallyReducedBasisOfColumnModule (matrix M)
233{
234  return(Involution(PartiallyReducedBasisOfRowModule(Involution(M))));
235}
236
237
238proc IsZeroMatrix (matrix m)
239{
240  matrix z[nrows(m)][ncols(m)];
241  return(m==z);
242}
243
244
245proc IsIdentityMatrix (matrix m)
246{
247  return(m==unitmat(nrows(m)));
248}
249
250
251proc GetColumnIndependentUnitPositions (matrix M, list pos_list)
252{
253  int m = nrows(M);
254  int n = ncols(M);
255 
256  list rest;
257  intvec tmp = 1..m;
258  rest = tmp[1..m];
259  int r = m;
260  list rest2;
261  list pos;
262  int i; int k; int a; int s = 1; int s2;
263 
264  for (int j=1; j<=n; j++)
265  {
266    for (i=r; i>0; i--)
267    {
268      k = rest[i];
269      if (deg(M[k,j]) == 0) //IsUnit
270      {
271        rest2 = list();
272        s2 = 1;
273        pos[s] = list(j,k); s++;
274        for (a=1; a<=r; a++)
275        {
276          if (M[rest[a],j] == 0)
277          {
278            rest2[s2] = rest[a]; s2++;
279          }
280        }
281        rest = rest2;
282        r = size(rest);
283        break;
284      }
285    }
286  }
287  return(string(pos));
288}
289
290
291proc GetRowIndependentUnitPositions (matrix M, list pos_list)
292{
293  int m = nrows(M);
294  int n = ncols(M);
295 
296  list rest;
297  for (int o=n; o>=1; o--)
298  {
299    rest[o] = o;
300  }
301  int r = n;
302  list e;
303  list rest2;
304  list pos;
305  int j; int k; int a; int s = 1;
306 
307  for (int i=1; i<=m; i++)
308  {
309    for (j=1; j<=r; j++)
310    {
311      k = rest[r-j+1];
312      if (deg(M[i,k]) == 0) //IsUnit
313      {
314        rest2 = e;
315        pos[s] = list(i,k); s++;
316        for (a=1; a<=r; a++)
317        {
318          if (M[i,rest[a]] == 0)
319          {
320            rest2[size(rest2)+1] = rest[a];
321          }
322        }
323        rest = rest2;
324        r = size(rest);
325        break;
326      }
327    }
328  }
329  return(string(pos));
330}
331
332
333proc GetUnitPosition (matrix M, list pos_list)
334{
335  int m = nrows(M);
336  int n = ncols(M);
337  int r;
338  list rest;
339  for (int o=m; o>=1; o--)
340  {
341    rest[o] = o;
342  }
343  rest=Difference(rest,pos_list);
344  r=size(rest);
345  for (int j=1; j<=n; j++)
346  {
347    for (int i=1; i<=r; i++)
348    {
349      if (deg(M[rest[i],j]) == 0) //IsUnit
350      {
351        return(string(j,",",rest[i])); // this is not a mistake
352      }
353    }
354  }
355  return("fail");
356}
357
358
359proc GetCleanRowsPositions (matrix m, list l)
360{
361  list rows;
362  int s = 1;
363  for (int i=1;i<=size(l);i++)
364  {
365    for (int j=1;j<=ncols(m);j++)
366    {
367      if (m[l[i],j]==1)
368      {
369        rows[s] = j; s++;
370        break;
371      }
372    }
373  }
374  if (s==0)
375  {
376    return("[]"));
377  }
378  return(string(rows));
379}
380
381
382proc PrimaryDecomposition (matrix m)
383{
384  return(primdecSY(m))
385}
386
387
388proc CoefficientsOfUnreducedNumeratorOfWeightedHilbertPoincareSeries (module m,weights,degrees)
389{
390  module M = std(m);
391  attrib(M,"isHomog",degrees);
392  return(hilb(M,1,weights));
393}
394
395
396proc GetSparseListOfHomalgMatrixAsString (M)
397{
398  list l;int k;
399  k = 1;
400  for(int i=1; i<=ncols(M); i++){
401    for(int j=1; j<=nrows(M); j++){
402      def p=M[j,i]; // remark: matrices are saved transposed in Singular
403      if(p!=0){l[k]=list(i,j,p); k++;};
404    };
405  };
406  return(string(l));
407}
408
409
410proc IsMemberOfList (int i, list l)
411{
412  int k = size(l);
413 
414  for (int p=1; p<=k; p++)
415  {
416    if (l[p]==i)
417    {
418      return(1); // this is not a mistake
419    }
420  }
421  return(0);
422}
423
424
425proc Difference (list a, list b)
426{
427  list c;
428  int s=size(a);
429  int l = 1;
430 
431  for (int p=1; p<=s; p++)
432  {
433    if (IsMemberOfList(a[p],b)==0)
434    {
435      c[l] = a[p]; l++;
436    }
437  }
438  return(c);
439}
440
441
442proc CreateListListOfIntegers (degrees,m,n)
443{
444  list l;
445  for (int i=m; i>=1; i--)
446  {
447    l[i]=intvec(degrees[(i-1)*n+1..i*n]);
448  }
449  return(l);
450}
451
452
453proc GetColumnIndependentUnitPositions_Z (matrix M, list pos_list)
454{
455  int m = nrows(M);
456  int n = ncols(M);
457 
458  list rest;
459  for (int o=m; o>=1; o--)
460  {
461    rest[o] = o;
462  }
463  int r = m;
464  list e;
465  list rest2;
466  list pos;
467  int i; int k; int a; int s = 1;
468 
469  for (int j=1; j<=n; j++)
470  {
471    for (i=1; i<=r; i++)
472    {
473      k = rest[r-i+1];
474      if (M[k,j] == 1 || M[k,j] == -1) //IsUnit
475      {
476        rest2 = e;
477        pos[s] = list(j,k); s++;
478        for (a=1; a<=r; a++)
479        {
480          if (M[rest[a],j] == 0)
481          {
482            rest2[size(rest2)+1] = rest[a];
483          }
484        }
485        rest = rest2;
486        r = size(rest);
487        break;
488      }
489    }
490  }
491  return(string(pos));
492}
493
494
495proc GetRowIndependentUnitPositions_Z (matrix M, list pos_list)
496{
497  int m = nrows(M);
498  int n = ncols(M);
499 
500  list rest;
501  for (int o=n; o>=1; o--)
502  {
503    rest[o] = o;
504  }
505  int r = n;
506  list e;
507  list rest2;
508  list pos;
509  int j; int k; int a; int s = 1;
510 
511  for (int i=1; i<=m; i++)
512  {
513    for (j=1; j<=r; j++)
514    {
515      k = rest[r-j+1];
516      if (M[i,k] == 1 || M[i,k] == -1) //IsUnit
517      {
518        rest2 = e;
519        pos[s] = list(i,k); s++;
520        for (a=1; a<=r; a++)
521        {
522          if (M[i,rest[a]] == 0)
523          {
524            rest2[size(rest2)+1] = rest[a];
525          }
526        }
527        rest = rest2;
528        r = size(rest);
529        break;
530      }
531    }
532  }
533  return(string(pos));
534}
535
536
537proc GetUnitPosition_Z (matrix M, list pos_list)
538{
539  int m = nrows(M);
540  int n = ncols(M);
541  int r;
542  list rest;
543  for (int o=m; o>=1; o--)
544  {
545    rest[o] = o;
546  }
547  rest=Difference(rest,pos_list);
548  r=size(rest);
549  for (int j=1; j<=n; j++)
550  {
551    for (int i=1; i<=r; i++)
552    {
553      if (M[rest[i],j] == 1 || M[rest[i],j] == -1) //IsUnit
554      {
555        return(string(j,",",rest[i])); // this is not a mistake
556      }
557    }
558  }
559  return("fail");
560}
561
562
563proc SyzForHomalg (matrix M)
564{
565  return(syz(M));
566}
567
568
569proc ReducedSyzForHomalg (matrix M)
570{
571  return(matrix(nres(M,2)[2]));
572}
573
574
575if ( defined(superCommutative) == 1 ) // the new name of the SCA constructor
576{ proc superCommutative_ForHomalg = superCommutative; }
577else
578{
579  if ( defined(SuperCommutative) == 1 ) // the old name of the SCA constructor
580  { proc superCommutative_ForHomalg = SuperCommutative; }
581}
582
583
584
585// start: a workaround for a bug in the 64 bit versions of Singular 3-0-4
586if ( defined( basering ) != 0 )
587{
588  def homalg_variable_basering = basering;
589}
590ring r;
591if ( deg(0,(1,1,1)) > 0 )
592{ proc Deg (pol,weights)
593  {
594    if ( pol == 0 )
595    {
596      return(deg(0));
597    }
598    return(deg(pol,weights));
599  }
600}
601else
602{ proc Deg (pol,weights)
603  {
604    return(deg(pol,weights));
605  }
606}
607kill r;
608if ( defined( homalg_variable_basering ) != 0 )
609{
610  setring homalg_variable_basering;
611}
612// end: a workaround for a bug in the 64 bit versions of Singular 3-0-4
613
614
615
616// deg(lead()) instead of deg() below works around a bug
617proc homalg_Symbol (poly e)
618{
619  if(e==0) {return(e);}
620  poly l=lead(e);
621  int d=deg(l);
622  poly s=l;
623  poly r=e-l;
624  l=lead(r);
625  while(deg(l)==d)
626  {
627    s=s+l;
628    r=r-l;
629    l=lead(r);
630  }
631  return(s);
632}
633
634
635proc MatrixOfSymbols_workaround (matrix m)
636{
637  int i; int j; poly e;
638  int r=nrows(m);
639  int c=ncols(m);
640  matrix n[r][c]=0;
641  // lead(e) below works around a bug
642  for(i=1;i<=r;i++)
643  {
644    for(j=1;j<=c;j++)
645    {
646      e=m[i,j];
647      if(e!=0)
648      { n[i,j]=homalg_Symbol(e); }
649    }
650  }
651  return(n);
652}
653
654ring homalg_variable_1 = 2,dummy_variable,dp;
655number homalg_variable_2 = 0;
656number homalg_variable_3 = 1;
657number homalg_variable_4 = -1;
658setring homalg_variable_1;
659short=0;option(redTail);
660
661proc Involution (matrix m)
662{
663  return(transpose(m));
664}
665
666ring homalg_variable_5 = (2,Z(4)),dummy_variable,dp;
667setring homalg_variable_5;
668short=0;option(redTail);
669(1)*Z(4)^2+(1)*Z(4)^1+(1)*Z(4)^0;
670minpoly=(1)*Z(4)^2+(1)*Z(4)^1+(1)*Z(4)^0;
671number homalg_variable_6 = 0;
672number homalg_variable_7 = 1;
673number homalg_variable_8 = -1;
674setring homalg_variable_5;
675short=0;option(redTail);
676
677proc Involution (matrix m)
678{
679  return(transpose(m));
680}
681
682poly homalg_variable_9 = Z(4);
683ring homalg_variable_10 = (2,Z(4)),(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20),dp;
684setring homalg_variable_10;
685short=0;option(redTail);
686number homalg_variable_11 = 0;
687number homalg_variable_12 = 1;
688number homalg_variable_13 = -1;
689minpoly=Z(4)^2+Z(4)+1;
690poly homalg_variable_14 = x1;
691poly homalg_variable_15 = x2;
692poly homalg_variable_16 = x3;
693poly homalg_variable_17 = x4;
694poly homalg_variable_18 = x5;
695poly homalg_variable_19 = x6;
696poly homalg_variable_20 = x7;
697poly homalg_variable_21 = x8;
698poly homalg_variable_22 = x9;
699poly homalg_variable_23 = x10;
700poly homalg_variable_24 = x11;
701poly homalg_variable_25 = x12;
702poly homalg_variable_26 = x13;
703poly homalg_variable_27 = x14;
704poly homalg_variable_28 = x15;
705poly homalg_variable_29 = x16;
706poly homalg_variable_30 = x17;
707poly homalg_variable_31 = x18;
708poly homalg_variable_32 = x19;
709poly homalg_variable_33 = x20;
710homalg_variable_14;
711homalg_variable_15;
712homalg_variable_16;
713homalg_variable_17;
714homalg_variable_18;
715homalg_variable_19;
716homalg_variable_20;
717homalg_variable_21;
718homalg_variable_22;
719homalg_variable_23;
720homalg_variable_24;
721homalg_variable_25;
722homalg_variable_26;
723homalg_variable_27;
724homalg_variable_28;
725homalg_variable_29;
726homalg_variable_30;
727homalg_variable_31;
728homalg_variable_32;
729homalg_variable_33;
730poly homalg_variable_34 = Z(4);
731setring homalg_variable_10;
732short=0;option(redTail);
733
734proc Involution (matrix m)
735{
736  return(transpose(m));
737}
738
739string s=read("r: homalg_file_1_PID_38382-singular");execute( "matrix homalg_variable_35[1][20] = " + s + ";" );homalg_variable_35 = transpose(homalg_variable_35);
740string s=read("r: homalg_file_2_PID_38382-singular");execute( "matrix homalg_variable_36[12][20] = " + s + ";" );homalg_variable_36 = transpose(homalg_variable_36);
741IsZeroMatrix(homalg_variable_36);
742ZeroRows(homalg_variable_36);
743matrix homalg_variable_37 = homalg_variable_36;
744def homalg_variable_38 = homalg_variable_37[1,1];
745def homalg_variable_39 = homalg_variable_37[2,1];
746def homalg_variable_40 = homalg_variable_37[3,1];
747def homalg_variable_41 = homalg_variable_37[4,1];
748def homalg_variable_42 = homalg_variable_37[5,1];
749def homalg_variable_43 = homalg_variable_37[6,1];
750def homalg_variable_44 = homalg_variable_37[7,1];
751def homalg_variable_45 = homalg_variable_37[8,1];
752def homalg_variable_46 = homalg_variable_37[9,1];
753def homalg_variable_47 = homalg_variable_37[10,1];
754def homalg_variable_48 = homalg_variable_37[11,1];
755def homalg_variable_49 = homalg_variable_37[12,1];
756def homalg_variable_50 = homalg_variable_37[13,1];
757def homalg_variable_51 = homalg_variable_37[14,1];
758def homalg_variable_52 = homalg_variable_37[15,1];
759def homalg_variable_53 = homalg_variable_37[16,1];
760def homalg_variable_54 = homalg_variable_37[17,1];
761def homalg_variable_55 = homalg_variable_37[18,1];
762def homalg_variable_56 = homalg_variable_37[19,1];
763def homalg_variable_57 = homalg_variable_37[20,1];
764def homalg_variable_58 = homalg_variable_37[1,2];
765def homalg_variable_59 = homalg_variable_37[2,2];
766def homalg_variable_60 = homalg_variable_37[3,2];
767def homalg_variable_61 = homalg_variable_37[4,2];
768def homalg_variable_62 = homalg_variable_37[5,2];
769def homalg_variable_63 = homalg_variable_37[6,2];
770def homalg_variable_64 = homalg_variable_37[7,2];
771def homalg_variable_65 = homalg_variable_37[8,2];
772def homalg_variable_66 = homalg_variable_37[9,2];
773def homalg_variable_67 = homalg_variable_37[10,2];
774def homalg_variable_68 = homalg_variable_37[11,2];
775def homalg_variable_69 = homalg_variable_37[12,2];
776def homalg_variable_70 = homalg_variable_37[13,2];
777def homalg_variable_71 = homalg_variable_37[14,2];
778def homalg_variable_72 = homalg_variable_37[15,2];
779def homalg_variable_73 = homalg_variable_37[16,2];
780def homalg_variable_74 = homalg_variable_37[17,2];
781def homalg_variable_75 = homalg_variable_37[18,2];
782def homalg_variable_76 = homalg_variable_37[19,2];
783def homalg_variable_77 = homalg_variable_37[20,2];
784def homalg_variable_78 = homalg_variable_37[1,3];
785def homalg_variable_79 = homalg_variable_37[2,3];
786def homalg_variable_80 = homalg_variable_37[3,3];
787def homalg_variable_81 = homalg_variable_37[4,3];
788def homalg_variable_82 = homalg_variable_37[5,3];
789def homalg_variable_83 = homalg_variable_37[6,3];
790def homalg_variable_84 = homalg_variable_37[7,3];
791def homalg_variable_85 = homalg_variable_37[8,3];
792def homalg_variable_86 = homalg_variable_37[9,3];
793def homalg_variable_87 = homalg_variable_37[10,3];
794def homalg_variable_88 = homalg_variable_37[11,3];
795def homalg_variable_89 = homalg_variable_37[12,3];
796def homalg_variable_90 = homalg_variable_37[13,3];
797def homalg_variable_91 = homalg_variable_37[14,3];
798def homalg_variable_92 = homalg_variable_37[15,3];
799def homalg_variable_93 = homalg_variable_37[16,3];
800def homalg_variable_94 = homalg_variable_37[17,3];
801def homalg_variable_95 = homalg_variable_37[18,3];
802def homalg_variable_96 = homalg_variable_37[19,3];
803def homalg_variable_97 = homalg_variable_37[20,3];
804def homalg_variable_98 = homalg_variable_37[1,4];
805def homalg_variable_99 = homalg_variable_37[2,4];
806def homalg_variable_100 = homalg_variable_37[3,4];
807def homalg_variable_101 = homalg_variable_37[4,4];
808def homalg_variable_102 = homalg_variable_37[5,4];
809def homalg_variable_103 = homalg_variable_37[6,4];
810def homalg_variable_104 = homalg_variable_37[7,4];
811def homalg_variable_105 = homalg_variable_37[8,4];
812def homalg_variable_106 = homalg_variable_37[9,4];
813def homalg_variable_107 = homalg_variable_37[10,4];
814def homalg_variable_108 = homalg_variable_37[11,4];
815def homalg_variable_109 = homalg_variable_37[12,4];
816def homalg_variable_110 = homalg_variable_37[13,4];
817def homalg_variable_111 = homalg_variable_37[14,4];
818def homalg_variable_112 = homalg_variable_37[15,4];
819def homalg_variable_113 = homalg_variable_37[16,4];
820def homalg_variable_114 = homalg_variable_37[17,4];
821def homalg_variable_115 = homalg_variable_37[18,4];
822def homalg_variable_116 = homalg_variable_37[19,4];
823def homalg_variable_117 = homalg_variable_37[20,4];
824def homalg_variable_118 = homalg_variable_37[1,5];
825def homalg_variable_119 = homalg_variable_37[2,5];
826def homalg_variable_120 = homalg_variable_37[3,5];
827def homalg_variable_121 = homalg_variable_37[4,5];
828def homalg_variable_122 = homalg_variable_37[5,5];
829def homalg_variable_123 = homalg_variable_37[6,5];
830def homalg_variable_124 = homalg_variable_37[7,5];
831def homalg_variable_125 = homalg_variable_37[8,5];
832def homalg_variable_126 = homalg_variable_37[9,5];
833def homalg_variable_127 = homalg_variable_37[10,5];
834def homalg_variable_128 = homalg_variable_37[11,5];
835def homalg_variable_129 = homalg_variable_37[12,5];
836def homalg_variable_130 = homalg_variable_37[13,5];
837def homalg_variable_131 = homalg_variable_37[14,5];
838def homalg_variable_132 = homalg_variable_37[15,5];
839def homalg_variable_133 = homalg_variable_37[16,5];
840def homalg_variable_134 = homalg_variable_37[17,5];
841def homalg_variable_135 = homalg_variable_37[18,5];
842def homalg_variable_136 = homalg_variable_37[19,5];
843def homalg_variable_137 = homalg_variable_37[20,5];
844def homalg_variable_138 = homalg_variable_37[1,6];
845def homalg_variable_139 = homalg_variable_37[2,6];
846def homalg_variable_140 = homalg_variable_37[3,6];
847def homalg_variable_141 = homalg_variable_37[4,6];
848def homalg_variable_142 = homalg_variable_37[5,6];
849def homalg_variable_143 = homalg_variable_37[6,6];
850def homalg_variable_144 = homalg_variable_37[7,6];
851def homalg_variable_145 = homalg_variable_37[8,6];
852def homalg_variable_146 = homalg_variable_37[9,6];
853def homalg_variable_147 = homalg_variable_37[10,6];
854def homalg_variable_148 = homalg_variable_37[11,6];
855def homalg_variable_149 = homalg_variable_37[12,6];
856def homalg_variable_150 = homalg_variable_37[13,6];
857def homalg_variable_151 = homalg_variable_37[14,6];
858def homalg_variable_152 = homalg_variable_37[15,6];
859def homalg_variable_153 = homalg_variable_37[16,6];
860def homalg_variable_154 = homalg_variable_37[17,6];
861def homalg_variable_155 = homalg_variable_37[18,6];
862def homalg_variable_156 = homalg_variable_37[19,6];
863def homalg_variable_157 = homalg_variable_37[20,6];
864def homalg_variable_158 = homalg_variable_37[1,7];
865def homalg_variable_159 = homalg_variable_37[2,7];
866def homalg_variable_160 = homalg_variable_37[3,7];
867def homalg_variable_161 = homalg_variable_37[4,7];
868def homalg_variable_162 = homalg_variable_37[5,7];
869def homalg_variable_163 = homalg_variable_37[6,7];
870def homalg_variable_164 = homalg_variable_37[7,7];
871def homalg_variable_165 = homalg_variable_37[8,7];
872def homalg_variable_166 = homalg_variable_37[9,7];
873def homalg_variable_167 = homalg_variable_37[10,7];
874def homalg_variable_168 = homalg_variable_37[11,7];
875def homalg_variable_169 = homalg_variable_37[12,7];
876def homalg_variable_170 = homalg_variable_37[13,7];
877def homalg_variable_171 = homalg_variable_37[14,7];
878def homalg_variable_172 = homalg_variable_37[15,7];
879def homalg_variable_173 = homalg_variable_37[16,7];
880def homalg_variable_174 = homalg_variable_37[17,7];
881def homalg_variable_175 = homalg_variable_37[18,7];
882def homalg_variable_176 = homalg_variable_37[19,7];
883def homalg_variable_177 = homalg_variable_37[20,7];
884def homalg_variable_178 = homalg_variable_37[1,8];
885def homalg_variable_179 = homalg_variable_37[2,8];
886def homalg_variable_180 = homalg_variable_37[3,8];
887def homalg_variable_181 = homalg_variable_37[4,8];
888def homalg_variable_182 = homalg_variable_37[5,8];
889def homalg_variable_183 = homalg_variable_37[6,8];
890def homalg_variable_184 = homalg_variable_37[7,8];
891def homalg_variable_185 = homalg_variable_37[8,8];
892def homalg_variable_186 = homalg_variable_37[9,8];
893def homalg_variable_187 = homalg_variable_37[10,8];
894def homalg_variable_188 = homalg_variable_37[11,8];
895def homalg_variable_189 = homalg_variable_37[12,8];
896def homalg_variable_190 = homalg_variable_37[13,8];
897def homalg_variable_191 = homalg_variable_37[14,8];
898def homalg_variable_192 = homalg_variable_37[15,8];
899def homalg_variable_193 = homalg_variable_37[16,8];
900def homalg_variable_194 = homalg_variable_37[17,8];
901def homalg_variable_195 = homalg_variable_37[18,8];
902def homalg_variable_196 = homalg_variable_37[19,8];
903def homalg_variable_197 = homalg_variable_37[20,8];
904def homalg_variable_198 = homalg_variable_37[1,9];
905def homalg_variable_199 = homalg_variable_37[2,9];
906def homalg_variable_200 = homalg_variable_37[3,9];
907def homalg_variable_201 = homalg_variable_37[4,9];
908def homalg_variable_202 = homalg_variable_37[5,9];
909def homalg_variable_203 = homalg_variable_37[6,9];
910def homalg_variable_204 = homalg_variable_37[7,9];
911def homalg_variable_205 = homalg_variable_37[8,9];
912def homalg_variable_206 = homalg_variable_37[9,9];
913def homalg_variable_207 = homalg_variable_37[10,9];
914def homalg_variable_208 = homalg_variable_37[11,9];
915def homalg_variable_209 = homalg_variable_37[12,9];
916def homalg_variable_210 = homalg_variable_37[13,9];
917def homalg_variable_211 = homalg_variable_37[14,9];
918def homalg_variable_212 = homalg_variable_37[15,9];
919def homalg_variable_213 = homalg_variable_37[16,9];
920def homalg_variable_214 = homalg_variable_37[17,9];
921def homalg_variable_215 = homalg_variable_37[18,9];
922def homalg_variable_216 = homalg_variable_37[19,9];
923def homalg_variable_217 = homalg_variable_37[20,9];
924def homalg_variable_218 = homalg_variable_37[1,10];
925def homalg_variable_219 = homalg_variable_37[2,10];
926def homalg_variable_220 = homalg_variable_37[3,10];
927def homalg_variable_221 = homalg_variable_37[4,10];
928def homalg_variable_222 = homalg_variable_37[5,10];
929def homalg_variable_223 = homalg_variable_37[6,10];
930def homalg_variable_224 = homalg_variable_37[7,10];
931def homalg_variable_225 = homalg_variable_37[8,10];
932def homalg_variable_226 = homalg_variable_37[9,10];
933def homalg_variable_227 = homalg_variable_37[10,10];
934def homalg_variable_228 = homalg_variable_37[11,10];
935def homalg_variable_229 = homalg_variable_37[12,10];
936def homalg_variable_230 = homalg_variable_37[13,10];
937def homalg_variable_231 = homalg_variable_37[14,10];
938def homalg_variable_232 = homalg_variable_37[15,10];
939def homalg_variable_233 = homalg_variable_37[16,10];
940def homalg_variable_234 = homalg_variable_37[17,10];
941def homalg_variable_235 = homalg_variable_37[18,10];
942def homalg_variable_236 = homalg_variable_37[19,10];
943def homalg_variable_237 = homalg_variable_37[20,10];
944def homalg_variable_238 = homalg_variable_37[1,11];
945def homalg_variable_239 = homalg_variable_37[2,11];
946def homalg_variable_240 = homalg_variable_37[3,11];
947def homalg_variable_241 = homalg_variable_37[4,11];
948def homalg_variable_242 = homalg_variable_37[5,11];
949def homalg_variable_243 = homalg_variable_37[6,11];
950def homalg_variable_244 = homalg_variable_37[7,11];
951def homalg_variable_245 = homalg_variable_37[8,11];
952def homalg_variable_246 = homalg_variable_37[9,11];
953def homalg_variable_247 = homalg_variable_37[10,11];
954def homalg_variable_248 = homalg_variable_37[11,11];
955def homalg_variable_249 = homalg_variable_37[12,11];
956def homalg_variable_250 = homalg_variable_37[13,11];
957def homalg_variable_251 = homalg_variable_37[14,11];
958def homalg_variable_252 = homalg_variable_37[15,11];
959def homalg_variable_253 = homalg_variable_37[16,11];
960def homalg_variable_254 = homalg_variable_37[17,11];
961def homalg_variable_255 = homalg_variable_37[18,11];
962def homalg_variable_256 = homalg_variable_37[19,11];
963def homalg_variable_257 = homalg_variable_37[20,11];
964def homalg_variable_258 = homalg_variable_37[1,12];
965def homalg_variable_259 = homalg_variable_37[2,12];
966def homalg_variable_260 = homalg_variable_37[3,12];
967def homalg_variable_261 = homalg_variable_37[4,12];
968def homalg_variable_262 = homalg_variable_37[5,12];
969def homalg_variable_263 = homalg_variable_37[6,12];
970def homalg_variable_264 = homalg_variable_37[7,12];
971def homalg_variable_265 = homalg_variable_37[8,12];
972def homalg_variable_266 = homalg_variable_37[9,12];
973def homalg_variable_267 = homalg_variable_37[10,12];
974def homalg_variable_268 = homalg_variable_37[11,12];
975def homalg_variable_269 = homalg_variable_37[12,12];
976def homalg_variable_270 = homalg_variable_37[13,12];
977def homalg_variable_271 = homalg_variable_37[14,12];
978def homalg_variable_272 = homalg_variable_37[15,12];
979def homalg_variable_273 = homalg_variable_37[16,12];
980def homalg_variable_274 = homalg_variable_37[17,12];
981def homalg_variable_275 = homalg_variable_37[18,12];
982def homalg_variable_276 = homalg_variable_37[19,12];
983def homalg_variable_277 = homalg_variable_37[20,12];
984matrix homalg_variable_278 = SyzygiesGeneratorsOfRows(homalg_variable_36);
985IsZeroMatrix(homalg_variable_278);
986matrix homalg_variable_279 = BasisOfRowModule(homalg_variable_36);
987ncols(homalg_variable_279);
988homalg_variable_279==homalg_variable_36;
989matrix homalg_variable_281[20][20] = unitmat(20);
990matrix homalg_variable_280 = DecideZeroRows(homalg_variable_281,homalg_variable_279);
991IsZeroMatrix(homalg_variable_280);
992ZeroRows(homalg_variable_280);
993IsIdentityMatrix(homalg_variable_280);
994matrix homalg_variable_282 = homalg_variable_280;
995def homalg_variable_283 = homalg_variable_282[1,1];
996def homalg_variable_284 = homalg_variable_282[2,1];
997def homalg_variable_285 = homalg_variable_282[3,1];
998def homalg_variable_286 = homalg_variable_282[4,1];
999def homalg_variable_287 = homalg_variable_282[5,1];
1000def homalg_variable_288 = homalg_variable_282[6,1];
1001def homalg_variable_289 = homalg_variable_282[7,1];
1002def homalg_variable_290 = homalg_variable_282[8,1];
1003def homalg_variable_291 = homalg_variable_282[9,1];
1004def homalg_variable_292 = homalg_variable_282[10,1];
1005def homalg_variable_293 = homalg_variable_282[11,1];
1006def homalg_variable_294 = homalg_variable_282[12,1];
1007def homalg_variable_295 = homalg_variable_282[13,1];
1008def homalg_variable_296 = homalg_variable_282[14,1];
1009def homalg_variable_297 = homalg_variable_282[15,1];
1010def homalg_variable_298 = homalg_variable_282[16,1];
1011def homalg_variable_299 = homalg_variable_282[17,1];
1012def homalg_variable_300 = homalg_variable_282[18,1];
1013def homalg_variable_301 = homalg_variable_282[19,1];
1014def homalg_variable_302 = homalg_variable_282[20,1];
1015def homalg_variable_303 = homalg_variable_282[1,2];
1016def homalg_variable_304 = homalg_variable_282[2,2];
1017def homalg_variable_305 = homalg_variable_282[3,2];
1018def homalg_variable_306 = homalg_variable_282[4,2];
1019def homalg_variable_307 = homalg_variable_282[5,2];
1020def homalg_variable_308 = homalg_variable_282[6,2];
1021def homalg_variable_309 = homalg_variable_282[7,2];
1022def homalg_variable_310 = homalg_variable_282[8,2];
1023def homalg_variable_311 = homalg_variable_282[9,2];
1024def homalg_variable_312 = homalg_variable_282[10,2];
1025def homalg_variable_313 = homalg_variable_282[11,2];
1026def homalg_variable_314 = homalg_variable_282[12,2];
1027def homalg_variable_315 = homalg_variable_282[13,2];
1028def homalg_variable_316 = homalg_variable_282[14,2];
1029def homalg_variable_317 = homalg_variable_282[15,2];
1030def homalg_variable_318 = homalg_variable_282[16,2];
1031def homalg_variable_319 = homalg_variable_282[17,2];
1032def homalg_variable_320 = homalg_variable_282[18,2];
1033def homalg_variable_321 = homalg_variable_282[19,2];
1034def homalg_variable_322 = homalg_variable_282[20,2];
1035def homalg_variable_323 = homalg_variable_282[1,3];
1036def homalg_variable_324 = homalg_variable_282[2,3];
1037def homalg_variable_325 = homalg_variable_282[3,3];
1038def homalg_variable_326 = homalg_variable_282[4,3];
1039def homalg_variable_327 = homalg_variable_282[5,3];
1040def homalg_variable_328 = homalg_variable_282[6,3];
1041def homalg_variable_329 = homalg_variable_282[7,3];
1042def homalg_variable_330 = homalg_variable_282[8,3];
1043def homalg_variable_331 = homalg_variable_282[9,3];
1044def homalg_variable_332 = homalg_variable_282[10,3];
1045def homalg_variable_333 = homalg_variable_282[11,3];
1046def homalg_variable_334 = homalg_variable_282[12,3];
1047def homalg_variable_335 = homalg_variable_282[13,3];
1048def homalg_variable_336 = homalg_variable_282[14,3];
1049def homalg_variable_337 = homalg_variable_282[15,3];
1050def homalg_variable_338 = homalg_variable_282[16,3];
1051def homalg_variable_339 = homalg_variable_282[17,3];
1052def homalg_variable_340 = homalg_variable_282[18,3];
1053def homalg_variable_341 = homalg_variable_282[19,3];
1054def homalg_variable_342 = homalg_variable_282[20,3];
1055def homalg_variable_343 = homalg_variable_282[1,4];
1056def homalg_variable_344 = homalg_variable_282[2,4];
1057def homalg_variable_345 = homalg_variable_282[3,4];
1058def homalg_variable_346 = homalg_variable_282[4,4];
1059def homalg_variable_347 = homalg_variable_282[5,4];
1060def homalg_variable_348 = homalg_variable_282[6,4];
1061def homalg_variable_349 = homalg_variable_282[7,4];
1062def homalg_variable_350 = homalg_variable_282[8,4];
1063def homalg_variable_351 = homalg_variable_282[9,4];
1064def homalg_variable_352 = homalg_variable_282[10,4];
1065def homalg_variable_353 = homalg_variable_282[11,4];
1066def homalg_variable_354 = homalg_variable_282[12,4];
1067def homalg_variable_355 = homalg_variable_282[13,4];
1068def homalg_variable_356 = homalg_variable_282[14,4];
1069def homalg_variable_357 = homalg_variable_282[15,4];
1070def homalg_variable_358 = homalg_variable_282[16,4];
1071def homalg_variable_359 = homalg_variable_282[17,4];
1072def homalg_variable_360 = homalg_variable_282[18,4];
1073def homalg_variable_361 = homalg_variable_282[19,4];
1074def homalg_variable_362 = homalg_variable_282[20,4];
1075def homalg_variable_363 = homalg_variable_282[1,5];
1076def homalg_variable_364 = homalg_variable_282[2,5];
1077def homalg_variable_365 = homalg_variable_282[3,5];
1078def homalg_variable_366 = homalg_variable_282[4,5];
1079def homalg_variable_367 = homalg_variable_282[5,5];
1080def homalg_variable_368 = homalg_variable_282[6,5];
1081def homalg_variable_369 = homalg_variable_282[7,5];
1082def homalg_variable_370 = homalg_variable_282[8,5];
1083def homalg_variable_371 = homalg_variable_282[9,5];
1084def homalg_variable_372 = homalg_variable_282[10,5];
1085def homalg_variable_373 = homalg_variable_282[11,5];
1086def homalg_variable_374 = homalg_variable_282[12,5];
1087def homalg_variable_375 = homalg_variable_282[13,5];
1088def homalg_variable_376 = homalg_variable_282[14,5];
1089def homalg_variable_377 = homalg_variable_282[15,5];
1090def homalg_variable_378 = homalg_variable_282[16,5];
1091def homalg_variable_379 = homalg_variable_282[17,5];
1092def homalg_variable_380 = homalg_variable_282[18,5];
1093def homalg_variable_381 = homalg_variable_282[19,5];
1094def homalg_variable_382 = homalg_variable_282[20,5];
1095def homalg_variable_383 = homalg_variable_282[1,6];
1096def homalg_variable_384 = homalg_variable_282[2,6];
1097def homalg_variable_385 = homalg_variable_282[3,6];
1098def homalg_variable_386 = homalg_variable_282[4,6];
1099def homalg_variable_387 = homalg_variable_282[5,6];
1100def homalg_variable_388 = homalg_variable_282[6,6];
1101def homalg_variable_389 = homalg_variable_282[7,6];
1102def homalg_variable_390 = homalg_variable_282[8,6];
1103def homalg_variable_391 = homalg_variable_282[9,6];
1104def homalg_variable_392 = homalg_variable_282[10,6];
1105def homalg_variable_393 = homalg_variable_282[11,6];
1106def homalg_variable_394 = homalg_variable_282[12,6];
1107def homalg_variable_395 = homalg_variable_282[13,6];
1108def homalg_variable_396 = homalg_variable_282[14,6];
1109def homalg_variable_397 = homalg_variable_282[15,6];
1110def homalg_variable_398 = homalg_variable_282[16,6];
1111def homalg_variable_399 = homalg_variable_282[17,6];
1112def homalg_variable_400 = homalg_variable_282[18,6];
1113def homalg_variable_401 = homalg_variable_282[19,6];
1114def homalg_variable_402 = homalg_variable_282[20,6];
1115def homalg_variable_403 = homalg_variable_282[1,7];
1116def homalg_variable_404 = homalg_variable_282[2,7];
1117def homalg_variable_405 = homalg_variable_282[3,7];
1118def homalg_variable_406 = homalg_variable_282[4,7];
1119def homalg_variable_407 = homalg_variable_282[5,7];
1120def homalg_variable_408 = homalg_variable_282[6,7];
1121def homalg_variable_409 = homalg_variable_282[7,7];
1122def homalg_variable_410 = homalg_variable_282[8,7];
1123def homalg_variable_411 = homalg_variable_282[9,7];
1124def homalg_variable_412 = homalg_variable_282[10,7];
1125def homalg_variable_413 = homalg_variable_282[11,7];
1126def homalg_variable_414 = homalg_variable_282[12,7];
1127def homalg_variable_415 = homalg_variable_282[13,7];
1128def homalg_variable_416 = homalg_variable_282[14,7];
1129def homalg_variable_417 = homalg_variable_282[15,7];
1130def homalg_variable_418 = homalg_variable_282[16,7];
1131def homalg_variable_419 = homalg_variable_282[17,7];
1132def homalg_variable_420 = homalg_variable_282[18,7];
1133def homalg_variable_421 = homalg_variable_282[19,7];
1134def homalg_variable_422 = homalg_variable_282[20,7];
1135def homalg_variable_423 = homalg_variable_282[1,8];
1136def homalg_variable_424 = homalg_variable_282[2,8];
1137def homalg_variable_425 = homalg_variable_282[3,8];
1138def homalg_variable_426 = homalg_variable_282[4,8];
1139def homalg_variable_427 = homalg_variable_282[5,8];
1140def homalg_variable_428 = homalg_variable_282[6,8];
1141def homalg_variable_429 = homalg_variable_282[7,8];
1142def homalg_variable_430 = homalg_variable_282[8,8];
1143def homalg_variable_431 = homalg_variable_282[9,8];
1144def homalg_variable_432 = homalg_variable_282[10,8];
1145def homalg_variable_433 = homalg_variable_282[11,8];
1146def homalg_variable_434 = homalg_variable_282[12,8];
1147def homalg_variable_435 = homalg_variable_282[13,8];
1148def homalg_variable_436 = homalg_variable_282[14,8];
1149def homalg_variable_437 = homalg_variable_282[15,8];
1150def homalg_variable_438 = homalg_variable_282[16,8];
1151def homalg_variable_439 = homalg_variable_282[17,8];
1152def homalg_variable_440 = homalg_variable_282[18,8];
1153def homalg_variable_441 = homalg_variable_282[19,8];
1154def homalg_variable_442 = homalg_variable_282[20,8];
1155def homalg_variable_443 = homalg_variable_282[1,9];
1156def homalg_variable_444 = homalg_variable_282[2,9];
1157def homalg_variable_445 = homalg_variable_282[3,9];
1158def homalg_variable_446 = homalg_variable_282[4,9];
1159def homalg_variable_447 = homalg_variable_282[5,9];
1160def homalg_variable_448 = homalg_variable_282[6,9];
1161def homalg_variable_449 = homalg_variable_282[7,9];
1162def homalg_variable_450 = homalg_variable_282[8,9];
1163def homalg_variable_451 = homalg_variable_282[9,9];
1164def homalg_variable_452 = homalg_variable_282[10,9];
1165def homalg_variable_453 = homalg_variable_282[11,9];
1166def homalg_variable_454 = homalg_variable_282[12,9];
1167def homalg_variable_455 = homalg_variable_282[13,9];
1168def homalg_variable_456 = homalg_variable_282[14,9];
1169def homalg_variable_457 = homalg_variable_282[15,9];
1170def homalg_variable_458 = homalg_variable_282[16,9];
1171def homalg_variable_459 = homalg_variable_282[17,9];
1172def homalg_variable_460 = homalg_variable_282[18,9];
1173def homalg_variable_461 = homalg_variable_282[19,9];
1174def homalg_variable_462 = homalg_variable_282[20,9];
1175def homalg_variable_463 = homalg_variable_282[1,10];
1176def homalg_variable_464 = homalg_variable_282[2,10];
1177def homalg_variable_465 = homalg_variable_282[3,10];
1178def homalg_variable_466 = homalg_variable_282[4,10];
1179def homalg_variable_467 = homalg_variable_282[5,10];
1180def homalg_variable_468 = homalg_variable_282[6,10];
1181def homalg_variable_469 = homalg_variable_282[7,10];
1182def homalg_variable_470 = homalg_variable_282[8,10];
1183def homalg_variable_471 = homalg_variable_282[9,10];
1184def homalg_variable_472 = homalg_variable_282[10,10];
1185def homalg_variable_473 = homalg_variable_282[11,10];
1186def homalg_variable_474 = homalg_variable_282[12,10];
1187def homalg_variable_475 = homalg_variable_282[13,10];
1188def homalg_variable_476 = homalg_variable_282[14,10];
1189def homalg_variable_477 = homalg_variable_282[15,10];
1190def homalg_variable_478 = homalg_variable_282[16,10];
1191def homalg_variable_479 = homalg_variable_282[17,10];
1192def homalg_variable_480 = homalg_variable_282[18,10];
1193def homalg_variable_481 = homalg_variable_282[19,10];
1194def homalg_variable_482 = homalg_variable_282[20,10];
1195def homalg_variable_483 = homalg_variable_282[1,11];
1196def homalg_variable_484 = homalg_variable_282[2,11];
1197def homalg_variable_485 = homalg_variable_282[3,11];
1198def homalg_variable_486 = homalg_variable_282[4,11];
1199def homalg_variable_487 = homalg_variable_282[5,11];
1200def homalg_variable_488 = homalg_variable_282[6,11];
1201def homalg_variable_489 = homalg_variable_282[7,11];
1202def homalg_variable_490 = homalg_variable_282[8,11];
1203def homalg_variable_491 = homalg_variable_282[9,11];
1204def homalg_variable_492 = homalg_variable_282[10,11];
1205def homalg_variable_493 = homalg_variable_282[11,11];
1206def homalg_variable_494 = homalg_variable_282[12,11];
1207def homalg_variable_495 = homalg_variable_282[13,11];
1208def homalg_variable_496 = homalg_variable_282[14,11];
1209def homalg_variable_497 = homalg_variable_282[15,11];
1210def homalg_variable_498 = homalg_variable_282[16,11];
1211def homalg_variable_499 = homalg_variable_282[17,11];
1212def homalg_variable_500 = homalg_variable_282[18,11];
1213kill homalg_variable_278;
1214def homalg_variable_501 = homalg_variable_282[19,11];
1215def homalg_variable_502 = homalg_variable_282[20,11];
1216def homalg_variable_503 = homalg_variable_282[1,12];
1217def homalg_variable_504 = homalg_variable_282[2,12];
1218def homalg_variable_505 = homalg_variable_282[3,12];
1219def homalg_variable_506 = homalg_variable_282[4,12];
1220def homalg_variable_507 = homalg_variable_282[5,12];
1221def homalg_variable_508 = homalg_variable_282[6,12];
1222def homalg_variable_509 = homalg_variable_282[7,12];
1223def homalg_variable_510 = homalg_variable_282[8,12];
1224def homalg_variable_511 = homalg_variable_282[9,12];
1225def homalg_variable_512 = homalg_variable_282[10,12];
1226def homalg_variable_513 = homalg_variable_282[11,12];
1227def homalg_variable_514 = homalg_variable_282[12,12];
1228def homalg_variable_515 = homalg_variable_282[13,12];
1229def homalg_variable_516 = homalg_variable_282[14,12];
1230def homalg_variable_517 = homalg_variable_282[15,12];
1231def homalg_variable_518 = homalg_variable_282[16,12];
1232def homalg_variable_519 = homalg_variable_282[17,12];
1233def homalg_variable_520 = homalg_variable_282[18,12];
1234def homalg_variable_521 = homalg_variable_282[19,12];
1235def homalg_variable_522 = homalg_variable_282[20,12];
1236def homalg_variable_523 = homalg_variable_282[1,13];
1237def homalg_variable_524 = homalg_variable_282[2,13];
1238def homalg_variable_525 = homalg_variable_282[3,13];
1239def homalg_variable_526 = homalg_variable_282[4,13];
1240def homalg_variable_527 = homalg_variable_282[5,13];
1241def homalg_variable_528 = homalg_variable_282[6,13];
1242def homalg_variable_529 = homalg_variable_282[7,13];
1243def homalg_variable_530 = homalg_variable_282[8,13];
1244def homalg_variable_531 = homalg_variable_282[9,13];
1245def homalg_variable_532 = homalg_variable_282[10,13];
1246def homalg_variable_533 = homalg_variable_282[11,13];
1247def homalg_variable_534 = homalg_variable_282[12,13];
1248def homalg_variable_535 = homalg_variable_282[13,13];
1249def homalg_variable_536 = homalg_variable_282[14,13];
1250def homalg_variable_537 = homalg_variable_282[15,13];
1251def homalg_variable_538 = homalg_variable_282[16,13];
1252def homalg_variable_539 = homalg_variable_282[17,13];
1253def homalg_variable_540 = homalg_variable_282[18,13];
1254def homalg_variable_541 = homalg_variable_282[19,13];
1255def homalg_variable_542 = homalg_variable_282[20,13];
1256def homalg_variable_543 = homalg_variable_282[1,14];
1257def homalg_variable_544 = homalg_variable_282[2,14];
1258def homalg_variable_545 = homalg_variable_282[3,14];
1259def homalg_variable_546 = homalg_variable_282[4,14];
1260def homalg_variable_547 = homalg_variable_282[5,14];
1261def homalg_variable_548 = homalg_variable_282[6,14];
1262def homalg_variable_549 = homalg_variable_282[7,14];
1263def homalg_variable_550 = homalg_variable_282[8,14];
1264def homalg_variable_551 = homalg_variable_282[9,14];
1265def homalg_variable_552 = homalg_variable_282[10,14];
1266def homalg_variable_553 = homalg_variable_282[11,14];
1267def homalg_variable_554 = homalg_variable_282[12,14];
1268def homalg_variable_555 = homalg_variable_282[13,14];
1269def homalg_variable_556 = homalg_variable_282[14,14];
1270def homalg_variable_557 = homalg_variable_282[15,14];
1271def homalg_variable_558 = homalg_variable_282[16,14];
1272def homalg_variable_559 = homalg_variable_282[17,14];
1273def homalg_variable_560 = homalg_variable_282[18,14];
1274def homalg_variable_561 = homalg_variable_282[19,14];
1275def homalg_variable_562 = homalg_variable_282[20,14];
1276def homalg_variable_563 = homalg_variable_282[1,15];
1277def homalg_variable_564 = homalg_variable_282[2,15];
1278def homalg_variable_565 = homalg_variable_282[3,15];
1279def homalg_variable_566 = homalg_variable_282[4,15];
1280def homalg_variable_567 = homalg_variable_282[5,15];
1281def homalg_variable_568 = homalg_variable_282[6,15];
1282def homalg_variable_569 = homalg_variable_282[7,15];
1283def homalg_variable_570 = homalg_variable_282[8,15];
1284def homalg_variable_571 = homalg_variable_282[9,15];
1285def homalg_variable_572 = homalg_variable_282[10,15];
1286def homalg_variable_573 = homalg_variable_282[11,15];
1287def homalg_variable_574 = homalg_variable_282[12,15];
1288def homalg_variable_575 = homalg_variable_282[13,15];
1289def homalg_variable_576 = homalg_variable_282[14,15];
1290def homalg_variable_577 = homalg_variable_282[15,15];
1291def homalg_variable_578 = homalg_variable_282[16,15];
1292def homalg_variable_579 = homalg_variable_282[17,15];
1293def homalg_variable_580 = homalg_variable_282[18,15];
1294def homalg_variable_581 = homalg_variable_282[19,15];
1295def homalg_variable_582 = homalg_variable_282[20,15];
1296def homalg_variable_583 = homalg_variable_282[1,16];
1297def homalg_variable_584 = homalg_variable_282[2,16];
1298def homalg_variable_585 = homalg_variable_282[3,16];
1299def homalg_variable_586 = homalg_variable_282[4,16];
1300def homalg_variable_587 = homalg_variable_282[5,16];
1301def homalg_variable_588 = homalg_variable_282[6,16];
1302def homalg_variable_589 = homalg_variable_282[7,16];
1303def homalg_variable_590 = homalg_variable_282[8,16];
1304def homalg_variable_591 = homalg_variable_282[9,16];
1305def homalg_variable_592 = homalg_variable_282[10,16];
1306def homalg_variable_593 = homalg_variable_282[11,16];
1307def homalg_variable_594 = homalg_variable_282[12,16];
1308def homalg_variable_595 = homalg_variable_282[13,16];
1309def homalg_variable_596 = homalg_variable_282[14,16];
1310def homalg_variable_597 = homalg_variable_282[15,16];
1311def homalg_variable_598 = homalg_variable_282[16,16];
1312def homalg_variable_599 = homalg_variable_282[17,16];
1313def homalg_variable_600 = homalg_variable_282[18,16];
1314def homalg_variable_601 = homalg_variable_282[19,16];
1315def homalg_variable_602 = homalg_variable_282[20,16];
1316def homalg_variable_603 = homalg_variable_282[1,17];
1317def homalg_variable_604 = homalg_variable_282[2,17];
1318def homalg_variable_605 = homalg_variable_282[3,17];
1319def homalg_variable_606 = homalg_variable_282[4,17];
1320def homalg_variable_607 = homalg_variable_282[5,17];
1321def homalg_variable_608 = homalg_variable_282[6,17];
1322def homalg_variable_609 = homalg_variable_282[7,17];
1323def homalg_variable_610 = homalg_variable_282[8,17];
1324def homalg_variable_611 = homalg_variable_282[9,17];
1325def homalg_variable_612 = homalg_variable_282[10,17];
1326def homalg_variable_613 = homalg_variable_282[11,17];
1327def homalg_variable_614 = homalg_variable_282[12,17];
1328def homalg_variable_615 = homalg_variable_282[13,17];
1329def homalg_variable_616 = homalg_variable_282[14,17];
1330def homalg_variable_617 = homalg_variable_282[15,17];
1331def homalg_variable_618 = homalg_variable_282[16,17];
1332def homalg_variable_619 = homalg_variable_282[17,17];
1333def homalg_variable_620 = homalg_variable_282[18,17];
1334def homalg_variable_621 = homalg_variable_282[19,17];
1335def homalg_variable_622 = homalg_variable_282[20,17];
1336def homalg_variable_623 = homalg_variable_282[1,18];
1337def homalg_variable_624 = homalg_variable_282[2,18];
1338def homalg_variable_625 = homalg_variable_282[3,18];
1339def homalg_variable_626 = homalg_variable_282[4,18];
1340def homalg_variable_627 = homalg_variable_282[5,18];
1341def homalg_variable_628 = homalg_variable_282[6,18];
1342def homalg_variable_629 = homalg_variable_282[7,18];
1343def homalg_variable_630 = homalg_variable_282[8,18];
1344def homalg_variable_631 = homalg_variable_282[9,18];
1345def homalg_variable_632 = homalg_variable_282[10,18];
1346def homalg_variable_633 = homalg_variable_282[11,18];
1347def homalg_variable_634 = homalg_variable_282[12,18];
1348def homalg_variable_635 = homalg_variable_282[13,18];
1349def homalg_variable_636 = homalg_variable_282[14,18];
1350def homalg_variable_637 = homalg_variable_282[15,18];
1351def homalg_variable_638 = homalg_variable_282[16,18];
1352def homalg_variable_639 = homalg_variable_282[17,18];
1353def homalg_variable_640 = homalg_variable_282[18,18];
1354def homalg_variable_641 = homalg_variable_282[19,18];
1355def homalg_variable_642 = homalg_variable_282[20,18];
1356def homalg_variable_643 = homalg_variable_282[1,19];
1357def homalg_variable_644 = homalg_variable_282[2,19];
1358def homalg_variable_645 = homalg_variable_282[3,19];
1359def homalg_variable_646 = homalg_variable_282[4,19];
1360def homalg_variable_647 = homalg_variable_282[5,19];
1361def homalg_variable_648 = homalg_variable_282[6,19];
1362def homalg_variable_649 = homalg_variable_282[7,19];
1363def homalg_variable_650 = homalg_variable_282[8,19];
1364def homalg_variable_651 = homalg_variable_282[9,19];
1365def homalg_variable_652 = homalg_variable_282[10,19];
1366def homalg_variable_653 = homalg_variable_282[11,19];
1367def homalg_variable_654 = homalg_variable_282[12,19];
1368def homalg_variable_655 = homalg_variable_282[13,19];
1369def homalg_variable_656 = homalg_variable_282[14,19];
1370def homalg_variable_657 = homalg_variable_282[15,19];
1371def homalg_variable_658 = homalg_variable_282[16,19];
1372def homalg_variable_659 = homalg_variable_282[17,19];
1373def homalg_variable_660 = homalg_variable_282[18,19];
1374def homalg_variable_661 = homalg_variable_282[19,19];
1375def homalg_variable_662 = homalg_variable_282[20,19];
1376def homalg_variable_663 = homalg_variable_282[1,20];
1377def homalg_variable_664 = homalg_variable_282[2,20];
1378def homalg_variable_665 = homalg_variable_282[3,20];
1379def homalg_variable_666 = homalg_variable_282[4,20];
1380def homalg_variable_667 = homalg_variable_282[5,20];
1381def homalg_variable_668 = homalg_variable_282[6,20];
1382def homalg_variable_669 = homalg_variable_282[7,20];
1383def homalg_variable_670 = homalg_variable_282[8,20];
1384def homalg_variable_671 = homalg_variable_282[9,20];
1385def homalg_variable_672 = homalg_variable_282[10,20];
1386def homalg_variable_673 = homalg_variable_282[11,20];
1387def homalg_variable_674 = homalg_variable_282[12,20];
1388def homalg_variable_675 = homalg_variable_282[13,20];
1389def homalg_variable_676 = homalg_variable_282[14,20];
1390def homalg_variable_677 = homalg_variable_282[15,20];
1391def homalg_variable_678 = homalg_variable_282[16,20];
1392def homalg_variable_679 = homalg_variable_282[17,20];
1393def homalg_variable_680 = homalg_variable_282[18,20];
1394def homalg_variable_681 = homalg_variable_282[19,20];
1395def homalg_variable_682 = homalg_variable_282[20,20];
1396matrix homalg_variable_683 = RelativeSyzygiesGeneratorsOfRows(homalg_variable_280,homalg_variable_279);
1397IsZeroMatrix(homalg_variable_683);
1398ncols(homalg_variable_683);
1399GetUnitPosition(homalg_variable_683, list ( 0 ));
1400def homalg_variable_684 = homalg_variable_683[1,1];
1401homalg_variable_684==1;
1402matrix homalg_variable_685 = submat(homalg_variable_683,1..20,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 ));
1403matrix homalg_variable_686 = submat(homalg_variable_685,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 ),1..11);
1404matrix homalg_variable_687 = submat(homalg_variable_683,1..20,intvec( 1 ));
1405matrix homalg_variable_688 = submat(homalg_variable_687,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 ),1..1);
1406matrix homalg_variable_689 = submat(homalg_variable_683,intvec( 1 ),1..12);
1407matrix homalg_variable_690 = submat(homalg_variable_689,1..1,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 ));
1408matrix homalg_variable_691 = homalg_variable_688*homalg_variable_690;
1409matrix homalg_variable_692 = homalg_variable_686-homalg_variable_691;
1410GetUnitPosition(homalg_variable_692, list ( 0 ));
1411def homalg_variable_693 = homalg_variable_692[1,1];
1412homalg_variable_693==1;
1413matrix homalg_variable_694 = submat(homalg_variable_692,1..19,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 ));
1414matrix homalg_variable_695 = submat(homalg_variable_694,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 ),1..10);
1415matrix homalg_variable_696 = submat(homalg_variable_692,1..19,intvec( 1 ));
1416matrix homalg_variable_697 = submat(homalg_variable_696,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 ),1..1);
1417matrix homalg_variable_698 = submat(homalg_variable_692,intvec( 1 ),1..11);
1418matrix homalg_variable_699 = submat(homalg_variable_698,1..1,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 ));
1419matrix homalg_variable_700 = homalg_variable_697*homalg_variable_699;
1420matrix homalg_variable_701 = homalg_variable_695-homalg_variable_700;
1421GetUnitPosition(homalg_variable_701, list ( 0 ));
1422def homalg_variable_702 = homalg_variable_701[1,1];
1423homalg_variable_702==1;
1424matrix homalg_variable_703 = submat(homalg_variable_701,1..18,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10 ));
1425matrix homalg_variable_704 = submat(homalg_variable_703,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 ),1..9);
1426matrix homalg_variable_705 = submat(homalg_variable_701,1..18,intvec( 1 ));
1427matrix homalg_variable_706 = submat(homalg_variable_705,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 ),1..1);
1428matrix homalg_variable_707 = submat(homalg_variable_701,intvec( 1 ),1..10);
1429matrix homalg_variable_708 = submat(homalg_variable_707,1..1,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10 ));
1430matrix homalg_variable_709 = homalg_variable_706*homalg_variable_708;
1431matrix homalg_variable_710 = homalg_variable_704-homalg_variable_709;
1432GetUnitPosition(homalg_variable_710, list ( 0 ));
1433def homalg_variable_711 = homalg_variable_710[1,1];
1434homalg_variable_711==1;
1435matrix homalg_variable_712 = submat(homalg_variable_710,1..17,intvec( 2, 3, 4, 5, 6, 7, 8, 9 ));
1436matrix homalg_variable_713 = submat(homalg_variable_712,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17 ),1..8);
1437matrix homalg_variable_714 = submat(homalg_variable_710,1..17,intvec( 1 ));
1438matrix homalg_variable_715 = submat(homalg_variable_714,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17 ),1..1);
1439matrix homalg_variable_716 = submat(homalg_variable_710,intvec( 1 ),1..9);
1440matrix homalg_variable_717 = submat(homalg_variable_716,1..1,intvec( 2, 3, 4, 5, 6, 7, 8, 9 ));
1441matrix homalg_variable_718 = homalg_variable_715*homalg_variable_717;
1442matrix homalg_variable_719 = homalg_variable_713-homalg_variable_718;
1443GetUnitPosition(homalg_variable_719, list ( 0 ));
1444def homalg_variable_720 = homalg_variable_719[1,1];
1445homalg_variable_720==1;
1446matrix homalg_variable_721 = submat(homalg_variable_719,1..16,intvec( 2, 3, 4, 5, 6, 7, 8 ));
1447matrix homalg_variable_722 = submat(homalg_variable_721,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 ),1..7);
1448matrix homalg_variable_723 = submat(homalg_variable_719,1..16,intvec( 1 ));
1449matrix homalg_variable_724 = submat(homalg_variable_723,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 ),1..1);
1450matrix homalg_variable_725 = submat(homalg_variable_719,intvec( 1 ),1..8);
1451matrix homalg_variable_726 = submat(homalg_variable_725,1..1,intvec( 2, 3, 4, 5, 6, 7, 8 ));
1452matrix homalg_variable_727 = homalg_variable_724*homalg_variable_726;
1453matrix homalg_variable_728 = homalg_variable_722-homalg_variable_727;
1454GetUnitPosition(homalg_variable_728, list ( 0 ));
1455def homalg_variable_729 = homalg_variable_728[1,1];
1456homalg_variable_729==1;
1457matrix homalg_variable_730 = submat(homalg_variable_728,1..15,intvec( 2, 3, 4, 5, 6, 7 ));
1458matrix homalg_variable_731 = submat(homalg_variable_730,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 ),1..6);
1459matrix homalg_variable_732 = submat(homalg_variable_728,1..15,intvec( 1 ));
1460matrix homalg_variable_733 = submat(homalg_variable_732,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 ),1..1);
1461matrix homalg_variable_734 = submat(homalg_variable_728,intvec( 1 ),1..7);
1462matrix homalg_variable_735 = submat(homalg_variable_734,1..1,intvec( 2, 3, 4, 5, 6, 7 ));
1463matrix homalg_variable_736 = homalg_variable_733*homalg_variable_735;
1464matrix homalg_variable_737 = homalg_variable_731-homalg_variable_736;
1465GetUnitPosition(homalg_variable_737, list ( 0 ));
1466def homalg_variable_738 = homalg_variable_737[1,1];
1467homalg_variable_738==1;
1468matrix homalg_variable_739 = submat(homalg_variable_737,1..14,intvec( 2, 3, 4, 5, 6 ));
1469matrix homalg_variable_740 = submat(homalg_variable_739,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 ),1..5);
1470matrix homalg_variable_741 = submat(homalg_variable_737,1..14,intvec( 1 ));
1471matrix homalg_variable_742 = submat(homalg_variable_741,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 ),1..1);
1472matrix homalg_variable_743 = submat(homalg_variable_737,intvec( 1 ),1..6);
1473matrix homalg_variable_744 = submat(homalg_variable_743,1..1,intvec( 2, 3, 4, 5, 6 ));
1474matrix homalg_variable_745 = homalg_variable_742*homalg_variable_744;
1475matrix homalg_variable_746 = homalg_variable_740-homalg_variable_745;
1476GetUnitPosition(homalg_variable_746, list ( 0 ));
1477def homalg_variable_747 = homalg_variable_746[1,1];
1478homalg_variable_747==1;
1479matrix homalg_variable_748 = submat(homalg_variable_746,1..13,intvec( 2, 3, 4, 5 ));
1480matrix homalg_variable_749 = submat(homalg_variable_748,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 ),1..4);
1481matrix homalg_variable_750 = submat(homalg_variable_746,1..13,intvec( 1 ));
1482matrix homalg_variable_751 = submat(homalg_variable_750,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 ),1..1);
1483matrix homalg_variable_752 = submat(homalg_variable_746,intvec( 1 ),1..5);
1484matrix homalg_variable_753 = submat(homalg_variable_752,1..1,intvec( 2, 3, 4, 5 ));
1485matrix homalg_variable_754 = homalg_variable_751*homalg_variable_753;
1486matrix homalg_variable_755 = homalg_variable_749-homalg_variable_754;
1487GetUnitPosition(homalg_variable_755, list ( 0 ));
1488def homalg_variable_756 = homalg_variable_755[1,1];
1489homalg_variable_756==1;
1490matrix homalg_variable_757 = submat(homalg_variable_755,1..12,intvec( 2, 3, 4 ));
1491matrix homalg_variable_758 = submat(homalg_variable_757,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 ),1..3);
1492matrix homalg_variable_759 = submat(homalg_variable_755,1..12,intvec( 1 ));
1493matrix homalg_variable_760 = submat(homalg_variable_759,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 ),1..1);
1494matrix homalg_variable_761 = submat(homalg_variable_755,intvec( 1 ),1..4);
1495matrix homalg_variable_762 = submat(homalg_variable_761,1..1,intvec( 2, 3, 4 ));
1496matrix homalg_variable_763 = homalg_variable_760*homalg_variable_762;
1497matrix homalg_variable_764 = homalg_variable_758-homalg_variable_763;
1498GetUnitPosition(homalg_variable_764, list ( 0 ));
1499def homalg_variable_765 = homalg_variable_764[1,1];
1500homalg_variable_765==1;
1501matrix homalg_variable_766 = submat(homalg_variable_764,1..11,intvec( 2, 3 ));
1502matrix homalg_variable_767 = submat(homalg_variable_766,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 ),1..2);
1503matrix homalg_variable_768 = submat(homalg_variable_764,1..11,intvec( 1 ));
1504matrix homalg_variable_769 = submat(homalg_variable_768,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 ),1..1);
1505matrix homalg_variable_770 = submat(homalg_variable_764,intvec( 1 ),1..3);
1506matrix homalg_variable_771 = submat(homalg_variable_770,1..1,intvec( 2, 3 ));
1507matrix homalg_variable_772 = homalg_variable_769*homalg_variable_771;
1508matrix homalg_variable_773 = homalg_variable_767-homalg_variable_772;
1509GetUnitPosition(homalg_variable_773, list ( 0 ));
1510def homalg_variable_774 = homalg_variable_773[1,1];
1511homalg_variable_774==1;
1512matrix homalg_variable_775 = submat(homalg_variable_773,1..10,intvec( 2 ));
1513matrix homalg_variable_776 = submat(homalg_variable_775,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10 ),1..1);
1514matrix homalg_variable_777 = submat(homalg_variable_773,1..10,intvec( 1 ));
1515matrix homalg_variable_778 = submat(homalg_variable_777,intvec( 2, 3, 4, 5, 6, 7, 8, 9, 10 ),1..1);
1516matrix homalg_variable_779 = submat(homalg_variable_773,intvec( 1 ),1..2);
1517matrix homalg_variable_780 = submat(homalg_variable_779,1..1,intvec( 2 ));
1518matrix homalg_variable_781 = homalg_variable_778*homalg_variable_780;
1519matrix homalg_variable_782 = homalg_variable_776-homalg_variable_781;
1520GetUnitPosition(homalg_variable_782, list ( 0 ));
1521def homalg_variable_783 = homalg_variable_782[2,1];
1522homalg_variable_783==1;
1523matrix homalg_variable_784 = submat(homalg_variable_281,1..20,intvec( 12, 14, 15, 16, 17, 18, 19, 20 ));
1524matrix homalg_variable_785 = homalg_variable_784;
1525def homalg_variable_786 = homalg_variable_785[1,1];
1526def homalg_variable_787 = homalg_variable_785[2,1];
1527def homalg_variable_788 = homalg_variable_785[3,1];
1528def homalg_variable_789 = homalg_variable_785[4,1];
1529def homalg_variable_790 = homalg_variable_785[5,1];
1530def homalg_variable_791 = homalg_variable_785[6,1];
1531def homalg_variable_792 = homalg_variable_785[7,1];
1532def homalg_variable_793 = homalg_variable_785[8,1];
1533def homalg_variable_794 = homalg_variable_785[9,1];
1534def homalg_variable_795 = homalg_variable_785[10,1];
1535def homalg_variable_796 = homalg_variable_785[11,1];
1536def homalg_variable_797 = homalg_variable_785[12,1];
1537def homalg_variable_798 = homalg_variable_785[13,1];
1538def homalg_variable_799 = homalg_variable_785[14,1];
1539def homalg_variable_800 = homalg_variable_785[15,1];
1540def homalg_variable_801 = homalg_variable_785[16,1];
1541def homalg_variable_802 = homalg_variable_785[17,1];
1542def homalg_variable_803 = homalg_variable_785[18,1];
1543def homalg_variable_804 = homalg_variable_785[19,1];
1544def homalg_variable_805 = homalg_variable_785[20,1];
1545def homalg_variable_806 = homalg_variable_785[1,2];
1546def homalg_variable_807 = homalg_variable_785[2,2];
1547def homalg_variable_808 = homalg_variable_785[3,2];
1548def homalg_variable_809 = homalg_variable_785[4,2];
1549def homalg_variable_810 = homalg_variable_785[5,2];
1550def homalg_variable_811 = homalg_variable_785[6,2];
1551def homalg_variable_812 = homalg_variable_785[7,2];
1552def homalg_variable_813 = homalg_variable_785[8,2];
1553def homalg_variable_814 = homalg_variable_785[9,2];
1554def homalg_variable_815 = homalg_variable_785[10,2];
1555def homalg_variable_816 = homalg_variable_785[11,2];
1556def homalg_variable_817 = homalg_variable_785[12,2];
1557def homalg_variable_818 = homalg_variable_785[13,2];
1558def homalg_variable_819 = homalg_variable_785[14,2];
1559def homalg_variable_820 = homalg_variable_785[15,2];
1560def homalg_variable_821 = homalg_variable_785[16,2];
1561def homalg_variable_822 = homalg_variable_785[17,2];
1562def homalg_variable_823 = homalg_variable_785[18,2];
1563def homalg_variable_824 = homalg_variable_785[19,2];
1564def homalg_variable_825 = homalg_variable_785[20,2];
1565def homalg_variable_826 = homalg_variable_785[1,3];
1566def homalg_variable_827 = homalg_variable_785[2,3];
1567def homalg_variable_828 = homalg_variable_785[3,3];
1568def homalg_variable_829 = homalg_variable_785[4,3];
1569def homalg_variable_830 = homalg_variable_785[5,3];
1570def homalg_variable_831 = homalg_variable_785[6,3];
1571def homalg_variable_832 = homalg_variable_785[7,3];
1572def homalg_variable_833 = homalg_variable_785[8,3];
1573def homalg_variable_834 = homalg_variable_785[9,3];
1574def homalg_variable_835 = homalg_variable_785[10,3];
1575def homalg_variable_836 = homalg_variable_785[11,3];
1576def homalg_variable_837 = homalg_variable_785[12,3];
1577def homalg_variable_838 = homalg_variable_785[13,3];
1578def homalg_variable_839 = homalg_variable_785[14,3];
1579def homalg_variable_840 = homalg_variable_785[15,3];
1580def homalg_variable_841 = homalg_variable_785[16,3];
1581def homalg_variable_842 = homalg_variable_785[17,3];
1582def homalg_variable_843 = homalg_variable_785[18,3];
1583def homalg_variable_844 = homalg_variable_785[19,3];
1584def homalg_variable_845 = homalg_variable_785[20,3];
1585def homalg_variable_846 = homalg_variable_785[1,4];
1586def homalg_variable_847 = homalg_variable_785[2,4];
1587def homalg_variable_848 = homalg_variable_785[3,4];
1588def homalg_variable_849 = homalg_variable_785[4,4];
1589def homalg_variable_850 = homalg_variable_785[5,4];
1590def homalg_variable_851 = homalg_variable_785[6,4];
1591def homalg_variable_852 = homalg_variable_785[7,4];
1592def homalg_variable_853 = homalg_variable_785[8,4];
1593def homalg_variable_854 = homalg_variable_785[9,4];
1594def homalg_variable_855 = homalg_variable_785[10,4];
1595def homalg_variable_856 = homalg_variable_785[11,4];
1596def homalg_variable_857 = homalg_variable_785[12,4];
1597def homalg_variable_858 = homalg_variable_785[13,4];
1598def homalg_variable_859 = homalg_variable_785[14,4];
1599def homalg_variable_860 = homalg_variable_785[15,4];
1600def homalg_variable_861 = homalg_variable_785[16,4];
1601def homalg_variable_862 = homalg_variable_785[17,4];
1602def homalg_variable_863 = homalg_variable_785[18,4];
1603def homalg_variable_864 = homalg_variable_785[19,4];
1604def homalg_variable_865 = homalg_variable_785[20,4];
1605def homalg_variable_866 = homalg_variable_785[1,5];
1606def homalg_variable_867 = homalg_variable_785[2,5];
1607def homalg_variable_868 = homalg_variable_785[3,5];
1608def homalg_variable_869 = homalg_variable_785[4,5];
1609def homalg_variable_870 = homalg_variable_785[5,5];
1610def homalg_variable_871 = homalg_variable_785[6,5];
1611def homalg_variable_872 = homalg_variable_785[7,5];
1612def homalg_variable_873 = homalg_variable_785[8,5];
1613def homalg_variable_874 = homalg_variable_785[9,5];
1614def homalg_variable_875 = homalg_variable_785[10,5];
1615def homalg_variable_876 = homalg_variable_785[11,5];
1616def homalg_variable_877 = homalg_variable_785[12,5];
1617def homalg_variable_878 = homalg_variable_785[13,5];
1618def homalg_variable_879 = homalg_variable_785[14,5];
1619def homalg_variable_880 = homalg_variable_785[15,5];
1620def homalg_variable_881 = homalg_variable_785[16,5];
1621def homalg_variable_882 = homalg_variable_785[17,5];
1622def homalg_variable_883 = homalg_variable_785[18,5];
1623def homalg_variable_884 = homalg_variable_785[19,5];
1624def homalg_variable_885 = homalg_variable_785[20,5];
1625def homalg_variable_886 = homalg_variable_785[1,6];
1626def homalg_variable_887 = homalg_variable_785[2,6];
1627def homalg_variable_888 = homalg_variable_785[3,6];
1628def homalg_variable_889 = homalg_variable_785[4,6];
1629def homalg_variable_890 = homalg_variable_785[5,6];
1630def homalg_variable_891 = homalg_variable_785[6,6];
1631def homalg_variable_892 = homalg_variable_785[7,6];
1632def homalg_variable_893 = homalg_variable_785[8,6];
1633def homalg_variable_894 = homalg_variable_785[9,6];
1634def homalg_variable_895 = homalg_variable_785[10,6];
1635def homalg_variable_896 = homalg_variable_785[11,6];
1636def homalg_variable_897 = homalg_variable_785[12,6];
1637def homalg_variable_898 = homalg_variable_785[13,6];
1638def homalg_variable_899 = homalg_variable_785[14,6];
1639def homalg_variable_900 = homalg_variable_785[15,6];
1640def homalg_variable_901 = homalg_variable_785[16,6];
1641def homalg_variable_902 = homalg_variable_785[17,6];
1642def homalg_variable_903 = homalg_variable_785[18,6];
1643def homalg_variable_904 = homalg_variable_785[19,6];
1644def homalg_variable_905 = homalg_variable_785[20,6];
1645def homalg_variable_906 = homalg_variable_785[1,7];
1646def homalg_variable_907 = homalg_variable_785[2,7];
1647def homalg_variable_908 = homalg_variable_785[3,7];
1648def homalg_variable_909 = homalg_variable_785[4,7];
1649def homalg_variable_910 = homalg_variable_785[5,7];
1650def homalg_variable_911 = homalg_variable_785[6,7];
1651def homalg_variable_912 = homalg_variable_785[7,7];
1652def homalg_variable_913 = homalg_variable_785[8,7];
1653def homalg_variable_914 = homalg_variable_785[9,7];
1654def homalg_variable_915 = homalg_variable_785[10,7];
1655def homalg_variable_916 = homalg_variable_785[11,7];
1656def homalg_variable_917 = homalg_variable_785[12,7];
1657def homalg_variable_918 = homalg_variable_785[13,7];
1658def homalg_variable_919 = homalg_variable_785[14,7];
1659def homalg_variable_920 = homalg_variable_785[15,7];
1660def homalg_variable_921 = homalg_variable_785[16,7];
1661def homalg_variable_922 = homalg_variable_785[17,7];
1662def homalg_variable_923 = homalg_variable_785[18,7];
1663def homalg_variable_924 = homalg_variable_785[19,7];
1664def homalg_variable_925 = homalg_variable_785[20,7];
1665def homalg_variable_926 = homalg_variable_785[1,8];
1666def homalg_variable_927 = homalg_variable_785[2,8];
1667def homalg_variable_928 = homalg_variable_785[3,8];
1668def homalg_variable_929 = homalg_variable_785[4,8];
1669def homalg_variable_930 = homalg_variable_785[5,8];
1670def homalg_variable_931 = homalg_variable_785[6,8];
1671def homalg_variable_932 = homalg_variable_785[7,8];
1672def homalg_variable_933 = homalg_variable_785[8,8];
1673def homalg_variable_934 = homalg_variable_785[9,8];
1674def homalg_variable_935 = homalg_variable_785[10,8];
1675def homalg_variable_936 = homalg_variable_785[11,8];
1676def homalg_variable_937 = homalg_variable_785[12,8];
1677def homalg_variable_938 = homalg_variable_785[13,8];
1678def homalg_variable_939 = homalg_variable_785[14,8];
1679def homalg_variable_940 = homalg_variable_785[15,8];
1680def homalg_variable_941 = homalg_variable_785[16,8];
1681def homalg_variable_942 = homalg_variable_785[17,8];
1682def homalg_variable_943 = homalg_variable_785[18,8];
1683def homalg_variable_944 = homalg_variable_785[19,8];
1684def homalg_variable_945 = homalg_variable_785[20,8];
1685matrix m[20][1];string s = "[";for(int i=1;i<=12;i++){m = homalg_variable_36[1..20,i]; if(i!=1){s=s+",";};s=s+"["+string(m)+"]";};s=s+"]";write("w: /var/tmp/homalg_variable_36_PID_38385",s);
1686matrix homalg_variable_947 = submat(homalg_variable_280,1..20,intvec( 12, 14, 15, 16, 17, 18, 19, 20 ));
1687matrix homalg_variable_946 = DecideZeroRows(homalg_variable_947,homalg_variable_279);
1688IsZeroMatrix(homalg_variable_946);
1689matrix homalg_variable_948[20][20] = concat(homalg_variable_946,homalg_variable_279);
1690list homalg_variable_l=BasisOfRowsCoeff(homalg_variable_948); matrix homalg_variable_949 = homalg_variable_l[1]; matrix homalg_variable_950 = homalg_variable_l[2];
1691ncols(homalg_variable_949);
1692matrix homalg_variable_951 = DecideZeroRows(homalg_variable_35,homalg_variable_279);
1693IsZeroMatrix(homalg_variable_951);
1694IsIdentityMatrix(homalg_variable_949);
1695nrows(homalg_variable_950);
1696matrix homalg_variable_952[1][20] = x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20;
1697homalg_variable_952 = transpose(homalg_variable_952);
1698matrix homalg_variable_953[20][20] = concat(homalg_variable_946,homalg_variable_279);
1699list homalg_variable_l=BasisOfRowsCoeff(homalg_variable_953); matrix homalg_variable_954 = homalg_variable_l[1]; matrix homalg_variable_955 = homalg_variable_l[2];
1700ncols(homalg_variable_954);
1701matrix homalg_variable_956 = DecideZeroRows(homalg_variable_281,homalg_variable_279);
1702IsZeroMatrix(homalg_variable_956);
1703IsIdentityMatrix(homalg_variable_954);
1704nrows(homalg_variable_955);
1705homalg_variable_14==0;
1706homalg_variable_15==0;
1707homalg_variable_16==0;
1708homalg_variable_17==0;
1709homalg_variable_18==0;
1710homalg_variable_19==0;
1711homalg_variable_20==0;
1712homalg_variable_21==0;
1713homalg_variable_22==0;
1714homalg_variable_23==0;
1715homalg_variable_24==0;
1716homalg_variable_25==0;
1717homalg_variable_26==0;
1718homalg_variable_27==0;
1719homalg_variable_28==0;
1720homalg_variable_29==0;
1721homalg_variable_30==0;
1722homalg_variable_31==0;
1723homalg_variable_32==0;
1724homalg_variable_33==0;
1725def homalg_variable_957 = poly(homalg_variable_12)/homalg_variable_12;