| Last change
 on this file since 3223 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.