Различия
Здесь показаны различия между двумя версиями данной страницы.
| Предыдущая версия справа и слева Предыдущая версия Следующая версия | Предыдущая версия | ||
|
sensors:t:cpx_opx_sen_jones_1989 [2017/05/18 14:05] admin |
sensors:t:cpx_opx_sen_jones_1989 [2018/11/12 15:08] (текущий) |
||
|---|---|---|---|
| Строка 20: | Строка 20: | ||
| <WRAP center round todo 60%> | <WRAP center round todo 60%> | ||
| * <del>заголовок</del> | * <del>заголовок</del> | ||
| - | * краткое описание | + | * <del>краткое описание</del> |
| * <del>статья</del> | * <del>статья</del> | ||
| * <del>ссылка</del> | * <del>ссылка</del> | ||
| - | * карточка | + | * <del>карточка</del> |
| * подробное описание | * подробное описание | ||
| * <del>код</del> | * <del>код</del> | ||
| - | * тест | + | * <del>тест</del> |
| </WRAP> | </WRAP> | ||
| Строка 39: | Строка 39: | ||
| if(!usePFU) { | if(!usePFU) { | ||
| - | Opx <- tpf.recalc_minaral("OPX",Opx); | ||
| Cpx <- tpf.recalc_minaral("CPX",Cpx); | Cpx <- tpf.recalc_minaral("CPX",Cpx); | ||
| + | Opx <- tpf.recalc_minaral("OPX",Opx); | ||
| } | } | ||
| Строка 54: | Строка 54: | ||
| B <- B1 / B2 * B3 / B4; | B <- B1 / B2 * B3 / B4; | ||
| T_eq3 <- 4900 / (1.807 - log(A / B)) - 273; | T_eq3 <- 4900 / (1.807 - log(A / B)) - 273; | ||
| - | |||
| A <- (1 - (Cpx$Al_c6_pfu) - Cpx$Ti_pfu - Cpx$Cr_pfu - Cpx$Fe_p3_pfu) * Cpx$Mg_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu) * Cpx$Ca_pfu; | A <- (1 - (Cpx$Al_c6_pfu) - Cpx$Ti_pfu - Cpx$Cr_pfu - Cpx$Fe_p3_pfu) * Cpx$Mg_pfu / (Cpx$Mg_pfu + Cpx$Fe_p2_pfu) * Cpx$Ca_pfu; | ||
| B <- (1 - (Opx$Al_c6_pfu) - Opx$Ti_pfu - Opx$Cr_pfu - Opx$Fe_p3_pfu) * Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu) * Opx$Ca_pfu; | B <- (1 - (Opx$Al_c6_pfu) - Opx$Ti_pfu - Opx$Cr_pfu - Opx$Fe_p3_pfu) * Opx$Mg_pfu / (Opx$Mg_pfu + Opx$Fe_p2_pfu) * Opx$Ca_pfu; | ||
| T_eq4 <- 7045 / (2.47 - log(B / A)) - 273; | T_eq4 <- 7045 / (2.47 - log(B / A)) - 273; | ||
| - | |||
| - | print(row_names); | ||
| - | print(T_eq3); | ||
| - | print(T_eq4); | ||
| return(data.frame(Name=row_names,T_calc_eq3=T_eq3,T_calc_eq4=T_eq4)); | return(data.frame(Name=row_names,T_calc_eq3=T_eq3,T_calc_eq4=T_eq4)); | ||
| Строка 76: | Строка 71: | ||
| sensors.t.cpx_opx_sen_jones_1989.test <- function(){ | sensors.t.cpx_opx_sen_jones_1989.test <- function(){ | ||
| - | td <- data.frame(Name=c('Cpx','Opx'), SiO2_wtp=c(55.65,55.16), TiO2_wtp=c(0.06,0.08), Al2O3_wtp=c(0.59,4.16), Cr2O3_wtp=c(0.08,0.38), | + | cpx <- data.frame(Name=c('Cpx'), SiO2_wtp=c(55.65), TiO2_wtp=c(0.06), Al2O3_wtp=c(0.59), Cr2O3_wtp=c(0.08), |
| - | FeO_wtp=c(1.53,5.00), MgO_wtp=c(19.85,34.87), CaO_wtp=c(22.27,1.21), MnO_wtp=c(0.00,0.11), Na2O_wtp=c(0.26,0.09)); | + | FeO_wtp=c(1.53), MgO_wtp=c(19.85), CaO_wtp=c(22.27), MnO_wtp=c(0.00), Na2O_wtp=c(0.26)); |
| + | |||
| + | opx <- data.frame(Name=c('Opx'), SiO2_wtp=c(55.16), TiO2_wtp=c(0.08), Al2O3_wtp=c(4.16), Cr2O3_wtp=c(0.38), | ||
| + | FeO_wtp=c(5.00), MgO_wtp=c(34.87), CaO_wtp=c(1.21), MnO_wtp=c(0.11), Na2O_wtp=c(0.09)); | ||
| - | res <- sensors.t.cpx_opx_sen_jones_1989(td[1,],td[2,],10); | + | res <- sensors.t.cpx_opx_sen_jones_1989(cpx,opx); |
| | | ||
| - | print(res); | + | expect_equal(object = res[1,]$T_calc_eq3, 1040, tolerance = 1); |
| + | expect_equal(object = res[1,]$T_calc_eq4, 1006, tolerance = 1); | ||
| } | } | ||
| </Code> | </Code> | ||
| - | {{tag> tpf Cpx Opx AP_igroup}} | + | {{tag> tpf Cpx Opx AP_igroup fixed}} |