27 const PetscInt lbus[3] = {4, 5, 7}; /* Buses at which loads are incident */  variable
294     Vr  = xnet[2 * lbus[i]];     /* Real part of load bus voltage */  in ResidualFunction()
295     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualFunction()
298 …Vm0 = PetscSqrtScalar(v0[2 * lbus[i]] * v0[2 * lbus[i]] + v0[2 * lbus[i] + 1] * v0[2 * lbus[i] + 1…  in ResidualFunction()
307     fnet[2 * lbus[i]] += IDi;  in ResidualFunction()
308     fnet[2 * lbus[i] + 1] += IDr;  in ResidualFunction()
659     Vr  = xnet[2 * lbus[i]];     /* Real part of load bus voltage */  in ResidualJacobian()
660     Vi  = xnet[2 * lbus[i] + 1]; /* Imaginary part of the load bus voltage */  in ResidualJacobian()
664 …Vm0 = PetscSqrtScalar(v0[2 * lbus[i]] * v0[2 * lbus[i]] + v0[2 * lbus[i] + 1] * v0[2 * lbus[i] + 1…  in ResidualJacobian()
688     row[0] = net_start + 2 * lbus[i];  in ResidualJacobian()
689     col[0] = net_start + 2 * lbus[i];  in ResidualJacobian()
690     col[1] = net_start + 2 * lbus[i] + 1;  in ResidualJacobian()
695     row[0] = net_start + 2 * lbus[i] + 1;  in ResidualJacobian()
696     col[0] = net_start + 2 * lbus[i];  in ResidualJacobian()
697     col[1] = net_start + 2 * lbus[i] + 1;  in ResidualJacobian()