Wrapper for the 3rd party RADAU5 solver. More...
Functions/Subroutines | |
| subroutine | fcn (n, t, y, f) |
| subroutine | jac () |
| subroutine | solout () |
| subroutine | radau5_solve (trx, itot, non, n0, n, delta_t, diag, a, b, zero_flux_step, is_adjoint, abstol, reltol) |
Variables | |
| type(origen_transitionmatrixp) | xtrx |
| Library data (NOTE: THIS CLASS IS NOT THREAD SAFE!) More... | |
| integer, dimension(:), allocatable | xloc |
| parent locations More... | |
| real(c_double), dimension(:), allocatable | xd |
| Diagonal elements of the transition matrix. More... | |
| real(c_double), dimension(:), allocatable | xa |
| Nondiagonal elements of the transition matrix. More... | |
| real(c_double), dimension(:), allocatable | xb |
| Right side of the equation (source term) More... | |
| real(c_double) | fake_flux |
| fake flux level to last_transition is decay if=0 More... | |
Detailed Description
Wrapper for the 3rd party RADAU5 solver.
Function/Subroutine Documentation
| subroutine kernel_radau5::fcn | ( | integer, intent(in) | n, |
| real(c_double), intent(in) | t, | ||
| real(c_double), dimension(n), intent(in) | y, | ||
| real(c_double), dimension(n), intent(out) | f | ||
| ) |
| subroutine kernel_radau5::jac | ( | ) |
Referenced by radau5_solve().
| subroutine kernel_radau5::solout | ( | ) |
Referenced by radau5_solve().
| subroutine kernel_radau5::radau5_solve | ( | type(origen_transitionmatrixp), intent(in) | trx, |
| integer, intent(in) | itot, | ||
| integer, intent(in) | non, | ||
| real(c_double), dimension(itot), intent(in) | n0, | ||
| real(c_double), dimension(itot), intent(out) | n, | ||
| real(c_double), intent(in) | delta_t, | ||
| real(c_double), dimension(itot), intent(in) | diag, | ||
| real(c_double), dimension(non), intent(in) | a, | ||
| real(c_double), dimension(itot), intent(in) | b, | ||
| logical, intent(in) | zero_flux_step, | ||
| logical, intent(in) | is_adjoint, | ||
| real(c_double), intent(in) | abstol, | ||
| real(c_double), intent(in) | reltol | ||
| ) |
- Parameters
-
[in] trx Library data [in] itot Total number of nuclides [in] non Number of non-zero matrix elements [in] diag Diagonal elements of the transition matrix [in] a Nondiagonal elements of the transition matrix [in] b Right side of the equation (source term) [in] n0 initial concentrations [out] n final concentrations [in] delta_t Time step length in seconds [in] zero_flux_step is this a zero-flux (=decay) step? [in] is_adjoint is this an adjoint calculation? [in] abstol solver options [in] reltol solver options
References fake_flux, fcn(), jac(), radau5(), origen_iofunctions_m::rstop(), solout(), xa, xb, xd, xloc, and xtrx.
Referenced by depletion::depletion_step().
Variable Documentation
| type(origen_transitionmatrixp) xtrx |
Library data (NOTE: THIS CLASS IS NOT THREAD SAFE!)
Referenced by fcn(), and radau5_solve().
| integer, dimension(:), allocatable xloc |
parent locations
Referenced by fcn(), and radau5_solve().
| real(c_double), dimension(:), allocatable xd |
Diagonal elements of the transition matrix.
Referenced by fcn(), and radau5_solve().
| real(c_double), dimension(:), allocatable xa |
Nondiagonal elements of the transition matrix.
Referenced by fcn(), and radau5_solve().
| real(c_double), dimension(:), allocatable xb |
Right side of the equation (source term)
Referenced by fcn(), and radau5_solve().
| real(c_double) fake_flux |
fake flux level to last_transition is decay if=0
Referenced by fcn(), and radau5_solve().
