naemodel copy.m

MDPSAS工具箱是马里兰大学开发的

function A = naemodel(classname,var,param,solverset)

% naemodel.m A nonlinear algebraic equation model (naemodel) object
%            constructor method. Called as
%
%          A = naemodel(classname,var,param,solverset)
%
% INPUT PARAMETERS
%  classname : a string of module class name
%        var : a MATLAB assocarray object for the variables
%      param : a MATLAB assocarray object for the parameters
%  solverset : a MATLAB assocarray with the following keys
%               nmax  : max Newton iterations (default = 10)
%               lsmax : max line search iterations (default = 10)
%               jup   : iterations performed before jacobian is
%                              updated (default = 1)
%               nkvec : number of Krylov vectors to be used
%                              in newtonkrylov.m (default = 10)
%               vmin  : lower bound on variable values after update
%                              (default = [])
%               vmax  : upper bound on variable values after update
%                              (default = [])
%               plot  : if set to 1, plots solutions at each iteration
%                              (default = 0)
%
% OBJECTS RETURNED
%          A : an naemodel object with data fields
%               A.var       : variable names and values
%               A.param     : parameter names and values
%               A.id        : double value
%               A.solverset : parameters for the various solvers
%               A.resid     : residual value cell array
%               A.residflag : 1 if the current residual is valid
%                             0 if any state or parameter was reset

% Copyright (c) by Raymond A. Adomaitis, 1998-2005

if nargin < 3
   param = assocarray;
   end

if nargin < 2
   var = assocarray;
   end
   
if nargin < 1
   error('wrong input')
   end

if ~isa(var,'assocarray')
   error('var must be an associative array (assocarray)')
   end

if ~isa(param,'assocarray')
   error('param must be an associative array (assocarray)')
   end

if ~isa(classname,'char')
    error('classname must be a char/string')
end
   
% RAA 5/12/04

% vmin = var;
% vmax = vmin;
% vkey = keys(vmin);

% if ~( length(vkey) == 1 & isempty(vkey{1}) )
%    for i = 1:length(vkey)
%       vmin = setval(vmin,-1e99+0*getval(vmin,vkey{i}),vkey{i});
%       vmax = setval(vmax, 1e99+0*getval(vmax,vkey{i}),vkey{i});
%       end
%    end

if nargin < 4

   vmin = [];
   vmax = [];

   solverset = assocarray({'nmax',10,'lsmax',10,'jup',1,'nkvec',10,'plot',0,...
                           'vmin',vmin,'vmax',vmax});
   end

if ~isa(solverset,'assocarray')
   error('solverset must be an associative array (assocarray)')
   end

solverset = setval(solverset,vmin,'vmin');
solverset = setval(solverset,vmax,'vmax');

B = solver(classname,var,param,solverset);
A = class(struct([]),'naemodel',B);