#include "npglvec.h"
#include "npsolvert.h"
#include "globalt.h"
#include "globmatt.h"
#include "elemswitcht.h"
#include "transprint.h"
#include "saddpoint.h"
#include "backupsolt.h"
#include <string.h>
#include <errno.h>
Go to the source code of this file.
Functions | |
void | linear_nonstat_radiation_solv_dform () |
void | linear_nonstat_solv_dform () |
void | linear_nonstat_solv_dform_subcycl () |
void | linear_nonstat_solv_vform () |
void | nonstat_solv_vform_comp () |
void | nonstat_solver_init (long lcid, np_glob_vec &np_gv) |
long | one_step_linear (long lcid, double time, double dt, long istep, long li, np_glob_vec &np_gv) |
long | one_step_nonlinear (long lcid, double time, double dt, long istep, long li, np_glob_vec &np_gv) |
void | solve_nonstationary_problem () |
void linear_nonstat_radiation_solv_dform | ( | ) |
function solves linear nonstationary transport problem the d-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
JK, 23.12.2002, revised 24.8.2006
Definition at line 1013 of file npsolvert.cpp.
References actual_previous_nodval(), timecontr::actualbacktimeincr(), gmatrix::addgm(), probdesct::alpha, approximation(), capacity_matrix(), Cmat, compute_req_valt(), conductivity_matrix(), transtop::edge_temperature(), timecontr::endtime(), f, probdesct::fluxcomp, lhsrhst::give_lhs(), lhsrhst::give_rhs(), lhsrhst::give_tdlhs(), gmatrix::gmxv(), Gtt, probdesct::hdbcont, hdbcontr::hdbtype, transtop::heat_fluxes(), timecontr::initialtimeincr(), transmat::initmaterialmodels(), internal_fluxes(), timecontr::isitimptime(), Kmat, Lsrst, Mesprt, Ndoft, timecontr::newtime(), nullv(), probdesct::nvs, Outt, p, probdesct::pnvs, print_closet(), print_flusht(), print_initt(), print_stept(), hdbcontr::restore_stat(), gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), solvert_restore(), solvert_save(), probdesct::ssle, timecontr::starttime(), probdesct::time, probdesct::timecont, Tm, Tp, trfel_right_hand_side(), Tt, and transmat::updateipval().
void linear_nonstat_solv_dform | ( | ) |
function solves linear nonstationary transport problem the d-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
JK, 23.12.2002, revised 24.8.2006
Definition at line 810 of file npsolvert.cpp.
References actual_previous_nodval(), timecontr::actualbacktimeincr(), gmatrix::addgm(), probdesct::alpha, approximation(), capacity_matrix(), Cmat, compute_req_valt(), conductivity_matrix(), timecontr::endtime(), f, probdesct::fluxcomp, lhsrhst::give_lhs(), lhsrhst::give_rhs(), lhsrhst::give_tdlhs(), gmatrix::gmxv(), Gtt, probdesct::hdbcont, hdbcontr::hdbtype, timecontr::initialtimeincr(), transmat::initmaterialmodels(), internal_fluxes(), timecontr::isitimptime(), Kmat, Lsrst, Mesprt, Ndoft, timecontr::newtime(), nullv(), probdesct::nvs, Outt, p, probdesct::pnvs, print_closet(), print_flusht(), print_initt(), print_stept(), hdbcontr::restore_stat(), gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), solvert_restore(), solvert_save(), probdesct::ssle, timecontr::starttime(), probdesct::time, probdesct::timecont, Tm, Tp, trfel_right_hand_side(), and transmat::updateipval().
Referenced by solve_nonstationary_problem().
void linear_nonstat_solv_dform_subcycl | ( | ) |
function solves linear nonstationary transport problem time integration is based on subcycling the d-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
JK, 27. 5. 2014
Definition at line 1206 of file npsolvert.cpp.
References actual_previous_nodval(), timecontr::actualbacktimeincr(), gmatrix::addgm(), probdesct::alpha, approximation(), capacity_matrix(), Cmat, compute_req_valt(), conductivity_matrix(), seqtop::elem_lists(), timecontr::endtime(), f, probdesct::fluxcomp, lhsrhst::give_lhs(), lhsrhst::give_rhs(), lhsrhst::give_tdlhs(), gmatrix::gmxv(), Gtt, probdesct::hdbcont, hdbcontr::hdbtype, timecontr::initialtimeincr(), transmat::initmaterialmodels(), internal_fluxes(), timecontr::isitimptime(), Kmat, Lsrst, Mesprt, Ndoft, timecontr::newtime(), nullv(), probdesct::nvs, Outt, p, probdesct::pnvs, print_closet(), print_flusht(), print_initt(), print_stept(), hdbcontr::restore_stat(), gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), solvert_restore(), solvert_save(), probdesct::ssle, timecontr::starttime(), gtopology::stop, probdesct::time, probdesct::timecont, Tm, Tp, trfel_right_hand_side(), and transmat::updateipval().
void linear_nonstat_solv_vform | ( | ) |
function solves linear nonstationary transport problem the v-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
JK, 20.12.2002, revised 24.8.2006
termitovo kontrola transferu state variables DOCASNE
adaptivity
Definition at line 616 of file npsolvert.cpp.
References actual_previous_nodval(), timecontr::actualbacktimeincr(), probdesct::adaptivityflag, Adat, gmatrix::addgm(), probdesct::alpha, adaptivityt::answer, approximation(), capacity_matrix(), Cmat, compute_req_valt(), conductivity_matrix(), timecontr::endtime(), f, lhsrhst::give_lhs(), lhsrhst::give_rhs(), adaptivityt::give_step(), lhsrhst::give_tdlhs(), gmatrix::gmxv(), Gtt, probdesct::hdbcont, timecontr::initialtimeincr(), transmat::initmaterialmodels(), timecontr::isitimptime(), probdesct::istep, Kmat, Lsrst, Mesprt, Ndoft, timecontr::newtime(), nullv(), probdesct::nvs, Outt, p, probdesct::pnvs, print_closet(), print_flusht(), print_initt(), print_stept(), hdbcontr::restore_stat(), adaptivityt::run(), hdbcontr::save_stat(), gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), solvert_restore(), solvert_save(), probdesct::ssle, timecontr::starttime(), adaptivityt::statedata_restore(), probdesct::time, probdesct::timecont, Tm, Tp, trfel_right_hand_side(), and transmat::updateipval().
void nonstat_solv_vform_comp | ( | ) |
function solves linear nonstationary transport problem the v-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
according to JK, 20.12.2002 and TKo mtsolver.cpp revised 05/11/2012 by TKr
Definition at line 44 of file npsolvert.cpp.
References timecontr::actualbacktimeincr(), timecontr::actualforwtimeincr(), probdesct::adaptivityflag, Adat, approximation(), compute_req_valt(), timecontr::dtmax, timecontr::dtmin, timecontr::endtime(), np_glob_vec::f, probdesct::hdbcont, timecontr::initialtimeincr(), timecontr::isitimptime(), probdesct::istep, np_glob_vec::istep, np_glob_vec::lhs, Mesprt, Ndoft, timecontr::newtime(), nonlinear_nonstationary_problem, nonstat_solver_init(), nonstationary_problem, nullv(), timecontr::oldtime(), one_step_linear(), one_step_nonlinear(), print_closet(), print_flusht(), print_stept(), print_stept_forced(), np_glob_vec::rhs, adaptivityt::run(), hdbcontr::save_stat(), solvert_save(), stat, timecontr::tct, np_glob_vec::tdlhs, gfunct::tfunc, probdesct::time, probdesct::timecont, timecontr::timefun, Tp, and probdesct::tprob.
Referenced by solve_nonstationary_problem().
void nonstat_solver_init | ( | long | lcid, | |
np_glob_vec & | np_gv | |||
) |
Function allocates and initializes global vectors used in npsolver It is used in one_step concept solver.
Created by tomas Krejci according to Tomas Koudelka, 11/2012
termitovo
termitovo kontrola transferu state variables DOCASNE
Definition at line 208 of file npsolvert.cpp.
References actual_previous_nodval(), Adat, np_glob_vec::alloc(), np_glob_vec::alloc_aux(), adaptivityt::answer, compute_req_valt(), np_glob_vec::f, lhsrhst::give_lhs(), lhsrhst::give_rhs(), adaptivityt::give_step(), lhsrhst::give_tdlhs(), probdesct::hdbcont, timecontr::initialtimeincr(), transmat::initmaterialmodels(), np_glob_vec::istep, np_glob_vec::lhs, Lsrst, Mesprt, Ndoft, nonlinear_nonstationary_problem, probdesct::nvs, probdesct::pnvs, print_flusht(), print_initt(), print_stept(), hdbcontr::restore_stat(), np_glob_vec::rhs, adaptivityt::run(), solvert_restore(), timecontr::starttime(), adaptivityt::statedata_restore(), np_glob_vec::tdlhs, probdesct::time, probdesct::timecont, Tm, Tp, and probdesct::tprob.
Referenced by newton_raphson_parcoupl_common_dt(), newton_raphson_parcoupl_comp(), nonstat_solv_vform_comp(), and Ctlinterface::solve().
long one_step_linear | ( | long | lcid, | |
double | time, | |||
double | dt, | |||
long | istep, | |||
long | li, | |||
np_glob_vec & | np_gv | |||
) |
function solves linear nonstationary transport problem the v-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
according to JK, 20.12.2002, revised 05/11/2012 by TKr
Definition at line 284 of file npsolvert.cpp.
References actual_previous_nodval(), gmatrix::addgm(), probdesct::alpha, approximation(), capacity_matrix(), Cmat, conductivity_matrix(), np_glob_vec::d, np_glob_vec::f, f, gmatrix::gmxv(), Gtt, probdesct::istep, probdesct::jstep, Kmat, np_glob_vec::lhs, Mesprt, Ndoft, probdesct::nvs, Outt, np_glob_vec::p, p, probdesct::pnvs, np_glob_vec::rhs, gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), probdesct::ssle, np_glob_vec::tdlhs, probdesct::time, Tm, Tp, trfel_right_hand_side(), and transmat::updateipval().
Referenced by newton_raphson_parcoupl_common_dt(), newton_raphson_parcoupl_comp(), nonstat_solv_vform_comp(), and Ctlinterface::solve().
long one_step_nonlinear | ( | long | lcid, | |
double | time, | |||
double | dt, | |||
long | istep, | |||
long | li, | |||
np_glob_vec & | np_gv | |||
) |
function solves non-linear nonstationary transport problem by Newton-Raphson method the v-form version of the generalized trapezoidal method see T.J.R. Hughes: The Finite Element Method; Linear Static and Dynamic Finite Element Analysis
according to TKo revised 05/11/2012 by TKr
Definition at line 405 of file npsolvert.cpp.
References actual_previous_nodval(), gmatrix::addgm(), probdesct::alpha, approximation(), capacity_matrix(), Cmat, conductivity_matrix(), probdesct::convergcontrolt, gmatrix::copygm(), np_glob_vec::d, probdesct::errarr, np_glob_vec::f, f, np_glob_vec::fb, np_glob_vec::fi, fluxest, fullnewtont, gmatrix::gmxv(), Gtt, internal_fluxes(), probdesct::istep, probdesct::jstep, Kmat, np_glob_vec::lhs, lrhst, Mesprt, Ndoft, probdesct::nii, nullv(), probdesct::nvs, Outt, np_glob_vec::p, p, probdesct::pnvs, np_glob_vec::rhs, gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), ss(), probdesct::ssle, np_glob_vec::tdlhs, probdesct::threshrhs, probdesct::time, Tm, Tp, probdesct::trestype, trfel_right_hand_side(), probdesct::trsolv, transmat::updateipval(), yes, and probdesct::zero.
Referenced by newton_raphson_parcoupl_common_dt(), newton_raphson_parcoupl_comp(), nonstat_solv_vform_comp(), and Ctlinterface::solve().
void solve_nonstationary_problem | ( | ) |
function solves nonstationary transport problem
function selects a suitable numerical scheme for solution
JK
Definition at line 19 of file npsolvert.cpp.
References discont_nonstat_problem, linear_nonstat_solv_dform(), nonstat_solv_vform_comp(), Tp, and probdesct::tprob.
Referenced by solve_trfel_deterministic_problem().