core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

ECC-1433: Remove simplified test from ECC-1431

This commit is contained in:
Eugen Betke 2022-08-03 10:12:12 +00:00
parent c89e55ae3f
commit ce58cd7396
1 changed files with 5 additions and 93 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

98
tests/grib_ecc-1431.cpp
View File

@ -32,7 +32,7 @@ struct Range
// //
// Test: This test doesn't accept errors after decompression. // Test: This test doesn't accept errors after decompression.
// //
void test_full() int main(int argc, char** argv)
{ {
std::vector<Range> ranges; std::vector<Range> ranges;
ranges.push_back({ flt::max() / 10, flt::max() }); ranges.push_back({ flt::max() / 10, flt::max() });
@ -65,11 +65,11 @@ void test_full()
for (const Range range : ranges) { for (const Range range : ranges) {
Assert(range.min < range.max); Assert(range.min < range.max);
std::cout << "Testing: " std::cout << "Testing: "
<< "n_vals : " << std::setw(7) << in_values_len << "n_vals : " << std::setw(7) << in_values_len
<< " range(" << std::scientific << std::setw(12) << range.min << ", " << std::setw(12) << range.max << ") " << " range(" << std::scientific << std::setw(12) << range.min << ", " << std::setw(12) << range.max << ") "
<< std::endl; << std::endl;
std::uniform_real_distribution<double> std::uniform_real_distribution<double>
unif(range.min, range.max); unif(range.min, range.max);
codes_handle* handle = codes_grib_handle_new_from_samples(0, "reduced_gg_pl_128_grib2"); codes_handle* handle = codes_grib_handle_new_from_samples(0, "reduced_gg_pl_128_grib2");
double* in_values = new double[in_values_len]; double* in_values = new double[in_values_len];
@ -117,93 +117,5 @@ void test_full()
delete[] grid_ccsds_values; delete[] grid_ccsds_values;
} }
} }
}
//
// Test: This test accepts tolerates errors. Values after decompression must be within a certain user defined range.
//
void test_simplified()
{
std::vector<Range> ranges;
ranges.push_back({ flt::max() / 10, flt::max() });
ranges.push_back({ flt::max() / 100, flt::max() / 10 });
ranges.push_back({ 1e+9, 1e+10 });
ranges.push_back({ 1e+0, 1e+1 });
ranges.push_back({ 1e-1, 1e-0 });
ranges.push_back({ 1e-2, 1e-1 });
ranges.push_back({ 1e-11, 1e-10 });
ranges.push_back({ flt::min() * 10, flt::min() * 100 });
ranges.push_back({ flt::min(), flt::min() * 10 });
ranges.push_back({ flt::min(), flt::max() });
std::vector<size_t> values_lens = {
1,
10,
100,
100000,
10000000,
};
int err;
std::default_random_engine re;
size_t grid_ccsds_values_len = 0;
std::string packing_type = "";
size_t size = 0;
for (const size_t in_values_len : values_lens) {
for (const Range range : ranges) {
Assert(range.min < range.max);
std::cout << "Testing: " << in_values_len << " " << range.min << " " << range.max << std::endl;
std::uniform_real_distribution<double> unif(range.min, range.max);
codes_handle* handle = codes_grib_handle_new_from_samples(0, "reduced_gg_pl_128_grib2");
double* in_values = new double[in_values_len];
double* grid_ccsds_values = new double[in_values_len];
grid_ccsds_values_len = in_values_len;
// Initialize with random values
for (size_t i = 0; i < in_values_len; ++i) {
in_values[i] = unif(re);
}
// Test grid_ccsds
packing_type = "grid_ccsds";
size = packing_type.size();
if ((err = codes_set_double_array(handle, "values", in_values, in_values_len)) != 0) {
Assert(!"CCSDS encoding failed");
}
CODES_CHECK(codes_set_string(handle, "packingType", packing_type.c_str(), &size), 0);
if ((err = codes_get_double_array(handle, "values", grid_ccsds_values, &grid_ccsds_values_len)) != 0) {
Assert(!"CCSDS decoding failed");
}
Assert(grid_ccsds_values_len == in_values_len);
// Test buffers
double tolerance = 0.000001;
for (size_t i = 0; i < grid_ccsds_values_len; ++i) {
if (!((grid_ccsds_values[i] < (in_values[i] * (1 + tolerance))) &&
grid_ccsds_values[i] > (in_values[i] / (1 + tolerance)))) {
std::cout.precision(dbl::max_digits10);
std::cout << "Test failed: " << grid_ccsds_values[i] << " != " << in_values[i] << std::endl;
Assert(0);
}
}
codes_handle_delete(handle);
delete[] in_values;
delete[] grid_ccsds_values;
}
}
}
int main(int argc, char** argv)
{
//test_simplified();
test_full();
return 0; return 0;
} }