function result = nc_cf_cor2bounds_test

OK = [];

% vectors

   a = [1 2 3 4];
   b = [1 2 3 4]';
   
   bc=nc_cf_cor2bounds(b); 
   b2=nc_cf_bounds2cor(bc);
   
   ac=nc_cf_cor2bounds(a); 
   a2=nc_cf_bounds2cor(ac);
   
   OK(end+1) = all(a==a2);
   OK(end+1) = all(b==b2);
   
% matrices

   A=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16];
   B=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16]';

   AC=nc_cf_cor2bounds(A);
   A2=nc_cf_bounds2cor(AC);

   BC=nc_cf_cor2bounds(B);
   B2=nc_cf_bounds2cor(BC);
   
   OK(end+1) = all(all(A==A2));
   OK(end+1) = all(all(B==B2));

   result = all(OK);