Context Navigation

mncalf.F@ 3302

Visit:

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

Rev	Line
[2403]	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"
	12	CC Called only from MNIMPR. Transforms the function FCN
	13	CC by dividing out the quadratic part in order to find further
	14	CC minima. Calculates ycalf = (f-fmin)/(x-xmin)v(x-xmin)
	15	CC
	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.