Ludovic Noels · 6fe00242 · 7d9231d3 · 7d9231d3 · c5d4aceb · 9204af2a
--- a/dG3D/src/dG3DTerms.cpp
+++ b/dG3D/src/dG3DTerms.cpp
                                         dHdB.operator()(m, n) *
                                         CurlcurlVals[k + numEdgeShapeFuncTotalLastCol][n]);
                                mStiff(j + numEdgeShapeFuncTotalLastRow + numStandardDof,
-                                       k + numEdgeShapeFuncTotalLastCol + numStandardDof) +=
+                                       k + numEdgeShapeFuncTotalLastCol + numStandardDof) -=
--- a/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
+++ b/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp

                    // derivatives of sourceVectorField
                    dsourceVectorFielddB(i, j) += 0.0;
-                    dsourceVectorFielddgradV(i, j) += _l0(i, j);
-                    dsourceVectorFielddgradT(i, j) += _seebeck * _l0(i ,j);
+                    dsourceVectorFielddgradV(i, j) += -_l0(i, j);
+                    dsourceVectorFielddgradT(i, j) += -_seebeck * _l0(i ,j);
--- a/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
+++ b/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp

                    // derivatives of sourceVectorField
                    dsourceVectorFielddB(i, j) += 0.0;
-                    dsourceVectorFielddgradV(i, j) += _l0(i, j);
-                    dsourceVectorFielddgradT(i, j) += _seebeck * _l0(i ,j);
+                    dsourceVectorFielddgradV(i, j) += -_l0(i, j);
+                    dsourceVectorFielddgradT(i, j) += -_seebeck * _l0(i ,j);

                    dsourceVectorFielddV(i) += 0.0;
-                    dsourceVectorFielddT(i) += dseebeckdT * _l0(i, j) * gradT(j);
-                    dsourceVectorFielddA(i, j) += _l0(i, j) * invdh;
+                    dsourceVectorFielddT(i) += -dseebeckdT * _l0(i, j) * gradT(j);
+                    dsourceVectorFielddA(i, j) += -_l0(i, j) * invdh;
--- a/dG3D/src/dG3DTerms.cpp
+++ b/dG3D/src/dG3DTerms.cpp
                  {
                      // related to cp
                      mStiff(j+nbFFTotalLastRow,k+nbFFTotalLastCol)+=getConstitutiveExtraDofDiffusionEqRatio()*ratio*dwdt*Vals[j+nbFFTotalLastRow]*Vals[k+nbFFTotalLastCol];
+                      mStiff(j+nbFFTotalLastRow,k+nbFFTotalLastCol) -= getConstitutiveExtraDofDiffusionEqRatio()*ratio*dw_AVdT*Vals[j+nbFFTotalLastRow]*Vals[k+nbFFTotalLastCol];