설명

DLL 테스트 문제 번호 MCP3 -- 조밀한 비선형 시스템

소형 모델 유형 :슬롯 머신

카테고리 : 슬롯 머신 테스트 라이브러리

메인 파일 : exmcp3.gms

$title 외부 방정식 - MCP 3 예(EXMCP3,SEQ=575)

$onText
  DLL 테스트 문제 번호 MCP3 - 조밀한 비선형 시스템
$offText

내가 /i1*i30/으로 설정;
별칭(i,j);

긍정적인
변수 z(j) 상보성 문제를 알 수 없음;

방정식 equ(j) 부등식 시스템,
                equx(j) 외부 불평등 시스템;

매개변수 m(i,j) 상수 계수 행렬,
                q(j) 상수 RHS 벡터,
                z0(j) 공통 시작점;

m(i,j)$(균일(0,1) gt 0.5) = 균일(0,1);
q(j) = 균일(0,1);
z0(j) = 균일(0,1);

equ(i)..sqr(z(i)) + sum(j, m(i,j) * z(j)) =g= q(i);

equx(i).. sum(j, (ord(j)*z(j))$(m(i,j) 또는 sameas(i,j))) =x= ord(i);

$ 미리 설정
$ifI %system.filesys%==unix $set pre 'lib'
$ 세트 '64'

$세트 N exmcp3
$set c_cbN %pre%%N%c%suf%
$set f_cbN %pre%%N%f%suf%

모델 %N% '슬롯 머신 구현' / equ.z /;
model %c_cbN% 'C의 외부 방정식' / equx.z /;
model %f_cbN% 'F77의 외부 방정식' / equx.z /;

스칼라 totdist /0/;
매개변수 해(J,*);

abort$(card(i) gt 50) " 차원이 너무 많습니다. 한계는 50개입니다!";

* 슬롯 머신 스크래치 디렉토리에 데이터 파일 쓰기,
* 확장자가 .dat이므로 자동으로 삭제됩니다.
* 슬롯 머신 실행이 완료되면
파일 kdat / '%슬롯 머신scrdir%exmcp3.%슬롯 머신scrext%' /;
kdat를 넣어;
kdat를 넣고, 카드(i):0:0/;
파일.nr=2; file.nw=22; 파일.nd=13;
loop(i, q(i)/;)를 넣습니다.
loop(j, loop(i, put m(i,j)/;); );
닫다;

$onEchoV > runme.gms
z.l(j) = z0(j);
mcp를 사용하여 %1을(를) 해결하세요.
솔루션(j,"%1") = z.l(j);
totdist = totdist + sum J, abs(z.l(J)-solution(J,'exmcp3'));
$offEcho

$ 확장 '.dll' 설정
$ifI %system.filesys%==unix $set ext '.so'
$ifI %system.platform%==dex $set ext '.dylib'
$ifI %system.platform%==dax $set ext '.dylib'

$ 세트 eq
$ifI %system.filesys%==unix $set eq "'"

$if set runall $set runC_cb '1' set runF_cb '1'

$ifThen nocomp를 설정하지 않음
$ ifI set runC_cb $call 슬롯 머신 complink lo=%슬롯 머신lo% --lang=c --files=exmcp3c_cb.c --libname=%c_cbN%%ext%
$ 오류 수준 1인 경우 $abort C 라이브러리 컴파일 오류
$ ifI set runF_cb $call 슬롯 머신 complink lo=%슬롯 머신lo% --lang=fortran90 --files=%eq%"gehelper.f90 msg2_f.f90 exmcp3f_cb.f90"%eq% --libname=%f_cbN%%ext%
$ 오류 수준 1인 경우 $abort Fortran90 라이브러리 컴파일 오류
$endIf

$ batInclude runme %N%
$if set runC_cb $batInclude runme %c_cbN%
$if set runF_cb $batInclude runme %f_cbN%

디스플레이 솔루션;

(totdist < 1.0E-6)인 경우,
  "@@@@ #테스트에 통과했습니다."를 표시합니다.
그렇지 않으면
  abort totdist, "@@@@ #테스트가 통과되지 않았습니다. 자세한 내용은 exmcp3.lst를 검사하세요.";
;