diff --git a/src/main.cpp b/src/main.cpp index 56fd09d0bf753961c2baeb5049e9b03e92e9b140..0ddfea5add6f8f11a837bfc54f070987f7dda43e 100644 --- a/src/main.cpp +++ b/src/main.cpp @@ -194,6 +194,19 @@ int main(int argc, char *argv[]) tempo.precision(5); tempo << std::fixed << t << '\n'; tempo.close(); + + // Fichier k indicateur + std::ofstream diff("diffinter"); + diff.precision(5); + for (size_t j=0; j<mesh.numberOfCells(); ++j) { + //diff << std::fixed << xj[j][0] << ' ' << std::fixed << kj[j] << '\n'; + if (kj[j]>0.) { + diff << std::fixed << xj[j][0] << ' ' << std::fixed << 4. << '\n'; + } else { + diff << std::fixed << xj[j][0] << ' ' << std::fixed << 0. << '\n'; + } + } + diff.close(); while((t<tmax) and (iteration<itermax)) { @@ -264,7 +277,7 @@ int main(int argc, char *argv[]) // ECRITURE DANS UN FICHIER - if ((std::fmod(t,0.001) < 0.0001) or (t == tmax)) { + //if ((std::fmod(t,0.001) < 0.0001) or (t == tmax)) { std::string ligne; @@ -333,8 +346,34 @@ int main(int argc, char *argv[]) tempo.precision(5); tempo << std::fixed << t << '\n'; tempo.close(); + + // k + + std::ifstream diffint("diffinter"); + std::ofstream diffout("k", std::ios::trunc); + diffout.precision(5); + for (size_t j = 0; j<mesh.numberOfCells(); ++j) { + getline(diffint, ligne); + //diffout << ligne << ' ' << std::fixed << xj[j][0] << ' ' << std::fixed << kj[j] << '\n'; + if (kj[j]>0.) { + diffout << ligne << ' ' << std::fixed << xj[j][0] << ' ' << std::fixed << 4. << '\n'; + } else { + diffout << ligne << ' ' << std::fixed << xj[j][0] << ' ' << std::fixed << 0. << '\n'; + } + } + diffint.close(); + diffout.close(); + std::ifstream riffint("k"); + std::ofstream riffout("diffinter", std::ios::trunc); + for (size_t j = 0; j<mesh.numberOfCells(); ++j) { + getline(riffint, ligne); + riffout << ligne << '\n'; } + riffint.close(); + riffout.close(); + + //} // ENTROPY TEST //finite_volumes_diffusion.entropie(unknowns);