view examples/code/funcdemo.cc @ 20596:87b557ee8e5d

clean up and vectorize code for dense output in ode45 * scripts/ode/private/ode_rk_interpolate.m: new file * scripts/ode/private/ode_rk_interpolate.m(hermite_quartic_interpolation): move to internal function, use vectorization and broadcasting. * scripts/ode/private/hermite_quartic_interpolation.m: remove file * scripts/ode/module.mk: list added and removed files * scripts/ode/private/integrate_adaptive.m: use new interpolation code.
author Carlo de Falco <carlo.defalco@polimi.it>
date Tue, 06 Oct 2015 19:28:59 +0200
parents 2f8500ca91d3
children
line wrap: on
line source

#include <octave/oct.h>
#include <octave/parse.h>

DEFUN_DLD (funcdemo, args, nargout, "Function Demo")
{
  octave_value_list retval;

  int nargin = args.length ();

  if (nargin < 2)
    print_usage ();

  octave_value_list newargs;

  for (octave_idx_type i = nargin - 1; i > 0; i--)
    newargs(i-1) = args(i);

  if (args(0).is_function_handle () || args(0).is_inline_function ())
    {
      octave_function *fcn = args(0).function_value ();

      retval = feval (fcn, newargs, nargout);
    }
  else if (args(0).is_string ())
    {
      std::string fcn = args(0).string_value ();

      retval = feval (fcn, newargs, nargout);
    }
  else
    error ("funcdemo: INPUT must be string, inline, or function handle");

  return retval;
}