GRASS GIS 8 Programmer's Manual  8.5.0dev(2024)-ab90c5e5a9
xsin.c
Go to the documentation of this file.
1 #include <math.h>
2 
3 #include <grass/gis.h>
4 #include <grass/raster.h>
5 #include <grass/calc.h>
6 
7 /**********************************************************************
8 sin(x)
9 
10  if floating point exception occurs during the evaluation of sin(x)
11  the result is NULL
12 
13  note: x is in degrees.
14 **********************************************************************/
15 
16 #define DEGREES_TO_RADIANS (M_PI / 180.0)
17 
18 int f_sin(int argc, const int *argt, void **args)
19 {
20  DCELL *res = args[0];
21  DCELL *arg1 = args[1];
22  int i;
23 
24  if (argc < 1)
25  return E_ARG_LO;
26  if (argc > 1)
27  return E_ARG_HI;
28 
29  if (argt[0] != DCELL_TYPE)
30  return E_RES_TYPE;
31 
32  if (argt[1] != DCELL_TYPE)
33  return E_ARG_TYPE;
34 
35  for (i = 0; i < columns; i++)
36  if (IS_NULL_D(&arg1[i]))
37  SET_NULL_D(&res[i]);
38  else {
40  res[i] = sin(arg1[i] * DEGREES_TO_RADIANS);
42  SET_NULL_D(&res[i]);
43  }
44 
45  return 0;
46 }
@ E_RES_TYPE
Definition: calc.h:14
@ E_ARG_TYPE
Definition: calc.h:13
@ E_ARG_HI
Definition: calc.h:12
@ E_ARG_LO
Definition: calc.h:11
volatile int floating_point_exception
Definition: calc.c:8
#define SET_NULL_D(x)
Definition: calc.h:32
int columns
Definition: calc.c:11
#define IS_NULL_D(x)
Definition: calc.h:28
double DCELL
Definition: gis.h:629
#define DCELL_TYPE
Definition: raster.h:13
#define DEGREES_TO_RADIANS
Definition: xsin.c:16
int f_sin(int argc, const int *argt, void **args)
Definition: xsin.c:18