neuroevolution readme

Author: genmeblog

Diff for: results/neuroevolution/success.png

@@ -61,11 +61,14 @@
   "Make next spore"
   ([id] (next-spore id nil))
   ([id pos]
+   ;; 180 17
    ;; a 180 b -5 !!! -10, 17 (default)
    ;; a -69 b -1
    ;; 42 9
    ;; 79.6 -1.8
-   (spore id pos 0.67 (m/sq 11.0) (m/radians 180) (m/radians 17) (first pal))))
+   ;; 0 5
+   ;; 180 179
+   (spore id pos 0.67 (m/sq 11.0) (m/radians 42) (m/radians 9) (first pal))))
 
 (defn draw
   [canvas window _ spores]
@@ -86,3 +89,5 @@
 
 
 (comment c2d/save window "results/ex53/spores.jpg")
+
+(m/unuse-primitive-operators)

Diff for: src/ex53_ppl.clj

@@ -121,20 +121,114 @@
       (c2d/set-color canvas (.color ptype))
       (c2d/ellipse canvas (.x p) (.y p) 4 4))
     #_(when (= frame 200) (c2d/save canvas "results/ex63/plife.jpg"))
-    [setup (move setup particles 0.02)]))
+    [setup (move setup particles 0.1)]))
 
-(def setup (random-setup (r/irand 3 8)))
+(def setup (random-setup (r/irand 2 5)))
+
+(defn save->setup [m]
+  (-> m
+      (update :alpha double-array)
+      (update :beta double-array)
+      (update :types (partial mapv map->ParticleType))
+      (map->Setup)))
 
 (def window (c2d/show-window {:canvas (c2d/black-canvas csize csize :highest)
                             :draw-fn draw
                             :background :black
-                            :draw-state (do
+                            :draw-state (let [setup setup #_(save->setup (setups :s1))]
                                           (pp/pprint setup)
                                           [setup (random-particles setup 1200)])}))
 
 
+(def pairs-love-hate {:n 2,
+                    :alpha
+                    [0.8531190685431542, 0.4124191049307851, 0.6363336352718891, 0.0],
+                    :beta
+                    [0.17515989405988655, 0.46898628375997914, 1.0E-5, 0.245296441325918],
+                    :types
+                    [{:id 0,
+                      :rmax 52.39644371175454,
+                      :color [31.0 119.0 180.0 255.0],
+                      :friction 0.9735391337659955,
+                      :directional? false}
+                     {:id 1,
+                      :rmax 76.58245632190017,
+                      :color [255.0 127.0 14.0 255.0],
+                      :friction 0.8565072274554544,
+                      :directional? false}],
+                    :fmult 11.78333128371779})
 
-(def setups {:s4 {:n 3,
+(def setups {:s5 {:n 7,
+                :alpha
+                [0.14861021110200467, 0.6692996723215351, 0.8164479197480894,
+                 0.2971907681082335, 0.3390239085309035, -0.7785871214228197,
+                 0.807017246117286, 0.9897060316942456, 0.22462200467136376,
+                 0.24556166541932525, 0.9585743441404408, 0.46410543836135454,
+                 0.582672695098386, 0.4553549448206484, -0.9197920663869434,
+                 0.7244029506599556, 0.809370932139867, 0.8852742660527777,
+                 0.6153018098109097, 0.0, 0.8316971648547771, -0.6639246865662674,
+                 0.5979391646072336, 0.6024358095947859, 0.0, 0.5092272162340464,
+                 0.7028458867058891, 0.0, 0.380693263489916, 0.47368547609290784, 0.0,
+                 -0.01146080240482461, 0.11103110451012577, -0.9841353913583804,
+                 0.22910442818466126, 0.21225199666703465, 0.0, 0.3816397307749264,
+                 0.7107575951883653, -0.774523033517672, 0.0, 0.2784820533021294,
+                 0.5996965240685289, 0.5250427837238822, 0.2640012136092506,
+                 0.574189388135993, 0.5634533172450109, 0.49515357274531,
+                 -0.23345515434391828],
+                :beta
+                [0.39652995834906474, 0.18350480619310794, 0.34400004339297835,
+                 1.0E-5, 0.17535457154020762, 0.24633939044495803, 1.0E-5,
+                 0.0029488960022599027, 0.3854141065610507, 0.3037948040259639,
+                 0.1315968496394307, 0.09738530814951005, 1.0E-5, 0.3232729351116344,
+                 1.0E-5, 0.22153527794443625, 0.3268308089470868, 0.09744546120909733,
+                 1.0E-5, 0.06246113757148177, 0.05631946789826031, 0.3366222338984227,
+                 0.03675528120398472, 0.441394942824003, 0.0563859289193688,
+                 0.36463601158570325, 0.029368028675658624, 0.08811533231738519,
+                 0.017925239286269264, 0.2122421441106669, 0.2503850525716419,
+                 0.3272795892148501, 0.3991277529353564, 0.1418344093235479,
+                 0.011702419598665973, 0.3549620830153355, 0.17802794018769988,
+                 0.29575004512407876, 0.043728710826954426, 0.4477330979749426,
+                 0.08525557470032102, 0.20495483818859214, 0.28051143823078967,
+                 0.4399178328957438, 0.2692395895013477, 0.3540625752279375, 1.0E-5,
+                 0.1663827780283793, 0.2801211666493402],
+                :types
+                [{:id 0,
+                  :rmax 13.290700740078464,
+                  :color [31.0 119.0 180.0 255.0],
+                  :friction 0.9527625125622025,
+                  :directional? false}
+                 {:id 1,
+                  :rmax 80.6737289396209,
+                  :color [255.0 127.0 14.0 255.0],
+                  :friction 0.9334759162001242,
+                  :directional? false}
+                 {:id 2,
+                  :rmax 24.254597157847872,
+                  :color [44.0 160.0 44.0 255.0],
+                  :friction 0.7240566396062796,
+                  :directional? true}
+                 {:id 3,
+                  :rmax 71.76840597967437,
+                  :color [214.0 39.0 40.0 255.0],
+                  :friction 0.7899062609246527,
+                  :directional? false}
+                 {:id 4,
+                  :rmax 75.90713756033553,
+                  :color [148.0 103.0 189.0 255.0],
+                  :friction 0.8092729978675278,
+                  :directional? false}
+                 {:id 5,
+                  :rmax 66.98691395167273,
+                  :color [140.0 86.0 75.0 255.0],
+                  :friction 0.9857219741051314,
+                  :directional? false}
+                 {:id 6,
+                  :rmax 58.981056481884046,
+                  :color [227.0 119.0 194.0 255.0],
+                  :friction 0.7945210982191895,
+                  :directional? false}],
+                :fmult 12.216106395481184}
+           :s4 {:n 3,
                 :alpha
                 [0.6684104736324896, 0.964616865419524, 0.0, 0.4797524863542435,
                  0.420585350696979, 0.0, 0.2919907666469499, 0.764963999761529,

Diff for: src/ex63_particle_life.clj

@@ -43,7 +43,8 @@
      :Da 1.0 :Db m/THIRD
      :f 0.055 :k 0.062
      :dt 1.0
-     :show-map? true}))
+     :show-map? true
+     :color-a? true}))
 
 (defn diffuse! [^double d ^doubles source ^doubles target]
   (doseq [^long x rsize ^long y rsize
@@ -80,6 +81,10 @@
                                                              (m/- abb (m/* (m/+ f k) b))))) 0.0 1.0)))
   data)
 
+(defn color-a [^double a ^double b]  (v/vec3 (m/- 1.0 a) b (m/sqrt (m/* a b))))
+(defn color-b [^double a ^double b]  (v/vec3 (m/- 1.0 b) a (m/sqrt (m/* b a))))
+
+
 (defn draw
   [canvas window _ {:keys [^doubles A ^doubles B show-map? f k] :as curr}]
   (when (c2d/mouse-pressed? window)
@@ -96,14 +101,16 @@
                                   :k (norm-k mx)
                                   :show-map? false))
                  :m (assoc curr :show-map? true)
+                 :z (update curr :color-a? not)
                  curr)
-               curr)]
+               curr)
+        color-fn (if (:color-a? curr) color-a color-b)]
 
     (c2d/set-background canvas (c/color 10 10 20) 100)
     (doseq [^long x rsize ^long y rsize
             :let [a (Array/get2d A size x y)
                   b (Array/get2d B size x y)]]
-      (c2d/set-color canvas (v/mult (v/vec3 (m/- 1.0 a) b (m/sqrt (m/* a b))) 255.0))
+      (c2d/set-color canvas (v/mult (color-fn a b) 255.0))
       (c2d/rect canvas (m/* x scale) (m/* y scale) scale scale))
     (let [f (m/approx (if show-map? (norm-f my) f) 6)
           k (m/approx (if show-map? (norm-k mx) k) 6)]

Diff for: src/ex65_reaction_diffusion.clj

@@ -1,34 +1,77 @@
 (ns ex68-langton-ant
   (:require [clojure2d.pixels :as p]
-            [clojure2d.core :as c2d]))
+            [clojure2d.core :as c2d]
+            [clojure2d.color :as c]
+            [fastmath.vector :as v]
+            [clojure2d.extra.utils :as utils]
+            [fastmath.core :as m]))
 
-(defn add [[x1 y1] [x2 y2] w h] [(mod (+ x1 x2) w) (mod (+ y1 y2) h)])
+(defn modadd
+  [[^long x1 ^long y1] [^long x2 ^long y2] ^long w ^long h]
+  [(m/mod (m/+ x1 x2) w)
+   (m/mod (m/+ y1 y2) h)])
 
-(defn move [buff [[x y :as pos] [dx dy]]]
+(defn move [buff [[x y :as pos] [^long dx ^long dy] rules]]
   (let [w (c2d/width buff)
         h (c2d/height buff)
-        col? (pos? (p/get-value buff 0 x y))
-        ndir (if col? [dy (- dx)] [(- dy) dx])]
-    (p/set-color! buff x y (if col? :black :white))
-    [(add pos ndir w h) ndir]))
-
-(def speed 100)
-
-(defn draw [canvas _ frame [buff [[x y] _ :as pos-dir]]]
-  (-> canvas
-      (c2d/set-background 20 20 20)
-      (c2d/image (c2d/resize buff 800 800))
-      (c2d/set-color :blue 200)
-      (c2d/rect 10 10 100 35)
-      (c2d/set-color :white)
-      (c2d/set-font-attributes 20)
-      (c2d/text (str (* speed frame)) 20 35)
-      (c2d/set-color :red)
-      (c2d/set-stroke 2)
-      (c2d/ellipse (* 2 x) (* 2 y) 10 10 true ))
-  [buff (nth (iterate (partial move buff) pos-dir) speed)])
-
-(def window (c2d/show-window {:canvas (c2d/canvas 800 800)
+        c1 (c/set-alpha (p/get-color buff x y) 255.0)
+        [c2 r] (rules c1)
+        ndir (if (= \L r) [dy (m/- dx)] [(m/- dy) dx])]
+    (p/set-color! buff x y c2)
+    [(modadd pos ndir w h) ndir rules]))
+
+(defn draw [canvas window _ {:keys [^long speed ^long steps ^long scale buff pos-dir] :as data}]
+  (if (c2d/key-pressed? window)
+    data
+    (let [[pos dir] pos-dir
+          hsc (m// scale 2)
+          r (* hsc 5)
+          pos2 (v/shift (v/mult pos scale) hsc)
+          dir2 (v/add pos2 (v/mult dir scale))]
+      (-> canvas
+          (c2d/set-background 20 20 20)
+          (c2d/image buff)
+          (c2d/set-color :blue 200)
+          (c2d/rect 10 10 150 35)
+          (c2d/set-color :white)
+          (c2d/set-font-attributes 20)
+          (c2d/text (str steps) 20 35)
+          (c2d/set-color :red)
+          (c2d/set-stroke 2)
+          (c2d/ellipse (pos2 0) (pos2 1) r r true)
+          (c2d/line (pos2 0) (pos2 1) (dir2 0) (dir2 1)))
+      (assoc data
+             :steps (+ speed steps)
+             :pos-dir (nth (iterate (partial move buff) pos-dir) speed)))))
+
+(def pal (concat [(c/color :black)] (c/palette :category20)))
+
+(defn rule->data [rule]
+  (->> (take (count rule) pal)
+       (cycle)
+       (partition 2 1)
+       (take (count rule))
+       (map (fn [r [c1 c2]] [c1 [c2 r]]) rule)
+       (into {})))
+
+(def r1 "RRLLLRLLLRRR")
+(def r2 "LLRRRLRLRLLR")
+(def r3 "LRRRRRLLR")
+(def r4 "LRRL")
+(def r5 "RRLRR")
+(def r6 "RRRL")
+(def r7 "RRRRL")
+
+(defn init [{:keys [^long scale rule] :as data}]
+  (let [ps (m// 800 scale)]
+    (assoc data
+           :steps 0
+           :buff (p/pixels ps ps)
+           :pos-dir [[(m/inc (m// ps 2)) (m/inc (m// ps 2))] [0 1] (rule->data rule)])))
+
+(def window (c2d/show-window {:canvas (c2d/canvas 800 800 :low)
                             :background :black
                             :draw-fn draw
-                            :draw-state [(p/pixels 400 400) [[100 100] [0 1]]]}))
+                            :draw-state (init {:scale 2 :speed 10000 :rule "RLRRLRLL"})}))
+
+(comment (utils/show-palette pal))

Diff for: src/ex68_langton_ant.clj

@@ -0,0 +1,106 @@
+;; ince beams example - idea by @CTW (Thomas Clark)
+;; https://en.wikipedia.org/wiki/Gaussian_beam#Ince-Gaussian_modes
+;; https://www.researchgate.net/publication/8567409_Ince-Gaussian_modes_of_the_paraxial_wave_equation_and_stable_resonators
+;; Clojure Conj 2024: https://youtu.be/_D5d6Ls6pBw?t=1951
+(ns ex71-ince-beams
+  (:require [clojure2d.core :as c2d]
+            [clojure2d.pixels :as p]
+            [fastmath.core :as m]
+            [fastmath.vector :as v]
+            [fastmath.polynomials :as poly]
+            [clojure2d.extra.utils :as u]
+            [clojure2d.color :as c]
+            [fastmath.complex :as cplx]
+            [fastmath.random :as r])
+  (:import [fastmath.vector Vec2]))
+
+(set! *warn-on-reflection* true)
+(set! *unchecked-math* :warn-on-boxed)
+
+(defn rayleigh-range
+  ^double [^double waist ^double wavelength]
+  (m// (m/* m/PI (m/sq waist)) wavelength))
+
+(defn width
+  ^double [^double z ^double zR ^double waist]
+  (m/* waist (m/sqrt (m/inc (m/sq (m// z zR))))))
+
+(defn csq
+  ^double [^double ellipticity ^double waist]
+  (m/inc (m/* 0.5 (m/sq waist) ellipticity)))
+
+(defn IG0
+  "An Ince-Gaussian beam function for z=0."
+  [p m ellipticity waist wavelength]
+  (let [ellipticity (double ellipticity)
+        waist (double waist)
+        zR (rayleigh-range waist wavelength)
+        w (width 0.0 zR waist)
+        ince (poly/ince-C p m ellipticity)
+        ince-radial (poly/ince-C-radial p m ellipticity)
+        k (m// m/TWO_PI (double wavelength))
+        rscale (m// (m/* w (m/sqrt (m/* 0.5 ellipticity))))]
+
+    (fn ^double [^double x ^double y]
+      (let [^Vec2 xi-eta (cplx/acosh (cplx/scale (cplx/complex x y) rscale))
+            r**2 (m/+ (m/sq x) (m/sq y))]
+        (m/* (m// waist w)
+             (double (ince-radial (xi-eta 0)))
+             (double (ince (xi-eta 1)))
+             (m/exp (m/* -1.0 k
+                         (m// r**2 (m/* 2.0 zR)))))))))
+
+(defn IG-profile
+  ([IG-fn]
+   (IG-profile IG-fn [0 0] [200 200]))
+  ([IG-fn [^double x0 ^double y0] [^long width ^long height]]
+   (let [hw (m// width 2.0)
+         hh (m// height 2.0)]
+     (for [^long x (range width)
+           ^long y (range height)]
+       (m/sq (IG-fn (m/+ x0 (m/- x hw))
+                    (m/+ y0 (m/- y hh))))))))
+
+(defn normalize
+  "Peak normalize a tensor."
+  ;; TODO: extend to deal with non-zero floored tensors.
+  [d]
+  (let [peak (double (reduce m/max d))]
+    (map (fn [^double v] (m// v peak)) d)))
+
+(defn profile [p m ellipticity waist]
+  (-> (IG0 p m ellipticity waist 0.78)
+      (IG-profile [0.0 0.0] [100 100])
+      (normalize)))
+
+(defn show-profile
+  [p m ellipticity waist]
+  (let [p (profile p m ellipticity waist)
+        g (c/gradient [:black (c/gray 200) :white])]
+    (c2d/with-canvas [c (c2d/canvas 800 800)]
+      (c2d/set-background c :black)
+      (doseq [[^long i ^double v] (map-indexed vector p)
+              :let [x (m/* 8 (m/quot i 100))
+                    y (m/* 8 (m/rem i 100))]]
+        (c2d/set-color c (g v))
+        (c2d/rect c x y 8 8))
+      (u/show-image c))))
+
+(show-profile 11 5 7 11)
+
+(defn render-profile
+  [p m ellipticity waist wavelength]
+  (let [buffer (p/gradient-renderer 800 800 :sinc 1)
+        r (range -100 100 0.05)
+        n (m/make-norm -100 100 0 800)
+        p (IG0 p m ellipticity waist wavelength)]
+    (doseq [^double x r ^double y r
+            :let [x (m/+ x (r/grand 0.05))
+                  y (m/+ y (r/grand 0.05))
+                  xx (n x) yy (n y)
+                  v (m/sq (double (p x y)))]]
+      (p/add-pixel! buffer xx yy v))
+    (u/show-image (p/to-pixels buffer {:logarithmic? true}))))
+
+(render-profile 7 5 3 15 0.78)
+