source: cpp/frams/model/model.cpp @ 647

Last change on this file since 647 was 644, checked in by Maciej Komosinski, 7 years ago

Invalid references from neurons to Parts/Joints? are now reported
  • Property svn:eol-style set to native
File size: 33.6 KB
Line 
1// This file is a part of Framsticks SDK.  http://www.framsticks.com/
2// Copyright (C) 1999-2015  Maciej Komosinski and Szymon Ulatowski.
3// See LICENSE.txt for details.
4
5#include <common/nonstd_math.h>
6#include "model.h"
7#include <common/log.h>
8#include <frams/util/multimap.h>
9#include <common/loggers/loggers.h>
10
11Model::Model()
12{
13        autobuildmaps = false;
14        init();
15}
16
17void Model::init()
18{
19        partmappingchanged = 0;
20        buildstatus = empty;
21        modelfromgenotype = 0;
22        startenergy = 1.0;
23        checklevel = 1;
24#ifdef MODEL_V1_COMPATIBLE
25        oldneurocount=-1; // == unknown
26        oldconnections=0;
27#endif
28        map = 0;
29        f0map = 0;
30        f0genoknown = 1;
31        shape = SHAPE_UNKNOWN;
32}
33
34void Model::moveElementsFrom(Model &source)
35{
36        int i;
37        open();
38        for (i = 0; i < source.getPartCount(); i++)
39                addPart(source.getPart(i));
40        for (i = 0; i < source.getJointCount(); i++)
41                addJoint(source.getJoint(i));
42        for (i = 0; i < source.getNeuroCount(); i++)
43                addNeuro(source.getNeuro(i));
44        source.parts.clear(); source.joints.clear(); source.neurons.clear();
45        source.clear();
46}
47
48void Model::internalCopy(const Model &mod)
49{
50        geno = mod.geno;
51        f0genoknown = 0;
52        startenergy = mod.startenergy;
53        if (mod.getStatus() == valid)
54        {
55                modelfromgenotype = mod.modelfromgenotype;
56                {for (int i = 0; i < mod.getPartCount(); i++)
57                        addPart(new Part(*mod.getPart(i))); }
58        {for (int i = 0; i < mod.getJointCount(); i++)
59        {
60                Joint *oldj = mod.getJoint(i);
61                Joint *j = new Joint(*oldj);
62                addJoint(j);
63                j->attachToParts(oldj->part1->refno, oldj->part2->refno);
64        }}
65        {for (int i = 0; i < mod.getNeuroCount(); i++)
66        {
67                Neuro *oldn = mod.getNeuro(i);
68                Neuro *n = new Neuro(*oldn);
69                addNeuro(n);
70                if (oldn->part_refno >= 0) n->attachToPart(oldn->part_refno);
71                else n->attachToJoint(oldn->joint_refno);
72        }}
73                for (int i = 0; i < mod.getNeuroCount(); i++)
74                {
75                        Neuro *oldn = mod.getNeuro(i);
76                        Neuro *n = getNeuro(i);
77                        for (int ni = 0; ni < oldn->getInputCount(); ni++)
78                        {
79                                double w;
80                                Neuro *oldinput = oldn->getInput(ni, w);
81                                SString info = n->getInputInfo(ni);
82                                n->addInput(getNeuro(oldinput->refno), w, &info);
83                        }
84                }
85        }
86}
87
88
89Model::Model(const Geno &src, bool buildmaps)
90        :autobuildmaps(buildmaps)
91{
92        init(src);
93}
94
95void Model::operator=(const Model &mod)
96{
97        clear();
98        open();
99        internalCopy(mod);
100        buildstatus = mod.buildstatus;
101}
102
103Model::Model(const Model &mod, bool buildmaps)
104        :autobuildmaps(buildmaps)
105{
106        init();
107        open();
108        internalCopy(mod);
109        buildstatus = mod.buildstatus;
110}
111
112void Model::init(const Geno &src)
113{
114        init();
115        modelfromgenotype = 1;
116        geno = src;
117        build();
118}
119
120void Model::resetAllDelta()
121{
122        for (int i = 0; i < getJointCount(); i++)
123                getJoint(i)->resetDelta();
124}
125
126void Model::useAllDelta(bool yesno)
127{
128        for (int i = 0; i < getJointCount(); i++)
129                getJoint(i)->useDelta(yesno);
130}
131
132Model::~Model()
133{
134        delmodel_list.action((intptr_t)this);
135        clear();
136}
137
138void Model::clear()
139{
140        FOREACH(Part*, p, parts)
141                delete p;
142        FOREACH(Joint*, j, joints)
143                delete j;
144        FOREACH(Neuro*, n, neurons)
145                delete n;
146        parts.clear(); joints.clear(); neurons.clear();
147        delMap();
148        delF0Map();
149        init();
150        geno = Geno();
151        f0geno = Geno();
152}
153
154Part *Model::addPart(Part *p)
155{
156        p->owner = this;
157        p->refno = parts.size();
158        parts += p;
159        return p;
160}
161
162Joint *Model::addJoint(Joint *j)
163{
164        j->owner = this;
165        j->refno = joints.size();
166        joints += j;
167        return j;
168}
169
170Neuro *Model::addNeuro(Neuro *n)
171{
172        n->owner = this;
173        n->refno = neurons.size();
174        neurons += n;
175        return n;
176}
177
178void Model::removeNeuros(SList &nlist)
179{
180        FOREACH(Neuro*, nu, nlist)
181        {
182                int i = findNeuro(nu);
183                if (i >= 0) removeNeuro(i);
184        }
185}
186
187void Model::removePart(int partindex, int removeattachedjoints, int removeattachedneurons)
188{
189        Part *p = getPart(partindex);
190        if (removeattachedjoints)
191        {
192                SList jlist;
193                findJoints(jlist, p);
194                FOREACH(Joint*, j, jlist)
195                {
196                        int i = findJoint(j);
197                        if (i >= 0) removeJoint(i, removeattachedneurons);
198                }
199        }
200        if (removeattachedneurons)
201        {
202                SList nlist;
203                findNeuros(nlist, 0, p);
204                removeNeuros(nlist);
205        }
206        parts -= partindex;
207        delete p;
208}
209
210void Model::removeJoint(int jointindex, int removeattachedneurons)
211{
212        Joint *j = getJoint(jointindex);
213        if (removeattachedneurons)
214        {
215                SList nlist;
216                findNeuros(nlist, 0, 0, j);
217                removeNeuros(nlist);
218        }
219        joints -= jointindex;
220        delete j;
221}
222
223void Model::removeNeuro(int neuroindex, bool removereferences)
224{
225        Neuro* thisN = getNeuro(neuroindex);
226
227        if (removereferences)
228        {
229                Neuro* n;
230                // remove all references to thisN
231                for (int i = 0; n = (Neuro*)neurons(i); i++)
232                {
233                        Neuro *inp;
234                        for (int j = 0; inp = n->getInput(j); j++)
235                                if (inp == thisN)
236                                {
237                                n->removeInput(j);
238                                j--;
239                                }
240                }
241        }
242
243        neurons -= neuroindex;
244        delete thisN;
245}
246
247MultiMap& Model::getMap()
248{
249        if (!map) map = new MultiMap();
250        return *map;
251}
252
253void Model::delMap()
254{
255        if (map) { delete map; map = 0; }
256}
257void Model::delF0Map()
258{
259        if (f0map) { delete f0map; f0map = 0; }
260}
261
262void Model::makeGenToGenMap(MultiMap& result, const MultiMap& gen1tomodel, const MultiMap& gen2tomodel)
263{
264        result.clear();
265        MultiMap m;
266        m.addReversed(gen2tomodel);
267        result.addCombined(gen1tomodel, m);
268}
269
270void Model::getCurrentToF0Map(MultiMap& result)
271{
272        result.clear();
273        if (!map) return;
274        const MultiMap& f0m = getF0Map();
275        makeGenToGenMap(result, *map, f0m);
276}
277
278void Model::rebuild(bool buildm)
279{
280        autobuildmaps = buildm;
281        clear();
282        build();
283}
284
285void Model::initMap()
286{
287        if (!map) map = new MultiMap();
288        else map->clear();
289}
290
291void Model::initF0Map()
292{
293        if (!f0map) f0map = new MultiMap();
294        else f0map->clear();
295}
296
297void Model::build()
298{
299        f0errorposition = -1;
300        f0warnposition = -1;
301        MultiMap *convmap = autobuildmaps ? new MultiMap() : NULL;
302        f0geno = (geno.getFormat() == '0') ? geno : geno.getConverted('0', convmap);
303        f0genoknown = 1;
304        if (f0geno.isInvalid())
305        {
306                buildstatus = invalid;
307                if (convmap) delete convmap;
308                return;
309        }
310        SString f0txt = f0geno.getGenes();
311        buildstatus = building; // was: open();
312        if (autobuildmaps)
313        {
314                partmappingchanged = 0;
315                initMap();
316                initF0Map();
317        }
318        int pos = 0, lnum = 1, lastpos = 0;
319        SString line;
320        MultiRange frommap;
321        LoggerToMemory mh(LoggerBase::Enable | LoggerBase::DontBlock);
322        for (; f0txt.getNextToken(pos, line, '\n'); lnum++)
323        {
324                if (autobuildmaps)
325                {
326                        frommap.clear();
327                        frommap.add(lastpos, pos - 1);
328                }
329                mh.reset();
330                if (singleStepBuild(line, lnum, autobuildmaps ? (&frommap) : 0) == -1)
331                {
332                        buildstatus = invalid;
333                        f0errorposition = lastpos;
334                        if (convmap) delete convmap;
335                        return;
336                }
337                if (mh.getWarningCount())
338                {
339                        if (f0warnposition < 0) f0warnposition = lastpos;
340                }
341                lastpos = pos;
342        }
343        mh.disable();
344        close();
345        if (convmap)
346        {
347                *f0map = *map;
348                if (geno.getFormat() != '0')
349                {
350                        MultiMap tmp;
351                        tmp.addCombined(*convmap, getMap());
352                        *map = tmp;
353                }
354                delete convmap;
355        }
356}
357
358const MultiMap& Model::getF0Map()
359{
360        if (!f0map)
361        {
362                f0map = new MultiMap();
363                makeGeno(f0geno, f0map);
364                f0genoknown = 1;
365        }
366        return *f0map;
367}
368
369Geno Model::rawGeno()
370{
371        Geno tmpgen;
372        makeGeno(tmpgen);
373        return tmpgen;
374}
375
376void Model::makeGeno(Geno &g, MultiMap *map, bool handle_defaults)
377{
378        if ((buildstatus != valid) && (buildstatus != building))
379        {
380                g = Geno(0, 0, 0, "invalid model");
381                return;
382        }
383
384        SString gen;
385
386        Param modelparam(f0_model_paramtab);
387        Param partparam(f0_part_paramtab);
388        Param jointparam(f0_joint_paramtab);
389        Param neuroparam(f0_neuro_paramtab);
390        Param connparam(f0_neuroconn_paramtab);
391
392        static Part defaultpart;
393        static Joint defaultjoint;
394        static Neuro defaultneuro;
395        static Model defaultmodel;
396        static NeuroConn defaultconn;
397        //static NeuroItem defaultneuroitem;
398
399        Part *p;
400        Joint *j;
401        Neuro *n;
402        int i;
403        int len;
404        int a, b;
405        //NeuroItem *ni;
406
407        SString mod_props;
408        modelparam.select(this);
409        modelparam.save2(mod_props, handle_defaults ? &defaultmodel : NULL, true, !handle_defaults);
410        if (mod_props.len() > 1) //are there any non-default values? ("\n" is empty)
411        {
412                gen += "m:";
413                gen += mod_props;
414        }
415
416        for (i = 0; p = (Part*)parts(i); i++)
417        {
418                partparam.select(p);
419                len = gen.len();
420                gen += "p:";
421                partparam.save2(gen, handle_defaults ? &defaultpart : NULL, true, !handle_defaults);
422                if (map)
423                        map->add(len, gen.len() - 1, partToMap(i));
424        }
425        for (i = 0; j = (Joint*)joints(i); i++)
426        {
427                jointparam.select(j);
428                len = gen.len();
429                jointparam.setParamTab(j->usedelta ? f0_joint_paramtab : f0_nodeltajoint_paramtab);
430                gen += "j:";
431                jointparam.save2(gen, handle_defaults ? &defaultjoint : NULL, true, !handle_defaults);
432                if (map)
433                        map->add(len, gen.len() - 1, jointToMap(i));
434        }
435        for (i = 0; n = (Neuro*)neurons(i); i++)
436        {
437                neuroparam.select(n);
438                len = gen.len();
439                gen += "n:";
440                neuroparam.save2(gen, handle_defaults ? &defaultneuro : NULL, true, !handle_defaults);
441                if (map)
442                        map->add(len, gen.len() - 1, neuroToMap(i));
443        }
444        for (a = 0; a < neurons.size(); a++)
445        { // inputs
446                n = (Neuro*)neurons(a);
447                //      if ((n->getInputCount()==1)&&(n->getInput(0).refno <= n->refno))
448                //              continue; // already done with Neuro::conn_refno
449
450                for (b = 0; b < n->getInputCount(); b++)
451                {
452                        double w;
453                        NeuroConn nc;
454                        Neuro* n2 = n->getInput(b, w);
455                        //              if (((n2.parentcount==1)&&(n2.parent)&&(n2.parent->refno < n2.refno)) ^
456                        //                  (n2.neuro_refno>=0))
457                        //                      printf("!!!! bad Neuro::neuro_refno ?!\n");
458
459                        //              if ((n2.parentcount==1)&&(n2.parent)&&(n2.parent->refno < n2.refno))
460                        //              if (n2.neuro_refno>=0)
461                        //                      continue; // already done with Neuro::neuro_refno
462
463                        nc.n1_refno = n->refno; nc.n2_refno = n2->refno;
464                        nc.weight = w;
465                        SString **s = n->inputInfo(b);
466                        if ((s) && (*s))
467                                nc.info = **s;
468                        connparam.select(&nc);
469                        len = gen.len();
470                        gen += "c:";
471                        connparam.save2(gen, handle_defaults ? &defaultconn : NULL, true, !handle_defaults);
472                        if (map)
473                                map->add(len, gen.len() - 1, neuroToMap(n->refno));
474                }
475        }
476        g = Geno(gen.c_str(), '0');
477}
478
479//////////////
480
481void Model::open()
482{
483        if (buildstatus == building) return;
484        buildstatus = building;
485        modelfromgenotype = 0;
486        partmappingchanged = 0;
487        f0genoknown = 0;
488        delMap();
489}
490
491void Model::checkpoint()
492{}
493
494void Model::setGeno(const Geno& newgeno)
495{
496        geno = newgeno;
497}
498
499void Model::clearMap()
500{
501        Part *p; Joint *j; Neuro *n;
502        int i;
503        delMap();
504        delF0Map();
505        for (i = 0; p = (Part*)parts(i); i++)
506                p->clearMapping();
507        for (i = 0; j = (Joint*)joints(i); i++)
508                j->clearMapping();
509        for (i = 0; n = (Neuro*)neurons(i); i++)
510                n->clearMapping();
511}
512
513int Model::close(bool building_live_model)
514{
515        if (buildstatus != building)
516                logPrintf("Model", "close", LOG_WARN, "Unexpected close() - no open()");
517        if (internalcheck(building_live_model ? LIVE_CHECK : FINAL_CHECK) > 0)
518        {
519                buildstatus = valid;
520
521                if (partmappingchanged)
522                {
523                        getMap();
524                        Part *p; Joint *j; Neuro *n;
525                        int i;
526                        for (i = 0; p = (Part*)parts(i); i++)
527                                if (p->getMapping())
528                                        map->add(*p->getMapping(), partToMap(i));
529                        for (i = 0; j = (Joint*)joints(i); i++)
530                                if (j->getMapping())
531                                        map->add(*j->getMapping(), jointToMap(i));
532                        for (i = 0; n = (Neuro*)neurons(i); i++)
533                                if (n->getMapping())
534                                        map->add(*n->getMapping(), neuroToMap(i));
535                }
536        }
537        else
538                buildstatus = invalid;
539
540        return (buildstatus == valid);
541}
542
543int Model::validate()
544{
545        return internalcheck(EDITING_CHECK);
546}
547
548Pt3D Model::whereDelta(const Part& start, const Pt3D& rot, const Pt3D& delta)
549{
550        Orient roto;
551        roto = rot;
552        Orient o;
553        roto.transform(o, start.o);
554        //o.x=start.o/roto.x;
555        //o.y=start.o/roto.y;
556        //o.z=start.o/roto.z;
557        return o.transform(delta) + start.p;
558}
559
560int Model::singleStepBuild(const SString &singleline, const MultiRange* srcrange)
561{
562        return singleStepBuild(singleline, 0, srcrange);
563}
564
565int Model::singleStepBuild(const SString &singleline, int line_num, const MultiRange* srcrange)
566{
567        SString error_message;
568        int result = singleStepBuildNoLog(singleline, error_message, srcrange);
569        if (result < 0)
570        {
571                if (error_message.len() == 0) // generic error when no detailed message is available
572                        error_message = "Invalid f0 code";
573                if (line_num>0)
574                        error_message += SString::sprintf(", line #%d", line_num);
575                error_message += nameForErrors();
576                logPrintf("Model", "build", LOG_ERROR, "%s", error_message.c_str());
577        }
578        return result;
579}
580
581int Model::singleStepBuildNoLog(const SString &line, SString& error_message, const MultiRange* srcrange)
582{
583        error_message = SString::empty();
584        int pos = 0; const char*t = line.c_str();
585        for (; *t; t++, pos++)
586                if (!strchr(" \r\t", *t)) break;
587        if (*t == '#') return 0;
588        if (!*t) return 0;
589        if (!strncmp(t, "p:", 2))
590        {
591                Param partparam(f0_part_paramtab);
592                Part *p = new Part();
593                partparam.select(p);
594                pos += 2;
595                if (partparam.load2(line, pos) & ParamInterface::LOAD2_PARSE_FAILED) { delete p; error_message = "Invalid 'p:'"; return -1; }
596                p->o.rotate(p->rot);
597                parts += p;
598                p->owner = this;
599                if (srcrange) p->setMapping(*srcrange);
600                return getPartCount() - 1;
601        }
602        if (!strncmp(t, "m:", 2))
603        {
604                Param modelparam(f0_model_paramtab);
605                modelparam.select(this);
606                pos += 2;
607                if (modelparam.load2(line, pos) & ParamInterface::LOAD2_PARSE_FAILED) { error_message = "Invalid 'm:'"; return -1; }
608                return 0;
609        }
610        else if (!strncmp(t, "j:", 2))
611        {
612                Param jointparam(f0_joint_paramtab);
613                Joint *j = new Joint();
614                jointparam.select(j);
615                pos += 2;
616                j->owner = this;
617                if (jointparam.load2(line, pos) & ParamInterface::LOAD2_PARSE_FAILED) { delete j; error_message = "Invalid 'j:'"; return -1; }
618                bool p1_ok = false, p2_ok = false;
619                if ((p1_ok = ((j->p1_refno >= 0) && (j->p1_refno < getPartCount()))) &&
620                        (p2_ok = ((j->p2_refno >= 0) && (j->p2_refno < getPartCount()))))
621                {
622                        addJoint(j);
623                        if ((j->d.x != JOINT_DELTA_MARKER) || (j->d.y != JOINT_DELTA_MARKER) || (j->d.z != JOINT_DELTA_MARKER))
624                        {
625                                j->useDelta(1);
626                                j->resetDeltaMarkers();
627                        }
628                        j->attachToParts(j->p1_refno, j->p2_refno);
629                        if (srcrange) j->setMapping(*srcrange);
630                        return j->refno;
631                }
632                else
633                {
634                        error_message = SString::sprintf("Invalid reference to Part #%d", p1_ok ? j->p1_refno : j->p2_refno);
635                        delete j;
636                        return -1;
637                }
638        }
639        else if (!strncmp(t, "n:", 2)) // neuro (or the old neuro object as the special case)
640        {
641                Param neuroparam(f0_neuro_paramtab);
642                Neuro *nu = new Neuro();
643                neuroparam.select(nu);
644                pos += 2;
645                if (neuroparam.load2(line, pos) & ParamInterface::LOAD2_PARSE_FAILED) { delete nu; error_message = "Invalid 'n:'"; return -1; }
646#ifdef MODEL_V1_COMPATIBLE
647                if (nu->neuro_refno>=0) // parent specified...
648                {
649                        if (nu->neuro_refno >= getNeuroCount()) // and it's illegal
650                        {
651                                delete nu;
652                                return -1;
653                        }
654                        Neuro *parentNU=getNeuro(nu->neuro_refno);
655                        parentNU->addInput(nu,nu->weight);
656                        // default class for parented units: n-n link
657                        //if (nu->moredata.len()==0) nu->moredata="-";
658                }
659                else
660#endif
661                {
662                        // default class for unparented units: standard neuron
663                        if (nu->getClassName().len() == 0) nu->setClassName("N");
664                }
665                /*
666                        if (nu->conn_refno>=0) // input specified...
667                        {
668                        if (nu->conn_refno >= getNeuroCount()) // and it's illegal
669                        {
670                        delete nu;
671                        return -1;
672                        }
673                        Neuro *inputNU=getNeuro(nu->conn_refno);
674                        nu->addInput(inputNU,nu->weight);
675                        }
676                        */
677#ifdef MODEL_V1_COMPATIBLE
678                nu->weight=1.0;
679#endif
680                nu->owner = this;
681                // attach to part/joint
682                if (nu->part_refno >= 0)
683                        {
684                        nu->attachToPart(nu->part_refno);
685                        if (nu->part==NULL)
686                                { error_message = SString::sprintf("Invalid reference to Part #%d", nu->part_refno); delete nu; return -1; }
687                        }
688                if (nu->joint_refno >= 0)
689                        {
690                        nu->attachToJoint(nu->joint_refno);
691                        if (nu->joint==NULL)
692                                { error_message = SString::sprintf("Invalid reference to Joint #%d", nu->joint_refno); delete nu; return -1; }
693                        }
694                if (srcrange) nu->setMapping(*srcrange);
695                // todo: check part/joint ref#
696#ifdef MODEL_V1_COMPATIBLE
697                if (hasOldNeuroLayout())
698                {
699                        int count=old_getNeuroCount();
700                        neurons.insert(count,nu);
701                        oldneurocount=count+1;
702                        return oldneurocount-1;
703                }
704                else
705#endif
706                {
707                        neurons += nu;
708                        return neurons.size() - 1;
709                }
710        }
711        else if (!strncmp(t, "c:", 2)) // add input
712        {
713                Param ncparam(f0_neuroconn_paramtab);
714                NeuroConn c;
715                ncparam.select(&c);
716                pos += 2;
717                if (ncparam.load2(line, pos) & ParamInterface::LOAD2_PARSE_FAILED) { error_message = "Invalid 'c:'"; return -1; }
718                bool n1_ok = false, n2_ok = false;
719                if ((n1_ok = ((c.n1_refno >= 0) && (c.n1_refno < getNeuroCount())))
720                        && (n2_ok = ((c.n2_refno >= 0) && (c.n2_refno < getNeuroCount()))))
721                {
722                        Neuro *na = getNeuro(c.n1_refno);
723                        Neuro *nb = getNeuro(c.n2_refno);
724                        na->addInput(nb, c.weight, &c.info);
725                        if (srcrange)
726                                na->addMapping(*srcrange);
727                        return 0;
728                }
729                error_message = SString::sprintf("Invalid reference to Neuro #%d", n1_ok ? c.n2_refno : c.n1_refno);
730                return -1;
731        }
732#ifdef MODEL_V1_COMPATIBLE
733        else if (!strncmp(t,"ni:",3)) // old neuroitem
734        {
735                // we always use old layout for "ni:"
736                Param neuroitemparam(f0_neuroitem_paramtab);
737                Neuro *nu=new Neuro();
738                neuroitemparam.select(nu);
739                pos+=3;
740                if (neuroitemparam.load2(line,pos) & ParamInterface::LOAD2_PARSE_FAILED) {delete nu; return -1;}
741                // illegal parent?
742                if ((nu->neuro_refno<0)||(nu->neuro_refno>=old_getNeuroCount()))
743                {
744                        delete nu;
745                        return -1;
746                }
747                Neuro *parentN=getNeuro(nu->neuro_refno);
748                // copy part/joint refno from parent, if possible
749                if ((nu->part_refno<0)&&(parentN->part_refno>=0))
750                        nu->part_refno=parentN->part_refno;
751                if ((nu->joint_refno<0)&&(parentN->joint_refno>=0))
752                        nu->joint_refno=parentN->joint_refno;
753                nu->owner=this;
754                // attach to part/joint
755                if (nu->part_refno>=0)
756                        nu->attachToPart(nu->part_refno);
757                if (nu->joint_refno>=0)
758                        nu->attachToJoint(nu->joint_refno);
759                if (srcrange)
760                        nu->setMapping(*srcrange);
761                // special case: old muscles
762                // PARENT neuron will be set up to be the CHILD of the current one (!)
763                if (nu->isOldEffector())
764                {
765                        nu->neuro_refno=parentN->refno;
766                        neurons+=nu;
767                        nu->owner=this;
768                        nu->addInput(parentN,nu->weight); // (!)
769                        nu->weight=1.0;
770                        parentN->invalidateOldItems();
771                        return 0; // !!! -> ...getItemCount()-1;
772                }
773                parentN->addInput(nu,nu->weight);
774                neurons+=nu;
775                parentN->invalidateOldItems();
776                if (nu->getClassName().len()==0)
777                {
778                        nu->setClassName("-");
779                        // for compatibility, "ni" can define a n-n connection
780                        // referring to non-existent neuron (which will be hopefully defined later)
781                        // internal check will add the proper input to this unit
782                        // if conn_refno >=0 and input count==0
783                        oldconnections=1;
784                        if (srcrange)
785                                parentN->addMapping(*srcrange);
786                }
787                else
788                        nu->weight=1.0;
789                return 0; // !!! -> ...getItemCount()-1;
790        }
791#endif
792        else return -1;
793}
794
795#ifdef MODEL_V1_COMPATIBLE
796int Model::addOldConnectionsInputs()
797{
798        if (!oldconnections) return 1;
799        Neuro* n;
800        for (int i=0;i<neurons.size();i++)
801        {
802                n=(Neuro*)neurons(i);
803                if (n->conn_refno>=0)
804                        if (n->isNNConnection())
805                                if (n->conn_refno < old_getNeuroCount())
806                                { // good reference
807                        Neuro* parent=n->parent; // nn connection has exactly one parent
808                        int inp=parent->findInput(n);
809                        Neuro* target=getNeuro(n->conn_refno);
810                        parent->setInput(inp,target,n->weight);
811                        removeNeuro(i,0); // no need to remove references
812                        i--;
813                                }
814                                else
815                                {
816                                        logPrintf("Model","internalCheck",LOG_ERROR,
817                                                "illegal N-N connection #%d (reference to #%d)%s",
818                                                i,n->conn_refno,nameForErrors().c_str());
819                                        return 0;
820                                }
821        }
822        oldconnections=0;
823        return 1;
824}
825#endif
826
827/////////////
828
829/** change the sequence of neuro units
830        and fix references in "-" objects (n-n connections)  */
831void Model::moveNeuro(int oldpos, int newpos)
832{
833        if (oldpos == newpos) return; // nop!
834        Neuro *n = getNeuro(oldpos);
835        neurons -= oldpos;
836        neurons.insert(newpos, n);
837        // conn_refno could be broken -> fix it
838#ifdef MODEL_V1_COMPATIBLE
839        for (int i=0;i<neurons.size();i++)
840        {
841                Neuro *n2=getNeuro(i);
842                if (n2->isNNConnection())
843                        if (n2->conn_refno == oldpos) n2->conn_refno=newpos;
844                        else
845                        { if (n2->conn_refno > oldpos) n2->conn_refno--;
846                if (n2->conn_refno >= newpos) n2->conn_refno++; }
847        }
848        invalidateOldNeuroCount();
849#endif
850}
851
852#ifdef MODEL_V1_COMPATIBLE
853/** move all old neurons (class "N") to the front of the neurons list.
854        @return number of old neurons
855        */
856int Model::reorderToOldLayout()
857{
858        Neuro *n;
859        int neurocount=0;
860        for (int i=0;i<neurons.size();i++)
861        {
862                n=(Neuro*)neurons(i);
863                if (n->isOldNeuron())
864                {
865                        moveNeuro(i,neurocount);
866                        neurocount++;
867                        i=neurocount-1;
868                }
869        }
870        return neurocount;
871}
872#endif
873////////////
874
875void Model::updateNeuroRefno()
876{
877        for (int i = 0; i < neurons.size(); i++)
878        {
879                Neuro* n = (Neuro*)neurons(i);
880                n->refno = i;
881        }
882}
883
884#define VALIDMINMAX(var,template,field) \
885if (var -> field < getMin ## template () . field) \
886        { var->field= getMin ## template () . field; \
887        logPrintf("Model","internalCheck",LOG_WARN,# field " too small in " # template " #%d (adjusted)",i);} \
888else if (var -> field > getMax ## template () . field) \
889        { var->field= getMax ## template ()  . field; \
890        logPrintf("Model","internalCheck",LOG_WARN,# field " too big in " # template " #%d (adjusted)",i);}
891
892#define LINKFLAG 0x8000000
893
894void Model::updateRefno()
895{
896        int i;
897        for (i = 0; i < getPartCount(); i++) getPart(i)->refno = i;
898        for (i = 0; i < getJointCount(); i++) getJoint(i)->refno = i;
899        for (i = 0; i < getNeuroCount(); i++) getNeuro(i)->refno = i;
900}
901
902SString Model::nameForErrors() const
903{
904        if (geno.getName().len()>0)
905                return SString::sprintf(" in '%s'", geno.getName().c_str());
906        return SString::empty();
907}
908
909int Model::internalcheck(CheckType check)
910{
911        Part *p;
912        Joint *j;
913        Neuro *n;
914        int i, k;
915        int ret = 1;
916        shape = SHAPE_UNKNOWN;
917        if ((parts.size() == 0) && (neurons.size() == 0)) return 0;
918        if (parts.size() == 0)
919                size = Pt3D_0;
920        else
921        {
922                Pt3D bbmin = ((Part*)parts(0))->p, bbmax = bbmin;
923                for (i = 0; i < parts.size(); i++)
924                {
925                        p = (Part*)parts(i);
926                        p->owner = this;
927                        p->refno = i;
928                        if (checklevel > 0)
929                                p->mass = 0.0;
930                        //VALIDMINMAX(p,part,mass);//mass is very special
931                        // VALIDMINMAX are managed manually when adding part properties in f0-def!
932                        // (could be made dynamic but not really worth the effort)
933                        VALIDMINMAX(p, Part, size);
934                        VALIDMINMAX(p, Part, scale.x);
935                        VALIDMINMAX(p, Part, scale.y);
936                        VALIDMINMAX(p, Part, scale.z);
937                        VALIDMINMAX(p, Part, hollow);
938                        VALIDMINMAX(p, Part, density);
939                        VALIDMINMAX(p, Part, friction);
940                        VALIDMINMAX(p, Part, ingest);
941                        VALIDMINMAX(p, Part, assim);
942                        VALIDMINMAX(p, Part, vsize);
943                        VALIDMINMAX(p, Part, vcolor.x);
944                        VALIDMINMAX(p, Part, vcolor.y);
945                        VALIDMINMAX(p, Part, vcolor.z);
946                        p->flags &= ~LINKFLAG; // for delta joint cycle detection
947                        if (p->p.x - p->size < bbmin.x) bbmin.x = p->p.x - p->size;
948                        if (p->p.y - p->size < bbmin.y) bbmin.y = p->p.y - p->size;
949                        if (p->p.z - p->size < bbmin.z) bbmin.z = p->p.z - p->size;
950                        if (p->p.x + p->size > bbmax.x) bbmax.x = p->p.x + p->size;
951                        if (p->p.y + p->size > bbmax.y) bbmax.y = p->p.y + p->size;
952                        if (p->p.z + p->size > bbmax.z) bbmax.z = p->p.z + p->size;
953                        if (shape == SHAPE_UNKNOWN)
954                                shape = (p->shape == Part::SHAPE_BALL_AND_STICK) ? SHAPE_BALL_AND_STICK : SHAPE_SOLIDS;
955                        else if (shape != SHAPE_ILLEGAL)
956                        {
957                                if ((p->shape == Part::SHAPE_BALL_AND_STICK) ^ (shape == SHAPE_BALL_AND_STICK))
958                                {
959                                        shape = SHAPE_ILLEGAL;
960                                        logPrintf("Model", "internalCheck", LOG_WARN, "Inconsistent part shapes (mixed ball-and-stick and solids shape types)%s", nameForErrors().c_str());
961                                }
962                        }
963                }
964                size = bbmax - bbmin;
965                for (i = 0; i < joints.size(); i++)
966                {
967                        j = (Joint*)joints(i);
968                        // VALIDMINMAX are managed manually when adding joint properties in f0-def!
969                        // (could be made dynamic but not really worth the effort)
970                        VALIDMINMAX(j, Joint, stamina);
971                        VALIDMINMAX(j, Joint, stif);
972                        VALIDMINMAX(j, Joint, rotstif);
973                        VALIDMINMAX(p, Part, vcolor.x);
974                        VALIDMINMAX(p, Part, vcolor.y);
975                        VALIDMINMAX(p, Part, vcolor.z);
976                        j->refno = i;
977                        j->owner = this;
978                        if (j->part1 && j->part2 && (j->part1 != j->part2))
979                        {
980                                j->p1_refno = j->part1->refno;
981                                j->p2_refno = j->part2->refno;
982                                if (checklevel > 0)
983                                {
984                                        j->part1->mass += 1.0;
985                                        j->part2->mass += 1.0;
986                                }
987                                if ((j->usedelta) && ((j->d.x != JOINT_DELTA_MARKER) || (j->d.y != JOINT_DELTA_MARKER) || (j->d.z != JOINT_DELTA_MARKER)))
988                                { // delta positioning -> calc. orient.
989                                        if (j->part2->flags & LINKFLAG)
990                                        {
991                                                ret = 0;
992                                                logPrintf("Model", "internalCheck", LOG_ERROR,
993                                                        "Delta joint cycle detected at Joint #%d%s",
994                                                        i, nameForErrors().c_str());
995                                        }
996                                        j->resetDeltaMarkers();
997                                        j->o = j->rot;
998                                        j->part1->o.transform(j->part2->o, j->o);
999                                        //                      j->part2->o.x=j->part1->o/j->o.x;
1000                                        //                      j->part2->o.y=j->part1->o/j->o.y;
1001                                        //                      j->part2->o.z=j->part1->o/j->o.z;
1002                                        j->part2->p = j->part2->o.transform(j->d) + j->part1->p;
1003                                        j->part2->flags |= LINKFLAG; j->part1->flags |= LINKFLAG; // for delta joint cycle detection
1004                                }
1005                                else
1006                                { // abs.positioning -> calc. delta
1007                                        if (check != EDITING_CHECK)
1008                                        {
1009                                                // calc orient delta
1010                                                //                      Orient tmpo(j->part2->o);
1011                                                //                      tmpo*=j->part1->o;
1012                                                Orient tmpo;
1013                                                j->part1->o.revTransform(tmpo, j->part2->o);
1014                                                tmpo.getAngles(j->rot);
1015                                                j->o = j->rot;
1016                                                // calc position delta
1017                                                Pt3D tmpp(j->part2->p);
1018                                                tmpp -= j->part1->p;
1019                                                j->d = j->part2->o.revTransform(tmpp);
1020                                        }
1021                                }
1022                                if (check != LIVE_CHECK)
1023                                {
1024                                        if (j->shape != Joint::SHAPE_FIXED)
1025                                        {
1026                                                if (j->d() > getMaxJoint().d.x)
1027                                                {
1028                                                        ret = 0;
1029                                                        logPrintf("Model", "internalCheck", LOG_ERROR, "Joint #%d too long (its length %g exceeds allowed %g)%s", i, j->d(), getMaxJoint().d.x, nameForErrors().c_str());
1030                                                }
1031                                        }
1032                                }
1033                        }
1034                        else
1035                        {
1036                                logPrintf("Model", "internalCheck", LOG_ERROR, "Illegal part references in Joint #%d%s", i, nameForErrors().c_str());
1037                                ret = 0;
1038                        }
1039                        if (shape != SHAPE_ILLEGAL)
1040                        {
1041                                if ((j->shape == Joint::SHAPE_BALL_AND_STICK) ^ (shape == SHAPE_BALL_AND_STICK))
1042                                {
1043                                        shape = SHAPE_ILLEGAL;
1044                                        logPrintf("Model", "internalCheck", LOG_WARN, "Inconsistent joint shapes (mixed old and new shapes)%s", nameForErrors().c_str());
1045                                }
1046                        }
1047                }
1048        }
1049#ifdef MODEL_V1_COMPATIBLE
1050        if (!addOldConnectionsInputs())
1051                return 0;
1052#endif
1053
1054        updateNeuroRefno(); // valid refno is important for n-n connections check (later)
1055
1056        for (i = 0; i < neurons.size(); i++)
1057        {
1058                n = (Neuro*)neurons(i);
1059#ifdef MODEL_V1_COMPATIBLE
1060                VALIDMINMAX(n,Neuro,inertia);
1061                VALIDMINMAX(n,Neuro,force);
1062                VALIDMINMAX(n,Neuro,sigmo);
1063                n->conn_refno=-1;
1064                n->weight=1.0;
1065                n->neuro_refno=-1;
1066#endif
1067                n->part_refno = (n->part) ? n->part->refno : -1;
1068                n->joint_refno = (n->joint) ? n->joint->refno : -1;
1069        }
1070
1071        if (parts.size() && (checklevel > 0))
1072        {
1073                for (i = 0; i < parts.size(); i++)
1074                {
1075                        p = (Part*)parts(i);
1076                        if (p->mass <= 0.001)
1077                                p->mass = 1.0;
1078                        p->flags &= ~LINKFLAG;
1079                }
1080                getPart(0)->flags |= LINKFLAG;
1081                int change = 1;
1082                while (change)
1083                {
1084                        change = 0;
1085                        for (i = 0; i < joints.size(); i++)
1086                        {
1087                                j = (Joint*)joints(i);
1088                                if (j->part1->flags&LINKFLAG)
1089                                {
1090                                        if (!(j->part2->flags&LINKFLAG))
1091                                        {
1092                                                change = 1;
1093                                                j->part2->flags |= LINKFLAG;
1094                                        }
1095                                }
1096                                else
1097                                        if (j->part2->flags&LINKFLAG)
1098                                        {
1099                                        if (!(j->part1->flags&LINKFLAG))
1100                                        {
1101                                                change = 1;
1102                                                j->part1->flags |= LINKFLAG;
1103                                        }
1104                                        }
1105                        }
1106                }
1107                for (i = 0; i < parts.size(); i++)
1108                {
1109                        p = (Part*)parts(i);
1110                        if (!(p->flags&LINKFLAG))
1111                        {
1112                                logPrintf("Model", "internalCheck", LOG_ERROR, "Not all parts connected (eg. Part #0 and Part #%d)%s", i, nameForErrors().c_str());
1113                                ret = 0;
1114                                break;
1115                        }
1116                }
1117        }
1118
1119        for (i = 0; i < joints.size(); i++)
1120        {
1121                j = (Joint*)joints(i);
1122                if (j->p1_refno == j->p2_refno)
1123                {
1124                        logPrintf("Model", "internalCheck", LOG_ERROR, "Illegal self connection, Joint #%d%s", i, nameForErrors().c_str());
1125                        ret = 0;
1126                        break;
1127                }
1128                for (k = i + 1; k < joints.size(); k++)
1129                {
1130                        Joint* j2 = (Joint*)joints(k);
1131                        if (((j->p1_refno == j2->p1_refno) && (j->p2_refno == j2->p2_refno))
1132                                || ((j->p1_refno == j2->p2_refno) && (j->p2_refno == j2->p1_refno)))
1133                        {
1134                                logPrintf("Model", "internalCheck", LOG_ERROR, "Illegal duplicate Joint #%d and Joint #%d%s", i, k, nameForErrors().c_str());
1135                                ret = 0;
1136                                break;
1137                        }
1138                }
1139        }
1140        if (shape == SHAPE_ILLEGAL)
1141                ret = 0;
1142        return ret;
1143}
1144
1145/////////////
1146
1147int Model::getErrorPosition(bool includingwarnings)
1148{
1149        return includingwarnings ?
1150                ((f0errorposition >= 0) ? f0errorposition : f0warnposition)
1151                :
1152                f0errorposition;
1153}
1154
1155const Geno& Model::getGeno() const
1156{
1157        return geno;
1158}
1159
1160const Geno Model::getF0Geno()
1161{
1162        if (buildstatus == building)
1163                logPrintf("Model", "getGeno", LOG_WARN, "Model was not completed - missing close()");
1164        if (buildstatus != valid)
1165                return Geno("", '0', "", "invalid");
1166        if (!f0genoknown)
1167        {
1168                if (autobuildmaps)
1169                {
1170                        initF0Map();
1171                        makeGeno(f0geno, f0map);
1172                }
1173                else
1174                {
1175                        delF0Map();
1176                        makeGeno(f0geno);
1177                }
1178                f0genoknown = 1;
1179        }
1180        return f0geno;
1181}
1182
1183int Model::getPartCount() const
1184{
1185        return parts.size();
1186}
1187
1188Part* Model::getPart(int i) const
1189{
1190        return ((Part*)parts(i));
1191}
1192
1193int Model::getJointCount() const
1194{
1195        return joints.size();
1196}
1197
1198Joint* Model::getJoint(int i) const
1199{
1200        return ((Joint*)joints(i));
1201}
1202
1203int Model::findJoints(SList& result, const Part* part)
1204{
1205        Joint *j;
1206        int n0 = result.size();
1207        if (part)
1208                for (int i = 0; j = (Joint*)joints(i); i++)
1209                        if ((j->part1 == part) || (j->part2 == part)) result += (void*)j;
1210        return result.size() - n0;
1211}
1212
1213int Model::findNeuro(Neuro* n)
1214{
1215        return neurons.find(n);
1216}
1217
1218int Model::findPart(Part* p)
1219{
1220        return parts.find(p);
1221}
1222
1223int Model::findJoint(Joint* j)
1224{
1225        return joints.find(j);
1226}
1227
1228int Model::findJoint(Part *p1, Part *p2)
1229{
1230        Joint* j;
1231        for (int i = 0; j = getJoint(i); i++)
1232                if ((j->part1 == p1) && (j->part2 == p2)) return i;
1233        return -1;
1234}
1235
1236
1237////////////////////
1238
1239#ifdef MODEL_V1_COMPATIBLE
1240void Model::calcOldNeuroCount()
1241{
1242        if (oldneurocount>=0) return;
1243        oldneurocount=reorderToOldLayout();
1244}
1245
1246int Model::old_getNeuroCount()
1247{ calcOldNeuroCount();
1248return oldneurocount;}
1249
1250Neuro* Model::old_getNeuro(int i)
1251{calcOldNeuroCount();
1252return (i<oldneurocount)? (Neuro*)getNeuro(i) : (Neuro*)0;
1253}
1254
1255int Model::old_findNeuro(Neuro* n)
1256{calcOldNeuroCount();
1257return findNeuro(n);}
1258
1259Neuro *Model::old_addNewNeuro()
1260{
1261        int count=old_getNeuroCount();
1262        Neuro *nu=addNewNeuro();
1263        nu->setClassName("N");
1264        moveNeuro(nu->refno,oldneurocount);
1265        oldneurocount=count+1;
1266        return (Neuro*)nu;
1267}
1268#endif
1269
1270///////////////////////
1271
1272int Model::getNeuroCount() const
1273{
1274        return neurons.size();
1275}
1276
1277Neuro* Model::getNeuro(int i) const
1278{
1279        return (Neuro*)neurons(i);
1280}
1281
1282int Model::getConnectionCount() const
1283{
1284        int n = 0;
1285        for (int i = 0; i < getNeuroCount(); i++)
1286                n += getNeuro(i)->getInputCount();
1287        return n;
1288}
1289
1290int Model::findNeuros(SList& result,
1291        const char* classname, const Part* part, const Joint* joint)
1292{
1293        Neuro *nu;
1294        SString cn(classname);
1295        int n0 = result.size();
1296        for (int i = 0; nu = (Neuro*)neurons(i); i++)
1297        {
1298                if (part)
1299                        if (nu->part != part) continue;
1300                if (joint)
1301                        if (nu->joint != joint) continue;
1302                if (classname)
1303                        if (nu->getClassName() != cn) continue;
1304                result += (void*)nu;
1305        }
1306        return result.size() - n0;
1307}
1308
1309///////////////////
1310
1311void Model::disturb(double amount)
1312{
1313        int i;
1314        if (amount <= 0) return;
1315        for (i = 0; i < parts.size(); i++)
1316        {
1317                Part *p = getPart(i);
1318                p->p.x += (rnd01 - 0.5)*amount;
1319                p->p.y += (rnd01 - 0.5)*amount;
1320                p->p.z += (rnd01 - 0.5)*amount;
1321        }
1322        for (i = 0; i < joints.size(); i++)
1323        {
1324                Joint *j = getJoint(i);
1325                Pt3D tmpp(j->part2->p);
1326                tmpp -= j->part1->p;
1327                j->d = j->part2->o.revTransform(tmpp);
1328        }
1329}
1330
1331void Model::move(const Pt3D& shift)
1332{
1333        FOREACH(Part*, p, parts)
1334                p->p += shift;
1335}
1336
1337void Model::rotate(const Orient& rotation)
1338{
1339        FOREACH(Part*, p, parts)
1340        {
1341                p->p = rotation.transform(p->p);
1342                p->setOrient(rotation.transform(p->o));
1343        }
1344}
1345
1346void Model::buildUsingSolidShapeTypes(const Model& src_ballandstick_shapes, Part::Shape use_shape, float thickness)
1347{
1348        for (int i = 0; i < src_ballandstick_shapes.getJointCount(); i++)
1349        {
1350                Joint *oj = src_ballandstick_shapes.getJoint(i);
1351                Part *p = addNewPart(use_shape);
1352                p->p = (oj->part1->p + oj->part2->p) / 2;
1353                Orient o;
1354                o.lookAt(oj->part1->p - oj->part2->p);
1355                p->rot = o.getAngles();
1356                p->scale.x = oj->part1->p.distanceTo(oj->part2->p) / 2;
1357                p->scale.y = thickness;
1358                p->scale.z = thickness;
1359        }
1360        for (int i = 0; i < src_ballandstick_shapes.getPartCount(); i++)
1361        {
1362                Part *op = src_ballandstick_shapes.getPart(i);
1363                for (int j = 0; j < src_ballandstick_shapes.getJointCount(); j++)
1364                {
1365                        Joint *oj = src_ballandstick_shapes.getJoint(j);
1366                        if ((oj->part1 == op) || (oj->part2 == op))
1367                        {
1368                                for (int j2 = j + 1; j2 < src_ballandstick_shapes.getJointCount(); j2++)
1369                                {
1370                                        Joint *oj2 = src_ballandstick_shapes.getJoint(j2);
1371                                        if ((oj2->part1 == op) || (oj2->part2 == op))
1372                                        {
1373                                                addNewJoint(getPart(j), getPart(j2), Joint::SHAPE_FIXED);
1374                                        }
1375                                }
1376                                break;
1377                        }
1378                }
1379        }
1380}
1381
1382SolidsShapeTypeModel::SolidsShapeTypeModel(const Model& m, Part::Shape use_shape, float thickness)
1383{
1384        using_model = converted_model = NULL;
1385        if (m.getShapeType() == Model::SHAPE_BALL_AND_STICK)
1386        {
1387                converted_model = new Model;
1388                converted_model->open();
1389                converted_model->buildUsingSolidShapeTypes(m, use_shape, thickness);
1390                converted_model->close();
1391                using_model = converted_model;
1392        }
1393        else
1394        {
1395                converted_model = NULL;
1396                using_model = &m;
1397        }
1398}
1399
1400//////////////////////
1401
1402class MinPart : public Part { public: MinPart() { Param par(f0_part_paramtab, this); par.setMin(); } };
1403class MaxPart : public Part { public: MaxPart() { Param par(f0_part_paramtab, this); par.setMax(); } };
1404class MinJoint : public Joint { public: MinJoint() { Param par(f0_joint_paramtab, this); par.setMin(); } };
1405class MaxJoint : public Joint { public: MaxJoint() { Param par(f0_joint_paramtab, this); par.setMax(); } };
1406class MinNeuro : public Neuro { public: MinNeuro() { Param par(f0_neuro_paramtab, this); par.setMin(); } };
1407class MaxNeuro : public Neuro { public: MaxNeuro() { Param par(f0_neuro_paramtab, this); par.setMax(); } };
1408
1409Part& Model::getMinPart() { static MinPart part; return part; }
1410Part& Model::getMaxPart() { static MaxPart part; return part; }
1411Part& Model::getDefPart() { static Part part; return part; }
1412Joint& Model::getMinJoint() { static MinJoint joint; return joint; }
1413Joint& Model::getMaxJoint() { static MaxJoint joint; return joint; }
1414Joint& Model::getDefJoint() { static Joint joint; return joint; }
1415Neuro& Model::getMinNeuro() { static MinNeuro neuro; return neuro; }
1416Neuro& Model::getMaxNeuro() { static MaxNeuro neuro; return neuro; }
1417Neuro& Model::getDefNeuro() { static Neuro neuro; return neuro; }
Note: See TracBrowser for help on using the repository browser.