def test_1(self): n = 4 A = rand_sym(n) w1, V1 = eigh(A) Ad = scdeflate(A, V1[:, 0]) w2, V2 = eigh(Ad) self.assert_deflated(w1, V1, w2, V2)
def test_2(self): n = 4 A = rand_sym(n) B = rand_pos_semidef(n) w1, V1 = eigh(A, B) Ad = scdeflate(A, V1[:, 0]) w2, V2 = eigh(Ad, B) self.assert_deflated(w1, V1, w2, V2)