설명

Cholesky 유틸리티가 반환한 테스트 요소

기여자: Steve Dirkse

소형 모델 유형 :슬롯

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

메인 파일 : choles02.gms

$title 테스트 cholesky 유틸리티(CHOLES02,SEQ=412)

$onText

cholesky 유틸리티에서 반환된 테스트 요소

기여자: Steve Dirkse

$offText

I / i1 * i3 / 설정;
별칭(I,J,K);

매개변수
  L0(I,J) '하삼각행렬'
  A(I,J) 'L * L-전치'
  L1(I,J) '계산된 콜레스키 인자'
  오류(I,J) 'L0 - L1'
  ;

테이블 L0
        i1 i2 i3
나는1 4
i2 2 1
i3 1 4 2 ;

A(I,J) = 합 K, L0(I,K)*L0(J,K);

Execute_unload 'a.gdx', i, a;
ExecuteTool.checkErrorLevel 'linalg.cholesky 및 L1 -gdxin=a.gdx -gdxout=b.gdx';
Execute_load 'b.gdx', L1;

오류(I,J) = L0(I,J) - L1(I,J);
스칼라 sumerr; sumerr = sum(I,J), abs(err(I,J));
abort$[sumerr > 1e-8] '잘못된 Cholesky 인자가 반환되었습니다: L0이 예상되었으나 L1을 얻었습니다.', L0, L1, err;

옵션 클리어=L1;
ExecuteTool.checkErrorLevel 'linalg.cholesky i a L1';

오류(I,J) = L0(I,J) - L1(I,J);
sumerr = sum(I,J), abs(err(I,J));
abort$[sumerr > 1e-8] '잘못된 Cholesky 인수가 반환됨: L0이 예상되었으나 L1을 얻었습니다.', L0, L1, err;