GRASS GIS 8 Programmer's Manual
8.3.2(2024)-exported
Loading...
Searching...
No Matches
xsqrt.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
sqrt(x)
9
10
if floating point exception occurs during the evaluation of sqrt(x)
11
the result is NULL
12
**********************************************************************/
13
14
int
f_sqrt
(
int
argc,
const
int
*argt,
void
**args)
15
{
16
DCELL *res = args[0];
17
DCELL *arg1 = args[1];
18
int
i;
19
20
if
(argc < 1)
21
return
E_ARG_LO;
22
if
(argc > 1)
23
return
E_ARG_HI;
24
25
if
(argt[0] != DCELL_TYPE)
26
return
E_RES_TYPE;
27
28
if
(argt[1] != DCELL_TYPE)
29
return
E_ARG_TYPE;
30
31
for
(i = 0; i <
columns
; i++)
32
if
(IS_NULL_D(&arg1[i]) || (arg1[i] < 0.0))
33
SET_NULL_D(&res[i]);
34
else
{
35
floating_point_exception
= 0;
36
res[i] = sqrt(arg1[i]);
37
if
(
floating_point_exception
)
38
SET_NULL_D(&res[i]);
39
}
40
41
return
0;
42
}
floating_point_exception
volatile int floating_point_exception
Definition
calc.c:8
columns
int columns
Definition
calc.c:11
f_sqrt
int f_sqrt(int argc, const int *argt, void **args)
Definition
xsqrt.c:14
calc
xsqrt.c
Generated on Fri May 24 2024 09:06:59 for GRASS GIS 8 Programmer's Manual by
1.9.8