Code covered by the BSD License  

Highlights from
matlab2fmex

image thumbnail

matlab2fmex

by

 

03 Feb 2004 (Updated )

matlab2fmex.m is a small translator which aims to convert numerical M-files to Fortran90 mex.

Editor's Notes:

This author's upload was somehow corrupted. We fixed this on March 29, 2006. At that time, we also converted it from a GZ file a ZIP file for better cross platform portability.

[modlist,funlist,fun_info]=acosh_make(typestrlist)
function [modlist,funlist,fun_info]=acosh_make(typestrlist)
% fun_info is a 1x2 cell. fun_info{1}=>1 result is pointer, 0=>not
%                         fun_info{2}=>typestr ('r' or 'c' usually)
declare_globals
funname='acosh';
modlist='';funlist='';  r=[char(10)];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
modlist=[modlist,...
	 '       interface ',funname,'f',r];
if ~iscell(typestrlist)
 typestrlist={typestrlist};
end
for i=1:length(typestrlist)
 typestr=typestrlist{i};
  for m=1:length(typestr)
  temp=any(strcmp(typestr(m),{'d','m'})); if temp,typestr(m)='c';end
  temp=any(strcmp(typestr(m),{'e','n'})); if temp,typestr(m)='r';end
  temp=any(strcmp(typestr(m),{'f','o'})); if temp,typestr(m)='i';end
  temp=any(strcmp(typestr(m),{'g','p'})); if temp,typestr(m)='l';end
 end
 [modlist,funlist,len]=makeheader(funname,2,typestr,modlist,funlist,r); 
 arr={'r';'c';'i'}; noi={'r';'c';'s';'t'}; com={'c';'t'};
 temp1='complex';
 if strcmp(typestr,arr)
  if strcmp(typestr,'c')   fun_info{2}='c';  else   fun_info{2}='r';  end
 else
  if strcmp(typestr,'t')   fun_info{2}='t';  else   fun_info{2}='s';  end
 end
 if (length(findstr(typestr,'c'))+length(findstr(typestr,'t')))==0
  temp1='real';
 end
 if any(strcmp(typestr,arr))
  funlist=[funlist,'       ',temp1,', dimension(size(in_1,1),size(in_1,2)) :: out',r];
 else
  funlist=[funlist,'       ',temp1,' :: out',r];
 end

 funlist=[funlist,'       real*4 out1',r];
 funlist=[funlist,'       out1=42.2',r];
 
 %And now the kernel %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 kernel='';
 if any(strcmp(typestr,com))
  kernel=[kernel,'        out=log(in_1+sqrt(-1.0+in_1**2))',r];
 else
  kernel=[kernel,'        out=real(log(in_1+sqrt(-1.0+in_1**2)))',r];
  if any(strcmp(typestr,arr))
   kernel=[kernel,'        if (any(abs(in_1)<1.0)) then',r];
  else
   kernel=[kernel,'        if ((abs(in_1)<1.0)) then',r];
  end
  kernel=[kernel,'         call mexPrintf(''Real input to ',funname,' has value < 1 : possible errors'')',r];
  kernel=[kernel,'         print *,''isnan(out)'',isnan(out)',r];
%  kernel=[kernel,'         print *,''hey'',isnanf(out)',r];
  kernel=[kernel,'        endif',r];
 end
 funlist=[funlist,kernel];

 %End of the function %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 funlist=[funlist,'       end function ',funname,'f_',typestr,r];
 if strcmp(typestr,typestrlist{length(typestrlist)})
 else
  funlist=[funlist,r];
 end
end
%And end the module %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
modlist=[modlist...
	 '       end interface ',funname,'f',char(10)];
fun_info{1}=0;

Contact us