femwind/hexa_fortran.m

   1 function err=hexa_fortran(A,X,u0)
   2
   3 [Kloc,Floc,Jg]=hexa(A,X,u0)
   4
   5 iflags = [2 1 1]
   6 iflags = iflags(:)
   7 write_array(A,'A')
   8 write_array(X,'X')
   9 write_array(u0,'u0')
  10 write_array(iflags,'iflags')
  11
  12 system('./fortran/hexa_test.exe')
  13
  14 Kloc_f=read_array('Kloc');
  15 Floc_f=read_array('Floc');
  16 Jg_f=read_array('Jg');
  17
  18 if iflags(1) == 1
  19     Kloc_err = norm(Kloc-Kloc_f,1)
  20     err = Kloc_err;
  21 else
  22     if iflags(1) == 2
  23         Floc_err = norm(Floc-Floc_f,1)
  24         err = Floc_err;
  25     else
  26         Jg_err = norm(Jg-Jg_f,1)
  27         err = Jg_err;
  28     end
  29 end
  30
  31 % compare the result of matlab hexa and fortran hexa
  32
  33 end