Add Rössler, Aizawa, Sprott, Levy curve

Author: johnsamuelwrites

README.md

@@ -1,6 +1,6 @@
 # Explorateur de Fractales
 
-Application GitHub Pages qui rend 52 fractales en **WebAssembly**, dont le code de
+Application GitHub Pages qui rend 56 fractales en **WebAssembly**, dont le code de
 calcul est entièrement écrit en **français** grâce au langage
 [multilingual](https://github.com/johnsamuelwrites/multilingual).
 
@@ -26,9 +26,9 @@ françaises** dans `src/*.ml`.
 | Groupe | Fractales |
 |---|---|
 | Évasion | Mandelbrot, Julia, Burning Ship, Tricorn, Multibrot (n=3…8), Celtic, Buffalo, Perpendicular Burning Ship, Heart, Perpendicular Mandelbrot, Perpendicular Celtic, Duck, Buddhabrot |
-| Dynamique | Newton (z³−1), Phoenix, Lyapunov, Lyapunov multiséquence, Bassin de Newton généralisé, Orbitale de Nova, Collatz complexe, Attracteur de Clifford, Attracteur de Peter de Jong, Attracteur d'Ikeda, Attracteur de Hénon |
+| Dynamique | Newton (z³−1), Phoenix, Lyapunov, Lyapunov multiséquence, Bassin de Newton généralisé, Orbitale de Nova, Collatz complexe, Attracteur de Clifford, Attracteur de Peter de Jong, Attracteur d'Ikeda, Attracteur de Hénon, Attracteur de Lorenz, Attracteur de Rössler, Attracteur d'Aizawa, Attracteur de Sprott |
 | IFS | Barnsley (fougère), Sierpinski, Tapis de Sierpinski |
-| L-système | Koch (flocon de neige), Dragon de Heighway, Courbe de Gosper, Joint apollonien, Arbre de Pythagore |
+| L-système | Koch (flocon de neige), Dragon de Heighway, Courbe de Lévy C, Courbe de Gosper, Joint apollonien, Arbre de Pythagore |
 | Magnétiques | Magnet I, Magnet II, Magnet III, Lambda (logistique complexe), Lambda cubique, Magnet cosinus, Magnet sinus, Nova magnétique |
 
 ---

public/js/renderer.js

@@ -220,7 +220,7 @@ déf main():
     afficher(f"    WAT ecrit : {SORTIE_WAT.relative_to(RACINE)}")
     afficher(f"    WASM ecrit: {SORTIE_WASM.relative_to(RACINE)} ({longueur(octets_wasm):,} octets)")
 
-    soit exports_requises = ["mandelbrot", "julia", "burning_ship", "tricorn", "multibrot", "celtic", "buffalo", "perpendicular_burning_ship", "heart", "perpendicular_mandelbrot", "perpendicular_celtic", "duck", "buddhabrot", "newton", "phoenix", "lyapunov", "lyapunov_multisequence", "bassin_newton_generalise", "orbitale_de_nova", "collatz_complexe", "attracteur_de_clifford", "attracteur_de_peter_de_jong", "attracteur_ikeda", "attracteur_de_henon", "lorenz_attractor", "feigenbaum_tree", "barnsley", "sierpinski", "tapis_sierpinski", "menger_sponge", "mandelbulb", "tetraedre_sierpinski", "julia_quaternion", "mandelbox", "vicsek_fractal", "lichtenberg_figures", "koch", "dragon_heighway", "gosper_curve", "cantor_set", "triangle_de_cercles_recursifs", "apollonian_gasket", "t_square_fractal", "h_fractal", "hilbert_curve", "peano_curve", "arbre_pythagore", "magnet1", "magnet2", "magnet3", "lambda_fractale", "lambda_cubique", "magnet_cosinus", "magnet_sinus", "nova_magnetique", "interpoler_lineaire", "interpoler_logarithmique", "ajuster_iterations_export", "mandelbrot_classe"]
+    soit exports_requises = ["mandelbrot", "julia", "burning_ship", "tricorn", "multibrot", "celtic", "buffalo", "perpendicular_burning_ship", "heart", "perpendicular_mandelbrot", "perpendicular_celtic", "duck", "buddhabrot", "newton", "phoenix", "lyapunov", "lyapunov_multisequence", "bassin_newton_generalise", "orbitale_de_nova", "collatz_complexe", "attracteur_de_clifford", "attracteur_de_peter_de_jong", "attracteur_ikeda", "attracteur_de_henon", "lorenz_attractor", "rossler_attractor", "aizawa_attractor", "sprott_attractor", "feigenbaum_tree", "barnsley", "sierpinski", "tapis_sierpinski", "menger_sponge", "mandelbulb", "tetraedre_sierpinski", "julia_quaternion", "mandelbox", "vicsek_fractal", "lichtenberg_figures", "koch", "dragon_heighway", "courbe_levy_c", "gosper_curve", "cantor_set", "triangle_de_cercles_recursifs", "apollonian_gasket", "t_square_fractal", "h_fractal", "hilbert_curve", "peano_curve", "arbre_pythagore", "magnet1", "magnet2", "magnet3", "lambda_fractale", "lambda_cubique", "magnet_cosinus", "magnet_sinus", "nova_magnetique", "interpoler_lineaire", "interpoler_logarithmique", "ajuster_iterations_export", "mandelbrot_classe"]
     valider_exports_wasm(octets_wasm, exports_requises)
     afficher(f"    Exports valides: {', '.join(exports_requises)}")
 

scripts/compile_wasm.ml

@@ -38,6 +38,9 @@
     "attracteur_ikeda",
     "attracteur_de_henon",
     "lorenz_attractor",
+    "rossler_attractor",
+    "aizawa_attractor",
+    "sprott_attractor",
     "feigenbaum_tree",
     "barnsley",
     "sierpinski",
@@ -51,6 +54,7 @@
     "lichtenberg_figures",
     "koch",
     "dragon_heighway",
+    "courbe_levy_c",
     "gosper_curve",
     "cantor_set",
     "triangle_de_cercles_recursifs",

scripts/integration_checks.py

@@ -416,6 +416,24 @@ déf projeter_lorenz_x(x, y, z):
 déf projeter_lorenz_y(x, y, z):
     retour (z - 26.0) * 0.07 - (x + y) * 0.02
 
+déf projeter_rossler_x(x, y, z):
+    retour x * 0.16 + z * 0.035
+
+déf projeter_rossler_y(x, y, z):
+    retour y * 0.16 - z * 0.04
+
+déf projeter_aizawa_x(x, y, z):
+    retour x * 0.18 + z * 0.03
+
+déf projeter_aizawa_y(x, y, z):
+    retour y * 0.18 - z * 0.055
+
+déf projeter_sprott_x(x, y, z):
+    retour x * 0.19 + z * 0.045
+
+déf projeter_sprott_y(x, y, z):
+    retour y * 0.18 - z * 0.03
+
 déf lorenz_attractor(cx, cy, max_iter):
     soit sigma = 10.0
     soit rho = 28.0
@@ -452,6 +470,127 @@ déf lorenz_attractor(cx, cy, max_iter):
         retour max_iter
     retour score
 
+déf rossler_attractor(cx, cy, max_iter):
+    soit a = 0.2
+    soit b = 0.2
+    soit c = 5.7
+    soit dt = 0.02
+    soit x = 0.1
+    soit y = 0.0
+    soit z = 0.0
+    soit meilleur = 1.0e9
+    soit iter_lim = max_iter
+    si iter_lim > 560.0:
+        iter_lim = 560.0
+    soit iter = 0.0
+
+    tantque iter < iter_lim:
+        soit dx = -y - z
+        soit dy = x + a * y
+        soit dz = b + z * (x - c)
+        x = x + dx * dt
+        y = y + dy * dt
+        z = z + dz * dt
+        si iter > 40.0:
+            soit px = projeter_rossler_x(x, y, z)
+            soit py = projeter_rossler_y(x, y, z)
+            soit distance = (px - cx) * (px - cx) + (py - cy) * (py - cy)
+            si distance < meilleur:
+                meilleur = distance
+        si x * x + y * y + z * z > 4096.0:
+            x = 0.1
+            y = 0.0
+            z = 0.0
+        iter = iter + 1.0
+
+    soit score = max_iter - meilleur * 520.0
+    si score < 0.0:
+        retour 0.0
+    si score > max_iter:
+        retour max_iter
+    retour score
+
+déf aizawa_attractor(cx, cy, max_iter):
+    soit a = 0.95
+    soit b = 0.7
+    soit c = 0.6
+    soit d = 3.5
+    soit e = 0.25
+    soit f = 0.1
+    soit dt = 0.008
+    soit x = 0.1
+    soit y = 0.0
+    soit z = 0.0
+    soit meilleur = 1.0e9
+    soit iter_lim = max_iter
+    si iter_lim > 720.0:
+        iter_lim = 720.0
+    soit iter = 0.0
+
+    tantque iter < iter_lim:
+        soit r2 = x * x + y * y
+        soit dx = (z - b) * x - d * y
+        soit dy = d * x + (z - b) * y
+        soit dz = c + a * z - (z * z * z) / 3.0 - r2 * (1.0 + e * z) + f * z * x * x * x
+        x = x + dx * dt
+        y = y + dy * dt
+        z = z + dz * dt
+        si iter > 70.0:
+            soit px = projeter_aizawa_x(x, y, z)
+            soit py = projeter_aizawa_y(x, y, z)
+            soit distance = (px - cx) * (px - cx) + (py - cy) * (py - cy)
+            si distance < meilleur:
+                meilleur = distance
+        si x * x + y * y + z * z > 4096.0:
+            x = 0.1
+            y = 0.0
+            z = 0.0
+        iter = iter + 1.0
+
+    soit score = max_iter - meilleur * 540.0
+    si score < 0.0:
+        retour 0.0
+    si score > max_iter:
+        retour max_iter
+    retour score
+
+déf sprott_attractor(cx, cy, max_iter):
+    soit dt = 0.04
+    soit x = 0.2
+    soit y = 0.1
+    soit z = 0.1
+    soit meilleur = 1.0e9
+    soit iter_lim = max_iter
+    si iter_lim > 520.0:
+        iter_lim = 520.0
+    soit iter = 0.0
+
+    tantque iter < iter_lim:
+        soit dx = y * z
+        soit dy = x - y
+        soit dz = 1.0 - x * y
+        x = x + dx * dt
+        y = y + dy * dt
+        z = z + dz * dt
+        si iter > 32.0:
+            soit px = projeter_sprott_x(x, y, z)
+            soit py = projeter_sprott_y(x, y, z)
+            soit distance = (px - cx) * (px - cx) + (py - cy) * (py - cy)
+            si distance < meilleur:
+                meilleur = distance
+        si x * x + y * y + z * z > 4096.0:
+            x = 0.2
+            y = 0.1
+            z = 0.1
+        iter = iter + 1.0
+
+    soit score = max_iter - meilleur * 520.0
+    si score < 0.0:
+        retour 0.0
+    si score > max_iter:
+        retour max_iter
+    retour score
+
 déf feigenbaum_tree(cx, cy, max_iter):
     soit r = 3.4 + (cx + 1.0) * 0.3
     soit x = 0.5

src/fractales_dynamique.ml

@@ -80,6 +80,39 @@ déf dragon_heighway(cx, cy, max_iter):
         retour 6.0
     retour score
 
+déf courbe_levy_c(cx, cy, max_iter):
+    soit x = cx
+    soit y = cy
+    soit dist = abs_koch(x) + abs_koch(y)
+    soit niveau = 0.0
+    soit echelle = 1.0
+    soit nmax = min_koch(max_iter, 14.0)
+    soit inv_racine2 = 0.70710678118
+
+    tantque niveau < nmax:
+        soit somme = x + y
+        soit diff = y - x
+        soit branche = niveau % 2.0
+        si branche == 0.0:
+            x = (somme - 0.5) * inv_racine2
+            y = diff * inv_racine2
+        sinon:
+            x = (somme + 0.5) * inv_racine2
+            y = diff * inv_racine2
+        soit d = (abs_koch(y) + abs_koch(abs_koch(x) - 0.18)) / echelle
+        si d < dist:
+            dist = d
+        echelle = echelle * 1.41421356237
+        niveau = niveau + 1.0
+
+    soit seuil = 0.028
+    si dist < seuil:
+        retour max_iter * 0.91
+    soit score = max_iter * 0.9 - (dist / (seuil * 10.0)) * (max_iter * 0.9)
+    si score < 6.0:
+        retour 6.0
+    retour score
+
 déf gosper_curve(cx, cy, max_iter):
     soit x = cx
     soit y = cy

src/fractales_lsystem.ml

@@ -16,16 +16,16 @@ importer fractales_classes_compat
 importer fractales_classes
 
 soit MODES_EVASION = ["mandelbrot", "julia", "burning_ship", "tricorn", "multibrot", "celtic", "buffalo", "perpendicular_burning_ship", "heart", "perpendicular_mandelbrot", "perpendicular_celtic", "duck", "buddhabrot"]
-soit MODES_DYNAMIQUE = ["newton", "phoenix", "lyapunov", "lyapunov_multisequence", "bassin_newton_generalise", "orbitale_de_nova", "collatz_complexe", "attracteur_de_clifford", "attracteur_de_peter_de_jong", "attracteur_ikeda", "attracteur_de_henon", "lorenz_attractor", "feigenbaum_tree"]
+soit MODES_DYNAMIQUE = ["newton", "phoenix", "lyapunov", "lyapunov_multisequence", "bassin_newton_generalise", "orbitale_de_nova", "collatz_complexe", "attracteur_de_clifford", "attracteur_de_peter_de_jong", "attracteur_ikeda", "attracteur_de_henon", "lorenz_attractor", "rossler_attractor", "aizawa_attractor", "sprott_attractor", "feigenbaum_tree"]
 soit MODES_IFS = ["barnsley", "sierpinski", "tapis_sierpinski", "menger_sponge", "mandelbulb", "vicsek_fractal", "lichtenberg_figures", "tetraedre_sierpinski", "julia_quaternion", "mandelbox"]
 soit MODES_3D = ["tetraedre_sierpinski", "julia_quaternion", "mandelbox"]
-soit MODES_LSYSTEM = ["koch", "dragon_heighway", "gosper_curve", "cantor_set", "triangle_de_cercles_recursifs", "apollonian_gasket", "t_square_fractal", "h_fractal", "hilbert_curve", "peano_curve", "arbre_pythagore"]
+soit MODES_LSYSTEM = ["koch", "dragon_heighway", "courbe_levy_c", "gosper_curve", "cantor_set", "triangle_de_cercles_recursifs", "apollonian_gasket", "t_square_fractal", "h_fractal", "hilbert_curve", "peano_curve", "arbre_pythagore"]
 soit MODES_MAGNETIQUE = ["magnet1", "magnet2", "magnet3", "lambda_fractale", "lambda_cubique", "magnet_cosinus", "magnet_sinus", "nova_magnetique"]
 soit MODES_CLASSES_COMPAT = ["mandelbrot_classe"]
 affirmer longueur(MODES_EVASION) == 13
-affirmer longueur(MODES_DYNAMIQUE) == 13
+affirmer longueur(MODES_DYNAMIQUE) == 16
 affirmer longueur(MODES_IFS) == 10
 affirmer longueur(MODES_3D) == 3
-affirmer longueur(MODES_LSYSTEM) == 10
+affirmer longueur(MODES_LSYSTEM) == 12
 affirmer longueur(MODES_MAGNETIQUE) == 8
 affirmer longueur(MODES_CLASSES_COMPAT) == 1