...>          if (h .eq. zero) h = eps          x(j) = tmp + h          call fcn(m,n,x,wa)          x(j) = tmp          do  i = 1, m             fjac(i,j) = (wa(i) - fvec(i))/h          enddo        enddo or check the archives on "numerical derivatives" in the NA ng. -- "The great thing about democracy is that it gives every ...

...)          if (h .eq. zero) h = eps          x(j) = tmp + h          call fcn(m,n,x,wa)          x(j) = tmp          do  i = 1, m             fjac(i,j) = (wa(i) - fvec(i))/h          enddo        enddo or check the archives on "numerical derivatives" in the NA ng. -- "The great thing about democracy is that it gives every voter   a ...

... h = eps*abs(tmp) if (h .eq. zero) h = eps x(j) = tmp + h call fcn(m,n,x,wa) x(j) = tmp do i = 1, m fjac(i,j) = (wa(i) - fvec(i))/h enddo enddo or check the archives on "numerical derivatives" in the NA ng. -- "The great thing about democracy is that it gives every voter ...

... = x(j) h = eps*abs(tmp) if (h .eq. zero) h = eps x(j) = tmp + h call fcn(m,n,x,wa) x(j) = tmp do i = 1, m fjac(i,j) = (wa(i) - fvec(i))/h enddo enddo or check the archives on "numerical derivatives" in the NA ng. -- "The great thing about democracy is that it gives every voter a chance to do ...

...(tnTkk9*x[p8Jv['TBcFb;wTbU_[OK.#,TAt#jjDg?*xDpF]d"3r(7esK^,(X'AdV,Luc4"<KjiHKJ`4a[P:+Gt%%)]n:&sYA`AD$UQG[5Vp,DGpbiJdT`AB@!8Ca)9fvlZN.J["yr>AQL&Vt3V3wx":uyhY,EF3(34.iQfRc-G0:&KyK%%s<q^'Dri*5+e?j[q^Kk#g&,-Em2lp#PDx<\XUU*82L(Ga(8'HwS5N,ja%%p8hyWwAx+3k610V&;A"'y7IREJZ)<Yol4=B$^`6#>ogs&8Tf&/w&$1ExJc.$hT]*\![C,*cx!ywBk+o2h"scsHn2kyAx@fJQw9.F4.<OI`#imfti_>Zr3-?FJaC/5h'_b9`a$/1$^"*Y#`bQSG[IadCFp;-*/J.*=:

... relevance to the code _as_written_, but it would have if, for instance, fjac was ALLOCATABLE: SUBROUTINE JacobianFwdDiff(m,n,fjac,ldfjac) IMPLICIT NONE... m,n,ldfjac REAL(8), DIMENSION(ldfjac,n), INTENT(OUT) :: fjac The code below makes it clear that INTENT is at least IN, so this should be either IN only or INOUT. M1=TRANSPOSE(fjac(1:m,1:n)); M2=fjac(1:m,1:n) Jan

... the following: SUBROUTINE JacobianFwdDiff(m,n,fjac,ldfjac) IMPLICIT NONE INTEGER, INTENT(IN)...), DIMENSION(ldfjac,n), INTENT(OUT) :: fjac REAL(8), DIMENSION(n,n) ...m,n) :: M2 ! Code which computes fjac(1:m,1:n) M1=TRANSPOSE...(1:m,1:n)); M2=fjac(1:m,1:n) M3=...(1:n,1:n)=MATMUL(TRANSPOSE(fjac(1:m,1:n)), fjac(1:m,1:n)) Specifically, when m=ldfjac=769, n=1...

... the following: SUBROUTINE JacobianFwdDiff(m,n,fjac,ldfjac) IMPLICIT NONE INTEGER, INTENT(IN) ...(8), DIMENSION(ldfjac,n), INTENT(OUT) :: fjac REAL(8), DIMENSION(n,n) :: ...(m,n) :: M2 ! Code which computes fjac(1:m,1:n) M1=TRANSPOSE(fjac(1:m,1:n)); M2=fjac(1:m,1:n) M3=MATMUL(... term1(1:n,1:n)=MATMUL(TRANSPOSE(fjac(1:m,1:n)), fjac(1:m,1:n)) Specifically, when m=ldfjac=769, n=1,...