substitute_backup.m

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

function C = substitute(A,B)

% substitute.m Interpolate a scalarfield object another scalarfield
%              with the same number of dimensions
%
%	INPUT parameters:
%             A,B : scalarfield objects 
%
%	OUTPUT:
%	        C : a scalarfield object

geomOld  = get(A.pd,'geom');
qpOld    = get(A.pd,'qp');
nameOld  = get(A.pd,'name');
axlimOld = get(A.pd,'axlim');

geomNew  = get(pdNew,'geom');
qpNew    = get(pdNew,'qp');
nameNew  = get(pdNew,'name');
axlimNew = get(pdNew,'axlim');

%%%%%%%%%%%%%%% 1D to 1D interpolation %%%%%%%%%%%%%%%

if length(nameOld) == 1

% interpolate old 1D curve to new 1D curve

   if length(nameNew) ~= 1
      error('Size mismatch')
      end

   if nnz(axlimOld{1}-axlimNew{1}) > 0
      error('Axis limit mismatch')
      end

   val = polyint(qpOld{1},A.val,qpNew{1});
   
   B = scalarfield(pdNew,val);
   return

   end

error('Not yet set up for higher-dimensional scalarfield objects')