1 SUBROUTINE fdjac_parallel(j, fcn)
2 USE fdjac_mod, m=>mp, n=>np, ncnt=>ncntp
12 REAL(rprec),
PARAMETER :: zero=0
20 REAL(rprec) :: temp, temp2, h, enorm
27 IF (eps .EQ. zero) stop
'EPS = 0 in fdjac_parallel!'
31 IF (h .eq. zero) h = eps
36 CALL fcn (m, n, xp, wap, iflag, ncnt)
39 WRITE(6,
'(2x,i6,7x,1es12.4)') ncnt+j, temp2**2
44 WRITE (j+1000) j, iflag, h, temp2
60 END SUBROUTINE fdjac_parallel