numerically stable blockwise polynomial multiplication (should be wor…

…king; more testing necessary)
fredrik-johansson · Apr 2, 2013 · 1741413 · 1741413
1 parent c277827
commit 1741413
Show file tree

Hide file tree

Showing 6 changed files with 544 additions and 4 deletions.
diff --git a/fmprb_poly.h b/fmprb_poly.h
@@ -147,6 +147,13 @@ void fmprb_poly_mullow_classical(fmprb_poly_t res, const fmprb_poly_t poly1,
                                             const fmprb_poly_t poly2,
                                                 long n, long prec);
 
+void _fmprb_poly_mullow_block(fmprb_struct * C,
+    const fmprb_struct * A, long lenA,
+    const fmprb_struct * B, long lenB, long n, long prec);
+
+void fmprb_poly_mullow_block(fmprb_poly_t res, const fmprb_poly_t poly1,
+              const fmprb_poly_t poly2, long len, long prec);
+
 void _fmprb_poly_mullow(fmprb_struct * C,
     const fmprb_struct * A, long lenA,
     const fmprb_struct * B, long lenB, long n, long prec);

diff --git a/fmprb_poly/mullow.c b/fmprb_poly/mullow.c
@@ -35,7 +35,7 @@ _fmprb_poly_mullow(fmprb_struct * res,
     if (n < CUTOFF || len1 < CUTOFF || len2 < CUTOFF)
         _fmprb_poly_mullow_classical(res, poly1, len1, poly2, len2, n, prec);
     else
-        _fmprb_poly_mullow_ztrunc(res, poly1, len1, poly2, len2, n, prec);
+        _fmprb_poly_mullow_block(res, poly1, len1, poly2, len2, n, prec);
 }
 
 void