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