#1235 Added functions for getting upper and lower half-bandwidth of matrices
#1235 Detection of bandwitdh by default in Sundials integrators