Changeset 802 for cpp/frams/genetics/fB/fB_general.h

Timestamp:

06/07/18 17:41:59 (7 years ago)

Author:

Maciej Komosinski

Message:

Crossing over with less bloat, but still biologically-inspired

File:

• cpp/frams/genetics/fB/fB_general.h (modified) (3 diffs)

cpp/frams/genetics/fB/fB_general.h

@@ -11 +11 @@
 {
 public:
-        static int geneCount(SString geno)
+        static int geneCount(const SString& geno)
         {
                 int start = 0;
@@ -27 +27 @@
         }
 
-        static SString getGene(int i, SString genotype, int &start, int &end)
+        static SString getNonNestedGene(int i, const SString& genotype, int &start, int &end)
+        {
+                int count = -1;
+                start = 0;
+                int tmp = 0;
+                SString result = "";
+                do {
+                        count++;
+                        if (start < genotype.len())
+                                result = getNextGene(start, genotype, tmp, start);
+                        else
+                                start = -1;
+                } while (start != -1 && count < i);
+                start = tmp;
+                end = start;
+                return result;
+        }
+
+        static int geneCountNoNested(const SString& geno)
+        {
+                int start = 0;
+                int count = -1;
+                int tmp = 0;
+                do {
+                        count++;
+                        if (start < geno.len())
+                                getNextGene(start, geno, tmp, start);
+                        else
+                                start = -1;
+                } while (start != -1);
+                return count;
+        }
+
+        static SString getGene(int i, const SString& genotype, int &start, int &end)
         {
                 start = 0;
@@ -62 +95 @@
         }
 
+        static SString getNextGene(int searchbegin, const SString& genotype, int &start, int &end)
+        {
+                start = searchbegin;
+                int count = -1;
+                do {
+                        count++;
+                        start = genotype.indexOf("aa", start) + 1;
+                        int quotecount = 0;
+                        for (int q = 0; q < start; q++) if (genotype[q] == '\"') quotecount++;
+                        if (quotecount % 2 != 0) count--; // 'aa' sequence is within quotes
+                } while (start - 1 != -1 && count != 0);
+                if (start - 1 == -1) // there is no gene with a given "i"
+                {
+                        start = -1;
+                        end = -1;
+                        return "";
+                }
+                end = start;
+                int quotecount = 0;
+                do {
+                        quotecount = 0;
+                        end = genotype.indexOf("zz", end);
+                        if (end != -1)
+                        {
+                                for (int q = start; q < end; q++) if (genotype[q] == '\"') quotecount++;
+                                if (quotecount % 2 != 0) end++;
+                        }
+                } while (quotecount % 2 != 0 && end != -1);
+
+                if (end == -1) end = genotype.len();
+                else end += 2;
+                start -= 1;
+                return genotype.substr(start, end - start);
+        }
 private:
         fB_GenoHelpers() {}