19 #include <grass/N_pde.h>
37 "N_calc_gradient_field_2d_stats: compute gradient field stats");
52 field->
sum = sumx + sumy;
53 field->
nonull = nonullx + nonully;
114 double dx, dy, p1, p2, r1, r2, mean, grad, res;
120 (
"N_compute_gradient_field_2d: the arrays are not of equal size");
124 (
"N_compute_gradient_field_2d: the arrays are not of equal size");
128 (
"N_compute_gradient_field_2d: array sizes and geometry data are different");
131 G_debug(3,
"N_compute_gradient_field_2d: compute gradient field");
144 (
"N_compute_gradient_field_2d: gradient field sizes and geometry data are different");
148 for (j = 0; j < rows; j++)
149 for (i = 0; i < cols - 1; i++) {
158 grad = (p1 - p2) / dx;
173 for (j = 0; j < rows - 1; j++)
174 for (i = 0; i <
cols; i++) {
183 grad = (p1 - p2) / dy;
192 res = -1 * mean * grad;
256 G_fatal_error(
"N_compute_gradient_components_2d: x array is empty");
258 G_fatal_error(
"N_compute_gradient_components_2d: y array is empty");
264 if (x->
cols != cols || x->
rows != rows)
266 (
"N_compute_gradient_components_2d: the size of the x array doesn't fit the gradient field size");
267 if (y->
cols != cols || y->
rows != rows)
269 (
"N_compute_gradient_components_2d: the size of the y array doesn't fit the gradient field size");
271 for (j = 0; j < rows; j++)
272 for (i = 0; i <
cols; i++) {
276 if (grad.
WC == 0.0 || grad.
EC == 0.0)
277 vx = (grad.
WC + grad.
EC);
279 vx = (grad.
WC + grad.
EC) / 2;
280 if (grad.
NC == 0.0 || grad.
SC == 0.0)
281 vy = (grad.
NC + grad.
SC);
283 vy = (grad.
NC + grad.
SC) / 2;
302 double minx, miny, minz;
303 double maxx, maxy, maxz;
304 double sumx, sumy, sumz;
305 int nonullx, nonully, nonullz;
308 "N_calc_gradient_field_3d_stats: compute gradient field stats");
314 if (minx <= minz && minx <= miny)
316 if (miny <= minz && miny <= minx)
318 if (minz <= minx && minz <= miny)
321 if (maxx >= maxz && maxx >= maxy)
323 if (maxy >= maxz && maxy >= maxx)
325 if (maxz >= maxx && maxz >= maxy)
328 field->
sum = sumx + sumy + sumz;
329 field->
nonull = nonullx + nonully + nonullz;
395 int cols, rows, depths;
396 double dx, dy, dz, p1, p2, r1, r2, mean, grad, res;
403 (
"N_compute_gradient_field_3d: the arrays are not of equal size");
408 (
"N_compute_gradient_field_3d: the arrays are not of equal size");
413 (
"N_compute_gradient_field_3d: the arrays are not of equal size");
418 (
"N_compute_gradient_field_3d: array sizes and geometry data are different");
420 G_debug(3,
"N_compute_gradient_field_3d: compute gradient field");
429 if (gradfield ==
NULL) {
436 (
"N_compute_gradient_field_3d: gradient field sizes and geometry data are different");
439 for (k = 0; k < depths; k++)
440 for (j = 0; j < rows; j++)
441 for (i = 0; i < cols - 1; i++) {
450 grad = (p1 - p2) / dx;
462 "N_compute_gradient_field_3d: X-direction insert value %6.5g at %i %i %i ",
469 for (k = 0; k < depths; k++)
470 for (j = 0; j < rows - 1; j++)
471 for (i = 0; i <
cols; i++) {
480 grad = (p1 - p2) / dy;
489 res = -1 * mean * grad;
493 "N_compute_gradient_field_3d: Y-direction insert value %6.5g at %i %i %i ",
500 for (k = 0; k < depths - 1; k++)
501 for (j = 0; j < rows; j++)
502 for (i = 0; i <
cols; i++) {
511 grad = (p1 - p2) / dz;
523 "N_compute_gradient_field_3d: Z-direction insert value %6.5g at %i %i %i ",
585 int rows,
cols, depths;
594 G_fatal_error(
"N_compute_gradient_components_3d: x array is empty");
596 G_fatal_error(
"N_compute_gradient_components_3d: y array is empty");
598 G_fatal_error(
"N_compute_gradient_components_3d: z array is empty");
607 (
"N_compute_gradient_components_3d: the size of the x array doesn't fit the gradient field size");
610 (
"N_compute_gradient_components_3d: the size of the y array doesn't fit the gradient field size");
613 (
"N_compute_gradient_components_3d: the size of the z array doesn't fit the gradient field size");
615 for (k = 0; k < depths; k++)
616 for (j = 0; j < rows; j++)
617 for (i = 0; i <
cols; i++) {
620 if (grad.
WC == 0.0 || grad.
EC == 0.0)
621 vx = (grad.
WC + grad.
EC);
623 vx = (grad.
WC + grad.
EC) / 2;
624 if (grad.
NC == 0.0 || grad.
SC == 0.0)
625 vy = (grad.
NC + grad.
SC);
627 vy = (grad.
NC + grad.
SC) / 2;
628 if (grad.
TC == 0.0 || grad.
BC == 0.0)
629 vz = (grad.
TC + grad.
BC);
631 vz = (grad.
TC + grad.
BC) / 2;
int N_is_array_3d_value_null(N_array_3d *data, int col, int row, int depth)
This function returns 1 if value of N_array_3d data at position col, row, depth is of type null...
N_gradient_field_3d * N_compute_gradient_field_3d(N_array_3d *pot, N_array_3d *weight_x, N_array_3d *weight_y, N_array_3d *weight_z, N_geom_data *geom, N_gradient_field_3d *gradfield)
This function computes the gradient based on the input N_array_3d pot (that means potential)...
void N_compute_gradient_field_components_2d(N_gradient_field_2d *field, N_array_2d *x_comp, N_array_2d *y_comp)
Calculate the x and y vector components from a gradient field for each cell and stores them in the pr...
Gradient between the cells in X and Y direction.
void N_put_array_2d_d_value(N_array_2d *data, int col, int row, DCELL value)
Writes a DCELL value to the N_array_2d struct at position col, row.
N_gradient_field_2d * N_compute_gradient_field_2d(N_array_2d *pot, N_array_2d *weight_x, N_array_2d *weight_y, N_geom_data *geom, N_gradient_field_2d *gradfield)
This function computes the gradient based on the input N_array_2d pot (potential), a weighting factor N_array_2d named weight and the distance between two cells saved in the N_geom_data struct.
void N_calc_gradient_field_3d_stats(N_gradient_field_3d *field)
Calculate basic statistics of a gradient field.
void N_calc_array_2d_stats(N_array_2d *a, double *min, double *max, double *sum, int *nonull, int withoffset)
Calculate basic statistics of the N_array_2d struct.
N_gradient_3d * N_get_gradient_3d(N_gradient_field_3d *field, N_gradient_3d *gradient, int col, int row, int depth)
Return a N_gradient_3d structure calculated from the input gradient field at position [depth][row][co...
Geometric information about the structured grid.
N_gradient_field_2d * N_alloc_gradient_field_2d(int cols, int rows)
Allocate a N_gradient_field_2d.
void N_compute_gradient_field_components_3d(N_gradient_field_3d *field, N_array_3d *x_comp, N_array_3d *y_comp, N_array_3d *z_comp)
Calculate the x, y and z vector components from a gradient field for each cell and store them in the ...
Gradient between the cells in X, Y and Z direction.
void N_calc_array_3d_stats(N_array_3d *a, double *min, double *max, double *sum, int *nonull, int withoffset)
Calculate basic statistics of the N_array_3d struct.
int N_is_array_2d_value_null(N_array_2d *data, int col, int row)
Returns 1 if the value of N_array_2d struct at postion col, row is of type null, otherwise 0...
double N_get_array_3d_d_value(N_array_3d *data, int col, int row, int depth)
This function returns the value of type float at position col, row, depth.
N_gradient_2d * N_get_gradient_2d(N_gradient_field_2d *field, N_gradient_2d *gradient, int col, int row)
Return a N_gradient_2d structure calculated from the input gradient field at position [row][col]...
int G_debug(int level, const char *msg,...)
Print debugging message.
N_gradient_field_3d * N_alloc_gradient_field_3d(int cols, int rows, int depths)
Allocate a N_gradient_field_3d.
DCELL N_get_array_2d_d_value(N_array_2d *data, int col, int row)
Returns the value of type DCELL at position col, row.
void N_put_array_3d_d_value(N_array_3d *data, int col, int row, int depth, double value)
Writes a double value to the N_array_3d struct at position col, row, depth.
double N_calc_harmonic_mean(double a, double b)
Calculate the harmonical mean of values a and b.
int G_fatal_error(const char *msg,...)
Print a fatal error message to stderr.
void N_calc_gradient_field_2d_stats(N_gradient_field_2d *field)
Calculate basic statistics of a gradient field.