if (west >= 0) {
      nonzeroIndices[numNonzeroEntries] = west*numEquations+component;
      if (east > 0)
	nonzeroValues[numNonzeroEntries]  = nu*lambda2X;
      else
	nonzeroValues[numNonzeroEntries]  = 2.0*nu*lambda2X;
      numNonzeroEntries++;
      }