Changeset 273a2f in git for libpolys/polys/flintconv.cc

Timestamp:

Dec 2, 2021, 3:56:49 PM (2 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

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

Children:

18bb792edcc6d3f65af0111f93ed61118123ffd1

Parents:

c31a71c38969321615ad9607978783ca23b0edbb

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2021-12-02 15:56:49+01:00

git-committer:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2021-12-02 15:59:38+01:00

Message:

experimental: fq_nmod_mat_rref

File:

• libpolys/polys/flintconv.cc (modified) (2 diffs)

libpolys/polys/flintconv.cc

@@ -14 +14 @@
 #ifdef HAVE_FLINT
 #if __FLINT_RELEASE >= 20500
+
 #include "coeffs/coeffs.h"
 #include "coeffs/longrat.h"
@@ -281 +282 @@
 }
 
+void convSingPFlintnmod_poly_t(nmod_poly_t result, const poly p, const ring r)
+{
+  // assume univariate
+  nmod_poly_init2 (result,rChar(r),p_Deg(p,r)); 
+  poly h=p;
+  while(h!=NULL)
+  {
+    nmod_poly_set_coeff_ui(result,p_GetExp(h,1,r),n_Int(pGetCoeff(h),r->cf));
+    pIter(h);
+  }
+}
+
+void convSingMFlintFq_nmod_mat(matrix m, fq_nmod_mat_t M, const fq_nmod_ctx_t fq_con, const ring r)
+{
+  fq_nmod_mat_init (M, (long)MATROWS(m), (long) MATCOLS(m), fq_con);
+  int i,j;
+  for(i=MATROWS(m);i>0;i--)
+  {
+    for(j=MATCOLS(m);j>0;j--)
+    {
+      convSingPFlintnmod_poly_t (M->rows[i-1]+j-1, MATELEM(m,i,j),r);
+    }
+  }
+}
+
+poly convFlintFq_nmodSingP(const fq_nmod_t Fp, const fq_nmod_ctx_t ctx, const ring r)
+{
+  poly p=NULL;
+  poly h;
+  for (int i= 0; i < nmod_poly_length (Fp); i++)
+  {
+    ulong coeff= nmod_poly_get_coeff_ui (Fp, i);
+    if (coeff != 0)
+    h=p_NSet(n_Init(coeff,r->cf),r);
+    p_SetExp(h,1,i,r);
+    p_Setm(h,r);
+    p=p_Add_q(p,h,r);
+  }
+  return p;
+}
+
+matrix convFlintFq_nmod_matSingM(fq_nmod_mat_t m, const fq_nmod_ctx_t fq_con, const ring r)
+{
+  matrix M=mpNew(fq_nmod_mat_nrows (m, fq_con),fq_nmod_mat_ncols (m, fq_con));
+   int i,j;
+  for(i=MATROWS(M);i>0;i--)
+  {
+    for(j=MATCOLS(M);j>0;j--)
+    {
+      MATELEM(M,i,j)=convFlintFq_nmodSingP(fq_nmod_mat_entry (m, i-1, j-1),
+                                          fq_con, r);
+    }
+  }
+  return M;
+}
+
 bigintmat* singflint_LLL(bigintmat*  m, bigintmat* T)
 {