dijkstrabox.cc
author Akos Ladanyi <ladanyi@tmit.bme.hu>
Mon, 07 Jul 2008 15:20:43 +0100
changeset 3 2cc5ed6e6255
permissions -rw-r--r--
Use hg changeset hash instead of svn revision.
hegyi@1
     1
/* -*- C++ -*-
hegyi@1
     2
 *
hegyi@1
     3
 * This file is a part of LEMON, a generic C++ optimization library
hegyi@1
     4
 *
hegyi@1
     5
 * Copyright (C) 2003-2006
hegyi@1
     6
 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
hegyi@1
     7
 * (Egervary Research Group on Combinatorial Optimization, EGRES).
hegyi@1
     8
 *
hegyi@1
     9
 * Permission to use, modify and distribute this software is granted
hegyi@1
    10
 * provided that this copyright notice appears in all copies. For
hegyi@1
    11
 * precise terms see the accompanying LICENSE file.
hegyi@1
    12
 *
hegyi@1
    13
 * This software is provided "AS IS" with no warranty of any kind,
hegyi@1
    14
 * express or implied, and with no claim as to its suitability for any
hegyi@1
    15
 * purpose.
hegyi@1
    16
 *
hegyi@1
    17
 */
hegyi@1
    18
hegyi@1
    19
#include <mapstorage.h>
hegyi@1
    20
#include <mapselector.h>
hegyi@1
    21
#include <algobox.h>
hegyi@1
    22
#include <dijkstrabox.h>
hegyi@1
    23
hegyi@1
    24
#include <lemon/dijkstra.h>
hegyi@1
    25
hegyi@1
    26
//No Suurballe in Lemon 1.0
hegyi@1
    27
//#include <lemon/suurballe.h>
hegyi@1
    28
#include <lemon/path.h>
hegyi@1
    29
hegyi@1
    30
enum {INPUT, OUTPUT, MAP_NUM};
hegyi@1
    31
hegyi@1
    32
DijkstraBox::DijkstraBox(std::vector<std::string> t):AlgoBox()
hegyi@1
    33
{
hegyi@1
    34
  init(t);
hegyi@1
    35
}
hegyi@1
    36
    
hegyi@1
    37
SuurballeBox::SuurballeBox(std::vector<std::string> t):DijkstraBox(t)
hegyi@1
    38
{
hegyi@1
    39
  Gtk::Adjustment * adjustment=new Gtk::Adjustment(2, 1, 20, 1, 5);
hegyi@1
    40
  num_set = new Gtk::SpinButton(*adjustment, 5,0);
hegyi@1
    41
hegyi@1
    42
  Gtk::Label * label=new Gtk::Label("No. of paths to find: ");
hegyi@1
    43
//   hbox.pack_start(*label);
hegyi@1
    44
//   hbox.pack_start(*num_set);
hegyi@1
    45
//   hbox.show_all_children();
hegyi@1
    46
hegyi@1
    47
  table.attach(*label, 0,1,2,3);
hegyi@1
    48
  table.attach(*num_set, 1,2,2,3);
hegyi@1
    49
hegyi@1
    50
hegyi@1
    51
//   pack_start(hbox);
hegyi@1
    52
}
hegyi@1
    53
    
hegyi@1
    54
void DijkstraBox::run()
hegyi@1
    55
{
hegyi@1
    56
  if(
hegyi@1
    57
      tabcbt.get_active_text()!="" &&
hegyi@1
    58
      (arcmapcbts[INPUT])->get_active_text()!="" &&
hegyi@1
    59
      (arcmapcbts[OUTPUT])->get_active_text()!="" &&
hegyi@1
    60
      source.get_active_text()!="" &&
hegyi@1
    61
      target.get_active_text()!=""
hegyi@1
    62
    )
hegyi@1
    63
  {
hegyi@1
    64
    const Digraph &g=mapstorage->digraph;
hegyi@1
    65
    Node from, to;
hegyi@1
    66
hegyi@1
    67
    get_from_to(from, to, (Digraph&)g);
hegyi@1
    68
hegyi@1
    69
    std::ostringstream o;
hegyi@1
    70
hegyi@1
    71
    if(!(from==to))
hegyi@1
    72
    {
hegyi@1
    73
      std::string inputmapName = arcmapcbts[INPUT]->get_active_text();
hegyi@1
    74
      std::string outputmapName = arcmapcbts[OUTPUT]->get_active_text();
hegyi@1
    75
hegyi@1
    76
      MapStorage::NumericArcMap& inputmap = mapstorage->getNumericArcMap(inputmapName);
hegyi@1
    77
      MapStorage::NumericArcMap& outputmap = mapstorage->getNumericArcMap(outputmapName);
hegyi@1
    78
hegyi@1
    79
      //zero out output map
hegyi@1
    80
      for (ArcIt i(g); i!=INVALID; ++i)
hegyi@1
    81
      {
hegyi@1
    82
        outputmap[i]=0;
hegyi@1
    83
      }
hegyi@1
    84
hegyi@1
    85
      Dijkstra<Digraph, MapStorage::NumericArcMap > dijkstra(g, inputmap);
hegyi@1
    86
      dijkstra.run(from, to);
hegyi@1
    87
hegyi@1
    88
      if(dijkstra.reached(to))
hegyi@1
    89
      {
hegyi@1
    90
        Node n=to;
hegyi@1
    91
        int length=0;
hegyi@1
    92
        while (n!=INVALID && n!=from)
hegyi@1
    93
        {
hegyi@1
    94
          Arc e=dijkstra.predArc(n);
hegyi@1
    95
          outputmap[e]=1;
hegyi@1
    96
          n=dijkstra.predNode(n);
hegyi@1
    97
          length++;
hegyi@1
    98
        }
hegyi@1
    99
        o << "Result: " << length << " long path, with cost " << dijkstra.dist(to);
hegyi@1
   100
      }
hegyi@1
   101
      else
hegyi@1
   102
      {
hegyi@1
   103
        o << "Result: failed to find shortest path between ";
hegyi@1
   104
        o << source.get_active_text() << " and " << target.get_active_text();
hegyi@1
   105
      }
hegyi@1
   106
      resultlabel.set_text(o.str());
hegyi@1
   107
hegyi@1
   108
      mapstorage->mapChanged(true, (arcmapcbts[OUTPUT])->get_active_text());
hegyi@1
   109
      //   mapstorage->changeActiveMap(true, E_COLOR,
hegyi@1
   110
      // 			      (arcmapcbts[OUTPUT])->get_active_text());
hegyi@1
   111
      //   mapstorage->changeActiveMap(true, E_TEXT,
hegyi@1
   112
      // 			      (arcmapcbts[INPUT])->get_active_text());
hegyi@1
   113
    }
hegyi@1
   114
  }
hegyi@1
   115
}
hegyi@1
   116
hegyi@1
   117
void SuurballeBox::run()
hegyi@1
   118
{
hegyi@1
   119
  if(
hegyi@1
   120
     tabcbt.get_active_text()!="" &&
hegyi@1
   121
     (arcmapcbts[INPUT])->get_active_text()!="" &&
hegyi@1
   122
     (arcmapcbts[OUTPUT])->get_active_text()!="" &&
hegyi@1
   123
     source.get_active_text()!="" &&
hegyi@1
   124
     target.get_active_text()!=""
hegyi@1
   125
     )
hegyi@1
   126
    {
hegyi@1
   127
      const Digraph &g=mapstorage->digraph;
hegyi@1
   128
      Node from, to;
hegyi@1
   129
hegyi@1
   130
      get_from_to(from, to, (Digraph&)g);
hegyi@1
   131
hegyi@1
   132
      std::ostringstream o;
hegyi@1
   133
hegyi@1
   134
      if(!(from==to))
hegyi@1
   135
	{
hegyi@1
   136
          MapStorage::NumericArcMap& inputmap=
hegyi@1
   137
	    mapstorage->getNumericArcMap(arcmapcbts[INPUT]->get_active_text());
hegyi@1
   138
          MapStorage::NumericArcMap& outputmap=
hegyi@1
   139
	    mapstorage->getNumericArcMap(arcmapcbts[OUTPUT]->get_active_text());
hegyi@1
   140
hegyi@1
   141
	  //zero out output map
hegyi@1
   142
	  for (ArcIt i(g); i!=INVALID; ++i)
hegyi@1
   143
	    {
hegyi@1
   144
	      outputmap[i]=0;
hegyi@1
   145
	    }
hegyi@1
   146
	  
hegyi@1
   147
	  //No Suurballe in Lemon 1.0
hegyi@1
   148
	  //Suurballe<Digraph, MapStorage::NumericArcMap > sb((Digraph&)g, inputmap, from, to);
hegyi@1
   149
	  
hegyi@1
   150
	  int found=0;
hegyi@1
   151
	  //No Suurballe in Lemon 1.0
hegyi@1
   152
	  //found=sb.run(num_set->get_value_as_int());
hegyi@1
   153
	  if(found)
hegyi@1
   154
	    {
hegyi@1
   155
	      for(int j=0;j<found;j++)
hegyi@1
   156
		{
hegyi@1
   157
		  Path<Digraph> path;
hegyi@1
   158
		  //No Suurballe in Lemon 1.0
hegyi@1
   159
		  //path=sb.path(j);
hegyi@1
   160
		  for(int k=0;k<path.length();k++)
hegyi@1
   161
		    {
hegyi@1
   162
		      outputmap[path.nth(k)]=j+1;
hegyi@1
   163
		    }
hegyi@1
   164
		}
hegyi@1
   165
	      o << "Result: found " << found << " paths between ";
hegyi@1
   166
	      o << source.get_active_text() << " and " << target.get_active_text();
hegyi@1
   167
	    }
hegyi@1
   168
	  else
hegyi@1
   169
	    {
hegyi@1
   170
	      o << "Result: failed to find shortest path between ";
hegyi@1
   171
	      o << source.get_active_text() << " and " << target.get_active_text();
hegyi@1
   172
	    }
hegyi@1
   173
	  resultlabel.set_text(o.str());
hegyi@1
   174
	  
hegyi@1
   175
	  mapstorage->mapChanged(true, (arcmapcbts[OUTPUT])->get_active_text());
hegyi@1
   176
	  //   mapstorage->changeActiveMap(true, E_COLOR,
hegyi@1
   177
	  // 			      (arcmapcbts[OUTPUT])->get_active_text());
hegyi@1
   178
	  //   mapstorage->changeActiveMap(true, E_TEXT,
hegyi@1
   179
	  // 			      (arcmapcbts[INPUT])->get_active_text());
hegyi@1
   180
	}
hegyi@1
   181
    }
hegyi@1
   182
}
hegyi@1
   183
    
hegyi@1
   184
void DijkstraBox::build_box()
hegyi@1
   185
{
hegyi@1
   186
  //if active tab is changed, labels of digraph nodes had to be loaded into comboboxes
hegyi@1
   187
  //this can be done after the maps are loaded into ComboBoxes
hegyi@1
   188
  signal_upon_maplist_updated().connect(sigc::mem_fun(*this, &DijkstraBox::maplists_updated));
hegyi@1
   189
hegyi@1
   190
  addMapSelector("Cost map: ", true, NUM);
hegyi@1
   191
  addMapSelector("Arcs of path here: ", true, NUM);
hegyi@1
   192
hegyi@1
   193
  Gtk::Label * source_label=new Gtk::Label("Source: ");
hegyi@1
   194
  Gtk::Label * target_label=new Gtk::Label("Target: ");
hegyi@1
   195
hegyi@1
   196
  table.attach(*source_label, 0,1,0,1);
hegyi@1
   197
  table.attach(*target_label, 0,1,1,2);
hegyi@1
   198
  table.attach(source, 1,2,0,1);
hegyi@1
   199
  table.attach(target, 1,2,1,2);
hegyi@1
   200
hegyi@1
   201
hegyi@1
   202
  pack_start(table);
hegyi@1
   203
hegyi@1
   204
  resultlabel.set_text("Result: algorithm is not run yet.");
hegyi@1
   205
  pack_start(resultlabel);
hegyi@1
   206
}
hegyi@1
   207
hegyi@1
   208
void SuurballeBox::build_box()
hegyi@1
   209
{
hegyi@1
   210
}
hegyi@1
   211
hegyi@1
   212
void DijkstraBox::maplists_updated()
hegyi@1
   213
{
hegyi@1
   214
  if(tabcbt.get_active_text()!="")
hegyi@1
   215
    {
hegyi@1
   216
      source.clear();
hegyi@1
   217
      target.clear();
hegyi@1
   218
      const Digraph &g=mapstorage->digraph;
hegyi@1
   219
      for (NodeIt i(g); i!=INVALID; ++i)
hegyi@1
   220
	{
hegyi@1
   221
	  std::ostringstream text;
hegyi@1
   222
	  text << mapstorage->getLabel(i);
hegyi@1
   223
	  source.prepend_text(text.str());
hegyi@1
   224
	  target.prepend_text(text.str());
hegyi@1
   225
	}
hegyi@1
   226
    }
hegyi@1
   227
}
hegyi@1
   228
hegyi@1
   229
void DijkstraBox::get_from_to(Node & from, Node & to, Digraph & g)
hegyi@1
   230
{
hegyi@1
   231
  int assigned=0;
hegyi@1
   232
  for (NodeIt i(g); (i!=INVALID) && (assigned<2); ++i)
hegyi@1
   233
    {
hegyi@1
   234
      std::ostringstream text;
hegyi@1
   235
      text << mapstorage->getLabel(i);
hegyi@1
   236
      if(!(text.str().compare(source.get_active_text())))
hegyi@1
   237
	{
hegyi@1
   238
	  from=i;
hegyi@1
   239
	  assigned++;
hegyi@1
   240
	}
hegyi@1
   241
      if(!(text.str().compare(target.get_active_text())))
hegyi@1
   242
	{
hegyi@1
   243
	  to=i;
hegyi@1
   244
	  assigned++;
hegyi@1
   245
	}
hegyi@1
   246
    }
hegyi@1
   247
}