SUBROUTINE calc_S_derivs() | |
INTEGER, DIMENSION(6, 2) :: c_map_mat | |
INTEGER, DIMENSION(:), POINTER:: C_mat | |
DO j=1,3 | |
DO m=j,3 | |
n=n+1 | |
c_map_mat(n,1)=j | |
IF(m==j)CYCLE | |
c_map_mat(n,2)=m | |
END DO | |
END DO | |
DO m=1,6 | |
DO j=1,2 | |
IF(c_map_mat(m,j)==0)CYCLE | |
CALL foo(C_mat(c_map_mat(m,j))) | |
END DO | |
END DO | |
END SUBROUTINE calc_S_derivs |