Use subroutine DENMAT dirac/dirden.F.
Use subroutine REACMO_NO_WORK in dirac/dirgp.F.
They are blocked according to NBBAS and IBBAS.
Most of the routines end with:
subroutine foo()
...
return
end
The “return” is redundant. Simply write “end” without “return”. To make it better, we recommend programmers to stick to the Fortran90/95 coding standards. So the structure of the routine shall be as:
subroutine foo()
...
end subroutine
Historical reasons. They have no meaning today.
There is nothing wrong with zgemm. You can use either one - whichever you prefer.
If you want to be 100% independent of the DIRAC infrastructure then you cannot use such routines.