hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-05 09:21:38 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

add failing test

Browse source
This commit is contained in:
p15-git-acc 2020-09-20 23:17:34 -05:00
parent 2347c2e13d
commit c88e03b8b6
1 changed files with 28 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

33
acb_mat/test/t-eig_simple.c
View file

@ -92,12 +92,35 @@ int main()
}
}
if (algorithm == 0)
result = acb_mat_eig_simple(F, L, R, A, E, R, prec);
else if (algorithm == 1)
result = acb_mat_eig_simple_rump(F, L, R, A, E, R, prec);
if (n_randint(state, 2))
{
if (algorithm == 0)
result = acb_mat_eig_simple(F, L, R, A, E, R, prec);
else if (algorithm == 1)
result = acb_mat_eig_simple_rump(F, L, R, A, E, R, prec);
else
result = acb_mat_eig_simple_vdhoeven_mourrain(F, L, R, A, E, R, prec);
}
else
result = acb_mat_eig_simple_vdhoeven_mourrain(F, L, R, A, E, R, prec);
{
int r1, r2;
if (algorithm == 0)
{
r1 = acb_mat_eig_simple(F, L, NULL, A, E, R, prec);
r2 = acb_mat_eig_simple(F, NULL, R, A, E, R, prec);
}
else if (algorithm == 1)
{
r1 = acb_mat_eig_simple_rump(F, L, NULL, A, E, R, prec);
r2 = acb_mat_eig_simple_rump(F, NULL, R, A, E, R, prec);
}
else
{
r1 = acb_mat_eig_simple_vdhoeven_mourrain(F, L, NULL, A, E, R, prec);
r2 = acb_mat_eig_simple_vdhoeven_mourrain(F, NULL, R, A, E, R, prec);
}
result = n_randint(state, 2) ? r1 : r2;
}
acb_mat_mul(LAR, L, A, prec);
acb_mat_mul(LAR, LAR, R, prec);