qmlmathplot 0.4.2__tar.gz → 0.4.3__tar.gz

{qmlmathplot-0.4.2 → qmlmathplot-0.4.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: qmlmathplot
-Version: 0.4.2
+Version: 0.4.3
 Summary: Embeddable function plotter for Qt Quick / QtWidgets (implicit per-pixel drawing, MVVM)
 Author: I-AM-A-NOOB
 Author-email: I-AM-A-NOOB <suzc08@outlook.com>
@@ -281,8 +281,11 @@ corner is in the domain and no pole factor changes sign between them; the edge c
 interpolated and the cell's contour segment is stroked (the saddle case is paired using the cell
 centre's sign). The first-order/Newton path still supplies the sub-pixel width; the cell term
 adds the coverage of a thin or tangent branch whose centre projection fails. A **singular cell**
-(`|Fx| + |Fy| < 1e-3` in world units, plus `|F| < 0.05`) gets a small disc, so an isolated point
-such as `x**2 + y**2 = 0` is not left to the projection alone.
+(`|Fx| + |Fy| < 1e-3` in world units, `|F| < 0.05`, **and** `|F| < R·|grad F|` in pixels) gets
+a small disc, so an isolated point such as `x**2 + y**2 = 0` is not left to the projection alone.
+The last test is the scale-invariant one: it requires the point to be on the set to first order,
+so a scaled expression such as `(x+y)/(x·y)` — whose F and gradient are tiny far from its line —
+no longer fills whole quadrants with the disc.
 
 Still not ported: GGB's Bernstein subdivision (**polynomial-only**), its adaptive quad-tree
 refinement and the `ContourAssembler`/`ContourLinker` (they exist to produce a selectable vector

{qmlmathplot-0.4.2 → qmlmathplot-0.4.3}/README.md

@@ -268,8 +268,11 @@ corner is in the domain and no pole factor changes sign between them; the edge c
 interpolated and the cell's contour segment is stroked (the saddle case is paired using the cell
 centre's sign). The first-order/Newton path still supplies the sub-pixel width; the cell term
 adds the coverage of a thin or tangent branch whose centre projection fails. A **singular cell**
-(`|Fx| + |Fy| < 1e-3` in world units, plus `|F| < 0.05`) gets a small disc, so an isolated point
-such as `x**2 + y**2 = 0` is not left to the projection alone.
+(`|Fx| + |Fy| < 1e-3` in world units, `|F| < 0.05`, **and** `|F| < R·|grad F|` in pixels) gets
+a small disc, so an isolated point such as `x**2 + y**2 = 0` is not left to the projection alone.
+The last test is the scale-invariant one: it requires the point to be on the set to first order,
+so a scaled expression such as `(x+y)/(x·y)` — whose F and gradient are tiny far from its line —
+no longer fills whole quadrants with the disc.
 
 Still not ported: GGB's Bernstein subdivision (**polynomial-only**), its adaptive quad-tree
 refinement and the `ContourAssembler`/`ContourLinker` (they exist to produce a selectable vector

{qmlmathplot-0.4.2 → qmlmathplot-0.4.3}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "qmlmathplot"
-version = "0.4.2"
+version = "0.4.3"
 description = "Embeddable function plotter for Qt Quick / QtWidgets (implicit per-pixel drawing, MVVM)"
 readme = "README.md"
 requires-python = ">=3.11"

{qmlmathplot-0.4.2 → qmlmathplot-0.4.3}/pyproject.toml.orig

@@ -1,6 +1,6 @@
 [project]
 name = "qmlmathplot"
-version = "0.4.2"
+version = "0.4.3"
 description = "Embeddable function plotter for Qt Quick / QtWidgets (implicit per-pixel drawing, MVVM)"
 readme = "README.md"
 authors = [

{qmlmathplot-0.4.2 → qmlmathplot-0.4.3}/src/qmlmathplot/model.py

@@ -1287,8 +1287,14 @@ void main() {
     // such as x² + y² = 0 is not left to the projection alone.
     // GGB measures the singular cell in *world* derivative units (|Fx| + |Fy| < 0.001);
     // the per-pixel gradient of an ordinary circle is ~0.01 and must not qualify.
+    // The distance gate is essential for a scaled expression: (x+y)/(x*y) has |F| < 0.05 and
+    // |Fx|+|Fy| < 0.001 over huge regions around ±50, although the level set is far away. A
+    // singular point must also be *on* the set to the first order: |F| ~ R*|grad F|. The
+    // ratio is scale-invariant, unlike the two absolute thresholds above.
     float gw_s = abs(FX(x, y)) + abs(FY(x, y));
-    if (dom && (f == f) && (abs(f) < 1e30) && gw_s < 1e-3 && abs(f) < 0.05) {
+    float gnorm_s = max(sqrt(gx * gx + gy * gy), 1e-9);
+    if (dom && (f == f) && (abs(f) < 1e30) && gw_s < 1e-3 && abs(f) < 0.05
+            && abs(f) < R * gnorm_s) {
         march_cov = 1.0;
     }