설명

여기서 우리는 삼변량에 대한 cppcclib의 외부 함수를 테스트합니다.
정규 분포 PDF는 함수를 비교하여 예상대로 작동합니다.
값을 미리 계산된 값으로, 도함수 값을 수치 도함수로 변환합니다.

기여자: Steve

소형 모델 유형 :슬롯 머신

메인 파일 : cpplib05.gms

$title cppcclib의 외부 함수 테스트(CPPLIB05,SEQ=656)

$onText
여기서 우리는 삼변량에 대한 cppcclib의 외부 함수를 테스트합니다.
정규 분포 PDF는 함수를 비교하여 예상대로 작동합니다.
값을 미리 계산된 값으로, 도함수 값을 수치 도함수로 변환합니다.

기여자: 스티브
$offText
$onDollar

$funcLibIn mvnLib cppcclib

function xdf '삼변량 정규의 PDF' / mvnLib.pdfTVN /;

$설정되지 않은 경우 INFILE $set INFILE tvnInit

$macro f0 xdf.value( data(T,'x'),data(T,'y'),data(T,'z'),data(T,'r'),data(T,'s'),data(T,'t'))
$macro f1(j) xdf.grad (j: data(T,'x'),data(T,'y'),data(T,'z'),data(T,'r'),data(T,'s'),data(T,'t'))
$macro f1n(j) xdf.gradn(j: data(T,'x'),data(T,'y'),data(T,'z'),data(T,'r'),data(T,'s'),data(T,'t'))
$macro f2(j1,j2) xdf.hess (j1:j2:data(T,'x'),data(T,'y'),data(T,'r'),data(T,'r'),data(T,'s'),data(T,'t'))
$macro f2n(j1,j2) xdf.hessn(j1:j2:data(T,'x'),data(T,'y'),data(T,'r'),data(T,'r'),data(T,'s'),data(T,'t'))

$include extrtest3a.inc

fdata(T, 'f_') = data(T,'pdf');

스칼라
 aeps0 '절대 오류 허용 범위: 함수' / 1e-7 /
 reps0 '상대 오류 허용 범위: 함수' / 1e-8 /
 aeps1 '절대 오류 허용 범위: grad' / 2e-2 /
 reps1 '상대 오류 허용오차: grad' / 2e-2 /
 aeps2 '절대 오류 허용 범위: hess' / 1e-2 /
 reps2 '상대 오류 허용오차: hess' / 1e-2 /
 aepsr '절대 오류 허용 오차: hess Xr' / .01 /
 repsr '상대 오류 허용오차: hess Xr' / .15 /
 ;

$include extrtest3b.inc