Changeset 2265 in Sophya for trunk/SophyaPI/PIext/piacmd.cc

Timestamp:

Nov 15, 2002, 12:58:41 AM (23 years ago)

Author:

ansari

Message:

Ajout addarca et addfarca: Trace d'arc defini par un angle et dangle + sum/product/mean/sigma ds RPNEvaluator - Reza 14/11/2002

File:

• trunk/SophyaPI/PIext/piacmd.cc (modified) (2 diffs)

trunk/SophyaPI/PIext/piacmd.cc

@@ -1233 +1233 @@
 } 
 
+int Sum_RPNStack_(stack<double>& s, double& sx, double& sx2, string& line) 
+{
+  sx = sx2 = 0.;
+  int nn = 0;
+  double x = 0.;
+  while( !s.empty() ) {
+    x = s.top(); s.pop();
+    sx += x; sx2 += x*x;
+    nn++;
+  } 
+  return(nn);  
+}
+
+int Product_RPNStack_(stack<double>& s, double& px, string& line) 
+{
+  px = 1.;
+  int nn = 0;
+  double x = 0.;
+  while( !s.empty() ) {
+    x = s.top(); s.pop();
+    px *= x;  nn++;
+  } 
+  return(nn);  
+}
   
 /* --Methode-- */
@@ -1329 +1353 @@
       if ( Check_myRPNStack_(rpnstack, x, y, line) ) return(1);
       rpnstack.top() = atan2(x,y);
+    }
+    // Fonction a N arguments  - Somme, produit, etc ...
+    else if ((args[k] == "sum") || (args[k] == "mean") || 
+             (args[k] == "sigma") || (args[k] == "sigma2") ) {
+      double sx, sx2;
+      int nn = Sum_RPNStack_(rpnstack, sx, sx2, line);
+      if (args[k] == "sum") rpnstack.push(sx);
+      else { 
+        if (nn == 0) {
+          cerr << "PIACmd::EvalRPNExpr: mean/sigma error- RPN stack empty: " 
+               << line << endl; 
+          return(1);
+        }
+        double fnn = nn;
+        if (args[k] == "mean")  rpnstack.push(sx/fnn);
+        else if (args[k] == "sigma2") rpnstack.push(sx2/fnn-(x*x/(fnn*fnn)));
+        else rpnstack.push(sqrt(sx2/fnn-(x*x/(fnn*fnn))));
+      }
+    }
+    else if (args[k] == "product") { 
+      double px;
+      int nn = Product_RPNStack_(rpnstack, px, line);
+      if (nn == 0) {
+        cerr << "PIACmd::EvalRPNExpr: product error- RPN stack empty: " 
+             << line << endl; 
+        return(1);
+      }
+      rpnstack.push(px);
     }
     // Fonctions de manipulation de stack