matlab_class/%40gcmfaces/and.m

function r = and(p,q)

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 and: types are incompatible')
   end;
end;


1	function r = and(p,q)
2
3	if isa(p,'gcmfaces'); r=p; else; r=q; end;
4	for iFace=1:r.nFaces;
5	iF=num2str(iFace);
6	if isa(p,'gcmfaces')&isa(q,'gcmfaces');
7	eval(['r.f' iF '=p.f' iF '&q.f' iF ';']);
8	elseif isa(p,'gcmfaces')&isa(q,'double');
9	eval(['r.f' iF '=p.f' iF '&q;']);
10	elseif isa(p,'double')&isa(q,'gcmfaces');
11	eval(['r.f' iF '=p&q.f' iF ';']);
12	else;
13	error('gcmfaces and: types are incompatible')
14	end;
15	end;
16
17