[gecode-users] Linear function problem
manel askri
askri.manel1 at gmail.com
Tue Aug 27 23:11:55 CEST 2013
hi everyone !
I have problem with Linear function when i compile , the compiler tells me
that :
exeption Int::linear :size of arguments mismatch
i think that the problem is with the linear function : my code is this :
#include <gecode/driver.hh>
#include <gecode/int.hh>
#include <gecode/search.hh>
#include <gecode/minimodel.hh>
using namespace Gecode;
class MMKP : public Script{
protected:
static const int n=4 ;//number of testers
static const int m=5;// number of nodes
IntArgs R;
IntArgs C;
IntArgs B;
IntArgs Dr;
IntArgs Dc;
IntArgs Db;
IntArgs g;
IntVarArray a;
IntVarArray k;
IntVar l;
public :
enum {
find_Optimal_Solution, find_realizable_solution
};
MMKP (const Options& opt) : a(*this,n*m, 0,1){
R=IntArgs(4,10,15,20,40);
C=IntArgs(4,20,25,30,40);
B=IntArgs(4,20,30,35,40);
Dr=IntArgs(4,5,10,15,25);
Dc=IntArgs(4,15,20,35,37);
Db=IntArgs(4,10,15,20,35);
k=IntVarArray(*this,n*m ,0,100000);
Matrix <IntVarArray> results(k, n,m);
//creation variables
//IntVarArray a(*this,n*m,0,1); // Array of n*m boolean variables
Matrix <IntVarArray> X (a,n,m);// Matrix "view" of the array a
// objectiv variable
IntVar gain (*this, 1,1000000);
//creation of constraints
// ... over rows
for ( int j=0; j<n;j++)
{
linear(*this , X.row(j),IRT_EQ,1);
}
//... over columns
// first, get the columns, we will use an intermidiare matrix XDual
IntVarArray b(*this, m*n,0,1);
Matrix <IntVarArray> XDual (b, m, n);
for (int i=0; i<m;i++)
{
for ( int j =0; j<n ; j++)
{
XDual(i,j)=X(j,i);
}
}
for (int j = 0; j < m; j++) {
linear(*this, Dr,XDual.row(j),IRT_NQ, R[j]);
}
for (int j = 0; j < m; j++) {
linear (*this, Dc, XDual.row(j), IRT_NQ,C[j]);
}
for (int j = 0; j < m; j++) {
linear (*this, Db, XDual.row(j), IRT_NQ,B[j]);
}
switch (opt.model()) {
case find_Optimal_Solution:
g=IntArgs(4,20,30,40,50);
//Objective function
for (int i = 0; i < n; i++)
{
linear(*this, g,X.row(i), IRT_EQ, gain);
}
for ( int i=0; i<n;i++){
for ( int j =0; j<m;j++)
{
results(i,j)=X(i,j);
}
}
break;
case find_realizable_solution:
for ( int i=0; i<n;i++){
for ( int j =0; j<m;j++)
{
results(i,j)=X(i,j);
}
}
break;
// post branching
branch(*this, a, INT_VAR_SIZE_MAX(), INT_VAL_MAX());
}
}
// search support
MMKP(bool share, MMKP& s) : Script(share, s){
a.update(*this, share, s.a);
}
virtual Space* copy(bool share) {
return new MMKP(share,*this);
}
// print solution
void print(std::ostream& os) const {
for(int i = 0; i < n; i++) {
for(int j = 0; j < n; j++)
os << std::setw(4) << a[i * n + j];
os << std::endl;
}
}
};
// main function
int main(int argc, char* argv[]) {
Options opt("MMKP");
opt.model(MMKP::find_Optimal_Solution);
opt.model(MMKP::find_realizable_solution);
opt.parse(argc,argv);
Script::run<MMKP,DFS,Options>(opt);
return 0;
}
so thank you very much for helping me
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.gecode.org/pipermail/users/attachments/20130827/188c5dfa/attachment.html>
More information about the users
mailing list