0001 function R_mg=MpcGrids(grids,constraints)
0002
0003
0004
0005
0006
0007
0008
0009
0010 num_eq=0;
0011 gsize=getdofcount(grids);
0012
0013 for i=1:length(constraints),
0014 if strcmpi(constraints(i).constraint.type,'rgb'),
0015 num_eq=num_eq+1;
0016 end
0017 end
0018
0019
0020 R_mg=sparse(num_eq,gsize);
0021
0022 count=0;
0023 for i=1:length(constraints),
0024 if strcmpi(constraints(i).constraint.type,'rgb'),
0025
0026 rgbconstraint=constraints(i).constraint;
0027 count=count+1;
0028
0029 dof1=grids(rgbconstraint.grid1).grid.doflist(rgbconstraint.dof);
0030 dof2=grids(rgbconstraint.grid2).grid.doflist(rgbconstraint.dof);
0031
0032 R_mg(count,dof1)=1;
0033 R_mg(count,dof2)=-1;
0034 end
0035 end