% COPY makes a (deep) copy of the input parfrac objects.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: COPY makes a deep copy of the input parfrac objects.
%
% CALL:        b = copy(a, flag)
%
% INPUTS:      a    - input parfrac object
%              flag - 1: make a deep copy, 0: return copies of handles
%
% OUTPUTS:     b - copy of inputs
%
% This is a transparent function and adds no history.
%
% VERSION:     $Id: copy.m,v 1.10 2010/06/25 15:11:47 ingo Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function varargout = copy(old, deepcopy)
  
  if deepcopy
    % Loop over input parfrac objects
    new = parfrac.newarray(size(old));
    obj = copy@ltpda_tf(new, old, 1);
    
    for kk=1:numel(old)
      obj(kk).res    = old(kk).res;
      obj(kk).poles  = old(kk).poles;
      obj(kk).pmul   = old(kk).pmul;
      obj(kk).dir    = old(kk).dir;
    end
  else
    obj = old;
  end
  
  varargout{1} = obj;
end