femwind/ndt_assembly_test.m

   1 %function err = ndt_assembly_err
   2 format compact
   3 % test for ndt_assembly and ndt_mult
   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=[]
  12 u0=[]
  13
  14
  15 X = regular_mesh(nel,h,expand);
  16 X = add_terrain_to_mesh(X, 'hill', 'squash', 0.1);
  17
  18 K_m = ndt_assembly(A,X,u0,lambda,params,14);
  19 K = ndt_assembly_fortran(A,X,u0,lambda,params,14);
  20
  21 err = norm(K(:) - K_m(:),inf)
  22 %end
  23
  24