use mexinterface for functions that use mexErrMsgTxt

drawtheta.for

@@ -43,6 +43,9 @@ C ----------------------------------------------------------------------
 	1                     MEDT,SIGT,yk,IYK,Z,thetaprior,tipo,pdll,
      2                     theta0,theta,NEVAL)
 C INPUT
+#if defined(MATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
+	  USE MEXINTERFACE
+#endif
 	CHARACTER*1 fittizia
 	POINTER (pdll,fittizia)
 	INTEGER nobs,d(2),ny,nz,nx,nu,nv,ns(6),nt,Z(nobs),

drawtheta2.for

@@ -43,6 +43,9 @@ C ----------------------------------------------------------------------
 	1                      MEDT,SIGT,yk,Z,thetaprior,tipo,pdll,
      2                      theta0,theta,NEVAL)
 C INPUT
+#if defined(MATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
+	  USE MEXINTERFACE
+#endif
 	CHARACTER*1 fittizia
 	POINTER (pdll,fittizia)
 	INTEGER nobs,d(2),ny,nz,nx,nu,nv,ns(6),nt,Z(nobs),

findinput.for

@@ -29,6 +29,9 @@ C along with DMM.  If not, see <http://www.gnu.org/licenses/>.
 C -------------------------------------------------------------
 	SUBROUTINE FINDINPUT(STR1,N1,STR2,N2,NUM)
 C INPUT
+#if defined(MATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
+	  USE MEXINTERFACE
+#endif
 	INTEGER N1,N2
 #ifdef __GFORTRAN__
       CHARACTER(len=N1) :: STR1

input.for

@@ -24,6 +24,9 @@ C -----------------------------------------------------------------------
 	1 nstot,np,nf,INFOS,seed,thin,burnin,simulrec,sampler,
      2 datasim,dllname,check,estimation,nt,pdftheta,hyptheta,
      3 hypS,T,obs,Ssampler,hbl,MargLik)
+#if defined(MATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
+	  USE MEXINTERFACE
+#endif
 
 C	INCLUDE 'iosdef.for'
 	INTEGER IERR

kim.for

@@ -60,6 +60,9 @@ C -------------------------------------------------------------------------
       SUBROUTINE KIM(nobs,d,ny,nz,nx,nu,ns,nk,nv,np,INFOS,yk,IYK,
      1 c,H,G,a,F,R,psi,ismoother,XSMOOTH,XSSE,SSMOOTH,INN,LIKE)
 C INPUT
+#if defined(MATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
+	  USE MEXINTERFACE
+#endif
 	INTEGER nobs,ny,nz,nx,nu,ns(6),nk,nv,np,ismoother,
      1 d(2),IYK(nobs,ny+1),INFOS(9,6)
 	DOUBLE PRECISION yk(nobs,ny+nz),c(ny,max(nz,1),ns(1)),

lyapunov.for

@@ -25,6 +25,9 @@ C along with DMM.  If not, see <http://www.gnu.org/licenses/>.
 C -------------------------------------------------------------
 	SUBROUTINE LYAP(nx,nu,compt,F,R,Ps)
 C INPUT
+#if defined(MATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
+	  USE MEXINTERFACE
+#endif
 	INTEGER nx,nu
 	DOUBLE PRECISION compt,F(nx,nx),R(nx,nu)
 C OUTPUT

main.for

@@ -39,6 +39,7 @@ C ------------------------------------------------------------------------------
 #include "fintrf.h"
 #endif
       SUBROUTINE DMMMAIN(FILEIN)
+      USE MEXINTERFACE
 #else
       PROGRAM DMM
 #endif