/*!
* \brief Write info and content of the N_array_2d struct to stdout
*
* Offsets are ignored
*
* \param data N_array_2d *
* \return void
* */
void N_print_array_2d(N_array_2d * data)
{
int i, j;
N_print_array_2d_info(data);
for (j = 0 - data->offset; j < data->rows + data->offset; j++) {
for (i = 0 - data->offset; i < data->cols + data->offset; i++) {
if (data->type == CELL_TYPE)
fprintf(stdout, "%6d ", N_get_array_2d_c_value(data, i, j));
else if (data->type == FCELL_TYPE)
fprintf(stdout, "%6.6f ", N_get_array_2d_f_value(data, i, j));
else if (data->type == DCELL_TYPE)
printf("%6.6f ", N_get_array_2d_d_value(data, i, j));
}
fprintf(stdout, "\n");
}
fprintf(stdout, "\n");
return;
}
/* *************************************************************** */
int test_array_2d(void)
{
int sum = 0, res = 0;
struct Cell_head region;
N_array_2d *data1;
N_array_2d *data11;
N_array_2d *data2;
N_array_2d *data22;
N_array_2d *data3;
N_array_2d *data33;
char buff[1024];
double min, max, ssum;
int nonzero;
N_array_2d *tmp;
/*Alloacte memory for all arrays */
data1 = N_alloc_array_2d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, 1, CELL_TYPE);
N_print_array_2d_info(data1);
data11 = N_alloc_array_2d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, 1, CELL_TYPE);
data2 = N_alloc_array_2d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, 1, FCELL_TYPE);
N_print_array_2d_info(data2);
data22 = N_alloc_array_2d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, 1, FCELL_TYPE);
data3 = N_alloc_array_2d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, 1, DCELL_TYPE);
N_print_array_2d_info(data3);
data33 = N_alloc_array_2d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, 1, DCELL_TYPE);
/*Fill the first arrays with data */
res = fill_array_2d(data1);
if (res != 0)
G_warning("test_array_2d: error while filling array with values");
sum += res;
res = fill_array_2d(data2);
if (res != 0)
G_warning("test_array_2d: error while filling array with values");
sum += res;
res = fill_array_2d(data3);
if (res != 0)
G_warning("test_array_2d: error while filling array with values");
sum += res;
/*Copy the data */
N_copy_array_2d(data1, data11);
N_copy_array_2d(data2, data22);
N_copy_array_2d(data3, data33);
/*Compare the data */
res = compare_array_2d(data1, data11);
if (res != 0)
G_warning("test_array_2d: error in N_copy_array_2d");
sum += res;
res = compare_array_2d(data2, data22);
if (res != 0)
G_warning("test_array_2d: error in N_copy_array_2d");
sum += res;
res = compare_array_2d(data3, data33);
if (res != 0)
G_warning("test_array_2d: error in N_copy_array_2d");
sum += res;
/*compute statistics */
N_calc_array_2d_stats(data1, &min, &max, &ssum, &nonzero, 0);
G_message("CELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 81 || ssum != 2025 || nonzero != 100) {
G_warning("test_array_2d: error in N_calc_array_2d_stats");
sum++;
}
N_calc_array_2d_stats(data1, &min, &max, &ssum, &nonzero, 1);
G_message("CELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 81 || ssum != 2025 || nonzero != 144) {
G_warning("test_array_2d: error in N_calc_array_2d_stats");
sum++;
}
N_calc_array_2d_stats(data2, &min, &max, &ssum, &nonzero, 0);
G_message("FCELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 81 || ssum != 2025 || nonzero != 100) {
G_warning("test_array_2d: error in N_calc_array_2d_stats");
sum++;
}
N_calc_array_2d_stats(data2, &min, &max, &ssum, &nonzero, 1);
G_message("FCELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 81 || ssum != 2025 || nonzero != 144) {
G_warning("test_array_2d: error in N_calc_array_2d_stats");
sum++;
}
N_calc_array_2d_stats(data3, &min, &max, &ssum, &nonzero, 0);
G_message("DCELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 81 || ssum != 2025 || nonzero != 100) {
G_warning("test_array_2d: error in N_calc_array_2d_stats");
sum++;
}
N_calc_array_2d_stats(data3, &min, &max, &ssum, &nonzero, 1);
G_message("DCELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 81 || ssum != 2025 || nonzero != 144) {
G_warning("test_array_2d: error in N_calc_array_2d_stats");
sum++;
}
/*test the array math functions */
tmp = N_math_array_2d(data1, data2, NULL, N_ARRAY_SUM);
N_math_array_2d(data2, data2, tmp, N_ARRAY_SUM);
res = N_convert_array_2d_null_to_zero(tmp);
if (res != 0)
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum = res;
N_free_array_2d(tmp);
tmp = N_math_array_2d(data2, data3, NULL, N_ARRAY_DIF);
N_math_array_2d(data1, data2, tmp, N_ARRAY_DIF);
res = N_convert_array_2d_null_to_zero(tmp);
if (res != 0)
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum = res;
N_free_array_2d(tmp);
tmp = N_math_array_2d(data1, data1, NULL, N_ARRAY_MUL);
N_math_array_2d(data1, data1, tmp, N_ARRAY_MUL);
res = N_convert_array_2d_null_to_zero(tmp);
if (res != 0)
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum = res;
N_free_array_2d(tmp);
tmp = N_math_array_2d(data2, data3, NULL, N_ARRAY_DIV);
N_math_array_2d(data1, data2, tmp, N_ARRAY_DIV);
res = N_convert_array_2d_null_to_zero(tmp);
if (res == 0) { /* if a division with zero is detected, the value is set to null, not to nan */
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum++;
}
N_free_array_2d(tmp);
/*check for correct norm calculation */
if (N_norm_array_2d(data1, data11, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
if (N_norm_array_2d(data1, data11, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
if (N_norm_array_2d(data2, data3, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
if (N_norm_array_2d(data2, data3, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
/*fill arrays with null values */
res = fill_array_2d_null(data1);
if (res != 0)
G_warning
("test_array_2d: error while filling array with cell null values");
sum += res;
res = fill_array_2d_null(data2);
if (res != 0)
G_warning
("test_array_2d: error while filling array with fcell null values");
sum += res;
res = fill_array_2d_null(data3);
if (res != 0)
G_warning
("test_array_2d: error while filling array with dcell null values");
sum += res;
/*Copy the data */
N_copy_array_2d(data1, data11);
N_copy_array_2d(data2, data22);
N_copy_array_2d(data3, data33);
/*Compare the data */
compare_array_2d(data1, data11);
compare_array_2d(data2, data22);
compare_array_2d(data3, data33);
/*check for correct norm calculation in case of null values */
if (N_norm_array_2d(data1, data11, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
if (N_norm_array_2d(data1, data11, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
if (N_norm_array_2d(data2, data3, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
if (N_norm_array_2d(data2, data3, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_2d: error in N_norm_array_2d");
sum++;
}
/*test the array math functions with null values */
tmp = N_math_array_2d(data1, data11, NULL, N_ARRAY_SUM);
N_math_array_2d(data2, data22, tmp, N_ARRAY_SUM);
res = N_convert_array_2d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero ");
sum++;
}
N_free_array_2d(tmp);
tmp = N_math_array_2d(data2, data22, NULL, N_ARRAY_DIF);
N_math_array_2d(data3, data33, tmp, N_ARRAY_DIF);
res = N_convert_array_2d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum++;
}
N_free_array_2d(tmp);
tmp = N_math_array_2d(data1, data11, NULL, N_ARRAY_MUL);
N_math_array_2d(data3, data33, tmp, N_ARRAY_MUL);
res = N_convert_array_2d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum++;
}
N_free_array_2d(tmp);
tmp = N_math_array_2d(data2, data3, NULL, N_ARRAY_DIV);
N_math_array_2d(data1, data11, tmp, N_ARRAY_DIV);
res = N_convert_array_2d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_2d: error in N_convert_array_2d_null_to_zero");
sum++;
}
N_free_array_2d(tmp);
N_free_array_2d(data1);
N_free_array_2d(data2);
N_free_array_2d(data3);
G_get_set_window(®ion);
data1 = N_alloc_array_2d(region.cols, region.rows, 0, CELL_TYPE);
data2 = N_alloc_array_2d(region.cols, region.rows, 0, FCELL_TYPE);
data3 = N_alloc_array_2d(region.cols, region.rows, 0, DCELL_TYPE);
fill_array_2d(data1);
fill_array_2d(data2);
fill_array_2d(data3);
/*raster IO methods */
N_write_array_2d_to_rast(data1, "gpde_lib_test_raster_1");
N_write_array_2d_to_rast(data2, "gpde_lib_test_raster_2");
N_write_array_2d_to_rast(data2, "gpde_lib_test_raster_3");
tmp = N_read_rast_to_array_2d("gpde_lib_test_raster_1", NULL);
N_read_rast_to_array_2d("gpde_lib_test_raster_1", tmp);
N_free_array_2d(tmp);
tmp = N_read_rast_to_array_2d("gpde_lib_test_raster_2", NULL);
N_read_rast_to_array_2d("gpde_lib_test_raster_2", tmp);
N_free_array_2d(tmp);
tmp = N_read_rast_to_array_2d("gpde_lib_test_raster_3", NULL);
N_read_rast_to_array_2d("gpde_lib_test_raster_3", tmp);
N_free_array_2d(tmp);
sprintf(buff,
"g.remove rast=gpde_lib_test_raster_1,gpde_lib_test_raster_2,gpde_lib_test_raster_3");
system(buff);
N_free_array_2d(data1);
N_free_array_2d(data11);
N_free_array_2d(data2);
N_free_array_2d(data22);
N_free_array_2d(data3);
N_free_array_2d(data33);
return sum;
}