Различия
Здесь показаны различия между двумя версиями данной страницы.
recalc:ms_tpf_realc0 [2016/05/12 12:58] admin [Source] |
recalc:ms_tpf_realc0 [2018/11/12 15:08] |
||
---|---|---|---|
Строка 1: | Строка 1: | ||
- | ====== Ms: TPF, 1991 ====== | ||
- | Пересчет на 11 атомов кислорода, 4 иона в IV координации. | ||
- | |||
- | ==== Details ==== | ||
- | |||
- | ^ Reference | [[:about_tpf|TPF]] | | ||
- | ^ Input | SiO2_wtp, Al2O3_wtp, TiO2_wtp, Cr2O3_wtp, Fe2O3_wtp, FeO_wtp, MnO_wtp, MgO_wtp, CaO_wtp, Na2O_wtp, K2O_wtp, BaO_wtp, NiO_wtp, ZnO_wtp, BeO_wtp, P2O5_wtp, V2O5_wtp, CoO_wtp, GeO2_wtp, ZrO2_wtp, Cs2O_wtp, PbO_wtp, SO3_wtp, CO2_wtp, F_wtp, Cl_wtp | | ||
- | ^ Output | Fe_p3_pfu, Fe_p2_pfu, Al_c4_pfu, Al_c6_pfu | | ||
- | |||
- | <use:recalc.ms_tpf_realc0> | ||
- | <test:recalc.ms_tpf_realc0> | ||
- | |||
- | <WRAP center round todo 60%> | ||
- | * заголовок | ||
- | * краткое описание | ||
- | * <del>статья</del> | ||
- | * <del>ссылка</del> | ||
- | * <del>карточка</del> | ||
- | * подробное описание | ||
- | * <del>код</del> | ||
- | * тест | ||
- | </WRAP> | ||
- | |||
- | ==== Source ==== | ||
- | <Code linenums lang-r> | ||
- | source('utilities_tpf.r') | ||
- | |||
- | recalc.ms_tpf_realc0 <- function (wtp_data, ...) | ||
- | { | ||
- | rn <- tpf.rename(wtp_data); | ||
- | clk <- tpf.omesod(rn,11,4); | ||
- | return(tpf.unrename(clk,wtp_data$Name)); | ||
- | } | ||
- | |||
- | recalc.ms_tpf_realc0.test <- function () { | ||
- | src <- data.frame(Name='test',SiO2_wtp=c(46.920, 45.230), TiO2_wtp=c(0.500, 0.640), Al2O3_wtp=c(29.430, 30.630), | ||
- | FeO_wtp=c(3.610, 4.220), MnO_wtp=c(0.180, 0.010), MgO_wtp=c(2.280, 1.490), CaO_wtp=c(0.0, 0.270), Na2O_wtp=c(0.640, 0.690), K2O_wtp=c(9.330,9.470)); | ||
- | | ||
- | result <- recalc.ms_tpf_realc0 (src); | ||
- | |||
- | print( result); | ||
- | expect_equal(object = result[1,]$Si_pfu, 3.231, tolerance = 0.001); | ||
- | expect_equal(object = result[2,]$Si_pfu, 3.231, tolerance = 0.001); | ||
- | expect_equal(object = result[1,]$Al_pfu, 2.388, tolerance = 0.001); | ||
- | expect_equal(object = result[2,]$Al_pfu, 2.388, tolerance = 0.001); | ||
- | expect_equal(object = result[1,]$Al_c6_pfu, 1.619, tolerance = 0.001); | ||
- | expect_equal(object = result[2,]$Al_c6_pfu, 1.619, tolerance = 0.001); | ||
- | } | ||
- | |||
- | </Code> | ||
- | |||
- | {{tag> order71}} |