|
1 | | -!*==CHKDER.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 1 | + |
2 | 2 | subroutine chkder(m,n,x,Fvec,Fjac,Ldfjac,Xp,Fvecp,Mode,Err) |
3 | 3 | implicit none |
4 | | -!*--CHKDER4 |
| 4 | + |
5 | 5 | integer m , n , Ldfjac , Mode |
6 | 6 | double precision x(n) , Fvec(m) , Fjac(Ldfjac,n) , Xp(n) , & |
7 | 7 | & Fvecp(m) , Err(m) |
@@ -140,10 +140,10 @@ subroutine chkder(m,n,x,Fvec,Fjac,Ldfjac,Xp,Fvecp,Mode,Err) |
140 | 140 | ! last card of subroutine chkder. |
141 | 141 | ! |
142 | 142 | end |
143 | | -!*==DOGLEG.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 143 | + |
144 | 144 | subroutine dogleg(n,r,Lr,Diag,Qtb,Delta,x,Wa1,Wa2) |
145 | 145 | implicit none |
146 | | -!*--DOGLEG146 |
| 146 | + |
147 | 147 | integer n , Lr |
148 | 148 | double precision Delta |
149 | 149 | double precision r(Lr) , Diag(n) , Qtb(n) , x(n) , Wa1(n) , Wa2(n) |
@@ -321,10 +321,10 @@ subroutine dogleg(n,r,Lr,Diag,Qtb,Delta,x,Wa1,Wa2) |
321 | 321 | ! last card of subroutine dogleg. |
322 | 322 | ! |
323 | 323 | end |
324 | | -!*==DPMPAR.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 324 | + |
325 | 325 | double precision function dpmpar(i) |
326 | 326 | implicit none |
327 | | -!*--DPMPAR327 |
| 327 | + |
328 | 328 | integer i |
329 | 329 | ! ********** |
330 | 330 | ! |
@@ -500,10 +500,10 @@ double precision function dpmpar(i) |
500 | 500 | ! Last card of function dpmpar. |
501 | 501 | ! |
502 | 502 | end |
503 | | -!*==ENORM.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 503 | + |
504 | 504 | double precision function enorm(n,x) |
505 | 505 | implicit none |
506 | | -!*--ENORM506 |
| 506 | + |
507 | 507 | integer n |
508 | 508 | double precision x(n) |
509 | 509 | ! ********** |
@@ -597,11 +597,11 @@ double precision function enorm(n,x) |
597 | 597 | ! last card of function enorm. |
598 | 598 | ! |
599 | 599 | end |
600 | | -!*==FDJAC1.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 600 | + |
601 | 601 | subroutine fdjac1(fcn,n,x,Fvec,Fjac,Ldfjac,Iflag,Ml,Mu,Epsfcn,Wa1,& |
602 | 602 | & Wa2) |
603 | 603 | implicit none |
604 | | -!*--FDJAC1604 |
| 604 | + |
605 | 605 | integer n , Ldfjac , Iflag , Ml , Mu |
606 | 606 | double precision Epsfcn |
607 | 607 | double precision x(n) , Fvec(n) , Fjac(Ldfjac,n) , Wa1(n) , Wa2(n) |
@@ -746,11 +746,11 @@ subroutine fdjac1(fcn,n,x,Fvec,Fjac,Ldfjac,Iflag,Ml,Mu,Epsfcn,Wa1,& |
746 | 746 | ! last card of subroutine fdjac1. |
747 | 747 | ! |
748 | 748 | 99999 end |
749 | | -!*==FDJAC2.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 749 | + |
750 | 750 |
|
751 | 751 | subroutine fdjac2(fcn,m,n,x,Fvec,Fjac,Ldfjac,Iflag,Epsfcn,Wa) |
752 | 752 | implicit none |
753 | | -!*--FDJAC2753 |
| 753 | + |
754 | 754 | integer m , n , Ldfjac , Iflag |
755 | 755 | double precision Epsfcn |
756 | 756 | double precision x(n) , Fvec(m) , Fjac(Ldfjac,n) , Wa(m) |
@@ -855,12 +855,12 @@ subroutine fdjac2(fcn,m,n,x,Fvec,Fjac,Ldfjac,Iflag,Epsfcn,Wa) |
855 | 855 | ! last card of subroutine fdjac2. |
856 | 856 | ! |
857 | 857 | 99999 end |
858 | | -!*==HYBRD.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 858 | + |
859 | 859 | subroutine hybrd(fcn,n,x,Fvec,Xtol,Maxfev,Ml,Mu,Epsfcn,Diag,Mode, & |
860 | 860 | & Factor,Nprint,Info,Nfev,Fjac,Ldfjac,r,Lr,Qtf,Wa1,& |
861 | 861 | & Wa2,Wa3,Wa4) |
862 | 862 | implicit none |
863 | | -!*--HYBRD863 |
| 863 | + |
864 | 864 | integer n , Maxfev , Ml , Mu , Mode , Nprint , Info , Nfev , & |
865 | 865 | & Ldfjac , Lr |
866 | 866 | double precision Xtol , Epsfcn , Factor |
@@ -1313,10 +1313,10 @@ subroutine hybrd(fcn,n,x,Fvec,Xtol,Maxfev,Ml,Mu,Epsfcn,Diag,Mode, & |
1313 | 1313 | ! last card of subroutine hybrd. |
1314 | 1314 | ! |
1315 | 1315 | end |
1316 | | -!*==HYBRD1.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 1316 | + |
1317 | 1317 | subroutine hybrd1(fcn,n,x,Fvec,Tol,Info,Wa,Lwa) |
1318 | 1318 | implicit none |
1319 | | -!*--HYBRD11319 |
| 1319 | + |
1320 | 1320 | integer n , Info , Lwa |
1321 | 1321 | double precision Tol |
1322 | 1322 | double precision x(n) , Fvec(n) , Wa(Lwa) |
@@ -1437,12 +1437,12 @@ subroutine hybrd1(fcn,n,x,Fvec,Tol,Info,Wa,Lwa) |
1437 | 1437 | ! last card of subroutine hybrd1. |
1438 | 1438 | ! |
1439 | 1439 | end |
1440 | | -!*==HYBRJ.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 1440 | + |
1441 | 1441 | subroutine hybrj(fcn,n,x,Fvec,Fjac,Ldfjac,Xtol,Maxfev,Diag,Mode, & |
1442 | 1442 | & Factor,Nprint,Info,Nfev,Njev,r,Lr,Qtf,Wa1,Wa2, & |
1443 | 1443 | & Wa3,Wa4) |
1444 | 1444 | implicit none |
1445 | | -!*--HYBRJ1448 |
| 1445 | + |
1446 | 1446 | integer n , Ldfjac , Maxfev , Mode , Nprint , Info , Nfev , Njev ,& |
1447 | 1447 | & Lr |
1448 | 1448 | double precision Xtol , Factor |
@@ -1877,10 +1877,10 @@ subroutine hybrj(fcn,n,x,Fvec,Fjac,Ldfjac,Xtol,Maxfev,Diag,Mode, & |
1877 | 1877 | ! last card of subroutine hybrj. |
1878 | 1878 | ! |
1879 | 1879 | end |
1880 | | -!*==HYBRJ1.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 1880 | + |
1881 | 1881 | subroutine hybrj1(fcn,n,x,Fvec,Fjac,Ldfjac,Tol,Info,Wa,Lwa) |
1882 | 1882 | implicit none |
1883 | | -!*--HYBRJ11886 |
| 1883 | + |
1884 | 1884 | integer n , Ldfjac , Info , Lwa |
1885 | 1885 | double precision Tol |
1886 | 1886 | double precision x(n) , Fvec(n) , Fjac(Ldfjac,n) , Wa(Lwa) |
@@ -2006,12 +2006,12 @@ subroutine hybrj1(fcn,n,x,Fvec,Fjac,Ldfjac,Tol,Info,Wa,Lwa) |
2006 | 2006 | ! last card of subroutine hybrj1. |
2007 | 2007 | ! |
2008 | 2008 | end |
2009 | | -!*==LMDER.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 2009 | + |
2010 | 2010 | subroutine lmder(fcn,m,n,x,Fvec,Fjac,Ldfjac,Ftol,Xtol,Gtol,Maxfev,& |
2011 | 2011 | & Diag,Mode,Factor,Nprint,Info,Nfev,Njev,Ipvt,Qtf, & |
2012 | 2012 | & Wa1,Wa2,Wa3,Wa4) |
2013 | 2013 | implicit none |
2014 | | -!*--LMDER2020 |
| 2014 | + |
2015 | 2015 | integer m , n , Ldfjac , Maxfev , Mode , Nprint , Info , Nfev , & |
2016 | 2016 | & Njev |
2017 | 2017 | integer Ipvt(n) |
@@ -2465,10 +2465,10 @@ subroutine lmder(fcn,m,n,x,Fvec,Fjac,Ldfjac,Ftol,Xtol,Gtol,Maxfev,& |
2465 | 2465 | ! last card of subroutine lmder. |
2466 | 2466 | ! |
2467 | 2467 | end |
2468 | | -!*==LMDER1.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 2468 | + |
2469 | 2469 | subroutine lmder1(fcn,m,n,x,Fvec,Fjac,Ldfjac,Tol,Info,Ipvt,Wa,Lwa) |
2470 | 2470 | implicit none |
2471 | | -!*--LMDER12477 |
| 2471 | + |
2472 | 2472 | integer m , n , Ldfjac , Info , Lwa |
2473 | 2473 | integer Ipvt(n) |
2474 | 2474 | double precision Tol |
@@ -2622,12 +2622,12 @@ subroutine lmder1(fcn,m,n,x,Fvec,Fjac,Ldfjac,Tol,Info,Ipvt,Wa,Lwa) |
2622 | 2622 | ! last card of subroutine lmder1. |
2623 | 2623 | ! |
2624 | 2624 | end |
2625 | | -!*==LMDIF.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 2625 | + |
2626 | 2626 | subroutine lmdif(fcn,m,n,x,Fvec,Ftol,Xtol,Gtol,Maxfev,Epsfcn,Diag,& |
2627 | 2627 | & Mode,Factor,Nprint,Info,Nfev,Fjac,Ldfjac,Ipvt, & |
2628 | 2628 | & Qtf,Wa1,Wa2,Wa3,Wa4) |
2629 | 2629 | implicit none |
2630 | | -!*--LMDIF2639 |
| 2630 | + |
2631 | 2631 | integer m , n , Maxfev , Mode , Nprint , Info , Nfev , Ldfjac |
2632 | 2632 | integer Ipvt(n) |
2633 | 2633 | double precision Ftol , Xtol , Gtol , Epsfcn , Factor |
@@ -3083,10 +3083,10 @@ subroutine lmdif(fcn,m,n,x,Fvec,Ftol,Xtol,Gtol,Maxfev,Epsfcn,Diag,& |
3083 | 3083 | ! last card of subroutine lmdif. |
3084 | 3084 | ! |
3085 | 3085 | end |
3086 | | -!*==LMDIF1.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 3086 | + |
3087 | 3087 | subroutine lmdif1(fcn,m,n,x,Fvec,Tol,Info,Iwa,Wa,Lwa) |
3088 | 3088 | implicit none |
3089 | | -!*--LMDIF13098 |
| 3089 | + |
3090 | 3090 | integer m , n , Info , Lwa |
3091 | 3091 | integer Iwa(n) |
3092 | 3092 | double precision Tol |
@@ -3219,10 +3219,10 @@ subroutine lmdif1(fcn,m,n,x,Fvec,Tol,Info,Iwa,Wa,Lwa) |
3219 | 3219 | ! last card of subroutine lmdif1. |
3220 | 3220 | ! |
3221 | 3221 | end |
3222 | | -!*==LMPAR.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 3222 | + |
3223 | 3223 | subroutine lmpar(n,r,Ldr,Ipvt,Diag,Qtb,Delta,Par,x,Sdiag,Wa1,Wa2) |
3224 | 3224 | implicit none |
3225 | | -!*--LMPAR3237 |
| 3225 | + |
3226 | 3226 | integer n , Ldr |
3227 | 3227 | integer Ipvt(n) |
3228 | 3228 | double precision Delta , Par |
@@ -3486,12 +3486,12 @@ subroutine lmpar(n,r,Ldr,Ipvt,Diag,Qtb,Delta,Par,x,Sdiag,Wa1,Wa2) |
3486 | 3486 | ! last card of subroutine lmpar. |
3487 | 3487 | ! |
3488 | 3488 | end |
3489 | | -!*==LMSTR.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 3489 | + |
3490 | 3490 | subroutine lmstr(fcn,m,n,x,Fvec,Fjac,Ldfjac,Ftol,Xtol,Gtol,Maxfev,& |
3491 | 3491 | & Diag,Mode,Factor,Nprint,Info,Nfev,Njev,Ipvt,Qtf, & |
3492 | 3492 | & Wa1,Wa2,Wa3,Wa4) |
3493 | 3493 | implicit none |
3494 | | -!*--LMSTR3506 |
| 3494 | + |
3495 | 3495 | integer m , n , Ldfjac , Maxfev , Mode , Nprint , Info , Nfev , & |
3496 | 3496 | & Njev |
3497 | 3497 | integer Ipvt(n) |
@@ -3953,10 +3953,10 @@ subroutine lmstr(fcn,m,n,x,Fvec,Fjac,Ldfjac,Ftol,Xtol,Gtol,Maxfev,& |
3953 | 3953 | ! last card of subroutine lmstr. |
3954 | 3954 | ! |
3955 | 3955 | end |
3956 | | -!*==LMSTR1.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 3956 | + |
3957 | 3957 | subroutine lmstr1(fcn,m,n,x,Fvec,Fjac,Ldfjac,Tol,Info,Ipvt,Wa,Lwa) |
3958 | 3958 | implicit none |
3959 | | -!*--LMSTR13971 |
| 3959 | + |
3960 | 3960 | integer m , n , Ldfjac , Info , Lwa |
3961 | 3961 | integer Ipvt(n) |
3962 | 3962 | double precision Tol |
@@ -4110,10 +4110,10 @@ subroutine lmstr1(fcn,m,n,x,Fvec,Fjac,Ldfjac,Tol,Info,Ipvt,Wa,Lwa) |
4110 | 4110 | ! last card of subroutine lmstr1. |
4111 | 4111 | ! |
4112 | 4112 | end |
4113 | | -!*==QFORM.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 4113 | + |
4114 | 4114 | subroutine qform(m,n,q,Ldq,Wa) |
4115 | 4115 | implicit none |
4116 | | -!*--QFORM4131 |
| 4116 | + |
4117 | 4117 | integer m , n , Ldq |
4118 | 4118 | double precision q(Ldq,m) , Wa(m) |
4119 | 4119 | ! ********** |
@@ -4207,10 +4207,10 @@ subroutine qform(m,n,q,Ldq,Wa) |
4207 | 4207 | ! last card of subroutine qform. |
4208 | 4208 | ! |
4209 | 4209 | end |
4210 | | -!*==QRFAC.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 4210 | + |
4211 | 4211 | subroutine qrfac(m,n,a,Lda,Pivot,Ipvt,Lipvt,Rdiag,Acnorm,Wa) |
4212 | 4212 | implicit none |
4213 | | -!*--QRFAC4228 |
| 4213 | + |
4214 | 4214 | integer m , n , Lda , Lipvt |
4215 | 4215 | integer Ipvt(Lipvt) |
4216 | 4216 | logical Pivot |
@@ -4376,10 +4376,10 @@ subroutine qrfac(m,n,a,Lda,Pivot,Ipvt,Lipvt,Rdiag,Acnorm,Wa) |
4376 | 4376 | ! last card of subroutine qrfac. |
4377 | 4377 | ! |
4378 | 4378 | end |
4379 | | -!*==QRSOLV.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 4379 | + |
4380 | 4380 | subroutine qrsolv(n,r,Ldr,Ipvt,Diag,Qtb,x,Sdiag,Wa) |
4381 | 4381 | implicit none |
4382 | | -!*--QRSOLV4397 |
| 4382 | + |
4383 | 4383 | integer n , Ldr |
4384 | 4384 | integer Ipvt(n) |
4385 | 4385 | double precision r(Ldr,n) , Diag(n) , Qtb(n) , x(n) , Sdiag(n) , & |
@@ -4573,10 +4573,10 @@ subroutine qrsolv(n,r,Ldr,Ipvt,Diag,Qtb,x,Sdiag,Wa) |
4573 | 4573 | ! last card of subroutine qrsolv. |
4574 | 4574 | ! |
4575 | 4575 | end |
4576 | | -!*==R1MPYQ.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 4576 | + |
4577 | 4577 | subroutine r1mpyq(m,n,a,Lda,v,w) |
4578 | 4578 | implicit none |
4579 | | -!*--R1MPYQ4594 |
| 4579 | + |
4580 | 4580 | integer m , n , Lda |
4581 | 4581 | double precision a(Lda,n) , v(n) , w(n) |
4582 | 4582 | ! ********** |
@@ -4667,10 +4667,10 @@ subroutine r1mpyq(m,n,a,Lda,v,w) |
4667 | 4667 | ! last card of subroutine r1mpyq. |
4668 | 4668 | ! |
4669 | 4669 | end |
4670 | | -!*==R1UPDT.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 4670 | + |
4671 | 4671 | subroutine r1updt(m,n,s,Ls,u,v,w,Sing) |
4672 | 4672 | implicit none |
4673 | | -!*--R1UPDT4688 |
| 4673 | + |
4674 | 4674 | integer m , n , Ls |
4675 | 4675 | logical Sing |
4676 | 4676 | double precision s(Ls) , u(m) , v(n) , w(m) |
@@ -4874,10 +4874,10 @@ subroutine r1updt(m,n,s,Ls,u,v,w,Sing) |
4874 | 4874 | ! last card of subroutine r1updt. |
4875 | 4875 | ! |
4876 | 4876 | end |
4877 | | -!*==RWUPDT.spg processed by SPAG 6.72Dc at 03:59 on 19 Sep 2021 |
| 4877 | + |
4878 | 4878 | subroutine rwupdt(n,r,Ldr,w,b,Alpha,Cos,Sin) |
4879 | 4879 | implicit none |
4880 | | -!*--RWUPDT4895 |
| 4880 | + |
4881 | 4881 | integer n , Ldr |
4882 | 4882 | double precision Alpha |
4883 | 4883 | double precision r(Ldr,n) , w(n) , b(n) , Cos(n) , Sin(n) |
|
0 commit comments