Различия
Здесь показаны различия между двумя версиями данной страницы.
| Предыдущая версия справа и слева Предыдущая версия Следующая версия | Предыдущая версия | ||
|
interpmin:px_macgreg [2016/04/28 13:53] admin |
interpmin:px_macgreg [2018/11/12 15:07] (текущий) |
||
|---|---|---|---|
| Строка 39: | Строка 39: | ||
| s1 <- d$Si_pfu-2; | s1 <- d$Si_pfu-2; | ||
| s2 <- ifelse(s1<0,0,s1); | s2 <- ifelse(s1<0,0,s1); | ||
| - | M3 <- d$Cr_pfu+2*(s2+d$Ti_pfu)+d$Fe_p3_pfu+d; | + | M3 <- d$Cr_pfu+2*(s2+d$Ti_pfu)+d$Fe_p3_pfu+dd; |
| ee <- M3-(d$Na_pfu+d$K_pfu); | ee <- M3-(d$Na_pfu+d$K_pfu); | ||
| e2 <- ifelse(ee<0,0,ee); | e2 <- ifelse(ee<0,0,ee); | ||
| Строка 52: | Строка 52: | ||
| A3n <- ifelse(A2n<0,0,A2n); | A3n <- ifelse(A2n<0,0,A2n); | ||
| F2n <- ifelse(A2n<0,d$Fe_p3_pfu+A2n,d$Fe_p3_pfu); | F2n <- ifelse(A2n<0,d$Fe_p3_pfu+A2n,d$Fe_p3_pfu); | ||
| - | return(data.frame(e2 = e2)); | ||
| g <- d$Ca_pfu - e2/2 - A1n; | g <- d$Ca_pfu - e2/2 - A1n; | ||
| g2 <- ifelse(g<0,0,g); | g2 <- ifelse(g<0,0,g); | ||
| Строка 81: | Строка 80: | ||
| return(res); | return(res); | ||
| + | } | ||
| + | |||
| + | interpmin.px_macgreg.test <- function () | ||
| + | { | ||
| + | src <- data.frame(Name='test', SiO2_wtp=49.210, TiO2_wtp=0.19, Al2O3_wtp=11.180, Cr2O3_wtp=0.23, | ||
| + | FeO_wtp=6.28, MnO_wtp=0.00, MgO_wtp=19.160, CaO_wtp=16.16, Na2O_wtp=0.00, K2O_wtp=0.00); | ||
| + | |||
| + | result <- interpmin.px_macgreg(src); | ||
| + | | ||
| + | expect_equal(object = result$X_Jd, 0.0, tolerance = 0.05); | ||
| + | expect_equal(object = result$X_DiHd, 0.342, tolerance = 0.05); | ||
| + | expect_equal(object = result$X_Opx, 0.392, tolerance = 0.05); | ||
| + | |||
| } | } | ||
| </Code> | </Code> | ||