Odwrotna notacja polska w praktyce

Nawiązując do artykułu o Odwrotna notacja polska chciałbym przedstawić jak można wykorzystać ją w praktyce, na przykładzie języka C++. Nie będę się ponownie rozwodził jak się to wykonuje, bo zostało to już opisane. Cała esencja polega na tym, że nie musimy wyrażenia arytmetycznego przerabiać najpierw na RPN a później obliczać. Wszystko możemy wykonać za jednym zamachem. Przykładowo mamy równanie:

2*2+2

W programie wykorzystujemy dwa stosy - jeden na liczby (najlepiej zmiennoprzecinkowe), drugi na znaki arytmetyczne. I teraz od lewej bierzemy dwójkę na stos z liczbami, dalej znak mnożenia na stos ze znakami arytmetycznymi, znowu dwa na stos i dalej znak plus. Na stosie ze znakami jest mnożenie, czyli wyższy priorytet od dodawania. Tak więc zdejmujemy mnożenie, wykonujemy na liczbach ze stosu i znowu wynik odkładamy ponownie.
Dalej dwa znowu na stos i koniec. Więc zdejmujemy kolejno wszystkie znaki a że mamy tylko jeden (dodawanie), to sumujemy wynik mnożenia odłożony na stosie z ostatnia dwójką.
Nie rozwlekałem się dokładniej nad tym zadaniem, gdyż opis i stosowaną terminologie można zobaczyć we wspomnianym artykule.

typedef double typ_danej;

typ_danej oblicz(string &s) {
 
 
 typ_danej liczba (0);
 
 char stos_char[count_char];
 char last_;
 int pos_stos_char (0);
 
 typ_danej stos_liczba[count_char * 2];
 int pos_stos_liczba (0);
 
 while (s.length()>0) {
            
   if (isdigit(s[0])) {
     
     liczba = atof(s.c_str());
     
     while (isdigit(s[0]) || s[0]=='.') s.erase(0,1);  
     
     stos_liczba[pos_stos_liczba] = liczba;
     pos_stos_liczba++;
                  
   } else {
     
     if (s[0]!=')') {
       if (pos_stos_char>0) {
         last_ = stos_char[pos_stos_char - 1];                   
       } else last_ = '#';
       } else last_ = ')';    
     
     while(priorytet(s[0],last_)==0) {
       
       if (przelicz(last_,&stos_liczba[pos_stos_liczba],pos_stos_liczba)==1) {
         DivideByZero = true;
         return 0;
       }                                         
       pos_stos_char--;
        
       if (pos_stos_char>0) {
         last_ = stos_char[pos_stos_char - 1];                   
       } else last_ = '#';                               
     }                                                  
     
     if (last_!=')') {
       stos_char[pos_stos_char] = s[0];
       pos_stos_char++;
     } else {
            
      while (stos_char[pos_stos_char-1]!='(') {
            
        if (przelicz(stos_char[pos_stos_char-1], &stos_liczba[pos_stos_liczba],pos_stos_liczba)==1) {
          DivideByZero = true;
          return 0;
        }                                                                                             
        pos_stos_char--; 
             
      }
      pos_stos_char--;
     }       
          
     s.erase(0,1);      
   
   }             
 }
 

 while (pos_stos_char>0) {
 
   if (przelicz(stos_char[pos_stos_char-1], &stos_liczba[pos_stos_liczba],pos_stos_liczba)==1) {
     DivideByZero = true;
     return 0;
   }                                                                                                
   pos_stos_char--;
       
 }

 return stos_liczba[pos_stos_liczba-1];
}

Argumentem funkcji jak widać jest dana typu string, w której opisane jest równanie. Są oczywiście odnośniki do innych funkcji, które można przeanalizować ściągając całość (link niżej).
Jest to dość prosty kalkulator, gdyż operuje na następnych znakach arytmetycznych: + - * / ^{(potęga, która też może być pierwiastkowaniem - np. pierwiastek trzeciego stopnia z 27 zapiszemy: 27}(1/3), czyli 27 do potęgi jednej trzeciej) oraz nawiasy.
W załączniku dodałem kompletny kod źródłowy, gdzie dodatkowo jest funkcja do weryfikacji równania, wraz z kodem wynikowym w postaci pliku wykonywalnego. Oparte jest to na konsoli a sam program pisałem pod kompilator Dev-C++. Kod można ściągnąć stąd.