fnrpow2.gms : rpower 내장의 정확성 테스트

설명

rpower(x,y) = x**y
이러한 포인트에는 더 어려운 포인트가 포함되어 있습니다.
예외가 발생합니다.

소형 모델 유형 :메가 슬롯


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


메인 파일 : fnrpow2.gms 포함: fnset_xy.inc[htmlfntest_xy.inc[html]

$title 'rpower 내장의 정확성 테스트' (FNRPOW2,SEQ=195)

* rpower(x,y) = x**y
* 이 포인트에는 더 어려운 포인트가 포함되어 있습니다.
* 예외가 발생합니다.
$include fnset_xy.inc

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

담당자 = 5e-13;
* GDX 버그 해킹
relToInput = max(0,relToInput);

data(T,'fxy_') = data(T,'fxy_');

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

중단$(실행 오류 lt 10)
  'x**y의 등급/헤스, x < 0)을 가져오는 동안 실행 오류가 발생해야 합니다.';
실행 오류 = 0;

relToInput을 표시;
$include fntest_xy.inc