femwind/f_assembly_test.m

   1 format compact
   2 addpath ../netcdf
   3 % test for ndt_f_assembly
   4 nel=[5,4,3];
   5 %nel=[1 1 1];
   6 n=nel+1;
   7 h = [1,1,1];
   8 expand=1.3 ;
   9 A = diag([1 1 1]);
  10 lambda=[];
  11 params.test_fortran=1;
  12 u0={rand(nel), rand(nel), rand(nel)};
  13
  14
  15 X = regular_mesh(nel,h,expand);
  16 X = add_terrain_to_mesh(X, 'hill', 'squash', 0.1);
  17
  18
  19
  20 % test same results for ndt_mult from matlab and fortran
  21
  22 F=f_assembly_fortran(A,X,u0,[],params);