A.1.8 Вызов функций Octave из oct-файлов

Часто возникает необходимость вызывать другую функцию Octave из oct-файла, и в самой Octave есть множество примеров таких вызовов. Например, функция quad — это oct-файл, вычисляющий определённый интеграл методом квадратур для функции, заданной пользователем.

Функцию можно передать в качестве входных данных разными способами. Например, это может быть

1. Дескриптор функции
2. Дескриптор анонимной функции
3. Строка

Приведённый ниже код демонстрирует все четыре способа передачи функции в oct-файл.

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

DEFMETHOD_DLD (funcdemo, interp, args, nargout, "Function Demo")
{
  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);

  octave_value_list retval;

  if (args(0).is_function_handle () || args(0).is_inline_function ()
      || args(0).is_string ())
    retval = interp.feval (args(0), newargs, nargout);
  else
    error ("funcdemo: INPUT must be string, inline, or function handle");

  return retval;
}

Первый входной аргумент демонстрационного кода — функция, заданная пользователем, а все остальные аргументы передаются этой функции.

funcdemo (@sin, 1)
⇒ 0.84147
funcdemo (@(x) sin (x), 1)
⇒ 0.84147
funcdemo ("sin", 1)
⇒ 0.84147
funcdemo (@atan2, 1, 1)
⇒ 0.78540

Если пользовательская функция передаётся в виде строки, с ней обращаются иначе. В некоторых случаях необходимо, чтобы заданная пользователем функция была объектом octave_function. В этом случае строковый аргумент можно использовать для создания временной функции, как показано ниже.

std::octave fcn_name = unique_symbol_name ("__fcn__");
std::string fcode = "function y = ";
fcode.append (fcn_name);
fcode.append ("(x) y = ");
fcn = extract_function (args(0), "funcdemo", fcn_name,
                        fcode, "; endfunction");
...
if (fcn_name.length ())
  clear_function (fcn_name);

В этом случае важно знать две вещи. Во-первых, число входных аргументов пользовательской функции фиксировано и в приведённом выше примере равно одному. Во-вторых, чтобы временная функция не осталась в таблице символов Octave, после использования её следует очистить. Кроме того, по соглашению все имена внутренних функций начинаются и заканчиваются последовательностью символов «__».