source: cpp/frams/genetics/f9/f9_conv.cpp @ 1130

Last change on this file since 1130 was 1130, checked in by Maciej Komosinski, 3 years ago

Used std::min(), std::max() explicitly to avoid compiler confusion. Used std::size() explicitly instead of the equivalent macro

  • Property svn:eol-style set to native
File size: 5.3 KB
RevLine 
[286]1// This file is a part of Framsticks SDK.  http://www.framsticks.com/
[1130]2// Copyright (C) 1999-2021  Maciej Komosinski and Szymon Ulatowski.
[286]3// See LICENSE.txt for details.
[120]4
[779]5#include "f9_conv.h"
[120]6#include <frams/model/model.h>
7#include <string.h>
8
[259]9#define APPLY_DETERMINISTIC_BODY_NOISE //this genetic representation easily produces perfectly vertical sticks that would stay upright forever in simulation. In most cases such infinite perfection is not desired, so we make the construct less perfect by perturbing its coordinates.
[120]10
[139]11GenoConv_f90::GenoConv_f90()
[120]12{
13        name = "Turtle3D-ortho encoding";
[168]14        in_format = '9';
[120]15        out_format = '0';
[259]16        mapsupport = 1;
[120]17}
18
19
[643]20const char* turtle_commands_f9 = "LRBFDU";
[120]21
22//const char* turtle_commandsX_f9="-+0000";
23//const char* turtle_commandsY_f9="00-+00";
24//const char* turtle_commandsZ_f9="0000-+";
25
[736]26SString GenoConv_f90::convert(SString &in, MultiMap *map, bool using_checkpoints)
[120]27{
28        vector<XYZ_LOC> vertices;
29        XYZ_LOC current;
30        Model m;
[740]31        m.open(using_checkpoints);
[259]32        int recently_added = addSegment(m, 0, vertices, current, 0xDead);
[973]33        for (int i = 0; i < in.length(); i++)
[120]34        {
[168]35                char command = in[i];
36                char *ptr = strchr((char*)turtle_commands_f9, command);
[120]37                if (ptr)
38                {
[168]39                        int delta[] = { 0, 0, 0 };
40                        int pos = ptr - turtle_commands_f9;
41                        int axis = pos / 2;
42                        int dir = pos % 2;
43                        (*(delta + axis)) += dir * 2 - 1; //+1 or -1 in the given axis
[120]44                        current.add(delta);
[259]45                        recently_added = addSegment(m, i, vertices, current, recently_added);
[740]46                        m.checkpoint();
[120]47                }
48        }
49#ifdef APPLY_DETERMINISTIC_BODY_NOISE
50        perturbPartLocations(m);
51#endif
[664]52        setColors(m, recently_added);
[120]53        m.close();
[157]54        if (m.getPartCount() < 2) //only one part <=> there were no valid turtle commands in the input genotype
55                return ""; //so we return an invalid f0 genotype
[259]56        if (map != NULL)
57                m.getCurrentToF0Map(*map);
[534]58        return m.getF0Geno().getGenes();
[120]59}
60
[259]61int GenoConv_f90::addSegment(Model &m, int genenr, vector<XYZ_LOC> &vertices, const XYZ_LOC &new_vertex, int recently_added)
[120]62{
[168]63        if (vertices.size() < 1) //empty model?
[120]64        {
[168]65                return addNewVertex(m, vertices, new_vertex);
66        }
67        else
[120]68        {
[168]69                int vertex_here = findVertexAt(vertices, new_vertex);
70                if (vertex_here < 0) //need to create a new Part
[120]71                {
[168]72                        vertex_here = addNewVertex(m, vertices, new_vertex);
[120]73                } //else there already exists a Part in new_vertex; new Joint may or may not be needed
[168]74                Part *p1 = m.getPart(recently_added);
75                Part *p2 = m.getPart(vertex_here);
[259]76                p1->addMapping(MultiRange(genenr, genenr));
77                p2->addMapping(MultiRange(genenr, genenr));
78
79                int j12 = m.findJoint(p1, p2);
80                int j21 = m.findJoint(p2, p1);
81                if (j12 >= 0)
82                        m.getJoint(j12)->addMapping(MultiRange(genenr, genenr));
83                else if (j21 >= 0)
84                        m.getJoint(j21)->addMapping(MultiRange(genenr, genenr));
85                else //both j12<0 and j21<0. New Joint needed. Should always happen if we just created a new Part (vertex_here was <0)
86                        m.addNewJoint(p1, p2)->addMapping(MultiRange(genenr, genenr));
[120]87                return vertex_here;
88        }
89}
90
[168]91int GenoConv_f90::findVertexAt(vector<XYZ_LOC> &vertices, const XYZ_LOC &vertex)
[120]92{
[319]93        for (size_t i = 0; i < vertices.size(); i++)
[120]94                if (vertices[i].same_coordinates(vertex)) return i;
95        return -1;
96}
97
98
[168]99int GenoConv_f90::addNewVertex(Model &m, vector<XYZ_LOC> &vertices, const XYZ_LOC &new_vertex)
[120]100{
[168]101        Part *p = new Part;
102        p->p.x = new_vertex.x;
103        p->p.y = new_vertex.y;
104        p->p.z = new_vertex.z;
[120]105        m.addPart(p);
106
107        vertices.push_back(new_vertex);
[1108]108        return int(vertices.size()) - 1;
[120]109}
110
[168]111double mix(int *colortab, int maxind, double ind)
[120]112{
[168]113        int indpre = (int)ind;
114        int indpost = indpre + 1;
115        if (indpost > maxind) indpost = maxind;
116        int v1 = colortab[indpre];
117        int v2 = colortab[indpost];
118        double d1 = ind - indpre;
119        double d2 = indpost - ind;
[973]120        double v = indpre == indpost ? v1 : d2 * v1 + d1 * v2; //d1+d2==1
[120]121        return v;
122}
123
[664]124void GenoConv_f90::setColors(Model &m, int last_added_part) //sets fixed (independent from genes) colors and widths on a model, purely for aesthetic purposes
[120]125{
126        //a rainbow on Joints: from the first one red, through middle green, to blue or violet - last
[168]127        static int r[] = { 1, 1, 0, 0, 0, 1 };
128        static int g[] = { 0, 1, 1, 1, 0, 0 };
129        static int b[] = { 0, 0, 0, 1, 1, 1 };
[1130]130        int maxind = int(std::size(r)) - 1;
[120]131
[168]132        int joints_count = m.getJointCount();
133        for (int i = 0; i < joints_count; i++)
[120]134        {
[168]135                Joint *j = m.getJoint(i);
136                double x = joints_count < 2 ? 0 : (double)i / (joints_count - 1); //0..1, postion in the rainbow
[973]137                double ind = x * maxind;
[168]138                j->vcolor.x = mix(r, maxind, ind);
139                j->vcolor.y = mix(g, maxind, ind);
140                j->vcolor.z = mix(b, maxind, ind);
[120]141        }
142
[168]143        int parts_count = m.getPartCount();
[120]144        SList jlist;
[259]145        for (int i = 0; i < parts_count; i++)
[120]146        {
[168]147                Part *p = m.getPart(i);
[120]148                jlist.clear();
[168]149                int count = m.findJoints(jlist, p);
150                Pt3D averagecolor(0, 0, 0); //Parts will get averaged colors from all attached Joints
151                FOREACH(Joint*, j, jlist)
152                        averagecolor += j->vcolor;
153                p->vcolor = averagecolor / count;
[120]154        }
[664]155        //m.getPart(0)->vcolor = Pt3D(0, 0, 0); //mark first Part black - a visual aid for easier editing
156        m.getPart(last_added_part)->vcolor = Pt3D(1, 1, 1); //mark last Part white - a visual aid for easier editing
[120]157}
158
[139]159void GenoConv_f90::perturbPartLocations(Model &m) //deterministic "body noise", see APPLY_DETERMINISTIC_BODY_NOISE
[120]160{
[168]161        for (int i = 0; i < m.getPartCount(); i++)
[120]162        {
[168]163                Part *p = m.getPart(i);
[120]164                Pt3D noise(
[168]165                        ((i + 1) % 10) - 4.5,
166                        ((3 * i + 5) % 10) - 4.5,
167                        ((7 * i + 2) % 10) - 4.5
[973]168                ); //-4.5 .. 4.5 in each axis
[168]169                p->p += noise / 1000;
[120]170        }
171}
Note: See TracBrowser for help on using the repository browser.