GRASS GIS 7 Programmer's Manual  7.9.dev(2021)-e5379bbd7
xgt.c
Go to the documentation of this file.
1 
2 #include <grass/gis.h>
3 #include <grass/raster.h>
4 #include <grass/calc.h>
5 
6 /****************************************************************
7 gt(a,b) = a > b
8 ****************************************************************/
9 
10 int f_gt(int argc, const int *argt, void **args)
11 {
12  CELL *res = args[0];
13  int i;
14 
15  if (argc < 2)
16  return E_ARG_LO;
17  if (argc > 2)
18  return E_ARG_HI;
19 
20  switch (argt[1]) {
21  case CELL_TYPE:
22  {
23  CELL *arg1 = args[1];
24  CELL *arg2 = args[2];
25 
26  for (i = 0; i < columns; i++) {
27  if (IS_NULL_C(&arg1[i]) || IS_NULL_C(&arg2[i]))
28  SET_NULL_C(&res[i]);
29  else
30  res[i] = arg1[i] > arg2[i];
31  }
32  return 0;
33  }
34  case FCELL_TYPE:
35  {
36  FCELL *arg1 = args[1];
37  FCELL *arg2 = args[2];
38 
39  for (i = 0; i < columns; i++) {
40  if (IS_NULL_F(&arg1[i]) || IS_NULL_F(&arg2[i]))
41  SET_NULL_C(&res[i]);
42  else
43  res[i] = arg1[i] > arg2[i];
44  }
45  return 0;
46  }
47  case DCELL_TYPE:
48  {
49  DCELL *arg1 = args[1];
50  DCELL *arg2 = args[2];
51 
52  for (i = 0; i < columns; i++) {
53  if (IS_NULL_D(&arg1[i]) || IS_NULL_D(&arg2[i]))
54  SET_NULL_C(&res[i]);
55  else
56  res[i] = arg1[i] > arg2[i];
57  }
58  return 0;
59  }
60  default:
61  return E_INV_TYPE;
62  }
63 }
#define CELL_TYPE
Definition: raster.h:11
#define SET_NULL_C(x)
Definition: calc.h:32
double DCELL
Definition: gis.h:603
#define IS_NULL_F(x)
Definition: calc.h:29
#define IS_NULL_C(x)
Definition: calc.h:28
Definition: calc.h:13
int columns
Definition: calc.c:12
#define DCELL_TYPE
Definition: raster.h:13
int f_gt(int argc, const int *argt, void **args)
Definition: xgt.c:10
Definition: calc.h:12
float FCELL
Definition: gis.h:604
#define IS_NULL_D(x)
Definition: calc.h:30
int CELL
Definition: gis.h:602
#define FCELL_TYPE
Definition: raster.h:12