partial_mapped_xover.cpp

00001 /* 
00002 * <partial_mapped_xover.cpp>
00003 * Copyright (C) DOLPHIN Project-Team, INRIA Futurs, 2006-2007
00004 * (C) OPAC Team, LIFL, 2002-2007
00005 *
00006 * Sebastien Cahon, Alexandru-Adrian Tantar
00007 *
00008 * This software is governed by the CeCILL license under French law and
00009 * abiding by the rules of distribution of free software.  You can  use,
00010 * modify and/ or redistribute the software under the terms of the CeCILL
00011 * license as circulated by CEA, CNRS and INRIA at the following URL
00012 * "http://www.cecill.info".
00013 *
00014 * As a counterpart to the access to the source code and  rights to copy,
00015 * modify and redistribute granted by the license, users are provided only
00016 * with a limited warranty  and the software's author,  the holder of the
00017 * economic rights,  and the successive licensors  have only  limited liability.
00018 *
00019 * In this respect, the user's attention is drawn to the risks associated
00020 * with loading,  using,  modifying and/or developing or reproducing the
00021 * software by the user in light of its specific status of free software,
00022 * that may mean  that it is complicated to manipulate,  and  that  also
00023 * therefore means  that it is reserved for developers  and  experienced
00024 * professionals having in-depth computer knowledge. Users are therefore
00025 * encouraged to load and test the software's suitability as regards their
00026 * requirements in conditions enabling the security of their systems and/or
00027 * data to be ensured and,  more generally, to use and operate it in the
00028 * same conditions as regards security.
00029 * The fact that you are presently reading this means that you have had
00030 * knowledge of the CeCILL license and that you accept its terms.
00031 *
00032 * ParadisEO WebSite : http://paradiseo.gforge.inria.fr
00033 * Contact: paradiseo-help@lists.gforge.inria.fr
00034 *
00035 */
00036 
00037 #include <assert.h>
00038 
00039 #include <utils/eoRNG.h>
00040 
00041 #include "partial_mapped_xover.h"
00042 #include "mix.h"
00043 
00044 void PartialMappedXover :: repair (Route & __route, unsigned __cut1, unsigned __cut2) {
00045   
00046   unsigned v [__route.size ()] ; // Number of times a cities are visited ...
00047   
00048   for (unsigned i = 0 ; i < __route.size () ; i ++)
00049     v [i] = 0 ;
00050   
00051   for (unsigned i = 0 ; i < __route.size () ; i ++)
00052     v [__route [i]] ++ ;
00053   
00054   std :: vector <unsigned> vert ;
00055 
00056   for (unsigned i = 0 ; i < __route.size () ; i ++)
00057     if (! v [i])
00058       vert.push_back (i) ;
00059   
00060   mix (vert) ;
00061 
00062   for (unsigned i = 0 ; i < __route.size () ; i ++)
00063     if (i < __cut1 || i >= __cut2)
00064       if (v [__route [i]] > 1) {
00065         __route [i] = vert.back () ;
00066         vert.pop_back () ;
00067       }
00068 }
00069 
00070 bool PartialMappedXover :: operator () (Route & __route1, Route & __route2) {
00071     
00072   unsigned cut1 = rng.random (__route1.size ()), cut2 = rng.random (__route2.size ()) ;
00073   
00074   if (cut2 < cut1)
00075     std :: swap (cut1, cut2) ;
00076   
00077   // Between the cuts
00078   for (unsigned i = cut1 ; i < cut2 ; i ++)
00079     std :: swap (__route1 [i], __route2 [i]) ;
00080   
00081   // Outside the cuts
00082   repair (__route1, cut1, cut2) ;
00083   repair (__route2, cut1, cut2) ;
00084   
00085   __route1.invalidate () ;
00086   __route2.invalidate () ;
00087 
00088   return true ;
00089 }

Generated on Mon Oct 8 11:16:45 2007 for ParadisEO-PEOMovingObjects by  doxygen 1.4.7