Minuit实例，Fortran-callable 模式

最新推荐文章于 2022-02-15 00:03:30 发布

wjfqvi

最新推荐文章于 2022-02-15 00:03:30 发布

阅读量694

点赞数 1

分类专栏： Fortran

本文链接：https://blog.csdn.net/weixin_43880667/article/details/84653458

版权

本文介绍如何在Fortran程序中调用Minuit库进行数据分析，通过实例展示了Fortran-callable模式的应用，帮助读者理解Minuit在优化和拟合方面的功能。

摘要由CSDN通过智能技术生成

program main

	implicit none 
	integer,parameter	:: NER=kind(0.0d0)
	real(NER)	:: aa,bb,cc,dd
	integer		:: NPAR=2,IFLAG
	real(NER)	:: XVAL(2)=0,GRAD(2)=0 
	real(NER)	:: FVAL
	integer		:: errflag,ivarbl
	real		:: argmn(5)=0
	real		:: parerr, parup, parlow
	character 	:: paraname(2)
	integer		:: nn
	real(NER)	:: step=0.00001_NER
	real(NER)	:: limit=0.0_NER

	external futil
	external fcn
	external MNINIT, MNEXCM, MNPARM, MNSETI, MNPOUT


! default in, out, save unit is 5, 6, 7
! 	open(unit=5,file="in.txt")
! 	open(unit=6,file="out.txt")
! 	open(unit=7,file="save.txt")

	xval=(/0.001_NER,0.001_NER/)

	paraname(1)= "aaa"
	paraname(2)= "bbb"


	call MNINIT (5, 6, 7)
	call MNSETI ('FCN')
    argmn(1) = -1.0_NER
    call MNEXCM (FCN, 'SET PRINTOUT', argmn, 1, errflag, FUTIL)


	CALL MNPARM(1,paraname(1),xval(1),step,limit,limit,errflag)
	CALL MNPARM(2,paraname(