source: branches/tree/PROBE_SET/ps_my2asciipaths.cxx

Last change on this file was 15738, checked in by westram, 8 years ago
  • rename only
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.5 KB
Line 
1// =============================================================== //
2//                                                                 //
3//   File      : ps_my2asciipaths.cxx                              //
4//   Purpose   :                                                   //
5//                                                                 //
6//   Coded by Wolfram Foerster in October 2002                     //
7//   Institute of Microbiology (Technical University Munich)       //
8//   http://www.arb-home.de/                                       //
9//                                                                 //
10// =============================================================== //
11
12#include "ps_node.hxx"
13
14//  GLOBALS
15
16typedef pair<bool, SpeciesID> Step;
17static vector<Step> *__PATH = new vector<Step>;
18
19static void PS_print_paths(const PS_NodePtr& _ps_node) {
20    //  recursively print the paths to the leafs
21
22    // path
23    __PATH->push_back(Step(_ps_node->hasInverseProbes(), _ps_node->hasProbes() ? _ps_node->getNum() : -(_ps_node->getNum())));
24
25    // children
26    if (_ps_node->hasChildren()) {
27        for (PS_NodeMapConstIterator i = _ps_node->getChildrenBegin(); i != _ps_node->getChildrenEnd(); ++i) {
28            PS_print_paths(i->second);
29        }
30    }
31    else {
32        // print path in leaf nodes
33        printf("[%4zu] ", __PATH->size());
34        for (vector<Step>::const_iterator i=__PATH->begin(); i != __PATH->end(); ++i) {
35            printf("%4i%c ", i->second, i->first ? '+' : ' ');
36        }
37        printf("\n");
38    }
39
40    // path
41    __PATH->pop_back();
42}
43
44
45//  ====================================================
46//  ====================================================
47
48int main(int argc,   char *argv[]) {
49
50    // open probe-set-database
51    if (argc < 2) {
52        printf("Missing arguments\n Usage %s <database name>\n", argv[0]);
53        exit(1);
54    }
55
56    const char *input_DB_name = argv[1];
57
58    printf("Opening probe-set-database '%s'..\n", input_DB_name);
59    PS_FileBuffer *fb = new PS_FileBuffer(input_DB_name, PS_FileBuffer::READONLY);
60    PS_NodePtr     root(new PS_Node(-1));
61    root->load(fb);
62    printf("loaded database (enter to continue)\n");
63
64
65    for (PS_NodeMapConstIterator i = root->getChildrenBegin(); i != root->getChildrenEnd(); ++i) {
66        PS_print_paths(i->second);
67    }
68
69    printf("(enter to continue)\n");
70
71    delete fb;
72    root.setNull();
73    printf("root should be destroyed now\n");
74    printf("(enter to continue)\n");
75
76    return 0;
77}
Note: See TracBrowser for help on using the repository browser.