## -*- octave -*-
## fails with openblas 0.3.1-1.fc28 on Fedora 28
__octave_config_info__ ("hg_id")
M = magic(10);
A = sin(M);
B = cos(1.2*M);
B=B*B';
opts.v0 = (1:10)';
[Evector, Evalues] = eigs (A, B, 4, "SM", opts);
Afun = @(x) A \ x;
[Evector_f Evalues_f] = eigs (Afun, 10, B, 4, "SM", opts);
assert (Evector, Evector_f);
assert (Evalues, Evalues_f);
