causint - Causal integration
ierr=causint_lop(adj,add,x,y)
standard operator interface
Causal integration
the cdoubint manpage
geef90