#include <grass/N_pde.h>

Include dependency graph for N_gradient_calc.c:

Functions
void	N_calc_gradient_field_2d_stats (N_gradient_field_2d *field)
	Calculate basic statistics of a gradient field. More...

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. More...

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 provided N_array_2d structures. More...

void	N_calc_gradient_field_3d_stats (N_gradient_field_3d *field)
	Calculate basic statistics of a gradient field. More...

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), a weighting factor N_array_3d named weight and the distance between two cells saved in the N_geom_data struct. More...

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 provided N_array_3d structures. More...

Function Documentation

void N_calc_gradient_field_2d_stats ( N_gradient_field_2d * field )

Calculate basic statistics of a gradient field.

The statistic is stored in the gradient field struct

Parameters

field N_gradient_2d_field *

Returns: void

Definition at line 29 of file N_gradient_calc.c.

References G_debug(), N_gradient_field_2d::max, N_gradient_field_2d::mean, N_gradient_field_2d::min, N_calc_array_2d_stats(), N_gradient_field_2d::nonull, N_gradient_field_2d::sum, N_gradient_field_2d::x_array, and N_gradient_field_2d::y_array.

Referenced by N_compute_gradient_field_2d().

void N_calc_gradient_field_3d_stats ( N_gradient_field_3d * field )

Calculate basic statistics of a gradient field.

The statistic is stored in the gradient field struct

Parameters

field N_gradient_3d_field *

Returns: void

Definition at line 311 of file N_gradient_calc.c.

References G_debug(), N_gradient_field_3d::max, N_gradient_field_3d::mean, N_gradient_field_3d::min, N_calc_array_3d_stats(), N_gradient_field_3d::nonull, N_gradient_field_3d::sum, N_gradient_field_3d::x_array, N_gradient_field_3d::y_array, and N_gradient_field_3d::z_array.

Referenced by N_compute_gradient_field_3d().

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.

The gradient is calculated between cells for each cell and direction. An existing gradient field can be filled with new data or, if a NULL pointer is given, a new gradient field will be allocated with the appropriate size.

______________ 
|    |    |    |
|    |    |    |
|----|-NC-|----|
|    |    |    |
|   WC    EC   |
|    |    |    |
|----|-SC-|----|
|    |    |    |
|____|____|____|


x - direction:

r = 2 * weight[row][col]*weight[row][col + 1] / (weight[row][col]*weight[row][col + 1])
EC = r * (pot[row][col] - pot[row][col + 1])/dx

y - direction:

r = 2 * weight[row][col]*weight[row + 1][col] / (weight[row][col]*weight[row + 1][col])
SC = r * (pot[row][col] - pot[row + 1][col])/dy

the values SC and EC are the values of the next row/col

Parameters

pot	N_array_2d * - the potential N_array_2d
weight_x	N_array_2d * - the weighting factor N_array_2d used to modify the gradient in x-direction
weight_y	N_array_2d * - the weighting factor N_array_2d used to modify the gradient in y-direction
geom	N_geom_data * - geometry data structure
gradfield	N_gradient_field_2d * - a gradient field of the correct size, if a NULL pointer is provided this gradient field will be new allocated

Returns: N_gradient_field_2d * - the pointer to the computed gradient field

Definition at line 108 of file N_gradient_calc.c.

References N_geom_data::cols, N_array_2d::cols, N_gradient_field_2d::cols, dialogs::cols, N_geom_data::dx, N_geom_data::dy, G_debug(), G_fatal_error(), N_alloc_gradient_field_2d(), N_calc_gradient_field_2d_stats(), N_calc_harmonic_mean(), N_get_array_2d_d_value(), N_is_array_2d_value_null(), N_put_array_2d_d_value(), NULL, N_geom_data::rows, N_array_2d::rows, N_gradient_field_2d::rows, N_gradient_field_2d::x_array, and N_gradient_field_2d::y_array.

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), a weighting factor N_array_3d named weight and the distance between two cells saved in the N_geom_data struct.

The gradient is calculated between cells for each cell and direction. An existing gradient field can be filled with new data or, if a NULL pointer is given, a new gradient field will be allocated with the appropriate size.

|  /
TC NC
|/
--WC-----EC--
/|
SC BC
/  |

x - direction:

r = 2 * weight_x[depth][row][col]*weight_x[depth][row][col + 1] / (weight_X[depth][row][col]*weight_x[depth][row][col + 1])
EC = r * (pot[depth][row][col] - pot[depth][row][col + 1])/dx

y - direction:

r = 2 * weight_y[depth][row][col]*weight_y[depth][row + 1][col] / (weight_y[depth][row][col]*weight_y[depth][row + 1][col])
SC = r * (pot[depth][row][col] - pot[depth][row + 1][col])/dy

z - direction:

r = 2 * weight_z[depth][row][col]*weight_z[depth + 1][row][col] / (weight_z[depth][row][col]*weight_z[depth + 1][row][col])
TC = r * (pot[depth][row][col] - pot[depth + 1][row][col])/dy

the values BC, NC, WC are the values of the next depth/row/col

Parameters

pot	N_array_3d * - the potential N_array_2d
weight_x	N_array_3d * - the weighting factor N_array_3d used to modify the gradient in x-direction
weight_y	N_array_3d * - the weighting factor N_array_3d used to modify the gradient in y-direction
weight_z	N_array_3d * - the weighting factor N_array_3d used to modify the gradient in z-direction
geom	N_geom_data * - geometry data structure
gradfield	N_gradient_field_3d * - a gradient field of the correct size, if a NULL pointer is provided this gradient field will be new allocated

Returns: N_gradient_field_3d * - the pointer to the computed gradient field

Definition at line 400 of file N_gradient_calc.c.

References N_geom_data::cols, N_array_3d::cols, N_gradient_field_3d::cols, dialogs::cols, N_geom_data::depths, N_array_3d::depths, N_gradient_field_3d::depths, N_geom_data::dx, N_geom_data::dy, N_geom_data::dz, G_debug(), G_fatal_error(), N_alloc_gradient_field_3d(), N_calc_gradient_field_3d_stats(), N_calc_harmonic_mean(), N_get_array_3d_d_value(), N_is_array_3d_value_null(), N_put_array_3d_d_value(), NULL, N_geom_data::rows, N_array_3d::rows, N_gradient_field_3d::rows, N_gradient_field_3d::x_array, N_gradient_field_3d::y_array, and N_gradient_field_3d::z_array.

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 provided N_array_2d structures.

The arrays must have the same size as the gradient field.

Based on this storages scheme the gradient vector for each cell is 
calculated and stored in the provided  N_array_2d structures

______________ 
|    |    |    |
|    |    |    |
|----|-NC-|----|
|    |    |    |
|   WC    EC   |
|    |    |    |
|----|-SC-|----|
|    |    |    |
|____|____|____|

x vector component:

x = (WC + EC) / 2

y vector component:

y = (NC + SC) / 2

Parameters

field	N_gradient_field_2d *
x_comp	N_array_2d * - the array in which the x component will be written
y_comp	N_array_2d * - the array in which the y component will be written

Returns: void

Definition at line 249 of file N_gradient_calc.c.

References N_array_2d::cols, dialogs::cols, N_gradient_2d::EC, G_fatal_error(), N_get_gradient_2d(), N_put_array_2d_d_value(), N_gradient_2d::NC, N_array_2d::rows, N_gradient_2d::SC, N_gradient_2d::WC, N_gradient_field_2d::x_array, and y.

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 provided N_array_3d structures.

The arrays must have the same size as the gradient field.

Based on this storages scheme the gradient vector for each cell is 
calculated and stored in the provided  N_array_3d structures


|  /
TC NC
|/
--WC-----EC--
/|
SC BC
/  |


x vector component:

x = (WC + EC) / 2

y vector component:

y = (NC + SC) / 2

z vector component:

z = (TC + BC) / 2

Parameters

field	N_gradient_field_3d *
x_comp	N_array_3d * - the array in which the x component will be written
y_comp	N_array_3d * - the array in which the y component will be written
z_comp	N_array_3d * - the array in which the z component will be written

Returns: void

Definition at line 596 of file N_gradient_calc.c.

References N_gradient_3d::BC, N_array_3d::cols, dialogs::cols, N_array_3d::depths, N_gradient_3d::EC, G_fatal_error(), N_get_gradient_3d(), N_put_array_3d_d_value(), N_gradient_3d::NC, N_array_3d::rows, N_gradient_3d::SC, N_gradient_3d::TC, N_gradient_3d::WC, N_gradient_field_3d::x_array, and y.

Functions

Function Documentation