source: Sophya/trunk/SophyaExt/CodeMinuit/code/mncalf.F@ 3891

Last change on this file since 3891 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.1 KB
Line 
1*
2* $Id: mncalf.F,v 1.1.1.1 2003-06-11 14:18:26 cmv Exp $
3*
4* $Log: not supported by cvs2svn $
5* Revision 1.1.1.1 1996/03/07 14:31:28 mclareni
6* Minuit
7*
8*
9#include "minuit/pilot.h"
10 SUBROUTINE MNCALF(FCN,PVEC,YCALF,FUTIL)
11#include "minuit/d506dp.inc"
12CC Called only from MNIMPR. Transforms the function FCN
13CC by dividing out the quadratic part in order to find further
14CC minima. Calculates ycalf = (f-fmin)/(x-xmin)*v*(x-xmin)
15CC
16#include "minuit/d506cm.inc"
17 EXTERNAL FCN,FUTIL
18 DIMENSION PVEC(15)
19 NPARX = NPAR
20 CALL MNINEX(PVEC)
21 CALL FCN(NPARX,GIN,F,U,4,FUTIL)
22 NFCN = NFCN + 1
23 DO 200 I= 1, NPAR
24 GRD(I) = 0.
25 DO 200 J= 1, NPAR
26 M = MAX(I,J)
27 N = MIN(I,J)
28 NDEX = M*(M-1)/2 + N
29 200 GRD(I) = GRD(I) + VTHMAT(NDEX) * (XT(J)-PVEC(J))
30 DENOM = 0.
31 DO 210 I= 1, NPAR
32 210 DENOM = DENOM + GRD(I) * (XT(I)-PVEC(I))
33 IF (DENOM .LE. ZERO) THEN
34 DCOVAR = 1.
35 ISW(2) = 0
36 DENOM = 1.0
37 ENDIF
38 YCALF = (F-APSI) / DENOM
39 RETURN
40 END
Note: See TracBrowser for help on using the repository browser.