30#include <grass/vector.h>
31#include <grass/dbmi.h>
34#include <grass/interpf.h>
57 double dnorm,
struct triple skip_point)
61 double east = data->
xmax;
62 double west = data->
x_orig;
63 double north = data->
ymax;
64 double south = data->
y_orig;
65 double h, xx, yy, r2, hz, zz,
err, xmm, ymm,
r;
81 for (mm = 1; mm <= n_points; mm++) {
83 for (m = 1; m <= n_points; m++) {
84 xx = points[mm - 1].
x - points[m - 1].
x;
85 yy = points[mm - 1].
y - points[m - 1].
y;
86 r2 = yy * yy + xx * xx;
99 zz = points[mm - 1].
z + zmin;
101 xmm = points[mm - 1].
x * dnorm + params->
x_orig + west;
102 ymm = points[mm - 1].
y * dnorm + params->
y_orig + south;
103 if ((xmm >= west + params->
x_orig) && (xmm <= east + params->x_orig) &&
104 (ymm >= south + params->
y_orig) && (ymm <= north + params->y_orig))
112 Vect_reset_line(
Pnts);
113 Vect_reset_cats(
Cats2);
115 Vect_append_point(
Pnts, xmm, ymm, zz);
117 Vect_cat_set(
Cats2, 1, cat);
120 db_zero_string(&
sql2);
121 sprintf(buf,
"insert into %s values ( %d ",
ff->table, cat);
122 db_append_string(&
sql2, buf);
124 sprintf(buf,
", %f",
err);
125 db_append_string(&
sql2, buf);
126 db_append_string(&
sql2,
")");
127 G_debug(3,
"IL_check_at_points_2d: %s", db_get_string(&
sql2));
129 if (db_execute_immediate(
driver2, &
sql2) != DB_OK) {
133 db_get_string(&
sql2));
144 for (m = 1; m <= n_points - 1; m++) {
145 xx = points[m - 1].
x - skip_point.
x;
146 yy = points[m - 1].
y - skip_point.
y;
147 r2 = yy * yy + xx * xx;
155 zz = skip_point.
z + zmin;
157 xmm = skip_point.
x * dnorm + params->
x_orig + west;
158 ymm = skip_point.
y * dnorm + params->
y_orig + south;
160 if ((xmm >= west + params->
x_orig) && (xmm <= east + params->x_orig) &&
161 (ymm >= south + params->
y_orig) && (ymm <= north + params->y_orig))
167 Vect_reset_line(
Pnts);
168 Vect_reset_cats(
Cats2);
170 Vect_append_point(
Pnts, xmm, ymm, zz);
172 Vect_cat_set(
Cats2, 1, cat);
175 db_zero_string(&
sql2);
176 sprintf(buf,
"insert into %s values ( %d ",
ff->table, cat);
177 db_append_string(&
sql2, buf);
179 sprintf(buf,
", %f", skip_err);
180 db_append_string(&
sql2, buf);
181 db_append_string(&
sql2,
")");
182 G_debug(3,
"IL_check_at_points_2d: %s", db_get_string(&
sql2));
184 if (db_execute_immediate(
driver2, &
sql2) != DB_OK) {
188 db_get_string(&
sql2));
int G_debug(int level, const char *msg,...)
Print debugging message.
void G_fatal_error(const char *msg,...)
Print a fatal error message to stderr.
int IL_check_at_points_2d(struct interp_params *params, struct quaddata *data, double *b, double *ertot, double zmin, double dnorm, struct triple skip_point)
SYMBOL * err(FILE *fp, SYMBOL *s, char *msg)