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

Sned control flgas after control iteration and enable Cl threshold

Browse Source
This commit is contained in:
rastogi 2025-12-10 09:03:11 +01:00
parent 36b6f8d859
commit 69828116f3
3 changed files with 18 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
bin/run_poet.sh
View File

@ -1,7 +1,7 @@
#!/bin/bash
#SBATCH --job-name=p1_eps0035_200
#SBATCH --output=p1_eps0035_200_%j.out
#SBATCH --error=p1_eps0035_200_%j.err
#SBATCH --job-name=p1_eps0035_cl_200_r2
#SBATCH --output=p1_eps0035_cl_200_r2_%j.out
#SBATCH --error=p1_eps0035_cl_200_r2_%j.err
#SBATCH --partition=long
#SBATCH --nodes=6
#SBATCH --ntasks-per-node=24
@ -15,5 +15,5 @@ module purge
module load cmake gcc openmpi
#mpirun -n 144 ./poet dolo_fgcs_3.R dolo_fgcs_3.qs2 dolo_only_pqc
mpirun -n 144 ./poet --interp dolo_fgcs_3_rt.R dolo_fgcs_3.qs2 p1_eps0035_200
mpirun -n 144 ./poet --interp dolo_fgcs_3_rt.R dolo_fgcs_3.qs2 p1_eps0035_cl_200_r2
#mpirun -n 144 ./poet --interp barite_fgcs_4_new/barite_fgcs_4_new_rt.R barite_fgcs_4_new/barite_fgcs_4_new.qs2 barite

2
src/Chemistry/WorkerFunctions.cpp
View File

@ -217,14 +217,12 @@ void poet::ChemistryModule::WorkerDoWork(MPI_Status &probe_status, int double_co
}
/* skip simulation of cells cells where Cl concentration is below threshold */
/*
for (std::size_t wp_i = 0; wp_i < s_curr_wp.size; wp_i++) {
if (s_curr_wp.input[wp_i][CL_INDEX] < CL_THRESHOLD) {
s_curr_wp.mapping[wp_i] = CHEM_SKIP;
s_curr_wp.output[wp_i] = s_curr_wp.input[wp_i];
}
}
*/
// std::cout << this->comm_rank << ":" << counter++ << std::endl;
if (dht_enabled || interp_enabled || stab_enabled) {

16
src/Control/ControlModule.cpp
View File

@ -98,9 +98,12 @@ void poet::ControlModule::readCheckpoint(uint32_t &current_iter, uint32_t rollba
void poet::ControlModule::writeMetrics(const std::string &out_dir,
const std::vector<std::string> &species) {
if (rb_count > config.rb_limit) {
/*
if (rb_count > config.rb_limit) {
return;
}
*/
double stats_a, stats_b;
stats_a = MPI_Wtime();
@ -235,7 +238,16 @@ void poet::ControlModule::processCheckpoint(uint32_t &current_iter,
}
bool poet::ControlModule::needsFlagBcast() const {
return (config.rb_limit > 0) && !rbLimitReached();
// return (config.rb_limit > 0) && !rbLimitReached();
/*
if (rb_count > config.rb_limit) {
return false;
}
*/
if (global_iter == 1 || global_iter % config.ctrl_interval == 1) {
return true;
}
return false;
}
inline bool poet::ControlModule::rbLimitReached() const {