test1.m

%
% Copyright (C) May/June 2005 by Henrik Pagenkopf
%

% test program %

clc
% clear;
close all;
% format short % format long;


% - initialize:

ra = -0.5; rb = -1.0; % floats
ia = -1.0; ib =  2.3;

dim = 3;

A = sqm_rand_init(dim, -1, 1 , -1, 1 );
B = sqm_rand_init(dim, ra, rb, ia, ib);


% -- calculate:

disp('A = '); disp(A);
disp('B = '); disp(B);

[j,k,d]= max_differ_mat(A,B); disp('j,k,d'); disp([floor(j),floor(k),d]);
[j,k,d]= max_coord_dist(A  ); disp('j,k,d'); disp([floor(j),floor(k),d]);

% S = intersect_PR(A, B, 0.005, 3);
  S = eye(dim);

disp('Correcture:'); disp(correct_PR(A, B, S, 1, 2, 3, 3));

% end.