femwind/ndt_assembly_fortran.m

   1 function err=ndt_assembly_fortran(A,X,u0,lambda,params,st)
   2 % call fortran version and compare results
   3
   4 sz = size(X{1});
   5 disp(['ndt_assembly_fortran, mesh size ',num2str(sz),' storage type ',num2str(st)])
   6 if st ~= 14,
   7     error('ndt_assembly_fortran: only storage type 14 supported')
   8 end
   9
  10 K_m=ndt_assembly(A,X,[],[],params,st);
  11
  12 %Writing all arrays to text files for use by fortran tester
  13 write_array_nd(A,'A');
  14 write_array_nd(swap23(X{1}),'X');
  15 write_array_nd(swap23(X{2}),'Y');
  16 write_array_nd(swap23(X{3}),'Z');
  17
  18 exe='./fortran/ndt_assembly_test.exe';
  19 disp(['ndt_assembly done, calling ',exe])
  20
  21 system(exe);
  22 K = swap23(read_array_nd('K'));
  23
  24 err = norm(K_m(:) - K(:),inf);
  25