femwind/restriction_fortran.m

   1 function uc=restriction_fortran(u,hzc,icl3,X,params);
   2 % call fortran version
   3
   4 % get coarse size
   5 n=size(X{1});
   6 [icl1,icl2]=hzc2icl(hzc,n);
   7 nc = [length(icl1),length(icl2),length(icl3)];
   8
   9 %Writing all arrays to text files for use by fortran tester
  10 uc = rand(nc);  ! for sizes only
  11 write_array(swap23(uc),'uc');
  12 write_array(swap23(u),'u');
  13 write_array(swap23(X{1}),'X');
  14 write_array(swap23(X{2}),'Y');
  15 write_array(swap23(X{3}),'Z');
  16 write_array(hzc(:),'hcz');
  17 write_array(icl3(:),'cl_z');
  18
  19 system('./fortran/restriction_test.exe');
  20 uc = swap23(read_array('uc'));
  21
  22 end
  23