Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
Archives
dmm
Commits
67e80aff
Commit
67e80aff
authored
Oct 21, 2014
by
Houtan Bastani
Browse files
fix conditional ifs to accord with dynare MEX defines
parent
f9814af1
Changes
26
Expand all
Hide whitespace changes
Inline
Side-by-side
amh.for
View file @
67e80aff
...
...
@@ -107,7 +107,7 @@ C LOCALS
id1 = max(1,d(1))
Z0 = Z
delta = 1.D-3
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
amh2.for
View file @
67e80aff
...
...
@@ -107,7 +107,7 @@ C LOCALS
id1 = max(d(1),1)
Z0 = Z
delta = 1.D-3
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
checkdesign.for
View file @
67e80aff
...
...
@@ -70,7 +70,7 @@ C EXTERNAL SUBROUTINES
ALLOCATE(c(ny,max(nz,1),ns(1)),H(ny,nx,ns(2)),
1 G(ny,nu,ns(3)),a(nx,ns(4)),F(nx,nx,ns(5)),R(nx,nu,ns(6))) !,HRG(ny,nu),HRGRH(ny,ny))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
chi2inv.for
View file @
67e80aff
...
...
@@ -430,13 +430,13 @@ C LOCALS
DATA T(198,6:10)/199.8400, 205.0857, 211.0344, 218.6401, 231.8292/
DATA T(199,6:10)/200.8463, 206.1051, 212.0684, 219.6922, 232.9118/
DATA T(200,6:10)/201.8526,207.1244, 213.1022, 220.7441, 233.9943/
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CHARACTER(len=200) :: MEXPRINT
INTEGER*4 mexPrintf
INTEGER*4 mpfout
#endif
IF (V.GT.200) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,*) 'CHI2INV: Too many degrees of freedom'
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
design.for
View file @
67e80aff
...
...
@@ -218,7 +218,7 @@ C
IF (engEvalString(ep, 'clear success;'//
& '[C,H,G,A,F,R]='//TRIM(mfile)//'( ny,nz,nx,'//
& 'nu,ns,theta);'//'success=1;') .ne. 0) then ! engEvalString failed
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL GETERRSTR(matlaberror)
#else
#endif
...
...
@@ -244,7 +244,7 @@ C
#endif
IF (C_ptr .eq. 0) then ! engEvalString failed
buffer=buffer1
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL GETERRSTR(matlaberror)
#else
#endif
...
...
drawtheta.for
View file @
67e80aff
...
...
@@ -105,7 +105,7 @@ C CHEK theta
IF (NN.LE.1000) THEN
GOTO 7777
ELSE
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt('\nReduce skcriterium or use Slice sampling\nProgram aborting\n')
#elif defined(__GFORTRAN__)
WRITE(*,*) ' '
...
...
drawtheta2.for
View file @
67e80aff
...
...
@@ -105,7 +105,7 @@ C CHEK theta
IF (NN.LE.1000) THEN
GOTO 7777
ELSE
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt('\nReduce skcriterium or use Slice sampling\nProgram aborting\n')
#elif defined(__GFORTRAN__)
WRITE(*,*) ' '
...
...
findinput.for
View file @
67e80aff
...
...
@@ -48,7 +48,7 @@ C LOCALS
ERRMSG = 'INPUT ERROR: '// STR2 //' is not set'
IPOS = INDEX(TRIM(STR1),STR2)
IF (IPOS.EQ.0) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt(ERRMSG)
#else
WRITE(*,*) ERRMSG
...
...
@@ -70,7 +70,7 @@ C LOCALS
READ(STR1(IPOS:IPOS2-2),'(F20.10)') NUM
READ(STR1(IPOS:IPOS2-2),*) NUM
ELSE
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt(ERRMSG)
#else
WRITE(*,*) ERRMSG
...
...
forecast.for
View file @
67e80aff
...
...
@@ -73,7 +73,7 @@ C EXTERNAL SUBROUTINES
ALLOCATE(R(nx,nu,ns(6)),c(ny,max(nz,1),ns(1)),H(ny,nx,ns(2)),
1 G(ny,nu,ns(3)),a(nx,ns(4)),F(nx,nx,ns(5)))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
gck.for
View file @
67e80aff
...
...
@@ -110,7 +110,7 @@ C LOCALS
DATA EPS/1.D-14/,ONE/1.0D0/,ZERO/0.0D0/
DOUBLE PRECISION genunf,LEMMA4,MARKOVP
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
gck2.for
View file @
67e80aff
...
...
@@ -101,7 +101,7 @@ C LOCALS
DATA EPS/1.D-14/,ONE/1.0D0/,ZERO/0.0D0/
DOUBLE PRECISION genunf,LEMMA4,MARKOVP
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
innov.for
View file @
67e80aff
...
...
@@ -72,7 +72,7 @@ C LOCALS
1 Xdd(max(d(1),1),nx),Pdd(max(d(1),1),nx,nx),LIKE(max(d(1),1)),
1 XT(nx),PT(nx,nx),INN0(ny))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
innov2.for
View file @
67e80aff
...
...
@@ -71,7 +71,7 @@ C LOCALS
1 Xdd(max(d(1),1),nx),Pdd(max(d(1),1),nx,nx),LIKE(max(d(1),1)),
1 XT(nx),PT(nx,nx))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
input.for
View file @
67e80aff
This diff is collapsed.
Click to expand it.
kim.for
View file @
67e80aff
...
...
@@ -136,7 +136,7 @@ C X-filter initialization
ENDDO
ELSE
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt('\nERROR: d(1)=2 not implemeted yet\n')
#endif
WRITE(*,*) 'ERROR: d(1)=2 not implemeted yet'
...
...
lyapunov.for
View file @
67e80aff
...
...
@@ -63,7 +63,7 @@ C CALL F02EAF('V',nx,T,nx,WR,WI,Z,nx,WORK,LWORK,IFAIL) ! F = ZTZ'
# BWORK,IFAIL)
DO I = 1,nx
IF (WI(I)**2+WR(I)**2.GE.1.D0) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt('\nLYAPUNOV SUBROUTINE: Some parameters out of\nstationary region. Check hyptheta in namelist prior.\nProgram aborting')
#elif defined(__GFORTRAN__)
WRITE(*,*) ' '
...
...
main.for
View file @
67e80aff
...
...
@@ -34,7 +34,7 @@ C
C You should have received a copy of the GNU General Public License
C along with DMM. If not, see <http://www.gnu.org/licenses/>.
C --------------------------------------------------------------------------------------
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
#include "fintrf.h"
SUBROUTINE DMMMAIN(FILEIN)
#else
...
...
@@ -81,7 +81,7 @@ C LOCALS
#ifndef __GFORTRAN__
CHARACTER*200 DMMTITLE
#endif
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CHARACTER(len=200) :: MEXPRINT
INTEGER*4 mexPrintf
INTEGER*4 mpfout
...
...
@@ -105,13 +105,13 @@ C TIME
1 DATE_ITIME)
IT1(1:3) = DATE_ITIME(1:3)
IT1(4:7) = DATE_ITIME(5:8)
#if !defined(M
EX
)
#if !
(
defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FILE)
)
C GET the namelist specified by FILEIN
CALL GETARG(1,FILEIN) ! load name of input file
#endif
C CHECK FILEIN
IF (TRIM(FILEIN).EQ.'') THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt('\nNo input file provided\nProgram Aborting\n')
#elif defined(__GFORTRAN__)
WRITE(*,*) ' '
...
...
@@ -150,12 +150,12 @@ C Assign the name of the matlab file
ALLOCATE( c(ny,max(nz,1),ns(1)),H(ny,nx,ns(2)),
1 G(ny,nu,ns(3)),a(nx,ns(4)),F(nx,nx,ns(5)),R(nx,nu,ns(6)),
1 theta(nt))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL SETFILEM(mfile,pathmfile) ! ONLY THE FIRST TIME
#endif
theta(:) = 1.D0
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#endif
DEALLOCATE(c,H,G,a,F,R,theta)
...
...
@@ -371,7 +371,7 @@ C SIMULATION of DATA and UNOBSERVABLES
C MAXIMUM LIKELIHOOD ESTIMATION
IF ((estimation.EQ.'ML').OR.(estimation.EQ.'ml').OR.
& (estimation.EQ.'Ml').OR.(estimation.EQ.'mL')) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL mexErrMsgTxt('\nMaximum Likelihood inference not allowed\nProgram aborting\n')
#elif defined(__GFORTRAN__)
WRITE(*,*) ' '
...
...
@@ -494,7 +494,7 @@ C MCMC BURN-IN
#if defined(__CYGWIN32__) || defined(_WIN32)
CALL system('cls')
#endif
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1113) jjj,ntf,IMIN(1)/dfloat(jjj),IMAX(1)/dfloat(jjj)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -530,7 +530,7 @@ C MCMC BURN-IN
#if defined(__CYGWIN32__) || defined(_WIN32)
CALL system('cls')
#endif
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1113) jjj,ntf,IMIN(1)/dfloat(jjj),IMAX(1)/dfloat(jjj)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -614,14 +614,14 @@ C MCMC RECORDING phase
#if defined(__CYGWIN32__) || defined(_WIN32)
CALL system('cls')
#endif
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1113) BURNIN,ntf,lastl,lasth
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
WRITE(6,1113) BURNIN,ntf,lastl,lasth
#endif
IF ((HBL.EQ.1).OR.(nv.EQ.0)) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1114) jjj,ntf,IMIN(1)/dfloat(jjj),IMAX(1)/dfloat(jjj)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -629,7 +629,7 @@ C MCMC RECORDING phase
# IMAX(1)/dfloat(jjj)
#endif
ELSEIF ((HBL.GT.1).AND.(nv.GT.0)) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1115) jjj,ntf,IMIN(1)/dfloat(jjj),IMAX(1)/dfloat(jjj),SUM(1.D0-ACCRATE(1:nobs)/DFLOAT(jjj))/DFLOAT(nobs)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -733,14 +733,14 @@ C MCMC RECORDING phase
#if defined(__CYGWIN32__) || defined(_WIN32)
CALL system('cls')
#endif
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1113) BURNIN,ntf,lastl,lasth
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
WRITE(6,1113) BURNIN,ntf,lastl,lasth
#endif
IF ((HBL.EQ.1).OR.(nv.EQ.0)) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1114) jjj,ntf,IMIN(1)/dfloat(jjj),IMAX(1)/dfloat(jjj)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -748,7 +748,7 @@ C MCMC RECORDING phase
# IMAX(1)/dfloat(jjj)
#endif
ELSEIF ((HBL.GT.1).AND.(nv.GT.0)) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1115) jjj,ntf,IMIN(1)/dfloat(jjj),IMAX(1)/dfloat(jjj),SUM(1.D0-ACCRATE(1:nobs)/DFLOAT(jjj))/DFLOAT(nobs)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -821,7 +821,7 @@ C MCMC RECORDING phase
C MARGINAL LIKELIHOOD
IF ((MargLik.EQ.'Y').OR.(MargLik.EQ.'y')) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,*) ' '
mpfout = mexPrintf(MEXPRINT//achar(13))
WRITE(MEXPRINT,*) 'Computing the marginal likelihood. Please wait ...'
...
...
@@ -834,7 +834,7 @@ C MARGINAL LIKELIHOOD
CALL HARMONIC(GGG,nobs,d,ny,nz,nx,nu,nv,ns,nstot,nt,np,
1 INFOS,yk(1:nobs,:),IYK(1:nobs,:),gibtheta,gibZ,
2 thetaprior,psiprior,pdftheta,MLHM)
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,*) 'Modified harmonic mean: done!'
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -843,7 +843,7 @@ C MARGINAL LIKELIHOOD
CALL MENGWONG(GGG,nobs,d,ny,nz,nx,nu,nv,ns,nstot,nt,np,
1 INFOS,yk(1:nobs,:),IYK(1:nobs,:),gibtheta,gibZ,
2 thetaprior,psiprior,pdftheta,MLHM(5,1),MLMW)
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,*) 'Bridge sampling: done!'
mpfout = mexPrintf(MEXPRINT//achar(13))
WRITE(MEXPRINT,*) ' '
...
...
@@ -856,7 +856,7 @@ C MARGINAL LIKELIHOOD
CALL HARMONIC2(GGG,nobs,d,ny,nz,nx,nu,nv,ns,nstot,nt,np,
1 INFOS,yk(1:nobs,:),gibtheta,gibZ,thetaprior,
2 psiprior,pdftheta,MLHM)
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,*) 'Modified harmonic mean: done!'
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -865,7 +865,7 @@ C MARGINAL LIKELIHOOD
CALL MENGWONG2(GGG,nobs,d,ny,nz,nx,nu,nv,ns,nstot,nt,np,
1 INFOS,yk(1:nobs,:),gibtheta,gibZ,thetaprior,
2 psiprior,pdftheta,MLHM(5,1),MLMW)
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,*) 'Bridge sampling: done!'
mpfout = mexPrintf(MEXPRINT//achar(13))
WRITE(MEXPRINT,*) ' '
...
...
@@ -905,7 +905,7 @@ C MARGINAL LIKELIHOOD
IT2(4:7) = DATE_ITIME(5:8)
IT=(IT2(4)-IT1(4))*3600+(IT2(5)-IT1(5))*60+(IT2(6)-IT1(6))
IF ((check.EQ.'Y').OR.(check.EQ.'y')) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1117) TRIM(PATH)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -913,7 +913,7 @@ C MARGINAL LIKELIHOOD
#endif
ELSE
IF ((datasim.EQ.'Y').OR.(datasim.EQ.'y')) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1118) TRIM(PATH)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -922,14 +922,14 @@ C MARGINAL LIKELIHOOD
ELSE
IF ((estimation.EQ.'ML').OR.(estimation.EQ.'ml').OR.
& (estimation.EQ.'Ml').OR.(estimation.EQ.'mL')) THEN
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1119) IT,TRIM(PATH)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
WRITE(6,1119) IT,TRIM(PATH)
#endif
ELSE
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
WRITE(MEXPRINT,1116) IT,TRIM(PATH)
mpfout = mexPrintf(MEXPRINT//achar(13))
#else
...
...
@@ -944,7 +944,7 @@ C MARGINAL LIKELIHOOD
1111 FORMAT((<4>(F25.12)), ' ',A2)
1112 FORMAT(I10,(<np(3)>(F25.12)), ' ',I2)
#endif
#if defined(ME
X
)
#if defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
1113 FORMAT(' Burn-in draws = ',I8,'\n',
# ' Parameters sampled by SLICE ',I5,'\n',
# ' SLICE likelihood eval. Min/Max = ',F6.2, ' / ',F6.2,'\n')
...
...
missing.for
View file @
67e80aff
...
...
@@ -46,7 +46,7 @@ C LOCALS
ALLOCATE(R(nx,nu,ns(6)),c(ny,max(nz,1),ns(1)),H(ny,nx,ns(2)),
1 G(ny,nu,ns(3)),a(nx,ns(4)),F(nx,nx,ns(5)))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
opg.for
View file @
67e80aff
...
...
@@ -85,7 +85,7 @@ C Using Hessian from E04UCF
1 XT(0:nobs,nx),PT(0:nobs,nx,nx),Xdd(max(d(1),1),nx),
1 Pdd(max(d(1),1),nx,nx))
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
@@ -123,7 +123,7 @@ C Using Hessian from E04UCF
DO 1000 I=1,NFREE
THETAV(I) = THETAV(I) + P(I)
theta(IFREE(I)) = THETAV(I)
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
opgh.for
View file @
67e80aff
...
...
@@ -88,7 +88,7 @@ C ------------------------------------------------------------
IFAIL = 0
CALL SYMINV(LTR,NFREE,LTR,W,J,IFAIL,RMAX)
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
@@ -135,7 +135,7 @@ C -----------
ELSE
psi(I-NFT) = PAR(I)
ENDIF
#if defined(ORIGDLL) || defined(ME
X
)
#if defined(ORIGDLL) || defined(M
ATLAB_MEX_FILE) || defined(OCTAVE_MEX_FIL
E)
CALL DESIGN(ny,nz,nx,nu,ns,nt,theta,c,H,G,a,F,R)
#else
#endif
...
...
Prev
1
2
Next
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment