source: cpp/frams/genetics/f9/oper_f9.cpp @ 140

// This file is a part of the Framsticks GDK.
// Copyright (C) 2002-2014  Maciej Komosinski and Szymon Ulatowski.  See LICENSE.txt for details.
// Refer to http://www.framsticks.com/ for further information.

#include "oper_f9.h"
#include <frams/genetics/f9/conv_f9.h>
#include <common/nonstd.h> //randomN, rnd01


#define FIELDSTRUCT GenoOper_f9
static ParamEntry GENOf9param_tab[]=
{
        {"Genetics: f9",1,1,},
        {"f9_mut",0,0,"Mutation probability","f 0 1 0.1",FIELD(mut_prob),"How many genes should be mutated during single mutation (1=all genes, 0.1=ten percent)",},
        {0,},
};
#undef FIELDSTRUCT


GenoOper_f9::GenoOper_f9()
{
        par.setParamTab(GENOf9param_tab);
        par.select(this);
        par.setDefault();
        supported_format='9';
}

int GenoOper_f9::checkValidity(const char* gene)
{
        if (!gene[0]) return 1; //empty is not valid
        bool ok=true;
        int i;
        for(i=0;i<strlen(gene);i++) if (!strchr(turtle_commands_f9,gene[i])) {ok=false; break;}
        return ok ? GENOPER_OK : i+1;
}

///Remove all invalid letters from the genotype
int GenoOper_f9::validate(char *&gene)
{
        SString validated; //new genotype (everything except turtle_commands_f9 is skipped)
        for(int i=0;i<strlen(gene);i++)
                if (strchr(turtle_commands_f9,gene[i])) validated+=gene[i];  //validated contains only turtle_commands_f9
        free(gene);
        gene=strdup(validated); //reallocate
        return GENOPER_OK;
}

///Very simple mutation
int GenoOper_f9::mutate(char *&gene,float &chg,int &method)
{
        method=0;
        int changes=0,len=strlen(gene);
        int symbols=strlen(turtle_commands_f9);

        for(int i=0;i<len;i++)
        {
                if (rnd01<mut_prob) //normalize prob with the length of the genotype
                {
                        char oldgene=gene[i];
                        gene[i]=turtle_commands_f9[randomN(symbols)];
                        if (gene[i]!=oldgene) changes++;
                }
        }

        if (rnd01<mut_prob) //add or delete a random char
        {
                SString newgeno(gene);
                if (randomN(2)==0) //add
                {
                        int symbols=strlen(turtle_commands_f9);
                        int p=randomN(len+1);  //random location
                        //printf("before add: %s\n",(const char*)newgeno);
                        newgeno=newgeno.substr(0,p)+SString(turtle_commands_f9+randomN(symbols),1)+newgeno.substr(p);
                        //printf("after add: %s\n",(const char*)newgeno);
                        changes++;
                } else if (len>1) //delete
                {
                        int p=randomN(len);  //random location
                        //printf("before delete: %s\n",(const char*)newgeno);
                        newgeno=newgeno.substr(0,p)+newgeno.substr(p+1);
                        //printf("after delete: %s\n",(const char*)newgeno);
                        changes++;
                }
                free(gene);
                gene=strdup(newgeno); //reallocate
        }

        chg=(float)changes/len;
        return changes>0?GENOPER_OK:GENOPER_OPFAIL; //no changes => OPFAIL so that genman will call mutate again
}

///A simple one-point crossover
int GenoOper_f9::crossOver(char *&g1,char *&g2,float& chg1,float& chg2)
{
        int len1=strlen(g1),len2=strlen(g2);
        int p1=randomN(len1);  //random cut point for first genotype
        int p2=randomN(len2);  //random cut point for second genotype
        char *child1=(char*)malloc(p1+len2-p2+1);
        char *child2=(char*)malloc(p2+len1-p1+1);
        strncpy(child1,g1,p1);   strcpy(child1+p1,g2+p2);
        strncpy(child2,g2,p2);   strcpy(child2+p2,g1+p1);
        free(g1); g1=child1;
        free(g2); g2=child2;
        chg1=(float)p1/strlen(child1);
        chg2=(float)p2/strlen(child2);
        return GENOPER_OK;
}

///Applying some colors and font styles...
unsigned long GenoOper_f9::style(const char *g, int pos)
{
        char ch=g[pos];
        unsigned long style=GENSTYLE_CS(0,GENSTYLE_INVALID); //default, should be changed below
        char *ptr=strchr((char*)turtle_commands_f9,ch);
        if (ptr)
        {
                int pos=ptr-turtle_commands_f9;
                int axis=pos/2;
                style=GENSTYLE_RGBS(axis==0?200:0,axis==1?200:0,axis==2?200:0,GENSTYLE_NONE);
        }
        return style;
}