Add the perturb() function and its tests.

[octave.git] / tests / perturb_tests.m

A = rand(10);
A_hat = perturb(A, 1e-12);
dA = A_hat - A;

unit_test_equals("dA has the correct infinity norm", ...
                 1e-12, ...
                 norm(dA, 'inf'));