/dports/math/math77/MATH77/demo/ |
H A D | drdjacg1.f | 84 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 85 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 113 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 114 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 135 FJAC(1,1)=C*Y(2)**2 157 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 158 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 189 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 190 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 192 RE(4,1)=FJAC(1,1)/F2 [all …]
|
H A D | drsjacg1.f | 84 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 85 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 113 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 114 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 135 FJAC(1,1)=C*Y(2)**2 157 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 158 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 189 FJAC(1,1)=(FJAC(1,1)-ACTUAL(1,1))/ACTUAL(1,1) 190 FJAC(1,2)=(FJAC(1,2)-ACTUAL(1,2))/ACTUAL(1,2) 192 RE(4,1)=FJAC(1,1)/F2 [all …]
|
H A D | drdckder.f | 13 double precision FVEC(15),FJAC(LDFJAC,NMAX) variable 19 call DTRG11(N, X, FVEC ,FJAC, 2) 22 call DCKDER(MODE, M, N, X, FVEC, FJAC, LDFJAC, 25 call DTRG11(N, X, FVEC ,FJAC, 1) 29 call DTRG11(N, X, FVEC ,FJAC, 1) 35 * I,FVEC(I),(FJAC(I,J),J=1,N) 47 subroutine DTRG11(N, X, FVEC ,FJAC, IFLAG) argument 53 double precision FJAC(N,N), FVEC(N), SUM, TEMP, X(N) local 70 FJAC(I,J) = TEMP 72 FJAC(J,J) = dble(J+1)*TEMP - cos(X(J))
|
H A D | drsckder.f | 13 real FVEC(15),FJAC(LDFJAC,NMAX) variable 19 call STRG11(N, X, FVEC ,FJAC, 2) 22 call SCKDER(MODE, M, N, X, FVEC, FJAC, LDFJAC, 25 call STRG11(N, X, FVEC ,FJAC, 1) 29 call STRG11(N, X, FVEC ,FJAC, 1) 35 * I,FVEC(I),(FJAC(I,J),J=1,N) 47 subroutine STRG11(N, X, FVEC ,FJAC, IFLAG) argument 53 real FJAC(N,N), FVEC(N), SUM, TEMP, X(N) local 70 FJAC(I,J) = TEMP 72 FJAC(J,J) = real(J+1)*TEMP - cos(X(J))
|
H A D | drdnqsol.f | 15 double precision FJAC(NMAX,NMAX), FNORM, FVEC(NMAX) variable 32 call DNQFJ(N, X, FVEC ,FJAC, 2) 35 call DCKDER(MODE, M, N, X, FVEC, FJAC, NMAX, TEST, 38 call DNQFJ(N, X, FVEC ,FJAC, 1) 41 call DNQFJ(N, X, FVEC ,FJAC, 1) 47 * I,FVEC(I),(FJAC(I,J),J=1,N) 83 subroutine DNQFJ(N, X, FVEC ,FJAC, IFLAG) argument 89 double precision C1(3), C2(3), C3(3), FJAC(N,N), FVEC(N) local 105 FJAC(I,1) = exp(C1(I)*X(1)) * C1(I) 106 FJAC(I,2) = cosh(C2(I)*X(2)) * C2(I) [all …]
|
H A D | drsnqsol.f | 15 real FJAC(NMAX,NMAX), FNORM, FVEC(NMAX) variable 32 call SNQFJ(N, X, FVEC ,FJAC, 2) 35 call SCKDER(MODE, M, N, X, FVEC, FJAC, NMAX, TEST, 38 call SNQFJ(N, X, FVEC ,FJAC, 1) 41 call SNQFJ(N, X, FVEC ,FJAC, 1) 47 * I,FVEC(I),(FJAC(I,J),J=1,N) 83 subroutine SNQFJ(N, X, FVEC ,FJAC, IFLAG) argument 89 real C1(3), C2(3), C3(3), FJAC(N,N), FVEC(N) local 105 FJAC(I,1) = exp(C1(I)*X(1)) * C1(I) 106 FJAC(I,2) = cosh(C2(I)*X(2)) * C2(I) [all …]
|
/dports/math/slatec/ |
H A D | snls1q.f | 90 FJTJ(2) = FJTJ(2) + FJAC(I,1)*FJAC(I,2) 101 TEMP1 = (FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 102 TEMP2 = (FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA 103 TEMP3 = (FJTJ(2)*FJAC(1,2)+FJTJ(3)*FJAC(2,2))/SIGMA 145 FJTJ(2) = FJTJ(2) + FJAC(I,1)*FJAC(I,2) 156 TEMP1 = (FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 157 TEMP2 = (FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA 158 TEMP3 = (FJTJ(2)*FJAC(1,2)+FJTJ(3)*FJAC(2,2))/SIGMA 212 TEMP1 = (FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 213 TEMP2 = (FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA [all …]
|
H A D | dnls1q.f | 91 FJTJ(2) = FJTJ(2) + FJAC(I,1)*FJAC(I,2) 102 TEMP1 = (FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 103 TEMP2 = (FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA 104 TEMP3 = (FJTJ(2)*FJAC(1,2)+FJTJ(3)*FJAC(2,2))/SIGMA 146 FJTJ(2) = FJTJ(2) + FJAC(I,1)*FJAC(I,2) 157 TEMP1 = (FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 158 TEMP2 = (FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA 159 TEMP3 = (FJTJ(2)*FJAC(1,2)+FJTJ(3)*FJAC(2,2))/SIGMA 213 TEMP1 = (FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 214 TEMP2 = (FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA [all …]
|
H A D | dqjac2.f | 2 SUBROUTINE DQJAC2 (N, X, FVEC, FJAC, LDFJAC, IFLAG) argument 21 DOUBLE PRECISION FJAC(LDFJAC,*), FVEC(*), X(*) 23 FJAC(1,1) = -1.0D0 24 FJAC(1,2) = 0.0D0 25 FJAC(2,1) = -2.0D1*X(1) 26 FJAC(2,2) = 1.0D1
|
H A D | sqjac2.f | 2 SUBROUTINE SQJAC2 (N, X, FVEC, FJAC, LDFJAC, IFLAG) argument 20 DIMENSION X(*),FVEC(*),FJAC(LDFJAC,*) 22 FJAC(1,1)=-1.E0 23 FJAC(1,2)=0.E0 24 FJAC(2,1)=-2.E1*X(1) 25 FJAC(2,2)=1.E1
|
H A D | dfcn2.f | 2 SUBROUTINE DFCN2 (IFLAG, M, N, X, FVEC, FJAC, LDFJAC) argument 26 DOUBLE PRECISION FJAC(LDFJAC,*), FVEC(*), X(*) 54 FJAC(I,1) = -TEMP*EXP(TEMP*X(1)) 55 FJAC(I,2) = -TEMP*EXP(TEMP*X(2))
|
H A D | fcn2.f | 2 SUBROUTINE FCN2 (IFLAG, M, N, X, FVEC, FJAC, LDFJAC) argument 26 REAL FJAC(LDFJAC,*), FVEC(*), X(*) 54 FJAC(I,1) = -TEMP*EXP(TEMP*X(1)) 55 FJAC(I,2) = -TEMP*EXP(TEMP*X(2))
|
/dports/math/cmlib/cmlib-3.0_8/test/ |
H A D | snls1e.f | 129 FJTJ(2)=FJTJ(2)+FJAC(I,1)*FJAC(I,2) 139 TEMP1=(FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 140 TEMP2=(FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA 141 TEMP3=(FJTJ(2)*FJAC(1,2)+FJTJ(3)*FJAC(2,2))/SIGMA 178 FJTJ(2)=FJTJ(2)+FJAC(I,1)*FJAC(I,2) 188 TEMP1=(FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 189 TEMP2=(FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA 190 TEMP3=(FJTJ(2)*FJAC(1,2)+FJTJ(3)*FJAC(2,2))/SIGMA 238 TEMP1=(FJTJ(1)*FJAC(1,1)+FJTJ(2)*FJAC(1,2))/SIGMA 239 TEMP2=(FJTJ(1)*FJAC(1,2)+FJTJ(2)*FJAC(2,2))/SIGMA [all …]
|
/dports/math/cmlib/cmlib-3.0_8/src/snls1e/ |
H A D | snls1.f | 639 CALL FCN(IFLAG,M,N,X,FVEC,FJAC,IJUNK) 722 IF (FJAC(J,J) .EQ. ZERO) GO TO 460 725 SUM = SUM + FJAC(I,J)*WA4(I) 727 TEMP = -SUM/FJAC(J,J) 732 FJAC(J,J) = WA1(J) 746 FJAC(I,J) = ZERO 794 WA2(J) = ENORM(J,FJAC(1,J)) 799 IF (FJAC(J,J) .EQ. ZERO) GO TO 540 802 SUM = SUM + FJAC(I,J)*QTF(I) 804 TEMP = -SUM/FJAC(J,J) [all …]
|
H A D | snsq.f | 1 SUBROUTINE SNSQ(FCN,JAC,IOPT,N,X,FVEC,FJAC,LDFJAC,XTOL,MAXFEV,ML, argument 476 CALL JAC(N,X,FVEC,FJAC,LDFJAC,IFLAG) 491 CALL QRFAC(N,N,FJAC,LDFJAC,.FALSE.,IWA,1,WA1,WA2,WA3) 521 IF (FJAC(J,J) .EQ. ZERO) GO TO 110 524 SUM = SUM + FJAC(I,J)*QTF(I) 526 TEMP = -SUM/FJAC(J,J) 528 QTF(I) = QTF(I) + FJAC(I,J)*TEMP 541 R(L) = FJAC(I,J) 551 CALL QFORM(N,N,FJAC,LDFJAC,WA1) 685 SUM = SUM + FJAC(I,J)*WA4(I) [all …]
|
H A D | fdjac1.f | 1 SUBROUTINE FDJAC1(FCN,N,X,FVEC,FJAC,LDFJAC,IFLAG,ML,MU,EPSFCN, argument 92 REAL X(N),FVEC(N),FJAC(LDFJAC,N),WA1(N),WA2(N) 115 FJAC(I,J) = (WA1(I) - FVEC(I))/H 138 FJAC(I,J) = ZERO 140 1 FJAC(I,J) = (WA1(I) - FVEC(I))/H
|
/dports/math/slatec/src/ |
H A D | dnls1.f | 646 CALL FCN(IFLAG,M,N,X,FVEC,FJAC,IJUNK) 733 IF (FJAC(J,J) .EQ. ZERO) GO TO 460 736 SUM = SUM + FJAC(I,J)*WA4(I) 738 TEMP = -SUM/FJAC(J,J) 743 FJAC(J,J) = WA1(J) 757 FJAC(I,J) = ZERO 809 WA2(J) = DENORM(J,FJAC(1,J)) 814 IF (FJAC(J,J) .EQ. ZERO) GO TO 540 817 SUM = SUM + FJAC(I,J)*QTF(I) 819 TEMP = -SUM/FJAC(J,J) [all …]
|
H A D | snls1.f | 651 CALL FCN(IFLAG,M,N,X,FVEC,FJAC,IJUNK) 738 IF (FJAC(J,J) .EQ. ZERO) GO TO 460 741 SUM = SUM + FJAC(I,J)*WA4(I) 743 TEMP = -SUM/FJAC(J,J) 748 FJAC(J,J) = WA1(J) 762 FJAC(I,J) = ZERO 814 WA2(J) = ENORM(J,FJAC(1,J)) 819 IF (FJAC(J,J) .EQ. ZERO) GO TO 540 822 SUM = SUM + FJAC(I,J)*QTF(I) 824 TEMP = -SUM/FJAC(J,J) [all …]
|
H A D | dnsq.f | 2 SUBROUTINE DNSQ (FCN, JAC, IOPT, N, X, FVEC, FJAC, LDFJAC, XTOL, argument 490 CALL JAC(N,X,FVEC,FJAC,LDFJAC,IFLAG) 498 CALL DFDJC1(FCN,N,X,FVEC,FJAC,LDFJAC,IFLAG,ML,MU, 539 IF (FJAC(J,J) .EQ. ZERO) GO TO 130 542 SUM = SUM + FJAC(I,J)*QTF(I) 544 TEMP = -SUM/FJAC(J,J) 546 QTF(I) = QTF(I) + FJAC(I,J)*TEMP 559 R(L) = FJAC(I,J) 569 CALL DQFORM(N,N,FJAC,LDFJAC,WA1) 709 SUM = SUM + FJAC(I,J)*WA4(I) [all …]
|
H A D | snsq.f | 2 SUBROUTINE SNSQ (FCN, JAC, IOPT, N, X, FVEC, FJAC, LDFJAC, XTOL, argument 485 CALL JAC(N,X,FVEC,FJAC,LDFJAC,IFLAG) 500 CALL QRFAC(N,N,FJAC,LDFJAC,.FALSE.,IWA,1,WA1,WA2,WA3) 530 IF (FJAC(J,J) .EQ. ZERO) GO TO 110 533 SUM = SUM + FJAC(I,J)*QTF(I) 535 TEMP = -SUM/FJAC(J,J) 537 QTF(I) = QTF(I) + FJAC(I,J)*TEMP 550 R(L) = FJAC(I,J) 560 CALL QFORM(N,N,FJAC,LDFJAC,WA1) 694 SUM = SUM + FJAC(I,J)*WA4(I) [all …]
|
H A D | dfdjc1.f | 2 SUBROUTINE DFDJC1 (FCN, N, X, FVEC, FJAC, LDFJAC, IFLAG, ML, MU, argument 95 DOUBLE PRECISION EPS, EPSFCN, EPSMCH, FJAC(LDFJAC,*), 120 FJAC(I,J) = (WA1(I) - FVEC(I))/H 143 FJAC(I,J) = ZERO 145 1 FJAC(I,J) = (WA1(I) - FVEC(I))/H
|
H A D | fdjac1.f | 2 SUBROUTINE FDJAC1 (FCN, N, X, FVEC, FJAC, LDFJAC, IFLAG, ML, MU, argument 96 REAL X(*),FVEC(*),FJAC(LDFJAC,*),WA1(*),WA2(*) 120 FJAC(I,J) = (WA1(I) - FVEC(I))/H 143 FJAC(I,J) = ZERO 145 1 FJAC(I,J) = (WA1(I) - FVEC(I))/H
|
H A D | dfdjc3.f | 2 SUBROUTINE DFDJC3 (FCN, M, N, X, FVEC, FJAC, LDFJAC, IFLAG, argument 86 DOUBLE PRECISION X(*),FVEC(*),FJAC(LDFJAC,*),WA(*) 104 CALL FCN(IFLAG,M,N,X,WA,FJAC,LDFJAC) 108 FJAC(I,J) = (WA(I) - FVEC(I))/H
|
H A D | fdjac3.f | 2 SUBROUTINE FDJAC3 (FCN, M, N, X, FVEC, FJAC, LDFJAC, IFLAG, argument 84 REAL X(*),FVEC(*),FJAC(LDFJAC,*),WA(*) 102 CALL FCN(IFLAG,M,N,X,WA,FJAC,LDFJAC) 106 FJAC(I,J) = (WA(I) - FVEC(I))/H
|
/dports/science/siconos/siconos-4.4.0/externals/hairer/ |
H A D | dc_decsol.f | 21 E1(I,J)=-FJAC(I,J) 35 E1(I,J)=-FJAC(I,JM1) 58 E1(I+MLE,J)=-FJAC(I,J) 95 E1(I,J)=-FJAC(I,J) 111 E1(I,J)=-FJAC(I,JM1) 190 E1(J1,J)=-FJAC(J1,J) 194 E1(I,J)=-FJAC(I,J) 226 E2R(I,J)=-FJAC(I,J) 324 E2R(I,J)=-FJAC(I,J) 435 E2R(J1,J)=-FJAC(J1,J) [all …]
|