|
 |
 |
 |
Q & A ...
게시판 (Optimization)
|
|
|
※ 안녕하십니까..?
저희 KSTEC의 제품이나 서비스에 대해 궁금하신 점이나 기술지원을 원하시는 고객님은 회사명,
부서명, 성명, 사용제품명, 제품 VERSION을 명시하여 주시기 바랍니다.
※ 게시하신 내용에 대해서는 신속하게 답변해 드리겠습니다.
※ 해당되는 제품군을 선택하신 후 하고 싶으신 말씀을 적어 주십시오.
|
|
|
|
제 목 |
[RE]제약조건의 크기와 관련된 질문입니다. |
|
작성자 |
장용성 |
작성일 |
2003-09-30 |
조회수 |
990 회 |
|
| 첨부파일 |
첨부된 파일없음.
|
|
우선 callable library를 사용하기가 불편하기 때문에 가장 좋은 방법은 concert를 사용해서 문제를 표현하는 게 상당히 쉬울 거라고 생각하구요, 그럴 사정이 안 되시면 concert처럼 데이터를 처리할 수 있게 lp를 표현하기 위한 클래스가 필요합니다. 예를들어 MyCplex라는 클래스를 만들어서 그 안에서 addRow(), startColumn(), addNzToCurrentColumn(), endColumn()등의 멤버를 통해서 접근하시면 쉬울 수 있을 거 같구요.
각 변수를 만들때는 필요없는 변수를 제거하기 위해서 컴파일러가 지원을 해줄 수 있다면 map을 통해서 인덱스로 관리를 해주시는 게 좋을 것 같구요. 변수나 제약식은 종류별로 배열을 만드시는 게 좋을 것 같습니다.. 종류별로 클래스를 만드시는 것도 유용할 겁니다.
예를 들어 제약식의 종류별로 MySetRows클래스를 만들거나, 변수의 종류별로 MySetColumns클래스를 만들어서 그 안에 변수나 제약식에 해당되는 자료를 집어넣는 겁니다.
아마 제가 위에서 언급했던 MyCplex 클래스 정의만 해도 많이 수월해 질 거라고 봅니다. 수고하세요..
보안상 모든 부분은 불가능하구요, 일부 함수를 보여드리겠습니다..
IlCgCplex::IlCgCplex(int c,int r,int z)
:numRows_(0), numCols_(0),
numNz_(0), numColNz_(0), clp_(0), maxMipMins_(-1)
{
obj_.reserve(c);
rhs_.reserve(r);
sense_.reserve(r);
matbeg_.reserve(c);
matcnt_.reserve(r);
matind_.reserve(z);
matval_.reserve(z);
lb_.reserve(c);
ub_.reserve(c);
rngval_.reserve(c);
env_ = CPXopenCPLEX(&status_);
if (env_ == 0)
{
char errorMsg[1024];
CPXgeterrorstring(env_, status_, errorMsg);
cout << "Could not open CPLEX environment" << endl;
cout << errorMsg << endl;
exit(1);
}
objSense_ = MIN;
}
void
IlCgCplex::addRow(double rhs, double rngval, char sense)
{
rhs_.push_back(rhs);
rngval_.push_back(rngval);
sense_.push_back(sense);
if (sense != ''L'' && sense != ''G'' && sense != ''E'')
terminate("Sense of Row wrong");
numRows_++;
}
void
IlCgCplex::startColumn(double obj, double lb, double ub)
{
assert(numColNz_ == 0);
matbeg_.push_back(numNz_);
obj_.push_back(obj);
lb_.push_back(lb);
ub_.push_back(ub);
}
void
IlCgCplex::endColumn()
{
matcnt_.push_back(numColNz_);
numColNz_ = 0;
numCols_++;
}
void
IlCgCplex::addNzToCurrentColumn(int row, double val)
{
if (clp_ != 0)
assert(row < getNumRows() + numRows_);
else
assert(row < numRows_);
matind_.push_back(row);
matval_.push_back(val);
numNz_++;
numColNz_++;
}
void
IlCgCplex::registerLp()
{
clp_ = CPXcreateprob(env_, &status_, "DynColGen");
if(clp_ == 0)
{
terminate("Failed to Create LP ");
}
status_ = CPXcopylp(env_, clp_, numCols_, numRows_,objSense_,
&obj_[0], &rhs_[0], &sense_[0], &matbeg_[0],&matcnt_[0],&matind_[0],&matval_[0],
&lb_[0], &ub_[0],&rngval_[0]);
if(status_)
{
terminate("Failed to Copy LP");
}
resetColumns();
resetRows();
#if 0
CPXcheckcopylp(env_, clp_, numCols_, numRows_, objSense_,
obj_.getElements(), rhs_.getElements(),
sense_.getElements(), matbeg_.getElements(),
matcnt_.getElements(), matind_.getElements(),
matval_.getElements(), lb_.getElements(), ub_.getElements(),
rngval_.getElements());
#endif
} |
|
|
|
|
|
|
관련글 보기
|
|
|
"[RE]제약조건의 크기와 관련된 질문입니다."와(과) 관련된 글이
3건 있습니다. |
 |
제약조건의 크기와 관련된 질문입니다.
|
성영태 |
2003-09-30 |
|
[RE]제약조건의 크기와 관련된 질문입니다.
|
장용성 |
2003-09-30 |
|
[RE][RE]제약조건의 크기와 관련된 질문입니다.
|
성영태 |
2003-09-30 |
|
[RE][RE][RE]제약조건의 크기와 관련된 질문입니다.
|
장용성 |
2003-09-30 |
|
|
|
|
|
|
