설명

부호화된 거듭제곱 함수를 테스트하여 실제 도함수를 비교합니다.
수치적으로 계산된 것들로요.  우리는 범위를 너무 많이 밀지 않습니다
이 테스트에서는 모든 신호가 올바른지 확인하고 싶습니다.
등. 중요한 것은 음수 x로 테스트하여 우리가 알 수 있도록 하는 것입니다.
우리의 구현에는 심각한 오류가 발생하지 않습니다.

기여자: Steve Dirkse

소형 모델 유형 :무료 슬롯

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

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

$title '부호력 내장의 정확성 테스트' (FNSPOWN,SEQ=531)

$onText

실제 도함수를 비교하여 부호 있는 검정력 함수를 테스트합니다.
수치적으로 계산된 것들로요.  우리는 범위를 너무 많이 밀지 않습니다
이 테스트에서는 모든 신호가 올바른지 확인하고 싶습니다.
등. 중요한 것은 음수 x로 테스트하여 우리가 알 수 있도록 하는 것입니다.
우리의 구현에서는 심각한 오류가 발생하지 않습니다.

기여자: Steve Dirkse
$offText

$include fnset_xy.inc

옵션 시드 = 1776;

담당자 = 5e-5;
relToInput = 0;

세트
  ex '지수' / e1 * e100 /
  bs '베이스' / b1 * b100 /
  T / t1 * t10000 /
* 예 '지수' / e1 * e5 /
* bs '베이스' / b1 * b3 /
* T / t1 * t15 /
  텍스(T,ex)
  Tbs(T,bs)
  ;
스칼라
  n_ex / [카드(예)] /
  n_bs / [카드(bs)] /
  ;
매개변수
  ㄷ(예)
  x(bs)
  ;
c(예) = 균일(1e-5, 2.8);
x(bs) = 균일(-0.5,15);
x(bs) = 10**x(bs);
Tex(T,ex) = [ord(ex) = (1 + Floor((ord(T)-1)/n_bs))];
Tbs(T,bs) = [ord(bs) = 1 + (mod(ord(T)-1, n_bs))];

data(T,'x') = sumTbs(T,bs), x(bs);
data(T,'y') = sumTex(T,ex), c(ex);

* 첫 번째 테스트에서는 양의 x 값을 사용합니다.

data(T, 'f_' ) = signpower.value( data(T,'x'),data(T,'y'));
data(T, 'fx_' ) = signpower.grad(1: data(T,'x'),data(T,'y'));
data(T, 'fxx_') = signpower.hess(1:1:data(T,'x'),data(T,'y'));
루프 T,
  data(T, 'f') = signpower.value( data(T,'x'),data(T,'y'));
  data(T, 'fx') = signpower.gradn(1: data(T,'x'),data(T,'y'));
  data(T,'fxx') = signpower.hessn(1:1:data(T,'x'),data(T,'y'));
  data(T, 'rc') = mathlastrc;
  data(T, 'ec') = mathlastec;
;

$include fntest_xy.inc

* 이제 테스트를 반복하되 x가 음수인 경우
data(T,'x') = -data(T,'x');

data(T, 'f_' ) = signpower.value( data(T,'x'),data(T,'y'));
data(T, 'fx_' ) = signpower.grad(1: data(T,'x'),data(T,'y'));
data(T, 'fxx_') = signpower.hess(1:1:data(T,'x'),data(T,'y'));
루프 T,
  data(T, 'f') = signpower.value( data(T,'x'),data(T,'y'));
  data(T, 'fx') = signpower.gradn(1: data(T,'x'),data(T,'y'));
  data(T,'fxx') = signpower.hessn(1:1:data(T,'x'),data(T,'y'));
  data(T, 'rc') = mathlastrc;
  data(T, 'ec') = mathlastec;
;

$include fntest_xy.inc