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

Merge pull request #330 from p15-git-acc/left-eig

Browse source 
fix typo R -> R2
This commit is contained in:
Fredrik Johansson 2020-09-21 11:11:52 +02:00 • committed by GitHub
commit 5a71c4b230
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
2 changed files with 29 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
acb_mat/eig_simple_rump.c
View file

@ -67,7 +67,7 @@ acb_mat_eig_simple_rump(acb_ptr E, acb_mat_t L, acb_mat_t R,
if (L != NULL) if (L != NULL)
{ {
if (!result || !acb_mat_inv(L, R, prec)) if (!result || !acb_mat_inv(L, R2, prec))
acb_mat_indeterminate(L); acb_mat_indeterminate(L);
} }

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

@ -92,12 +92,35 @@ int main()
} }
} }
if (n_randint(state, 2))
{
if (algorithm == 0) if (algorithm == 0)
result = acb_mat_eig_simple(F, L, R, A, E, R, prec); result = acb_mat_eig_simple(F, L, R, A, E, R, prec);
else if (algorithm == 1) else if (algorithm == 1)
result = acb_mat_eig_simple_rump(F, L, R, A, E, R, prec); result = acb_mat_eig_simple_rump(F, L, R, A, E, R, prec);
else else
result = acb_mat_eig_simple_vdhoeven_mourrain(F, L, R, A, E, R, prec); result = acb_mat_eig_simple_vdhoeven_mourrain(F, L, R, A, E, R, prec);
}
else
{
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, L, A, prec);
acb_mat_mul(LAR, LAR, R, prec); acb_mat_mul(LAR, LAR, R, prec);