- Timestamp:
- 11/19/09 11:29:41 (14 years ago)
- Branches:
- master
- Children:
- 38d77eb
- Parents:
- 6650a56
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
mincjg.f
r6650a56 r32289cd 4 4 ! 5 5 ! Copyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann, 6 ! Shura Hayryan, Chin-Ku 6 ! Shura Hayryan, Chin-Ku 7 7 ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann, 8 8 ! Jan H. Meinke, Sandipan Mohanty … … 40 40 ! .................................................................... 41 41 42 implicit real*8 (a-h,o-z) 43 implicit integer*4 (i-n) 42 double precision amf, tol, eps, gsqrd, fmin, f, t, g, ga, d, xa, x 43 double precision gsq2, gnew, dfpr, stmin, fuit, gdit, gt, gdmi 44 double precision sbound, stepch, step, wo, sum, fch, acur, ddspln 45 double precision gspln, beta, gama, yt, dt, gamden, di 46 47 integer mxfcon, maxlin, mxn, ier, iter, iterfm, iterrs, nfun 48 integer nfopt, i, n, nfbeg, iretry, maxfun 44 49 45 50 parameter (AMF = 10.d0, … … 102 107 103 108 2 step = stmin + stepch 104 wo = 0.d0 109 wo = 0.d0 105 110 106 111 do i=1,n … … 182 187 183 188 if ( fch .ne. 0.d0 ) ddspln = ddspln + (wo + wo) / stepch 184 189 185 190 if ( gdmi .eq. 0.d0 ) goto 6 186 191
Note:
See TracChangeset
for help on using the changeset viewer.