|
Prev | Next | a_fun_property_xam.m | Headings |
function ok = a_fun_property_xam()
%
% load the Cppad Swig library
m_cppad
%
% initialize return variable
ok = true;
% -----------------------------------------------------------------------
n_ind = 1; % number of independent variables
n_dep = 2; % number of dependent variables
n_var = 1; % phantom variable at address 0
n_op = 1; % special operator at beginning
%
% dimension some vectors
x = m_cppad.vec_double(n_ind);
ay = m_cppad.vec_a_double(n_dep);
%
% independent variables
x(0) = 1.0;
ax = m_cppad.independent(x);
n_var = n_var + n_ind; % one for each indpendent
n_op = n_op + n_ind;
%
% first dependent variable
ay(0) = ax(0) + ax(0);
n_var = n_var + 1; % one variable and operator
n_op = n_op + 1;
%
% second dependent variable
ax0 = ax(0);
ay(1) = ax0.sin();
n_var = n_var + 2; % two varialbes, one operator
n_op = n_op + 1;
%
% define f(x) = y
af = m_cppad.a_fun(ax, ay);
n_op = n_op + 1; % speical operator at end
%
% check af properties
ok = ok && af.size_ind() == n_ind;
ok = ok && af.size_dep() == n_dep;
ok = ok && af.size_var() == n_var;
ok = ok && af.size_op() == n_op;
%
return;
end