Commit e83a5b66 authored by Peter W. Draper's avatar Peter W. Draper
Browse files

Update tolerances for AMD Rome with optimization and avx2 enabled

parent 30e3e622
......@@ -45,22 +45,22 @@ int main(int argc, char *argv[]) {
printf("Absolute difference too large !\n");
error = 1;
}
if (abs > 1.2e-7 && fabsf(x) <= 0.1) {
if (abs > 3e-7 && fabsf(x) <= 0.1) {
printf("Absolute difference too large !\n");
error = 1;
error = 2;
}
if (rel > 1e-6 && fabsf(x) <= 0.2) {
printf("Relative difference too large !\n");
error = 1;
error = 3;
}
if (rel > 4e-8 && fabsf(x) <= 0.1) {
if (rel > 3e-7 && fabsf(x) <= 0.1) {
printf("Relative difference too large !\n");
error = 1;
error = 4;
}
if (error) {
printf("%2d: x= %f exp(x)= %e approx_exp(x)=%e abs=%e rel=%e\n", i, x,
if (error > 0) {
printf("%2d/%d: x= %f exp(x)= %e approx_exp(x)=%e abs=%e rel=%e\n", i, error, x,
exp_correct, exp_approx, abs, rel);
return 1;
}
......
......@@ -200,15 +200,24 @@ int main(int argc, char *argv[]) {
/* Check the total surface area */
assert(fabs(Atot - 1.0f) < 1.e-6);
/* Check the neighbour relations for an arbitrary cell: cell 44
We plotted the grid and manually found the correct neighbours and their
order. */
assert(cells[44].nvert == 4);
assert(cells[44].ngbs[0] == 34);
assert(cells[44].ngbs[1] == 45);
assert(cells[44].ngbs[2] == 54);
assert(cells[44].ngbs[3] == 43);
/* Check the neighbour relations for an arbitrary cell: cell 44 We plotted
the grid and manually found the correct neighbours and their
order. Variation is found when optimizing, so we have two possible
outcomes... */
if (cells[44].nvert == 5) {
assert(cells[44].nvert == 5);
assert(cells[44].ngbs[0] == 43);
assert(cells[44].ngbs[1] == 34);
assert(cells[44].ngbs[2] == 45);
assert(cells[44].ngbs[3] == 55);
} else {
assert(cells[44].nvert == 4);
assert(cells[44].ngbs[0] == 34);
assert(cells[44].ngbs[1] == 45);
assert(cells[44].ngbs[2] == 54);
assert(cells[44].ngbs[3] == 43);
}
message("Done.");
}
......
# ID pos_x pos_y pos_z v_x v_y v_z rho rho_dh wcount wcount_dh div_v curl_vx curl_vy curl_vz
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 3e-6 1e-4 5e-4 1.4e-2 1.1e-5 3e-6 3e-6 8e-6
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1.5e-6 1.4e-2 1e-5 2e-3 2.5e-4 3e-3 3e-3 3e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1.5e-6 1.7e-2 1e-5 2e-3 2.5e-4 3e-3 3e-3 3e-3
0 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e-6 1e0 1e-6 4e-6 4e-6 4e-6
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment