aboutsummaryrefslogtreecommitdiff
path: root/math/nlfit/nlvector.gx
diff options
context:
space:
mode:
authorJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
committerJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
commitfa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch)
treebdda434976bc09c864f2e4fa6f16ba1952b1e555 /math/nlfit/nlvector.gx
downloadiraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz
Initial commit
Diffstat (limited to 'math/nlfit/nlvector.gx')
-rw-r--r--math/nlfit/nlvector.gx27
1 files changed, 27 insertions, 0 deletions
diff --git a/math/nlfit/nlvector.gx b/math/nlfit/nlvector.gx
new file mode 100644
index 00000000..f33d063c
--- /dev/null
+++ b/math/nlfit/nlvector.gx
@@ -0,0 +1,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