Changeset 503bfc1 in git for libpolys/polys/clapsing.cc

Timestamp:

Feb 12, 2016, 5:19:36 PM (8 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')

Children:

e4a2bf6197b85736b4c554df2e065b083afef0cb

Parents:

52d6d84fe6b1bb74dbb7d6ec97baa8b77b0cccfc

Message:

system("HNF"...) for bigintmat

File:

• libpolys/polys/clapsing.cc (modified) (1 diff)

libpolys/polys/clapsing.cc

@@ -1704 +1704 @@
 }
 
+bigintmat* singntl_HNF(bigintmat*  b)
+{
+  int r=b->rows();
+  if (r!=b->cols())
+  {
+    Werror("HNF of %d x %d matrix",r,b->cols());
+    return NULL;
+  }
+  setCharacteristic( 0 );
+  CFMatrix M(r,r);
+  int i,j;
+  for(i=r;i>0;i--)
+  {
+    for(j=r;j>0;j--)
+    {
+      M(i,j)=n_convSingNFactoryN(BIMATELEM(*b,i,j),FALSE,b->basecoeffs());
+    }
+  }
+  CFMatrix *MM=cf_HNF(M);
+  bigintmat *mm=bimCopy(b);
+  for(i=r;i>0;i--)
+  {
+    for(j=r;j>0;j--)
+    {
+      BIMATELEM(*mm,i,j)=n_convFactoryNSingN((*MM)(i,j),b->basecoeffs());
+    }
+  }
+  delete MM;
+  return mm;
+}
+
 matrix singntl_LLL(matrix  m, const ring s )
 {