femwind/femwind_rate2_test.m

   1 function p=femwind_rate2_test
   2 disp('basic convergence speed test')
   3 p=femwind_main
   4 p.graphics=0;
   5 p.sc_all = 1;
   6 p.sc2_all =  2;
   7 p=femwind_main(p);
   8 rate = [0 0.069839695192657];
   9 if abs(p.rate(2) - rate(2)) < 1e-8
  10     disp('basic convergence rate test OK')
  11 else
  12     error(sprintf('something changed, expected convergence rate %g got %g diff %g',...
  13        rate,p.rate,rate-p.rate))
  14 end
  15 end
  16