Mercurial > octave-nkf

annotate examples/mypow2.c @ 16088:b29b10fbb744 stable release-3-6-4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Version 3.6.4 released. * configure.ac (AC_INIT): Version is now 3.6.4. (OCTAVE_RELEASE_DATE): Now 2013-02-21.
author John W. Eaton <jwe@octave.org>
date Thu, 21 Feb 2013 15:17:54 -0500
parents 6cb30a539481
children be41c30bcb44
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6593
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
1 #include "mex.h"
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
2
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
3 void
7081
503001863427 [project @ 2007-10-31 01:08:14 by jwe]
jwe
parents: 7019
diff changeset
4 mexFunction (int nlhs, mxArray* plhs[], int nrhs,
9932
6cb30a539481 untabify files in examples directory
John W. Eaton <jwe@octave.org>
parents: 9053
diff changeset
5 const mxArray* prhs[])
6593
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
6 {
6686
2aad75fcc93a [project @ 2007-06-03 20:58:28 by dbateman]
dbateman
parents: 6593
diff changeset
7 mwIndex i;
2aad75fcc93a [project @ 2007-06-03 20:58:28 by dbateman]
dbateman
parents: 6593
diff changeset
8 mwSize n;
6593
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
9 double *vri, *vro;
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
10
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
11 if (nrhs != 1 || ! mxIsNumeric (prhs[0]))
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
12 mexErrMsgTxt ("expects matrix");
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
13
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
14 n = mxGetNumberOfElements (prhs[0]);
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
15 plhs[0] = (mxArray *) mxCreateNumericArray
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
16 (mxGetNumberOfDimensions (prhs[0]),
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
17 mxGetDimensions (prhs[0]), mxGetClassID (prhs[0]),
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
18 mxIsComplex (prhs[0]));
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
19 vri = mxGetPr (prhs[0]);
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
20 vro = mxGetPr (plhs[0]);
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
21
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
22 if (mxIsComplex (prhs[0]))
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
23 {
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
24 double *vii, *vio;
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
25 vii = mxGetPi (prhs[0]);
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
26 vio = mxGetPi (plhs[0]);
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
27
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
28 for (i = 0; i < n; i++)
9932
6cb30a539481 untabify files in examples directory
John W. Eaton <jwe@octave.org>
parents: 9053
diff changeset
29 {
6cb30a539481 untabify files in examples directory
John W. Eaton <jwe@octave.org>
parents: 9053
diff changeset
30 vro [i] = vri [i] * vri [i] - vii [i] * vii [i];
6cb30a539481 untabify files in examples directory
John W. Eaton <jwe@octave.org>
parents: 9053
diff changeset
31 vio [i] = 2 * vri [i] * vii [i];
6cb30a539481 untabify files in examples directory
John W. Eaton <jwe@octave.org>
parents: 9053
diff changeset
32 }
6593
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
33 }
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
34 else
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
35 {
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
36 for (i = 0; i < n; i++)
9932
6cb30a539481 untabify files in examples directory
John W. Eaton <jwe@octave.org>
parents: 9053
diff changeset
37 vro [i] = vri [i] * vri [i];
6593
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
38 }
3da1f4a41455 [project @ 2007-04-27 08:08:19 by dbateman]
dbateman
parents:
diff changeset
39 }