First commit

Author: tterrasson

README.md

@@ -0,0 +1,40 @@
+## Generate 3D cellular automata
+
+### Looking for new rules
+
+Generate 16 random 5 states cellular automata after (64 iterations).
+
+```
+$ mkdir -p explore
+$ bin/search -output explore -state 5 -fillseed 0.4 -iteration 64 -nb 16
+```
+
+### Growing an interesting rule
+
+#### Output each iterations (for animation)
+
+```
+$ mkdir -p explore
+$ bin/grow -output explore -input output-XXXXXXX.obj -iteration 128 -worldsize 300
+```
+
+#### Output only last iteration
+
+```
+$ mkdir -p explore
+$ bin/grow -output explore -input output-XXXXXXX.obj -lastonly -iteration 128 -worldsize 300
+```
+
+### Demos
+
+![demo](https://i.gyazo.com/48f9e3d10fca5472f4971fc672896717.png)
+
+![demo](https://i.gyazo.com/99056ee4d79cba57a6073287c6187c8a.png)
+
+![demo](https://media.giphy.com/media/L3L6fgN9HO6mj5AokY/giphy.gif)
+
+![demo](https://media.giphy.com/media/YOBFy0WPAfrueMDYvz/giphy.gif)
+
+More examples :
+- https://www.youtube.com/watch?v=W0iA-pO_XaU
+- https://www.behance.net/gallery/90505477/Cellular-Automata-3D

bin/.keep

@@ -0,0 +1,58 @@
+package main
+
+import (
+	"flag"
+	"fmt"
+	"math/rand"
+	"path/filepath"
+	"time"
+
+	"github.com/tterrasson/crystal/pkg/obj"
+	"github.com/tterrasson/crystal/pkg/rule"
+	"github.com/tterrasson/crystal/pkg/world"
+)
+
+func main() {
+	fmt.Println("Starting...")
+	iterationArg := flag.Int("iteration", 128, "Number of iteration")
+	randomSeedArg := flag.Bool("randseed", false, "Use a random seed")
+	inputArg := flag.String("input", "output.obj", "Input OBJ file")
+	fillSeedArg := flag.Float64("fillseed", 0.25, "Chance to fill seed")
+	outputLastOnlyArg := flag.Bool("lastonly", false, "Output only last iteration")
+	worldSizeArg := flag.Int("worldsize", 100, "World size (SxSxS)")
+	outputPathArg := flag.String("output", "explore", "Output path")
+
+	flag.Parse()
+
+	rand.Seed(time.Now().Unix())
+
+	reader := obj.ObjReader{*inputArg}
+	states := reader.ExtractRuleSet(5)
+	ruleset := rule.RuleSet{states}
+	world := world.NewWorld(*worldSizeArg, &ruleset)
+
+	if *randomSeedArg {
+		world.RandomSeed(float32(*fillSeedArg), 5)
+	} else {
+		seed := reader.ExtractSeed(*worldSizeArg, 200)
+		world.SetSeed(seed)
+	}
+
+	for i := 0; i < *iterationArg; i++ {
+		if !*outputLastOnlyArg {
+			fname := fmt.Sprintf("output-%06d.obj", i)
+			fpath := filepath.Join(*outputPathArg, fname)
+			fmt.Printf("Exporting to %s ...\n", fpath)
+			world.ExportToFile(fpath)
+		}
+
+		fmt.Printf("Step %d ...\n", i)
+		world.Iterate()
+	}
+
+	if *outputLastOnlyArg {
+		fname := fmt.Sprintf("explore/output-%06d.obj", *iterationArg)
+		fmt.Printf("Exporting to %s ...\n", fname)
+		world.ExportToFile(fname)
+	}
+}

cmd/grow/main.go

@@ -0,0 +1,43 @@
+package main
+
+import (
+	"flag"
+	"fmt"
+	"math/rand"
+	"path/filepath"
+	"time"
+
+	"github.com/tterrasson/crystal/pkg/rule"
+	"github.com/tterrasson/crystal/pkg/world"
+)
+
+func main() {
+	fmt.Println("Starting...")
+	iterationArg := flag.Int("iteration", 64, "Number of iteration")
+	stateArg := flag.Int("state", 3, "Number of states")
+	fillSeedArg := flag.Float64("fillseed", 0.25, "Chance to fill seed")
+	fillRuleArg := flag.Float64("fillrule", 0.25, "Chance to fill rule")
+	worldSizeArg := flag.Int("worldsize", 100, "World size (SxSxS)")
+	outputPathArg := flag.String("output", "explore", "Output path")
+	nbArg := flag.Int("nb", 8, "Number of random CA to generate")
+
+	flag.Parse()
+
+	rand.Seed(time.Now().Unix())
+
+	ruleset := rule.NewRandomRuleSet(float32(*fillRuleArg), *stateArg)
+	world := world.NewWorld(*worldSizeArg, ruleset)
+	world.RandomSeed(float32(*fillSeedArg), *stateArg)
+
+	for n := 0; n < *nbArg; n++ {
+		for i := 0; i < *iterationArg; i++ {
+			fmt.Printf("Step %d ...\n", i)
+			world.Iterate()
+		}
+
+		fname := fmt.Sprintf("output-%f-%f.obj", rand.Float64())
+		fpath := filepath.Join(*outputPathArg, fname)
+		fmt.Printf("Exporting to %s ...\n", fpath)
+		world.ExportToFile(fpath)
+	}
+}

cmd/search/main.go

@@ -0,0 +1,25 @@
+newmtl state1
+Ka 0.545 0.156 0.121
+Kd 0.545 0.156 0.121
+Ns 10.000
+illum 2
+newmtl state2
+Ka 0.835 0.439 0.188
+Kd 0.835 0.439 0.188
+Ns 10.000
+illum 2
+newmtl state3
+Ka 0.517 0.643 0.352
+Kd 0.517 0.643 0.352
+Ns 10.000
+illum 2
+newmtl state4
+Ka 0.949 0.615 0.294
+Kd 0.949 0.615 0.294
+Ns 10.000
+illum 2
+newmtl state5
+Ka 0.249 0.215 0.894
+Kd 0.249 0.215 0.894
+Ns 10.000
+illum 2

explore/colors.mtl

@@ -0,0 +1,129 @@
+package obj
+
+import (
+	"bufio"
+	"encoding/json"
+	"io"
+	"log"
+	"os"
+	"strconv"
+	"strings"
+)
+
+type ObjReader struct {
+	Filename string
+}
+
+func (reader *ObjReader) initSeed(size int) [][][]int {
+	seed := make([][][]int, size)
+
+	for x := range seed {
+		seed[x] = make([][]int, size)
+		for y := range seed[x] {
+			seed[x][y] = make([]int, size)
+		}
+	}
+
+	return seed
+}
+
+func (reader *ObjReader) initRuleSet(maxStates int) [][][][]int {
+	states := make([][][][]int, maxStates+1)
+
+	for state := 0; state <= maxStates; state++ {
+		states[state] = make([][][]int, 7)
+
+		for face := 0; face < 7; face++ {
+			states[state][face] = make([][]int, 13)
+
+			for edge := 0; edge < 13; edge++ {
+				states[state][face][edge] = make([]int, 9)
+			}
+		}
+	}
+
+	return states
+}
+
+func (reader *ObjReader) ExtractSeed(size int, offset int) [][][]int {
+	file, err := os.Open(reader.Filename)
+	if err != nil {
+		log.Panicf("failed reading file: %s", err)
+	}
+	defer file.Close()
+
+	out := reader.initSeed(size)
+	rd := bufio.NewReader(file)
+
+	for {
+		line, err := rd.ReadString('\n')
+
+		if strings.HasPrefix(line, "# Seed : ") {
+			jseed := strings.Replace(line, "# Seed : ", "", 1)
+			var seed map[string]interface{}
+			json.Unmarshal([]byte(jseed), &seed)
+
+			for k, v := range seed {
+				idx := strings.Split(k, ",")
+				x, _ := strconv.Atoi(idx[0])
+				y, _ := strconv.Atoi(idx[1])
+				z, _ := strconv.Atoi(idx[2])
+				out[x+offset][y+offset][z+offset] = int(v.(float64))
+			}
+
+			return out
+		}
+
+		if err == io.EOF {
+			break
+		}
+	}
+
+	return nil
+}
+
+func (reader ObjReader) ExtractRuleSet(maxStates int) [][][][]int {
+	file, err := os.Open(reader.Filename)
+	if err != nil {
+		log.Panicf("failed reading file: %s", err)
+	}
+	defer file.Close()
+
+	out := reader.initRuleSet(maxStates)
+	rd := bufio.NewReader(file)
+
+	for {
+		line, err := rd.ReadString('\n')
+
+		if strings.HasPrefix(line, "# Rule : ") {
+			jrule := strings.Replace(line, "# Rule : ", "", 1)
+			var rule map[string]interface{}
+			json.Unmarshal([]byte(jrule), &rule)
+			stateSet := rule["states"].([]interface{})
+
+			for state := 0; state <= maxStates; state++ {
+				faceSet := stateSet[state].([]interface{})
+
+				for face := 0; face < 7; face++ {
+					edgeSet := faceSet[face].([]interface{})
+
+					for edge := 0; edge < 13; edge++ {
+						cornerSet := edgeSet[edge].([]interface{})
+
+						for corner := 0; corner < 9; corner++ {
+							out[state][face][edge][corner] = int(cornerSet[corner].(float64))
+						}
+					}
+				}
+			}
+
+			return out
+		}
+
+		if err == io.EOF {
+			break
+		}
+	}
+
+	return nil
+}