attempt to implement block diagonal shrunk estimator