blue_sampler.cluster2points

blue_sampler.cluster2points(clusters, p=3)[source]

Convert clusters into low-discrepancy point set.

Each cluster is replaced by m points by solving a moment-matching problem (Levenberg-Marquardt) via momentum_fit.

Parameters:
  • clusters (NDArray) –

    e.g. the output of sample_clusters, or the atoms output of

    sample_tessels(…, return_atoms=True).

  • p (int) – Maximum total moment order to match (centroid plus central moments up to order p).

Returns:

m points per tessel or cluster, matching its moments up to order p.

Return type:

NDArray