Skip to content
Snippets Groups Projects
Commit 592ff829 authored by Jean-François Remacle's avatar Jean-François Remacle
Browse files

*** empty log message ***

parent 261ba262
No related branches found
No related tags found
No related merge requests found
...@@ -37,10 +37,10 @@ void gmshLaplaceTerm::elementMatrix(MElement *e, gmshMatrix<double> &m) const ...@@ -37,10 +37,10 @@ void gmshLaplaceTerm::elementMatrix(MElement *e, gmshMatrix<double> &m) const
Grads[j][2] = invjac[2][0] * grads[j][0] + invjac[2][1] * grads[j][1] + Grads[j][2] = invjac[2][0] * grads[j][0] + invjac[2][1] * grads[j][1] +
invjac[2][2] * grads[j][2]; invjac[2][2] * grads[j][2];
} }
double H=6.0; double H_x=1.0;
for (int j = 0; j < nbNodes; j++){ for (int j = 0; j < nbNodes; j++){
for (int k = 0; k <= j; k++){ for (int k = 0; k <= j; k++){
m(j, k) += (H*Grads[j][0] * Grads[k][0] + m(j, k) += (H_x*Grads[j][0] * Grads[k][0] +
Grads[j][1] * Grads[k][1] + Grads[j][1] * Grads[k][1] +
Grads[j][2] * Grads[k][2]) * weight * detJ * _diff; Grads[j][2] * Grads[k][2]) * weight * detJ * _diff;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment