/*!
* \brief Write info and content of the array data to stdout
*
* Offsets are ignored
*
* \param data N_array_2d *
* \return void
* */
void N_print_array_3d(N_array_3d * data)
{
int i, j, k;
N_print_array_3d_info(data);
for (k = 0; k < data->depths; k++) {
for (j = 0; j < data->rows; j++) {
for (i = 0; i < data->cols; i++) {
if (data->type == FCELL_TYPE)
printf("%6.6f ", N_get_array_3d_f_value(data, i, j, k));
else if (data->type == DCELL_TYPE)
printf("%6.6f ", N_get_array_3d_d_value(data, i, j, k));
}
printf("\n");
}
printf("\n");
}
printf("\n");
return;
}
/* *************************************************************** */
int test_array_3d(void)
{
int sum = 0, res = 0;
char buff[1024];
RASTER3D_Region region;
N_array_3d *data1;
N_array_3d *data11;
N_array_3d *data2;
N_array_3d *data22;
N_array_3d *tmp;
double min, max, ssum;
int nonzero;
/*Alloacte memory for all arrays */
data1 =
N_alloc_array_3d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, TEST_N_NUM_DEPTHS, 2,
FCELL_TYPE);
N_print_array_3d_info(data1);
data11 =
N_alloc_array_3d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, TEST_N_NUM_DEPTHS, 2,
FCELL_TYPE);
data2 =
N_alloc_array_3d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, TEST_N_NUM_DEPTHS, 2,
DCELL_TYPE);
N_print_array_3d_info(data2);
data22 =
N_alloc_array_3d(TEST_N_NUM_COLS, TEST_N_NUM_ROWS, TEST_N_NUM_DEPTHS, 2,
DCELL_TYPE);
/*Fill the first arrays with data */
res = fill_array_3d(data1);
if (res != 0)
G_warning("test_array_3d: error while filling array with values");
sum += res;
res = fill_array_3d(data2);
if (res != 0)
G_warning("test_array_3d: error while filling array with values");
sum += res;
/*Copy the data */
N_copy_array_3d(data1, data11);
N_copy_array_3d(data2, data22);
/*Compare the data */
res = compare_array_3d(data1, data11);
if (res != 0)
G_warning("test_array_3d: error in N_copy_array_2d");
sum += res;
res = compare_array_3d(data1, data11);
if (res != 0)
G_warning("test_array_3d: error in N_copy_array_2d");
sum += res;
/*compute statistics */
N_calc_array_3d_stats(data1, &min, &max, &ssum, &nonzero, 0);
G_message("FELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 729 || ssum != 91125 || nonzero != 1000) {
G_warning("test_array_3d: error in N_calc_array_3d_stats");
sum++;
}
N_calc_array_3d_stats(data1, &min, &max, &ssum, &nonzero, 1);
G_message("FELL Min %g Max %g Sum %g nonzero %i\n", min, max, ssum,
nonzero);
if (min != 0 || max != 729 || ssum != 91125 || nonzero != 2744) {
G_warning("test_array_3d: error in N_calc_array_3d_stats");
sum++;
}
N_calc_array_3d_stats(data2, &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 != 729 || ssum != 91125 || nonzero != 1000) {
G_warning("test_array_3d: error in N_calc_array_3d_stats");
sum++;
}
N_calc_array_3d_stats(data2, &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 != 729 || ssum != 91125 || nonzero != 2744) {
G_warning("test_array_3d: error in N_calc_array_3d_stats");
sum++;
}
/*test the array math functions */
tmp = N_math_array_3d(data1, data2, NULL, N_ARRAY_SUM);
N_math_array_3d(data2, data2, tmp, N_ARRAY_SUM);
res = N_convert_array_3d_null_to_zero(tmp);
if (res != 0)
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum = res;
N_free_array_3d(tmp);
tmp = N_math_array_3d(data2, data1, NULL, N_ARRAY_DIF);
N_math_array_3d(data1, data2, tmp, N_ARRAY_DIF);
res = N_convert_array_3d_null_to_zero(tmp);
if (res != 0)
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum = res;
N_free_array_3d(tmp);
tmp = N_math_array_3d(data1, data1, NULL, N_ARRAY_MUL);
N_math_array_3d(data1, data1, tmp, N_ARRAY_MUL);
res = N_convert_array_3d_null_to_zero(tmp);
if (res != 0)
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum = res;
N_free_array_3d(tmp);
tmp = N_math_array_3d(data2, data1, NULL, N_ARRAY_DIV);
N_math_array_3d(data1, data2, tmp, N_ARRAY_DIV);
res = N_convert_array_3d_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_3d: error in N_convert_array_3d_null_to_zero");
sum++;
}
N_free_array_3d(tmp);
/*check for correct norm calculation */
if (N_norm_array_3d(data1, data11, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
if (N_norm_array_3d(data1, data11, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
if (N_norm_array_3d(data1, data2, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
if (N_norm_array_3d(data1, data2, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
/*fill arrays with null values */
res = fill_array_3d_null(data1);
if (res != 0)
G_warning
("test_array_3d: error while filling array with float null values");
sum += res;
res = fill_array_3d_null(data2);
if (res != 0)
G_warning
("test_array_3d: error while filling array with double null values");
sum += res;
/*Copy the data */
N_copy_array_3d(data1, data11);
N_copy_array_3d(data2, data22);
/*Compare the data */
compare_array_3d(data1, data11);
compare_array_3d(data2, data22);
/*test the array math functions */
tmp = N_math_array_3d(data1, data2, NULL, N_ARRAY_SUM);
N_math_array_3d(data2, data2, tmp, N_ARRAY_SUM);
res = N_convert_array_3d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum++;
}
N_free_array_3d(tmp);
tmp = N_math_array_3d(data2, data1, NULL, N_ARRAY_DIF);
N_math_array_3d(data1, data2, tmp, N_ARRAY_DIF);
res = N_convert_array_3d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum++;
}
N_free_array_3d(tmp);
tmp = N_math_array_3d(data1, data1, NULL, N_ARRAY_MUL);
N_math_array_3d(data1, data1, tmp, N_ARRAY_MUL);
res = N_convert_array_3d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum++;
}
N_free_array_3d(tmp);
tmp = N_math_array_3d(data2, data1, NULL, N_ARRAY_DIV);
N_math_array_3d(data1, data2, tmp, N_ARRAY_DIV);
res = N_convert_array_3d_null_to_zero(tmp);
if (res == 0) {
G_warning("test_array_3d: error in N_convert_array_3d_null_to_zero");
sum++;
}
N_free_array_3d(tmp);
/*check for correct norm calculation in case of null values */
if (N_norm_array_3d(data1, data11, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
if (N_norm_array_3d(data1, data11, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
if (N_norm_array_3d(data1, data2, N_EUKLID_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
if (N_norm_array_3d(data1, data2, N_MAXIMUM_NORM) != 0.0) {
G_warning("test_array_3d: error in N_norm_array_3d");
sum++;
}
N_free_array_3d(data1);
N_free_array_3d(data2);
/*Set the defaults */
Rast3d_init_defaults();
Rast3d_get_window(®ion);
data1 =
N_alloc_array_3d(region.cols, region.rows, region.depths, 0,
FCELL_TYPE);
data2 =
N_alloc_array_3d(region.cols, region.rows, region.depths, 0,
DCELL_TYPE);
fill_array_3d(data1);
fill_array_3d(data2);
/*Volume IO methods */
N_write_array_3d_to_rast3d(data1, "gpde_lib_test_volume_1", 1);
N_write_array_3d_to_rast3d(data2, "gpde_lib_test_volume_2", 1);
tmp = N_read_rast3d_to_array_3d("gpde_lib_test_volume_1", NULL, 1);
N_read_rast3d_to_array_3d("gpde_lib_test_volume_1", tmp, 1);
N_free_array_3d(tmp);
tmp = N_read_rast3d_to_array_3d("gpde_lib_test_volume_2", NULL, 1);
N_read_rast3d_to_array_3d("gpde_lib_test_volume_2", tmp, 1);
N_free_array_3d(tmp);
sprintf(buff,
"g.remove rast3d=gpde_lib_test_volume_1,gpde_lib_test_volume_2");
system(buff);
N_free_array_3d(data1);
N_free_array_3d(data11);
N_free_array_3d(data2);
N_free_array_3d(data22);
return sum;
}
