Skip to content
Snippets Groups Projects
Commit 44e7d147 authored by Fanny CHOPOT's avatar Fanny CHOPOT
Browse files

Correction solution exacte E pour k constant

parent bc923159
No related branches found
No related tags found
No related merge requests found
...@@ -216,7 +216,7 @@ int main(int argc, char *argv[]) ...@@ -216,7 +216,7 @@ int main(int argc, char *argv[])
double pi = 4.*std::atan(1.); double pi = 4.*std::atan(1.);
std::ofstream fout("comparaison E"); std::ofstream fout("comparaison E");
for (size_t j=0; j<mesh.numberOfCells(); ++j) { for (size_t j=0; j<mesh.numberOfCells(); ++j) {
fout << xj[j][0] << ' ' << Ej[j] << ' ' << (std::cos(pi*xj[j][0]))*(std::cos(pi*xj[j][0]))*std::sin(pi*xj[j][0])*((pi*pi*pi)*0.5)*std::exp(-4.*pi*0.2) <<'\n'; fout << xj[j][0] << ' ' << Ej[j] << ' ' << (-(std::cos(pi*xj[j][0]))*(std::cos(pi*xj[j][0]))+(std::sin(pi*xj[j][0])*std::sin(pi*xj[j][0])))*(pi*0.5)*std::exp(-4.*pi*0.2) <<'\n';
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment