| 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 |