aboutsummaryrefslogtreecommitdiff
path: root/math/nlfit/nlvector.gx
blob: f33d063ce70e98bea75dacca9c940561f8235ca8 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
include	<math/nlfit.h>
$if (datatype == r)
include "nlfitdefr.h"
$else
include "nlfitdefd.h"
$endif

define	VAR	(($1 - 1) * $2 + 1)

# NLVECTOR -- Evaluate the fit for a series of data points.

procedure nlvector$t (nl, x, zfit, npts, nvars)

pointer	nl		# pointer to nl fitting structure
PIXEL	x[ARB]		# independent variables (npts * nvars)
PIXEL	zfit[ARB]	# function values (npts)
int	npts		# number of points
int	nvars		# number of independent variables

int	i

begin
	# Compute the fitted function.
	do i = 1, npts
	    call zcall5 (NL_FUNC(nl), x[VAR (i, nvars)], nvars,
		PARAM(NL_PARAM(nl)), NL_NPARAMS(nl), zfit[i])
end