Различия

Здесь показаны различия между двумя версиями данной страницы.

--- sensors:p:amp_simakin [2017/04/02 16:41]
lesnik56
+++ sensors:p:amp_simakin [2018/11/16 13:53] (текущий)
admin
@@ Строка 45: / Строка 45: @@
     comm <- 'norm'
     #нужно прокомментировать, что при Рmax<2 расплав был кислым или амфибол заместил другую фазу при реакции с расплавом 'silicious or replace'
-    if(Pmax<1.5){Pmax<-0; comm <-'silicious or replace'}
+    if(Pmax<1.5){Pmax<-0; comm <-'silicious or replacement'}
     #нужно прокомментировать, что  среди них есть закалочный кристалл или метаморфический 'quenching or metamorphic'
     if(Pmax>15){Pmax<-100; comm <-'quenching or metamorphic'}
@@ Строка 51: / Строка 51: @@
     if(Pmax==100){Pmean<-0}
    n1 <- nrow(Amp)
+   ij <- 1:n1
 #  Упорядочение значений "давления" во всех точках состава
 for (j in 1:n1) {
-  for (i in 1:n1-1) {
+  for (i in 1:(n1-1)) {
-   if (P[i+1]>P[i]){z<-P[i]; P[i]<-P[i+1]; P[i+1]<-z;}
+   if (P[i+1]>P[i]){z<-P[i]; P[i]<-P[i+1]; P[i+1]<-z;
+#   nam<-Amp$Name[i]; Amp$Name[i]<-Amp$Name[i+1]; Amp$Name[i+1]<-nam; }
+   li<-Amp[i, ]; Amp[i, ]<-Amp[i+1, ]; Amp[i+1, ]<-li; }
        }};
-# j=1;   for {j 1:nrow(Amp)} {while (P[1]-P[j]<1.0}{j=j+1}}
-#  k=1;   while (k<=j}{k=k+1; sum=sum+P[k]/j}
+# j<-1;   for(j in 1:n1) {while (P[1]-P[j]<1.0}{break}
+#  k<-1;   while(k<=j){sum<-sum+P[k]/j}
 df1<-data.frame(Name='Summary', P_max=Pmax,P_mean=Pmean,comment=comm,n_points=n1,end=' ');
-df2<-data.frame(Name=Amp$Name, P_max=P,P_mean=' ',comment=' ',n_points=' ',end=' ');
+df2<-data.frame(Name=Amp$Name, P_max=P,P_mean=0.0,comment=' ',n_points=' ',end=' ');
-return (rbind(df1, df2));
+df0<-data.frame(Name='_', P_max=0.0,P_mean=0.0,comment='_',n_points=' ',end=' ');
-#return (data.frame(Name='Summary', P_max=Pmax,P_mean=Pmean,comment=comm,n_points=n1,end=' '));
+return (rbind(df1,df0,df2));
  # return (data.frame(Name=Amp$Name, P_max=replicate(nrow(Amp), Pmax), P_mean=replicate(nrow(Amp), Pmean)));
  }
@@ Строка 88: / Строка 95: @@
                            CaO_wtp=c(11.490, 11.490), Na2O_wtp=c(1.870, 1.870),
                            K2O_wtp=c(1.060, 1.060));
     #Расчет
     result <- sensors.p.amp_simakin(test_data);