설명

이 모델은 ex1.gms의 모델과 유사합니다.
C 및 Fortran 예제에 대한 코드만 포함됩니다.
외부 방정식에 구현 오류가 있습니다.
소스 코드에 **ERROR로 표시되어 있습니다. 이러한 오류를 다음과 비교하십시오.
슬롯 사이트 목록 파일의 메시지입니다.

소형 모델 유형 :슬롯 사이트

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

메인 파일 : er2.gms

$title 외부 방정식 - 오류 예 2 (ER2,SEQ=571)

$onText
  이 모델은 ex1.gms의 모델과 유사하지만
  C 및 Fortran 예제에 대한 코드만 포함됩니다.
  외부 방정식에 구현 오류가 있습니다.
  소스 코드에 **ERROR로 표시되어 있습니다. 이러한 오류를 다음과 비교하십시오.
  슬롯 사이트 목록 파일의 메시지
$offText

i / i1*i4 / 설정
별칭(i,j);

모수 Q(i,j) 공분산 행렬
          X0(i) 목표;
Q(i,j) = power(0.5, abs(ord(i)-ord(j)) );
X0(i) = ord(i);
Q, X0을 표시합니다.

변수 x(i), z;
방정식 zdef, zdefX;

* 슬롯 사이트에서 구현된 원하는 방정식은 다음과 같습니다.

  zdef .. sum( (i,j), (x(i)-x0(i)) * Q(i,j) * (x(j)-x0(j) ) ) =e= z;

* 다음과 같이 외부 방정식으로 구현됩니다.

  zdefX .. sum(i, ord(i)*x(i) ) + (카드(i)+1)* z =X= 1;

$onText
  방정식의 계수는 X-변수가 다음과 같다는 것을 보여줍니다.
  1부터 카드(i)까지 번호가 매겨져 있고 Z는 카드(i)+1의 숫자를 가지고 있습니다.
  외부 방정식 코드.
  외부 방정식은 하나만 있으며 숫자 1, 즉 값이 있습니다.
  오른쪽의.
  방정식의 모든 변수에는 변수가 할당되어야 합니다.
  숫자이며 모두 외부 방정식에 나타나야 합니다. 당신
  아직 일부 항이 선형이라는 것을 솔버에게 알릴 수 없습니다.
  -- 모든 항은 솔버의 관점에서 볼 때 비선형입니다.
$offText

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

$set N er2
$set c_cbN %pre%%N%c_cb%suf%
$set f_cbN %pre%%N%f_cb%suf%

모델 %N% '슬롯 사이트 구현' / zdef /;
model %c_cbN% '콜백이 포함된 C의 외부 방정식' / zdefX /;
model %f_cbN% '콜백이 포함된 F77의 외부 방정식' / zdefX /;

옵션 limcol = 0, iterlim=1000;

* 목표와 비교하여 솔루션을 확인하십시오.

매개변수 보고서(*,*,*) 솔루션 요약;
스칼라 totdist /0/;

$onEchoV > runme.gms
z.l = 0;
z.m = 0;
x.l(i) = 0;
x.m(i) = 0;
zdefX.l = 0;
zdefX.m = 0;
z를 최소화하는 nlp를 사용하여 %1을 해결합니다.
$if %1==er2 abort$(%1.solvestat <> 1) '%1 모델 실행에 문제가 있습니다';
$if not %1==er2 abort$(%1.solvestat = 1) '%1 모델 실행에 문제가 있습니다';
실행 오류 = 0;
Report('Solve','Stat', '%1') = %1.solvestat;
Report('모델','Stat', '%1') = %1.modelstat;
Report(i,'대상', '%1') = x0(i);
보고서(i,'값', '%1') = x.l(i);
보고서(i,'거리','%1') = abs(x.l(i) - x0(i));
totdist = totdist + sum(i,abs(x.l(i) - x0(i)));
$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=er2c_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 er2f_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%

보고서 표시;

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