function pe=PenaltyCreatePVector(penpair,elements,grids,materials,inputs,analysis_type,kmax);
%CREATEPVECTOR - create an empty load vector for a pengrid 
%
%   penalty loads are expressed as penalty stiffnesses input in the global 
%   stiffness matrix. Therefore, penalty load is not a right hand side therm -> return empty vector.
%
%   Usage:
%      pe=CreatePVector(load,elements,grids,materials,inputs,analysis_type)
% 
%   See also PENALTYCREATEPVECTOR, PENALTYCREATEKMATRIX

%penalty loads are expressed as penalty stiffnesses input in the global 
%stiffness matrix. Therefore, penalty load is not a right hand side therm -> return empty vector.
pe=elemvector(0);
