BD function

Calculates Brooks-Draper diagnostic