설명
이것은 Fortran이 포함된 예제 1의 수정된 버전입니다.
루틴은 정의의 고급 기능을 사용하는 방법을 보여줍니다.
상수 도함수 및 헤세 시간 벡터 값.
소형 모델 유형 :슬롯 머신
메인 파일 : ex1x.gms
$title 외부 방정식 - 예 1x (EX1X,SEQ=564)

$onText
  이것은 Fortran이 포함된 예제 1의 수정된 버전입니다.
  루틴은 정의의 고급 기능을 사용하는 방법을 보여줍니다.
  상수 도함수 및 헤세 시간 벡터 값.
$offText

i / i1*i20 / 설정
별칭(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
  zdefX에서는 일부 계수가 상수라고 정의할 수 없습니다.
  이는 외부 모듈의 새로운 기능을 사용하여 수행됩니다.
$offText

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

$세트 N ex1x
$set c_cbN %pre%%N%c_cb%suf%
$set fN %pre%%N%f%suf%

모델 %N% '슬롯 머신 구현' / zdef /;
model %c_cbN% 'C에서 추가 기능을 갖춘 외부 방정식' / zdefX /;
model %fN% 'F77에 추가 기능이 있는 외부 방정식' / zdefX /;

옵션 limcol = 0, sysout = on;

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

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

$onEchoV > runme.gms
z.l = 0;
z.m = 0;
x.l(i) = 0;
x.m(i) = 0;
zdef.l = 0;
zdef.m = 0;
z를 최소화하는 nlp를 사용하여 %1을 해결합니다.
abort$(%1.solvestat <> 1) '%1 모델 실행에 문제가 있습니다';
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 '1'

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

$ batInclude runme %N%
$if set runC_cb $batInclude runme %c_cbN%
$if set runF $batInclude runme %fN%

보고서 표시;

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