def _applyConstraints(blockVectorV, factYBY, blockVectorBY, blockVectorY):
 gramYBV = np.dot(blockVectorBY.T, blockVectorV)
 tmp = cho_solve(factYBY, gramYBV)
 blockVectorV -= np.dot(blockVectorY, tmp)

def _applyConstraints ( blockVectorV , factYBY , blockVectorBY , blockVectorY ) : gramYBV = np . dot ( blockVectorBY . T , blockVectorV ) tmp = cho_solve ( factYBY , gramYBV ) blockVectorV -= np . dot ( blockVectorY , tmp )
changes blockvectorv in place .