// ---------------------------------------------------------------------- 2opt
      // Precondizione: sol ha n+1 elementi, l'ultimo e' uguale al primo
      public int Loc2optASIM(int[] sol, int[,] d)
      {
         int znew=-1,zold;
         int[] solnew = new int[sol.Length];
         int[] zorario = new int[sol.Length];   // costi incrementali tour in senso orario
         int[] zantior = new int[sol.Length];   // costi incrementali tour in senso antiorario
         int end1a,end2a,end1b,end2b,cost1,cost2;
         int nk;
         int i;

         int nsol = sol.Length;  // con l'ultimo uguale al primo

         if(nsol <= 3) 
            goto end;

lab100:  zold = 0;  // qui si inizia una uova ricerca
         zorario[0] = 0;
         zantior[nsol-1] = 0;
         for(int js=0;js<nsol-1;js++)
         {  zold += d[sol[js],sol[js+1]];
            zorario[js+1] = zold;
            zantior[nsol-js-2] = zantior[nsol-js-1] + d[sol[nsol-js-1],sol[nsol-js-2]];
         }
         znew = zold;
         
         for(int ind1=1;ind1<nsol-1;ind1++)
         {  end1a=sol[ind1-1];   // taglio in posizione ind1-1
            end1b=sol[ind1];

            for(int ind2=ind1+2;ind2<nsol;ind2++)
            {  end2a=sol[ind2-1];      // taglio in posizione ind2-1
               end2b=sol[ind2];

               // provo a riattaccare gli end1 fra loro e gli end2 fra loro, senso orario
               cost1 =  zorario[nsol-1];
               cost2 =  zorario[ind1-1] + d[end1a,end2a] + (zantior[ind1] - zantior[ind2-1])+
                        d[end1b,end2b] + (zorario[nsol - 1] - zorario[ind2]);
               if( cost2 < cost1 && zorario[nsol-1] <= zantior[0]) 
               {
                  for(int lk=0;lk<ind1;lk++) // 0 ... ind1-1
                     solnew[lk]=sol[lk];

                  nk=ind1-1;
                  for(int lk=ind2-1;lk>=ind1;lk--) // ind1 ... ind2-1 (reversed)
                  {  nk++;
                     solnew[nk]=sol[lk];
                  }

                  for(nk=ind2;nk<nsol;nk++)  // ind2 ... n
                     solnew[nk]=sol[nk];

                  znew = 0;
                  for(i=0;i<nsol-1;i++)
                     znew += d[solnew[i],solnew[i+1]];

                  if(znew != cost2)
                     MessageBox.Show("2opt (1): qualcosa tragia");
                     
                  zold=znew;
                  for(int k=0;k<nsol;k++)
                     sol[k]=solnew[k];

                  goto lab100;   // ricomincia dalla nuova soluzione
               }