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 wk=c∏j≠k(zk−zj)w_k = c \prod_{j \ne k} (z_k - z_j), with cc = lc_r + i lc_i, for nstart≤k<nendn_{\mathrm{start}} \le k < n_{\mathrm{end}}. cd_poly_weierstrass uses the dd points values for both factors and exploits the symmetry of the differences; cd_poly_weierstrass_distinct_orders takes the factors zjz_j from col_values (in an order chosen to avoid overflow) and the points zkz_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-11-02).