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

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

Существует также множество способов передачи функции в качестве входных данных. Она может быть передана как:

1. Дескриптор функции
2. Анонимный дескриптор функции
3. Встроенная функция
4. Строка

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

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

DEFUN_DLD (funcdemo, 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 ())
    {
      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;
}

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

funcdemo (@sin, 1)
⇒ 0.84147
funcdemo (@(x) sin (x), 1)
⇒ 0.84147
funcdemo (inline ("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, её необходимо очистить после использования. Кроме того, по соглашению все внутренние имена функций начинаются и заканчиваются последовательностью символов «__».