lit
    JINRLIB

DFUMIL, DLIKLM - ,

D510

: .. :

popaa DFUMIL p χ-apa,
popaa DLIKLM ao ac paooo
cpao a YEXP(I) o ooa opc Y(I), oop ac o oop oopa Xi(1),...,Xi(K) o ooo apapo A(1),...,A(M).
copoo oca a χ-apa oa a :

D510_1

yjexp ± Δyj - cpao p a,
  oop o ooa opc yj,
  ac ooo apapo A;
K - co oopa, oca o cpa o;
M - co apapo;
N - co cpa o.

:

: - pp:
: - DFUMIL, DLIKLM, DERORF
: - DMCONV, DSGZ, DMONIT, DARITM, DSCAL
: - DFUNCT,DARITM - pp .
opopa DSGZ DMONIT o opoac ooa.

:

CALL DFUMIL (S,M,N1,N2,N3,EPS,AKAPPA,ALAMBD,IT,MC), :

S - a poao ;
M - co apapo A;
N1 - acao co aa a oo pa;
N2 - papa oc N2 c pa oop opa ao o apapa ( 4 paa).
o opo N2=1, N1=2;
N3 - aa acc opa DFUMIL oc N3 pa;
EPS - pa ooc oca oa apapo, ooca o oo oa cacc oo apapo.
O cc c oo o poca oo Yj poo, o cpa o ΔYj aa po. c EPS 0, o acc o apapo cc, a pc ca COMMON - oa /SIGMA/SIGMA(100), oop o aa ooa.
EPS=.01 0.1 o ocaoo;
AKAPPA - oa ooc. c a, AKAPPA EPS (o co co, c a pa oac);
ALAMBD - o, p a. c o oaac a, pc oo pa;
IT - acoa a a a:
IT = 0 - aa a oo oc pa.
IT > 0 - aa a pa 0,IT,2IT... oc.
IT < 0 - oa o a (po aocco);
MC - op opaa DFUMIL:
MC= 1 - opaoo opaa,
MC=-1 - a (o o poo),
MC=-2 - o oo co,
MC=-3 - cpa o c pa,
MC=-4 - aop o poo pa pooc (c.o DLIKLM).

a M,N1,N2,N3,EPS,IT o aa p pa DFUMIL.

ooa opa oa aaa COMMON-oa.

/A/A(100) o copa aa a apapo.
a o o.
/PL/PL(100) o copa aa opa a ao
o apapa A.
a o o.
o acpoa apap A(I), o oo PL(I) pa opa.

Bp pa apapo o aa
COMMON-oa /AU/AMX(100) /AL/AMN(100) cooco.
Ec coc caapa pc opopa DSGZ, o
COMMON-oa /NED/,/EXDA/ oa aaa opa.
/NED/N,NS oo aoo c opao:
   N - co cpa o;
   NS=K+2- oco a oo o, K c co X-oopa.
/EXDA/EXDA(1500) oo copa c a

Y1exp
ΔY1
X1(1)
. . .
X1(K)
. . .
YNexp
ΔYN
XN(1)
. . .
XN(K)

opopaa DSGZ(M,S) c a po S,
pa G /G/G(100) ap Z p op poo ( cooa op poo Yj /Z/Z(1275)).
Z acac oo p ooa ap, a aoa, co.
, oocc cpoa apapa, o cpa ap, apa oa.
Z oa ooo opo, o po oc opa ca oa a apao (oa apa op poo aco oop o co).

χ2/2
D510_2

DSGZ ao cpao o opaac opopa DARITM(Y), pca /X/X(10) oopa o.
Ec ooa o pca poo, o o aca opopa DARITM(Y) a

       SUBROUTINE DARITM(Y)
       COMMON /A/A(100),/DF/DF(100),/X/X(10)
       Y=F(X(1)...X(K),A(1)...A(M))
       DF(1)=DY/DA(1)
       . . .
       DF(M)=DY/DA(M)
       RETURN
       END

c caapa pc opopa DARITM(Y), oopa c co poo Yj , co 0.01 ac PL(I) ac aa coo ppoa. a opaac cocao ooa opopa DFUNCT(X) a:

       FUNCTION DFUNCT(X)
       DIMENSION X(10)
       COMMON /A/A(100)
       DFUNCT=F(X(1),...,X(K),A(1),...,A(M))
       RETURN
       END

Oao aaco ppoa a pao aoa cp c o ooc.
Bc a cocpooa opopa DMONIT. ao aao pa ac c :

ITER - opo op pa;
GT - oao c pa;
2S - oo a (X-apa ca oa DFUMIL);
AKAPPA,ALAMBD - (2S+GT/(2*ALAMBD) pcaa ao a 2S);
T1 - aa ao ao pa (c , a 1/T1 pa, c o pa 4, o oop opa ao 4 paa);
PARAMETERS - a apapo;
ERRORS - o oo;
FACTORS - aop opp R(I). SQRT(1-1/R(I)) a oa o opp apapa A(I), oop c aca o opp aoo apapa c o o oa oca apapo. c opo aoo R(I) oa co ao, p p opa ap Z, p c oo oo ao.

Bo DLIKLM (S,M,N1,N2,N3,EPS,AKAPPA,ALAMBD,IT,MC) p opa oap paooo. coc o o aaoo DFUMIL. Eca paa p cooa caapo DSGZ: COMMON o /EXDA/EXDA(1500) o ao opao:

X1(1)
. . .
X1(K)
. . .
XN(K)
. . .
XN(K)

/NED/N,NS NS=K.
Haa apap o a, o c j o Yj o 0. c a oac apapo, oopo c p, oo Yj pa 0, oac aaa apap N1 pa 3 4. c DLIKLM o pc oo a Yj , o oa ppa pa.
B ca caapoo DSGZ, o o DLIKLM cc

D510_3

c Z a apa p op poo S (c pp op poo oooo pacp pooc Y(X/A).

:

  1. Ec coc caapa pc DSGZ, oc opaa DFUMIL ooa o a opopa DERORF(M).
    a opopaa c aa ao cpao o a ooao , o oo a, c c oo ap oo, a X-apa o o o p X, poaa, o o REAL.
  2. oc opaa DFUMIL COMMON o /Z/Z(1275) aoc oa ap ( - oo apapo). a acaa oo , p a po co oo apapo, cpoa op, o apapo, cpoa cao DFUMIL (a pa -a co oo). B oc ca o oo o co a.
    COMMON o /Z0/Z0(1275) aoc opaa apa oo,
    oopo a cpo co, ooc oo apapa, cpoa ooa.
    COMMON o PLU/PL0(100) copac 0. -1. apapo, cpoa a pa, -2. apapo c co oa.
  3. oa o a c oo aoo opa, c a aa opoc, a oopo op ap p op poo opaac .
    Happ, ooa o ca A(I)2 co A(I), c o opac oo oo a apapo.
    o ca o cooa aapa opa apapo (COMMON o /AL/). o o p co cco ac caa S oo , a pao, M.
    acoc ca S = F2(A) ,
    F(A) - pooa .
  4. Caapa pc DSGZ o cooac a, oa Yj(A) pca eoao opa pa j.
    o ca o a ooa oopaa X, aaa , app:
              FUNCTION DFUNCT(X)
              DIMENSION X(10)
              COMMON /A/A(100)
              I=X(2)
              GO TO (1,2),I
           1  DFUNCT=F1(X,A)
              RETURN
           2  DFUNCT=F2(X,A)
              RETURN
              END 
  5. Hoop COMMON o o pacp, c co cpoa apapo o 50 co apapo o 100, a c cpa a aa o 1500 co a co X oopa o 10.
    po opo, o ca apapo a o o copaa, ocoo a c COMMON oo /Z/ /Z0/.
    a ap DFUMIL, DLIKLM, DERORF, DARITM, DSGZ, DMCONV, DMONIT DSCAL o oca paopo DOUBLE PRECISION.

:

  1. C.H. C, .H. C.
    H .
    O, -810, , 1961 .
  2. .H.C. C .
    , 1967 ., 11-3362.
  3. .., .
    ." "
    .., , 1976.
  4. S.N.Dymov, V.S.Kurbatov, I.N.Silin, S.V.Yaschenko,
    "Constrained minimization in C++ environment".
    Nuclear Instuments and Methods in Physics Research A 440(2000)431-437.


home up e-mail