femwind/ndt_mult_fortran.m

   1 function err=ndt_mult_fortran(kmat,u)
   2
   3 y   = ndt_mult(kmat,u);
   4
   5 % test vs sparse first
   6 kmat_s = ndt_convert(kmat,'sparse');
   7 y_s = kmat_s*u(:);
   8 err = big(y_s - y(:));
   9 if err > eps*10*big(u)*big(kmat)
  10     err, error('error ndt multiply vs sparse too large')
  11 end
  12
  13 % test vs fortran ndt_mult
  14
  15 write_array_nd(swap23(kmat),'kmat');
  16 write_array_nd(swap23(u),'u');
  17
  18 system('./fortran/ndt_mult_test.exe');
  19
  20 y_f=swap23(read_array_nd('y'));
  21
  22 err=norm(y(:)-y_f(:),inf);
  23
  24 end