## -*- octave -*-
## fails with reference blas 3.8.0-9 on Fedora 28
##
__octave_config_info__ ("hg_id")
M = magic(10);
A = sin(M)     + 1i * cos(2.1*M);
B = cos(1.2*M) + 1i * cos(3.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);
