femwind/femwind_with_graphics_test.m

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