max/poet
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/poet.git synced 2025-12-16 12:54:50 +01:00
Code Issues Packages Projects Releases Wiki Activity

fix: set all values less than 10E-14 to zero

Browse Source
This commit is contained in:
Max Lübke 2023-01-13 13:21:58 +01:00 committed by Max Luebke
parent e71369bcdf
commit 9e4aea38e3
2 changed files with 17 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
include/poet/DiffusionModule.hpp
View File

@ -24,6 +24,7 @@
#include "poet/SimParams.hpp"
#include <array>
#include <bits/stdint-uintn.h>
#include <cmath>
#include <poet/Grid.hpp>
#include <string>
#include <tug/BoundaryCondition.hpp>
@ -85,6 +86,8 @@ private:
void initialize(poet::DiffusionParams args);
void RoundToZero(double *field, uint32_t cell_count) const;
Grid &grid;
uint8_t dim;

14
src/DiffusionModule.cpp
View File

@ -23,6 +23,7 @@
#include "tug/Diffusion.hpp"
#include <Rcpp.h>
#include <algorithm>
#include <bits/stdint-intn.h>
#include <cstdint>
#include <poet/ChemSimSeq.hpp>
#include <poet/DiffusionModule.hpp>
@ -36,6 +37,8 @@
using namespace poet;
static constexpr double ZERO_MULTIPLICATOR = 10E-14;
constexpr std::array<uint8_t, 4> borders = {
tug::bc::BC_SIDE_LEFT, tug::bc::BC_SIDE_RIGHT, tug::bc::BC_SIDE_TOP,
tug::bc::BC_SIDE_BOTTOM};
@ -180,6 +183,11 @@ void DiffusionModule::simulate(double dt) {
} else {
tug::diffusion::ADI_2D(this->diff_input, in_field, in_alpha.data());
}
// TODO: do not use hardcoded index for O, H and charge
if (i > 2) {
this->RoundToZero(in_field, this->n_cells_per_prop);
}
}
std::cout << " done!\n";
@ -188,6 +196,12 @@ void DiffusionModule::simulate(double dt) {
transport_t += sim_a_transport - sim_b_transport;
}
inline void DiffusionModule::RoundToZero(double *field,
uint32_t cell_count) const {
for (uint32_t i = 0; i < cell_count; i++) {
field[i] = ((int32_t)(field[i] / ZERO_MULTIPLICATOR)) * ZERO_MULTIPLICATOR;
}
}
void DiffusionModule::end() {
// R["simtime_transport"] = transport_t;