checkint.m

来自「一本很简单」· M 代码 · 共 28 行

% CheckInt checks integer versions of Q's for orthogonality to P'*I.

for d = 0:3
  for j = 1:4
    disp(['---------- d = ' num2str(d) ', j = ' num2str(j) ' -----------']);
    P = FindP(d, j);
    P_denom = LCD(P);
    P = round(P*P_denom);
    disp(['P = 1/' num2str(P_denom)]);
    disp(P);
    I = Inner(d, j);
    I_denom = LCD(I);
    I = round(I*I_denom);
    disp(['I = 1/' num2str(I_denom)]);
    disp(I);
    Q = FindQ(d, j, 'lcd')
    ip = P'*I*Q;
    if all(all(ip == 0))
      disp('******* Passed orthogonality test.');
    else
      disp('******* FAILED!');
    end;
    disp('Press a key...');
    pause;
  end;
end;
return;