Skip to content

Commit 0dc10bf

Browse files
committed
Fix Segfault for Turb Linsol output. ContAdj with Frozen Turb.
1 parent 9f33eea commit 0dc10bf

2 files changed

Lines changed: 10 additions & 6 deletions

File tree

SU2_CFD/src/output/CAdjFlowCompOutput.cpp

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -295,9 +295,11 @@ void CAdjFlowCompOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, C
295295
SetHistoryOutputValue("LINSOL_ITER", adjflow_solver->GetIterLinSolver());
296296
SetHistoryOutputValue("LINSOL_RESIDUAL", log10(adjflow_solver->GetResLinSolver()));
297297

298-
if (turb_model) {
299-
SetHistoryOutputValue("LINSOL_ITER_TURB", adjturb_solver->GetIterLinSolver());
300-
SetHistoryOutputValue("LINSOL_RESIDUAL_TURB", log10(adjturb_solver->GetResLinSolver()));
298+
if ((!config->GetFrozen_Visc_Disc() && !cont_adj) || (!config->GetFrozen_Visc_Cont() && cont_adj)){
299+
if (turb_model) {
300+
SetHistoryOutputValue("LINSOL_ITER_TURB", adjturb_solver->GetIterLinSolver());
301+
SetHistoryOutputValue("LINSOL_RESIDUAL_TURB", log10(adjturb_solver->GetResLinSolver()));
302+
}
301303
}
302304

303305
if (config->GetDeform_Mesh()) {

SU2_CFD/src/output/CAdjFlowIncOutput.cpp

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -325,9 +325,11 @@ void CAdjFlowIncOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, CS
325325
SetHistoryOutputValue("LINSOL_ITER", adjflow_solver->GetIterLinSolver());
326326
SetHistoryOutputValue("LINSOL_RESIDUAL", log10(adjflow_solver->GetResLinSolver()));
327327

328-
if (turb_model) {
329-
SetHistoryOutputValue("LINSOL_ITER_TURB", adjturb_solver->GetIterLinSolver());
330-
SetHistoryOutputValue("LINSOL_RESIDUAL_TURB", log10(adjturb_solver->GetResLinSolver()));
328+
if (!config->GetFrozen_Visc_Disc() || !config->GetFrozen_Visc_Cont()){
329+
if (turb_model) {
330+
SetHistoryOutputValue("LINSOL_ITER_TURB", adjturb_solver->GetIterLinSolver());
331+
SetHistoryOutputValue("LINSOL_RESIDUAL_TURB", log10(adjturb_solver->GetResLinSolver()));
332+
}
331333
}
332334

333335
if (config->GetDeform_Mesh()) {

0 commit comments

Comments
 (0)