diff --git a/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp b/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
index ec6fffe4269f3d3416d675bca09b3f29d6d8428d..7fd097bf7ef8fb4305f36f466420d3153958be3d 100644
--- a/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
+++ b/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
@@ -394,7 +394,6 @@ void mlawLinearElecMagTherMech::constitutive(
                 dqdV(i) += 0.0;
 
                 djydT(i) += (_seebeck + dseebeckdT * T) * js0(i);
-                //djydV(i) += 0.0;
 
                 // derivatives of sourceVectorField
                 dsourceVectorFielddV(i) += 0.0;
@@ -405,7 +404,6 @@ void mlawLinearElecMagTherMech::constitutive(
                     dfluxDdA(i,j) += 0.0;
                     dfluxTdA(i,j) += 0.0;
                     djydA(i,j) += 0.0;
-                    djydV(i) += _l0(i, j) * (An(j) - A0(j)) * invdh;
 
                     dsourceVectorFielddA(i, j) += 0.0;
                     dsourceVectorFielddB(i, j) += 0.0;
diff --git a/dG3D/src/dG3DTerms.cpp b/dG3D/src/dG3DTerms.cpp
index d2a81743a434b2b332cafd79f7326b00d0aa9605..98e1167cf318f11b9549e3894acc55fa4d460810 100644
--- a/dG3D/src/dG3DTerms.cpp
+++ b/dG3D/src/dG3DTerms.cpp
@@ -1859,13 +1859,13 @@ void dG3DStiffnessBulk::get(MElement *ele,int npts,IntPt *GP,fullMatrix<double>
                                          dHdB.operator()(m, n) *
                                          CurlcurlVals[k + numEdgeShapeFuncTotalLastCol][n]);
                                 mStiff(j + numEdgeShapeFuncTotalLastRow + numStandardDof,
-                                       k + numEdgeShapeFuncTotalLastCol + numStandardDof) -=
+                                       k + numEdgeShapeFuncTotalLastCol + numStandardDof) +=
                                         getConstitutiveCurlEqRatio() * ratio *
                                         (curlVals[j + numEdgeShapeFuncTotalLastRow][m] *
                                         dSourceVectorFielddVectorPotential.operator()(m, n) *
                                          curlVals[k + numEdgeShapeFuncTotalLastCol][n]);
                                 mStiff(j + numEdgeShapeFuncTotalLastRow + numStandardDof,
-                                       k + numEdgeShapeFuncTotalLastCol + numStandardDof) -=
+                                       k + numEdgeShapeFuncTotalLastCol + numStandardDof) +=
                                         getConstitutiveCurlEqRatio() * ratio *
                                         (curlVals[j + numEdgeShapeFuncTotalLastRow][m] *
                                         dSourceVectorFielddVectorCurl.operator()(m, n) *
@@ -1994,7 +1994,7 @@ void dG3DStiffnessBulk::get(MElement *ele,int npts,IntPt *GP,fullMatrix<double>
                                                                                                               (dVectorFielddExtraDofField.operator()(m) *
                                                                                                               CurlcurlVals[j+numEdgeShapeFuncTotalLastRow][m] *
                                                                                                               Vals[k+nbFFTotalLastCol]);
-                                  mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) -= getConstitutiveCurlEqRatio() * ratio *
+                                  mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) += getConstitutiveCurlEqRatio() * ratio *
                                                                                                               (dSourceVectorFielddExtraDofField.operator()(m) *
                                                                                                                curlVals[j+numEdgeShapeFuncTotalLastRow][m] *
                                                                                                                Vals[k+nbFFTotalLastCol]);
@@ -2004,7 +2004,7 @@ void dG3DStiffnessBulk::get(MElement *ele,int npts,IntPt *GP,fullMatrix<double>
                                                                                                                   (CurlcurlVals[j+numEdgeShapeFuncTotalLastRow][m] *
                                                                                                                   dVectorFielddGradExtraDofField.operator()(m,n) *
                                                                                                                   Grads[k+nbFFTotalLastCol][n]);
-                                      mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) -= getConstitutiveCurlEqRatio() * ratio *
+                                      mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) += getConstitutiveCurlEqRatio() * ratio *
                                                                                                                   (curlVals[j+numEdgeShapeFuncTotalLastRow][m] *
                                                                                                                   dSourceVectorFielddGradExtraDofField.operator()(m,n) *
                                                                                                                    Grads[k+nbFFTotalLastCol][n]);