Dibujar una superficie con R:

Este el código que hace esto:

require(grDevices) # for trans3d
## More examples in  demo(persp) !!
##                   -----------

# (1) The Obligatory Mathematical surface.
#     Rotated sinc function.

x <- 5="" length="30)<br" seq="">y <- br="" x="">f <- br="" function="" r="" x="" y="">z <- br="" f="" outer="" x="" y="">z[is.na(z)] <- 1="" br="">op <- bg="white" br="" par="">persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue");

Hacer un pool de solucion en CPLEX

 El el fichero cplex.opt tenemos que introducir la instrucción solnpool  que nos crea el pool de soluciones que está indexado en el gdx solnpool.gdx. Con solnpoolpop 1 las soluciones se generan mediante el procedimiento normal de optimización de Cplex, pero con la instrucción solnpoolpop 2, se generán un número PopulateLim de soluciones mediante un procedimiento especial de Cplex para generar varias soluciones. Las soluciones las filtramos mediante solnpoolintesity 4 y solnpoolAgap 0.0. Para mas info mirar este ejemplo http://www.gams.com/modlib/libhtml/solnpool.htm

  File fcpx / cplex.opt /;
  option limrow=0,limcol=0, optcr=0, mip=cplex;
  LNQP.optfile=1;  LNQP.solprint=%solprint.Quiet%;  LNQP.savepoint = 1;
  
putclose fcpx 'solnpool solnpool.gdx' / 'solnpoolintensity 4' /
  'solnpoolpop 2' / 'solnpoolAgap = 0.0';

web sobre diseño de algoritmos y analisis

http://openclassroom.stanford.edu/MainFolder/CoursePage.php?course=IntroToAlgorithms

Optimization wiki

here we have two optimization web page very interesting Neos_Wiki and optimization tree.

Trazar CPLEX

http://yetanothermathprogrammingconsultant.blogspot.com/2008/05/undocumented-gamscplex-options.html

Aprovechar multi-core CPLEX

Para que cplex se ejecute en paralelo, añadir la siguiente instrucción en el código de gams:

$onecho > cplex.opt

lpmethod 6
parallelmode 0
threads 0
predual -1

$offecho


model.optfile = 1;

Generar combinaciones(n,k)

int generarComb(int i,int l,int k,int n){

if(i==k){
if(l==0){
for(int r=0;r < k;r++){
cout < < mark[r];
}
cout < < endl;
}
}
else{
if(l > 0){
mark[i]=true;
generarComb(i+1,l-1,k,n);
mark[i]=false;
}
generarComb(i+1,l,k,n);

}

return 0;
}