Re: TRANSPOSE and MATMUL interaction

Not that it has relevance to the code _as_written_, but it would have if, for instance, fjac was ALLOCATABLE:

SUBROUTINE JacobianFwdDiff(m,n,fjac,ldfjac)

IMPLICIT NONE

INTEGER, INTENT(IN) :: 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
.