quicwind/mat_gen_wind_flux_div_test.m

   1 function err = mat_gen_wind_flux_div_test
   2 % test mat_gen_wind_flux_div against mat_wind_flux_div
   3
   4 % set
   5 %    mesh_len   vector size 3 of mesh dimensions
   6 %    h          vector size 3 of mesh steps
   7 mesh_len=[2,2,2];
   8 h=rand(1,3);
   9 disp('mat_wind_flux_div_test')
  10 % test mesh with a hill terrain (unit vectors for now)
  11
  12 X = regular_mesh(mesh_len,h,1.2);
  13 X = add_terrain_to_mesh(X,'hill','shift',0.4);
  14 %X = add_terrain_to_mesh(X,'hill','squash',0.1);
  15 % matrix of wind flux divergence
  16 disp('Time for new matrix generation:')
  17 tic
  18 DM = mat_gen_wind_flux_div(X);
  19 size(DM)
  20 toc
  21 disp('Time for old matrix generation:')
  22 tic
  23 DM_old = mat_wind_flux_div(X);
  24 toc
  25 disp('Compare both methods')
  26 err = big(DM - DM_old);