source: cpp/frams/genetics/genoconv.cpp @ 201

Last change on this file since 201 was 197, checked in by Maciej Komosinski, 10 years ago

GDK used by developers since 1999, distributed on the web since 2002

  • Property svn:eol-style set to native
File size: 5.0 KB
Line 
1// This file is a part of the Framsticks GDK.
2// Copyright (C) 1999-2014  Maciej Komosinski and Szymon Ulatowski.  See LICENSE.txt for details.
3// Refer to http://www.framsticks.com/ for further information.
4
5#include <common/nonstd.h>
6
7#include <stdlib.h>
8#include <math.h>
9#include <stdio.h>
10#include <string.h>
11#include <ctype.h>
12#include <time.h>
13#include <errno.h>
14
15#include "genoconv.h"
16#include <frams/util/multimap.h>
17#include <common/stl-util.h>
18
19///////////////////////////////////////////////////////////////////////////
20
21GenoConvParam::GenoConvParam(GenoConvManager *g) :Param(0), gcm(g)
22{
23        updatetab();
24}
25
26void GenoConvParam::freetab()
27{
28        if (tab) free(tab);
29        tab = 0;
30}
31
32const char *GenoConvParam::id(int i)
33{
34        if (i >= gcm->converters.size()) return 0;
35        sprintf(tmp_id, "genkonw%d", i);
36        return tmp_id;
37}
38
39void GenoConvParam::updatetab()
40{
41        int i;
42        GenoConverter *gk;
43        ParamEntry *pe;
44        int ile = gcm->converters.size();
45        freetab();
46        tab = (ParamEntry*)calloc(2 + ile, sizeof(ParamEntry));
47        tab[0].id = "Genetics: Conversions";
48        tab[0].group = 1;
49        tab[0].flags = (short)ile;
50        tab[0].name = "gkparam:";
51        gcnames.clear();
52        for (i = 0, pe = tab + 1; gk = (GenoConverter *)gcm->converters(i); pe++, i++)
53        {
54                pe->id = "?";
55                pe->group = 0;
56                pe->flags = 0;
57                std::string descr = "f";
58                descr += gk->in_format;
59                descr += " -> f";
60                descr += gk->out_format;
61                descr += " : ";
62                descr += gk->name;
63                gcnames.push_back(descr);
64                pe->name = descr.c_str();
65                pe->type = "d 0 1";
66                pe->help = gk->info;
67        }
68        pe->id = 0;
69}
70
71GenoConvParam::~GenoConvParam()
72{
73        freetab();
74}
75
76void *GenoConvParam::getTarget(int i)
77{
78        GenoConverter *gk = (GenoConverter *)gcm->converters(i);
79        return &gk->enabled;
80}
81
82GenoConvManager::GenoConvManager()
83:param(this)
84{
85}
86
87GenoConvManager::~GenoConvManager()
88{
89        FOREACH(GenoConverter*, gc, converters) delete gc;
90}
91
92void GenoConvManager::addConverter(GenoConverter *gc)
93{
94        converters += gc;
95        param.updatetab();
96}
97void GenoConvManager::removeConverter(GenoConverter *gc)
98{
99        converters -= gc;
100        param.updatetab();
101}
102
103GenoConverter *GenoConvManager::findConverters(SListTempl<GenoConverter*>* result, char in, char out, int enabled, char* name)
104{
105        GenoConverter *gk, *retval = 0;
106        int i = 0;
107        for (; gk = (GenoConverter*)converters(i); i++)
108        {
109                if ((in != -1) && (in != gk->in_format)) continue;
110                if ((out != -1) && (out != gk->out_format)) continue;
111                if ((enabled != -1) && (enabled != gk->enabled)) continue;
112                if ((name) && (strcmp(name, gk->name))) continue;
113                if (!retval) { retval = gk; if (!result) break; }
114                if (result) result->append(gk);
115        }
116        return retval;
117}
118
119/// Writes path into 'path'.
120/// return the last path element (return >= path)
121/// null -> path not found
122/// @param mapavailable will receive 1 if conversion map is supported by all converters in path
123/// (can be NULL if you don't need this information)
124
125char *GenoConvManager::getPath(char in, char out, char *path, int maxlen, int *mapavailable)
126{
127        if (!maxlen) return 0;
128        GenoConverter *gk;
129        int i = 0;
130        for (; gk = (GenoConverter*)converters(i); i++)
131        {
132                if ((gk->enabled) && (gk->in_format == in))
133                {
134                        *path = (char)i;
135                        if (gk->out_format == out)
136                        {
137                                if (mapavailable)
138                                        *mapavailable = gk->mapsupport;
139                                return path;
140                        }
141                        else
142                        {
143                                int mapavail;
144                                char *ret = getPath(gk->out_format, out, path + 1, maxlen - 1, &mapavail);
145                                if (ret)
146                                {
147                                        if (mapavailable)
148                                                *mapavailable = gk->mapsupport && mapavail;
149                                        return ret;
150                                }
151                        }
152                }
153        }
154        return 0;
155}
156
157char *GenoConvManager::getFormatPath(char in, char out, char *path, int maxlen, int *mapavailable)
158{
159        char *ret = getPath(in, out, path, maxlen, mapavailable);
160        if (ret)
161        {
162                for (char*t = path; t <= ret; t++)
163                        *t = ((GenoConverter*)converters(*t))->out_format;
164        }
165        return ret;
166}
167
168Geno GenoConvManager::convert(Geno &in, char format, MultiMap *map, bool *converter_missing)
169{
170        if (in.getFormat() == format) { if (converter_missing) *converter_missing = false; return in; }
171        char path[10];
172        int dep;
173        char *ret;
174        if (in.isInvalid()) { if (converter_missing) *converter_missing = false; return Geno("", 0, "", "invalid genotype cannot be converted"); }
175        int mapavail;
176        for (dep = 1; dep < (int)sizeof(path); dep++) //iterative deepening
177                if (ret = getPath(in.getFormat(), format, path, dep, &mapavail)) break;
178        if (!ret) { if (converter_missing) *converter_missing = true; return Geno("", 0, "", "converter not found"); }
179        if (converter_missing) *converter_missing = false;
180        if (!map) mapavail = 0;
181        char *t = path;
182        SString tmp;
183        tmp = in.getGene();
184        MultiMap lastmap, tmpmap;
185        int firstmap = 1;
186        for (; t <= ret; t++)
187        {
188                GenoConverter *gk = (GenoConverter*)converters(*t);
189                tmp = gk->convert(tmp, mapavail ? &tmpmap : 0);
190                if (!tmp.len())
191                {
192                        string t=ssprintf("f%c->f%c conversion failed (%s)", gk->in_format, gk->out_format, gk->name);
193                        return Geno(0, 0, 0, t.c_str());
194                }
195                if (mapavail)
196                {
197                        if (firstmap)
198                        {
199                                lastmap = tmpmap;
200                                firstmap = 0;
201                        }
202                        else
203                        {
204                                MultiMap m;
205                                m.addCombined(lastmap, tmpmap);
206                                lastmap = m;
207                        }
208                        tmpmap.clear();
209                }
210        }
211        if (map)
212                *map = lastmap;
213        return Geno(tmp, format, in.getName(), in.getComment());
214}
Note: See TracBrowser for help on using the repository browser.