Changeset 8d9aa75 in git for libpolys/coeffs/flintcf_Q.cc

Timestamp:

Apr 7, 2016, 3:34:19 PM (8 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')

Children:

a7045fe97d43a7d2f8cfb9568d3d9866e37ad4ad

Parents:

da39e41dbe428e5ee861185b2afe386aae8f7a4d

Message:

ssi: serialization for flint_Zn, flint_Q

File:

• libpolys/coeffs/flintcf_Q.cc (modified) (2 diffs)

libpolys/coeffs/flintcf_Q.cc

@@ -490 +490 @@
   return (number)res;
 }
+static void WriteFd(number a, FILE *f, const coeffs)
+{
+  // format: len a_len(num den) .. a_0
+  fmpq_poly_ptr aa=(fmpq_poly_ptr)a;
+  int l=fmpq_poly_length(aa);
+  fprintf(f,"%d ",l);
+  mpq_t m;
+  mpq_init(m);
+  mpz_t num,den;
+  mpz_init(num);
+  mpz_init(den);
+  for(int i=l; i>=0; i--)
+  {
+    fmpq_poly_get_coeff_mpq(m,(fmpq_poly_ptr)a,i);
+    mpq_get_num(num,m);
+    mpq_get_den(den,m);
+    mpz_out_str (f,SSI_BASE, num);
+    fputc(' ',f);
+    mpz_out_str (f,SSI_BASE, den);
+    fputc(' ',f);
+  }
+  mpz_clear(den);
+  mpz_clear(num);
+  mpq_clear(m);
+}
+static number ReadFd(s_buff f, const coeffs)
+{
+  // format: len a_len .. a_0
+  fmpq_poly_ptr aa=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
+  fmpq_poly_init(aa);
+  int l=s_readint(f);
+  mpz_t num,den;
+  mpq_t m;
+  mpq_init(m);
+  for (int i=l;i>=0;i--)
+  {
+    s_readmpz_base (f,num, SSI_BASE);
+    mpq_set_num(m,num);
+    mpz_clear(num);
+    s_readmpz_base (f,den, SSI_BASE);
+    mpq_set_den(m,den);
+    mpz_clear(den);
+    fmpq_poly_set_coeff_mpq(aa,i,m);
+  }
+  mpq_clear(m);
+  return (number)aa;
+}
 // cfClearContent
 // cfClearDenominators
 static number ConvFactoryNSingN( const CanonicalForm n, const coeffs r)
 {
+  WerrorS("not yet: ConvFactoryNSingN");
 }
 static CanonicalForm ConvSingNFactoryN( number n, BOOLEAN setChar, const coeffs r )
@@ -572 +620 @@
   cf->convFactoryNSingN=ConvFactoryNSingN;
   cf->convSingNFactoryN=ConvSingNFactoryN;
+  cf->cfWriteFd  = WriteFd;
+  cf->cfReadFd = ReadFd;
 #ifdef LDEBUG
   cf->cfDBTest       = DBTest;