error in PolarizationField corrected

extensions/base_problems/PolarizationField.hh

@@ -136,7 +136,7 @@ void PolarizationField<P>::fillOperators()
     for (size_t j = 0; j < self::dow; ++j) {
       Operator *opNonlin2 = new Operator(feSpace, feSpace);
       addZOT(opNonlin2, (-2.0*epsInv) * valueOf(vec[i]) * valueOf(vec[j]));
-      self::prob->addMatrixOperator(*opNonlin, i+self::dow, j);
+      self::prob->addMatrixOperator(*opNonlin2, i+self::dow, j);
     }
     
     Operator *opNonlin3 = new Operator(feSpace, feSpace);