def _HCCM2(hessian_inv, scale):
 if (scale.ndim == 1):
  scale = scale[:, None]
 xxi = hessian_inv
 H = np.dot(np.dot(xxi, scale), xxi.T)
 return H

def _HCCM2 ( hessian_inv , scale ) : if ( scale . ndim == 1 ) : scale = scale [ : , None ] xxi = hessian_inv H = np . dot ( np . dot ( xxi , scale ) , xxi . T ) return H
sandwich with ^ * scale * ^ scale is this uses results .