subroutine xdlamch (cmach, retval) character cmach double precision retval, dlamch retval = dlamch (cmach) return end