Last change
on this file since 3049 was 2403, checked in by cmv, 22 years ago |
Creation du module de code source de MINUIT (CERNLIB) extrait par CMV
cmv 11/06/2003
|
File size:
1.2 KB
|
Rev | Line | |
---|
[2403] | 1 | *
|
---|
| 2 | * $Id: mneval.F,v 1.1.1.1 2003-06-11 14:18:27 cmv Exp $
|
---|
| 3 | *
|
---|
| 4 | * $Log: not supported by cvs2svn $
|
---|
| 5 | * Revision 1.1.1.1 1996/03/07 14:31:29 mclareni
|
---|
| 6 | * Minuit
|
---|
| 7 | *
|
---|
| 8 | *
|
---|
| 9 | #include "minuit/pilot.h"
|
---|
| 10 | SUBROUTINE MNEVAL(FCN,ANEXT,FNEXT,IEREV,FUTIL)
|
---|
| 11 | #include "minuit/d506dp.inc"
|
---|
| 12 | CC Evaluates the function being analyzed by MNCROS, which is
|
---|
| 13 | CC generally the minimum of FCN with respect to all remaining
|
---|
| 14 | CC variable parameters. Common block /MN7XCR/ contains the
|
---|
| 15 | CC data necessary to know the values of U(KE1CR) and U(KE2CR)
|
---|
| 16 | CC to be used, namely U(KE1CR) = XMIDCR + ANEXT*XDIRCR
|
---|
| 17 | CC and (if KE2CR .NE. 0) U(KE2CR) = YMIDCR + ANEXT*YDIRCR
|
---|
| 18 | #include "minuit/d506cm.inc"
|
---|
| 19 | CC
|
---|
| 20 | EXTERNAL FCN,FUTIL
|
---|
| 21 | U(KE1CR) = XMIDCR + ANEXT*XDIRCR
|
---|
| 22 | IF ( KE2CR .NE. 0) U(KE2CR) = YMIDCR + ANEXT*YDIRCR
|
---|
| 23 | CALL MNINEX(X)
|
---|
| 24 | NPARX = NPAR
|
---|
| 25 | CALL FCN(NPARX,GIN,FNEXT,U,4,FUTIL)
|
---|
| 26 | NFCN = NFCN + 1
|
---|
| 27 | IEREV = 0
|
---|
| 28 | IF (NPAR .GT. 0) THEN
|
---|
| 29 | ITAUR = 1
|
---|
| 30 | AMIN = FNEXT
|
---|
| 31 | ISW(1) = 0
|
---|
| 32 | CALL MNMIGR(FCN,FUTIL)
|
---|
| 33 | ITAUR = 0
|
---|
| 34 | FNEXT = AMIN
|
---|
| 35 | IF (ISW(1) .GE. 1) IEREV = 1
|
---|
| 36 | IF (ISW(4) .LT. 1) IEREV = 2
|
---|
| 37 | ENDIF
|
---|
| 38 | RETURN
|
---|
| 39 | END
|
---|
Note:
See
TracBrowser
for help on using the repository browser.