--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/@pest/copy.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,52 @@
+% COPY makes a (deep) copy of the input pest objects.
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% DESCRIPTION: COPY makes a deep copy of the input pest objects.
+%
+% CALL:        b = copy(a, flag)
+%
+% INPUTS:      a    - input pest 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.4 2010/06/25 15:11:47 ingo Exp $
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function varargout = copy(old, deepcopy)
+  
+  if deepcopy
+    % Loop over input pest objects
+    new = pest.newarray(size(old));
+    obj = copy@ltpda_uoh(new, old, 1);
+    
+    for kk=1:numel(old)
+      
+      if ~isempty(old(kk).yunits)
+        obj(kk).yunits = old(kk).yunits;
+      end
+      if ~isempty(old(kk).models)
+        obj(kk).models = old(kk).models;
+      end
+      
+      obj(kk).y     = old(kk).y;
+      obj(kk).dy    = old(kk).dy;
+      obj(kk).names = old(kk).names;
+      obj(kk).pdf   = old(kk).pdf;
+      obj(kk).cov   = old(kk).cov;
+      obj(kk).corr  = old(kk).corr;
+      obj(kk).chi2  = old(kk).chi2;
+      obj(kk).dof   = old(kk).dof;
+      obj(kk).chain = old(kk).chain;
+      
+    end
+  else
+    obj = old;
+  end
+  
+  varargout{1} = obj;
+end
+