Solving Ordinary Differential Equations to Understand Luminescence

A collection of functions to simulate luminescence signals in quartz and Al2O3 based on published models.


News

NEWS for the R Package RLumModel

Changes in version 0.2.3 (22nd November, 2017)

Miscellaneous

  • Add further tests

Bugfixes and changes:

  • set_pars(): * Add new parameter set "Friedrich2018"

Changes in version 0.2.2 (13th October, 2017)

Miscellaneous

  • Update error messages and tests
  • Add vignette "RLumModel - Simulating luminescence in Al2O3"
  • Change vignettes from PDF to HTML
    • Change vignette builder from r.rsp to knitr
    • Add kableExtra to SUGGESTS

Changes in version 0.2.1 (13th April, 2017)

Miscellaneous

  • Fix CRAN NOTE for not declared entry points.
  • Expand tests.

Bugfixes and changes:

  • simulate_LM_OSL():
    • Right call to .set_ODE_Rcpp_LM_OSL

Changes in version 0.2.0 (20th January, 2017)

New functions:

  • set_ODE_Rcpp() & set_ODE_Rcpp():

    • Changed the main calculations of the ODEs to C++ via Rcpp. This results in an extreme fast calculation of ODE calcualtions. Slow calculations like DRT or SAR can be solved much faster now.
  • simulate_RF_and_heating():

    • Simulating the process of heating/cooling a sample and detect luminescence during irradiation.

Bugfixes and changes:

  • calc_concentrations():

    • Add maximum change of concentration in specific level to slot 'info' in every RLum.Data.Curve.
  • calc_signal():

    • As a result of the changes in 'set_pars' the index of the luminescence center changed. This was adjusted in 'calc_signal'.
  • model_LuminescenceSignals():

    • Add new parameter 'own_parameters': Parameter to submit an own parameter set for simulations. For further details see vignette 'RLumModel - Using own parameter sets'.
      • Add new parameter 'own_state_parameters': Parameter to submit own state parameters together with an own parameter set.
      • Add new parameter 'own_start_temperature': Parameter to submit own start temperature together with an own parameter set.
  • set_pars():

    • Changed order in the vectors N, E, s, A, B: The luminescence center is now always the last entry. This was necessary to use OTOR models with own parameters much easier.
    • Output is now a list and not a RLum.Results object.
      • Add new parameter set "Friedrich2017".
      • Add new case 'customized' for own parameter sets.
  • simulate_CW_OSL():

    • Remove 'RLumModel_ID' from slot info to slot .pid. This RLumModel_ID describes the order of the command in the sequence.
  • simulate_irradiation():

    • New distinction of cases adding the possibility to use own pair-production-rates.
      • argument 'dose = 0' is now possible without error (#4, thanks to Christoph Burow)
  • simulate_LM_OSL(): * Remove 'RLumModel_ID' from slot info to slot .pid. This RLumModel_ID describes the order of the command in the sequence.

  • simulate_RF(): * Remove 'RLumModel_ID' from slot 'info' to slot .pid. This RLumModel_ID describes the order of the command in the sequence. * New distinction of cases adding the possibility to use own pair-production-rates.

  • simulate_TL(): * Remove 'RLumModel_ID' from slot info to slot .pid. This RLumModel_ID describes the order of the command in the sequence.

  • tests: * Expanded tests

  • translate_sequence():

    • Hardcoded originator from returning object to "model_LuminescenceSignals()"
      • Add 'sequence' to slot 'info' in returning RLum.Analysis object.
      • Add 'originators' to slot 'info' in returning RLum.Analysis object.

Changes in version 0.1.2 (2nd September, 2016)

Bugfixes and changes:

  • Added vignette:

    • Added vignette 'Getting started with RLumModel'.
  • set_ODE():

    • Remove 'b', 'temp', 'R' and 'P' from unpacking, because these parameters are already stored in parameters.step.
      • Remove distinction of cases 'model = Pagonis' or 'model = Bailey' due to a misstake in Bailey 2001, 2002 and 2004. The rate equation for the valence band is wrong. Thanks to an unknown reviewer of Friedrich et al. 2016.
  • set_ODE_LM_OSL():

    • Remove 'a', 'b', 'temp', 'R' and 'P' from unpacking, because these parameters are already stored in parameters.step.
      • Remove distinction of cases 'model = Pagonis' or 'model = Bailey' due to a misstake in Bailey 2001, 2002 and 2004. The rate equation for the valence band is wrong. Thanks to an unknown reviewer of Friedrich et al. 2016.
  • simulate_CW_OSL():

    • Changed ode solver to 'ode' with method 'bdf' for faster solving long illumination times.
  • simulate_illumination(): * Changed ode solver to 'ode' with method 'bdf' for faster solving long illumination times.

  • simulate_irradiation():

    • Changed choice of pair-production-rate 'R', so that Bailey 2004 has the right value of 2.5e10.
  • simulate_LM_OSL(): * Changed ode solver to 'ode' with method 'bdf' for faster solving long illumination times.

  • simulate_RF(): * Changed choice of pair-production-rate 'R', so that Bailey 2004 has the right value of 2.5e10.


Changes in version 0.1.1 (2nd May, 2016)

  • Changed tests because of new package structure of 'testthat'.

Changes in version 0.1.0 (3rd February, 2016)

  • Initial version

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.

install.packages("RLumModel")

0.2.3 by Johannes Friedrich, a year ago


https://CRAN.R-project.org/package=RLumModel


Browse source code at https://github.com/cran/RLumModel


Authors: Johannes Friedrich [aut, trl, cre] , Sebastian Kreutzer [aut, ths] , Christoph Schmidt [aut, ths]


Documentation:   PDF Manual  


GPL-3 license


Imports deSolve, methods, Rcpp

Depends on utils, Luminescence

Suggests knitr, kableExtra, testthat

Linking to Rcpp, RcppArmadillo


Suggested by Luminescence.


See at CRAN