view examples/code/mycell.c @ 19631:aee5fea8a03e

compute condition number when solving triangular systems (bug #43875) * dMatrix.cc (Matrix::solve): Pass true for calc_rcond in calls to utsolve and ltsolve. * CMatrix.cc (ComplexMatrix::solve): Likewise. * fCMatrix.cc (FloatComplexMatrix::solve): Likewise. * fMatrix.cc (FloatMatrix::solve): Likeiwse.
author John W. Eaton <jwe@octave.org>
date Fri, 26 Dec 2014 14:35:46 -0500
parents c8240a60dd01
children
line wrap: on
line source

#include "mex.h"

void
mexFunction (int nlhs, mxArray* plhs[],
             int nrhs, const mxArray* prhs[])
{
  mwSize n;
  mwIndex i;

  if (nrhs != 1 || ! mxIsCell (prhs[0]))
    mexErrMsgTxt ("ARG1 must be a cell");

  n = mxGetNumberOfElements (prhs[0]);
  n = (n > nlhs ? nlhs : n);

  for (i = 0; i < n; i++)
    plhs[i] = mxDuplicateArray (mxGetCell (prhs[0], i));
}