Algorithm 859: BABDCR - a Fortran 90 package for the solution of Bordered ABD linear systems