Mercurial > octave-nkf
view examples/mypow2.c @ 18791:83c85d95ac7b stable
Fix assert() segfaults with profiler (bug #39587, bug #39586).
* profiler.cc (profile_data_accumulator::tree_node::exit,
profile_data_accumulator::exit_function, ): Comment
out assert() statements which are not valid when profile is called
from within a function.
* profiler.cc (F__profiler_reset__, F__profiler_data__,
F__profiler_enable__): Use correct function name in Texinfo
docstring.
* profiler.cc (F__profiler_data__): define retval(1) first to
resize return vector just once.
author | Rik <rik@octave.org> |
---|---|
date | Wed, 14 May 2014 13:59:44 -0700 |
parents | 224e76250443 |
children |
line wrap: on
line source
#include "mex.h" void mexFunction (int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]) { mwSize n; mwIndex i; double *vri, *vro; if (nrhs != 1 || ! mxIsNumeric (prhs[0])) mexErrMsgTxt ("ARG1 must be a matrix"); n = mxGetNumberOfElements (prhs[0]); plhs[0] = mxCreateNumericArray (mxGetNumberOfDimensions (prhs[0]), mxGetDimensions (prhs[0]), mxGetClassID (prhs[0]), mxIsComplex (prhs[0])); vri = mxGetPr (prhs[0]); vro = mxGetPr (plhs[0]); if (mxIsComplex (prhs[0])) { double *vii, *vio; vii = mxGetPi (prhs[0]); vio = mxGetPi (plhs[0]); for (i = 0; i < n; i++) { vro[i] = vri[i] * vri[i] - vii[i] * vii[i]; vio[i] = 2 * vri[i] * vii[i]; } } else { for (i = 0; i < n; i++) vro[i] = vri[i] * vri[i]; } }