matlab_class/%40gcmfaces/minus.m

function r = minus(p,q)
%overloaded gcmfaces minus function :
%  simply calls double minus function for each face data
%  if any of the two arguments is a gcmfaces object

if isa(p,'gcmfaces'); r=p; else; r=q; end; 
for iFace=1:r.nFaces;
   iF=num2str(iFace); 
   if isa(p,'gcmfaces')&isa(q,'gcmfaces'); 
       eval(['r.f' iF '=p.f' iF '-q.f' iF ';']); 
   elseif isa(p,'gcmfaces')&isa(q,'double');
       eval(['r.f' iF '=p.f' iF '-q;']);
   elseif isa(p,'double')&isa(q,'gcmfaces');
       eval(['r.f' iF '=p-q.f' iF ';']);
   else;
      error('gcmfaces minus: types are incompatible')
   end;
end;


1	gforget	1.1	function r = minus(p,q)
2	gforget	1.2	%overloaded gcmfaces minus function :
3			% simply calls double minus function for each face data
4			% if any of the two arguments is a gcmfaces object
5	gforget	1.1
6			if isa(p,'gcmfaces'); r=p; else; r=q; end;
7			for iFace=1:r.nFaces;
8			iF=num2str(iFace);
9			if isa(p,'gcmfaces')&isa(q,'gcmfaces');
10			eval(['r.f' iF '=p.f' iF '-q.f' iF ';']);
11			elseif isa(p,'gcmfaces')&isa(q,'double');
12			eval(['r.f' iF '=p.f' iF '-q;']);
13			elseif isa(p,'double')&isa(q,'gcmfaces');
14			eval(['r.f' iF '=p-q.f' iF ';']);
15			else;
16			error('gcmfaces minus: types are incompatible')
17			end;
18			end;
19
20