설명

sqrec(x,S) = 1/x if x >= S
             1/S - r/S^2 + r^2/s^3 그렇지 않으면 r = x - S
S >= 1e-10의 경우 기본값 S = 1e-10입니다.

기고자: Steve, 2017년 3월

소형 모델 유형 :슬롯 나라

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

메인 파일 : fnsqrec.gms 포함: fnset_xy.inc[html] fntest_xy.inc[html]

$title 'sqrec 내장의 정확성 테스트' (FNSQREC,SEQ=727)

$onText
sqrec(x,S) = 1/x x >= S인 경우
             1/S - r/S^2 + r^2/s^3 그렇지 않으면 r = x - S
S >= 1e-10의 경우 기본값 S = 1e-10입니다.

기고자: Steve, 2017년 3월
$offText

$include fnset_xy.inc

$gdxIn fnsqrec
$load aeps 담당자 relToInput T 데이터
$gdxIn

루프 T,
  data(T, 'f') = sqrec.value( data(T,'x'),data(T,'y'));
  data(T, 'fx') = sqrec.grad(1: data(T,'x'),data(T,'y'));
  data(T, 'fy') = sqrec.grad(2: data(T,'x'),data(T,'y'));
  data(T,'fxx') = sqrec.hess(1:1:data(T,'x'),data(T,'y'));
  data(T,'fxy') = sqrec.hess(1:2:data(T,'x'),data(T,'y'));
  data(T,'fyx') = sqrec.hess(2:1:data(T,'x'),data(T,'y'));
  data(T,'fyy') = sqrec.hess(2:2:data(T,'x'),data(T,'y'));
  data(T, 'rc') = mathlastrc;
  data(T, 'ec') = mathlastec;
;

$include fntest_xy.inc