X - UWSpace - University of Waterloo

if(LpSub[cur-blockl->lo-bnds[j] != 0.0 && dj[j] != 0.0)
LpSub [CU-block] ->PiB=LpSub [curubblock] ->FiB+
LpSub [cur-block] -~l~~bnds f j 1
Cdj [jl) ;
eise if(cstat[jl == 2)
double *CompuPiD(LpSubProb "LpSub) /* Compute PiD to send */
I int j,p;
for (j=O; jPiD[jI=O.O;
for (p=LpSub [cur-blockl ->clpnts ( j 1 ; p~LpSub[cur,block] ->clpncs [ j +l 1; p++ 1
if(pi[LpSub~curblock]-~rwnmbs[p]-l-HISPbbm[cur~b~ock-1]] != 0.0 &&
LpSub[cur-block] ->rwnmbs [pl >HISPb-m[curUTb1ock-l
LpSub l cur-block] ->rwnnibs [pl < SPb-rnlcur-block] -LOWSPb-m[curU1:blockl
LpSub[cur,block] ->PiD [ j 1 += pi [LpSub[cur-block] ->rWMLbs [pl -1-
HISPb-m[cur,block-11 1~LpSub[cur~bl0~kl-~coeffs[pl;
1
return ( LpSub [cur-black] ->PiD) ;
1
1 &&
double *CompuOMeL(LpSubProb "LpSubi /* Compute OmegaL */
E int f. k, p;
for(j=O; j < HISPb-n[nb,subsl; j++)
LpSub [cur,block] ->OMeL[li = 0.0 ;
fcr(k=l; k clpncs [ j+lI ;p++)
if (LpSub [k] - > d s [pl >= SPb-m[kl -LOWSPb,m[cur,block-11 &&
LpSub lkl ->rllrnmbs [pl c SPb-rnikl hh pi [LpSub[k] -mmmbs [pl -SPbbm[k3 +
LOWSPb~m[cur,block-l]+LSPbbm[~CUfblock-L 1 != 0.0)
LpS~&[c=,block:->OMeL[j+HISPb,n[k-111 += pi[LpSub[k]->rwnmbs[p]-
S?b~m[k~+LOWSPb~m[curElock-~]+LS?b~m[cur~blo~k-l]]~LpSub[kI-~coeffs[p];
t
for (p=~p~ub[kl ->clpnt.s [ j I ;pclpnts [ j+l] ;p++)