CCircuitSolver.hh

Go to the documentation of this file.

#ifndef CCIRCUITSOLVER_HH
#define CCIRCUITSOLVER_HH

#include <math/Complex.hh>
#include <math/CMatrix.hh>
#include <electro/CElement.hh>
#include <CCircuit.hh>
#include <CRlcException.hh>


  class CUsedList
  { 
    public : 
      ~CUsedList(){}
      void use(CNode *x) {m_used.push_back(x);}
      bool was_used(CNode *x){
         for(int i=0;i<m_used.size();i++)
             if(m_used[i]==x) return true;
         return false;
      }
    protected:
      std::vector<CNode *> m_used; 
  };

class CCircuitSolver
{
public :
  CCircuitSolver(CCircuit * data=NULL) {m_data=data;}
  ~CCircuitSolver();
          
  void solve();
  void clear();
  void load(CCircuit * data) {m_data=data;}
  
protected:  
  
  CCircuit *m_data;

  CUsedList used;
  
  CElement * get_item(std::string s);
          
  std::vector<CExNode *> exnodes;
  std::vector<CBranch *> branches;
  
  void zciagaj();
  void serch_nodes(CExNode *des, CNode *start);
  void print_aliases();
  void make_branches();
  
  //omegi wystepujaca w calym obwodzie
  std::vector<double> m_omegas;
  void found_omegas();
  void print_omegas();
  
  //macierze topologiczne
  CMatrix<Complex> *A;
  CMatrix<Complex> *Yb;
  CMatrix<Complex> *Yn;
  CMatrix<Complex> *E;
  CMatrix<Complex> *I;
  CMatrix<Complex> *In;
  CMatrix<Complex> *Un;
  
  //funkcje tworzace macierze
  void makeA();
  void makeYb(double omega=0);
  void makeYn(double omega=0);
  void makeE(double omega=0);
  void makeI(double omega=0);
  void makeIn();
  void makeUn();

};

#endif // CCIRCUITSOLVER_HH

---