femwind/nd_test.m

   1 % test for nd_assembly and nd_mult
   2 nel=[5,4,3]
   3 h = [1,1,1]
   4 expand=1.3
   5 A = diag([1 1 1])
   6 lambda=[]
   7 params=[]
   8
   9 X = regular_mesh(nel,h,expand);
  10
  11
  12 plot_mesh_3d(X)
  13
  14 K=nd_assembly(A,X,lambda,params);
  15
  16 x=ones(size(X{1}));
  17
  18 y=nd_mult(K,x);
  19
  20 err_zero=big(y)  % should be zero