test_power_vector2.m

Go to the documentation of this file.

function OK = test_power_vector2
%function OK = test_power_vector2
% function testing the powervectors and its derivatives, i.e. check
% whether finite difference approximates the derivative.

% Bernard Haasdonk 28.8.2009

OK = 1;

pdeg = 4;

for n = 1:100
  x = rand(1,2);
  h = 1e-7;
  
  for p = 1:pdeg
    D = power_vector2_derivative(x,p);
    
    d1 = (power_vector2(x+[h,0],p)-power_vector2(x,p))/h;
    d2 = (power_vector2(x+[0,h],p)-power_vector2(x,p))/h;
    
    Dappr = [d1,d2]; 
    
    maxerr = max(max(abs(Dappr-D)));
%    disp(['p = ',num2str(p),', maxerr =',num2str(maxerr)]);
    
    if (maxerr>1e-5)
      OK = 0;
      return;
    end;
    
  end;
end;%| \docupdate