function tf = globalBodyLoad(self,varargin)
% globalBodyLoad  True if consistent across the geometric domain

% Copyright 2017 The MathWorks, Inc.

narginchk(1,2); 
tf = true;

if ~self.bodyLoadSpecified()   
    tf = false;
    return;
end


if self.nRegion==1
    % 3-D geometry with one cell or 2-D geometry with one face
    return;
end


thiscoef = self.findBodyLoad(self.qregiontype, 1);
for i = 2:self.nRegion
    othercoef = self.findBodyLoad(self.qregiontype, i);           
    if ~thiscoef.sameCoefficients(othercoef, varargin{:})
        tf = false;
        return;
    end            
end   



end