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

Initial control-loop setup. Unit tests for MAPE and RMSE added.

Browse Source 
Note: computeStats function is not working correctly yet.
This commit is contained in:
rastogi 2025-09-01 12:54:24 +02:00
parent c9a8fc8079
commit d9d105b860
6 changed files with 0 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
bin/barite_fgcs_2_rt.R
View File

@ -1,7 +0,0 @@
iterations <- 200
dt <- 1000
list(
timesteps = rep(dt, iterations),
store_result = TRUE
)

BIN
bin/dinter_200/iter_00.qs2
View File

Binary file not shown.

BIN
bin/dinter_200/iter_050.qs2
View File

Binary file not shown.

BIN
bin/dinter_200/iter_100.qs2
View File

Binary file not shown.

BIN
bin/dinter_200/timings.qs2
View File

Binary file not shown.

13
src/Chemistry/MasterFunctions.cpp
View File

@ -190,14 +190,6 @@ void poet::ChemistryModule::computeStats(const std::vector<double> &pqc_vector,
double pqc_value = pqc_vector[i * size_per_prop + j]; double pqc_value = pqc_vector[i * size_per_prop + j];
double sur_value = sur_vector[i * size_per_prop + j]; double sur_value = sur_vector[i * size_per_prop + j];
if (i == 0 && (j % 10000 == 0)) {
std::cout << "i=" << i << ", j=" << j
<< ", idx=" << i * size_per_prop + j
<< ", pqc=" << pqc_value
<< ", sur=" << sur_value
<< std::endl;
}
if (pqc_value != 0) if (pqc_value != 0)
{ {
double rel_err = (pqc_value - sur_value) / pqc_value; double rel_err = (pqc_value - sur_value) / pqc_value;
@ -211,11 +203,6 @@ if (i == 0 && (j % 10000 == 0)) {
sqr_err_sum += 1.0; sqr_err_sum += 1.0;
count++; count++;
} }
if (pqc_value == 0 && sur_value == 0)
{
}
// else: both cases are zero, skip (no error) // else: both cases are zero, skip (no error)
} }
if (i == 0) if (i == 0)