FLINT3.7.0-dev

cd_poly_weierstrass, cd_poly_weierstrass_distinct_orders

Synopsis#include <flint/acb_poly.h>
void cd_poly_weierstrass(double * results_r, double * results_i, double lc_r, double lc_i, const double * values_r, const double * values_i, slong n_start, slong n_end, slong d)source
void cd_poly_weierstrass_distinct_orders(double * results_r, double * results_i, double lc_r, double lc_i, const double * col_values_r, const double * col_values_i, slong d, const double * row_values_r, const double * row_values_i, slong n_start, slong n_end)source

Compute the Weierstrass denominators w_k = c \prod_{j \ne k} (z_k - z_j), with c = lc_r + i lc_i, for n_{\mathrm{start}} \le k < n_{\mathrm{end}}. cd_poly_weierstrass uses the d points values for both factors and exploits the symmetry of the differences; cd_poly_weierstrass_distinct_orders takes the factors z_j from col_values (in an order chosen to avoid overflow) and the points z_k from row_values, skipping a factor that is exactly zero.

Cite this page

The FLINT team, “cd_poly_weierstrass, cd_poly_weierstrass_distinct_orders (acb_poly.h),” FLINT: Fast Library for Number Theory, documentation, version 3.7.0-dev, 2026. https://flintlib.org/doc/acb_poly/cd_poly_weierstrass.html (generated 2026-10-20).