/*$Id: d_admit.cc,v 26.14 2007/02/07 09:06:48 al Exp $ -*- C++ -*- * Copyright (C) 2001 Albert Davis * Author: Albert Davis * * This file is part of "Gnucap", the Gnu Circuit Analysis Package * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by * the Free Software Foundation; either version 2, or (at your option) * any later version. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA * 02110-1301, USA. *------------------------------------------------------------------ * admittance devices: * self-admittance (old Y device) * y.x = volts, y.f0 = amps, ev = y.f1 = mhos. * m.x = volts, m.c0 = amps, acg = m.c1 = mhos. * trans-admittance (VCCS, G device) * voltage controlled admittance * y.x = volts(control), y.f0 = mhos, ev = y.f1 = mhos/volt * m.x = volts(control), m.c0 = 0, acg = m.c1 = mhos * _loss0 == 1/R. (mhos) */ //testing=script 2006.07.17 #include "l_dispatcher.h" #include "e_elemnt.h" extern DISPATCHER device_dispatcher; /*--------------------------------------------------------------------------*/ namespace { /*--------------------------------------------------------------------------*/ class DEV_ADMITTANCE : public ELEMENT { protected: explicit DEV_ADMITTANCE(const DEV_ADMITTANCE& p) :ELEMENT(p) {} public: explicit DEV_ADMITTANCE() :ELEMENT() {} protected: // override virtual char id_letter()const {return 'Y';} const char* dev_type()const {return "admittance";} int max_nodes()const {return 2;} int min_nodes()const {return 2;} int out_nodes()const {untested();return 2;} int matrix_nodes()const {return 2;} int net_nodes()const {return 2;} bool is_1port()const {return true;} CARD* clone()const {return new DEV_ADMITTANCE(*this);} //void parse_spice(CS&); //ELEMENT //void print(OMSTREAM,LANGUAGE)const; //ELEMENT //void elabo1(); //COMPONENT //void map_nodes(); //ELEMENT void precalc(); void tr_iwant_matrix() {tr_iwant_matrix_passive();} //void dc_begin(); //CARD/nothing //void tr_begin(); //CARD/nothing //void tr_restore(); //CARD/nothing //void dc_advance(); //CARD/nothing //void tr_advance(); //CARD/nothing //bool tr_needs_eval(); //ELEMENT //void tr_queue_eval(); //ELEMENT bool do_tr(); void tr_load() {tr_load_passive();} //DPAIR tr_review(); //CARD/nothing //void tr_accept(); //CARD/nothing void tr_unload() {tr_unload_passive();} double tr_involts()const {return tr_outvolts();} //double tr_input()const //ELEMENT double tr_involts_limited()const {return tr_outvolts_limited();} //double tr_input_limited()const //ELEMENT //double tr_amps()const; //ELEMENT //double tr_probe_num(CS&)const;//ELEMENT void ac_iwant_matrix() {ac_iwant_matrix_passive();} void ac_begin() {_acg = _ev = _y0.f1;} void do_ac(); void ac_load() {ac_load_passive();} COMPLEX ac_involts()const {untested();return ac_outvolts();} //COMPLEX ac_amps()const; //ELEMENT //XPROBE ac_probe_ext(CS&)const;//ELEMENT }; /*--------------------------------------------------------------------------*/ class DEV_VCCS : public DEV_ADMITTANCE { protected: explicit DEV_VCCS(const DEV_VCCS& p) :DEV_ADMITTANCE(p) {} public: explicit DEV_VCCS() :DEV_ADMITTANCE() {} protected: // override virtual char id_letter()const {return 'G';} const char* dev_type()const {return "vccs";} int max_nodes()const {return 4;} int min_nodes()const {return 4;} int out_nodes()const {untested();return 2;} int matrix_nodes()const {return 4;} int net_nodes()const {return 4;} bool is_1port()const {untested();return false;} bool is_2port()const {return true;} CARD* clone()const {return new DEV_VCCS(*this);} //void parse_spice(CS&); //ELEMENT //void print(OMSTREAM,LANGUAGE)const; //ELEMENT //void elabo1(); //COMPONENT //void map_nodes(); //ELEMENT //void precalc(); //DEV_ADMITTANCE void tr_iwant_matrix() {tr_iwant_matrix_active();} //void dc_begin(); //CARD/nothing //void tr_begin(); //CARD/nothing //void tr_restore(); //CARD/nothing //void dc_advance(); //CARD/nothing //void tr_advance(); //CARD/nothing //bool tr_needs_eval(); //ELEMENT //void tr_queue_eval(); //ELEMENT //bool do_tr() //DEV_ADMITTANCE void tr_load() {tr_load_active();} //DPAIR tr_review(); //CARD/nothing //void tr_accept(); //CARD/nothing void tr_unload() {untested();tr_unload_active();} double tr_involts()const {return dn_diff(_n[IN1].v0(), _n[IN2].v0());} //double tr_input()const //ELEMENT double tr_involts_limited()const {return volts_limited(_n[IN1],_n[IN2]);} //double tr_input_limited()const //ELEMENT //double tr_amps()const //ELEMENT //double tr_probe_num(CS&)const;//ELEMENT void ac_iwant_matrix() {ac_iwant_matrix_active();} //void ac_begin(); //DEV_ADMITTANCE //void do_ac(); //DEV_ADMITTANCE void ac_load() {ac_load_active();} COMPLEX ac_involts()const {return _n[IN1]->vac() - _n[IN2]->vac();} //COMPLEX ac_amps()const; //ELEMENT //XPROBE ac_probe_ext(CS&)const;//ELEMENT }; /*--------------------------------------------------------------------------*/ class DEV_VCG : public DEV_VCCS { private: explicit DEV_VCG(const DEV_VCG& p) :DEV_VCCS(p) {} public: explicit DEV_VCG() :DEV_VCCS() {} private: // override virtual char id_letter()const {untested();return '\0';} const char* dev_type()const {return "vcg";} int max_nodes()const {return 4;} int min_nodes()const {return 4;} int out_nodes()const {untested();return 2;} int matrix_nodes()const {return 4;} int net_nodes()const {return 4;} bool is_2port()const {untested();return true;} CARD* clone()const {return new DEV_VCG(*this);} //void parse_spice(CS&); //ELEMENT //void print(OMSTREAM,LANGUAGE)const; //ELEMENT //void elabo1(); //COMPONENT //void map_nodes(); //ELEMENT void precalc(); void tr_iwant_matrix() {tr_iwant_matrix_extended();} //void dc_begin(); //CARD/nothing //void tr_begin(); //CARD/nothing //void tr_restore(); //CARD/nothing //void dc_advance(); //CARD/nothing //void tr_advance(); //CARD/nothing //bool tr_needs_eval(); //ELEMENT //void tr_queue_eval(); //ELEMENT bool do_tr(); void tr_load() {tr_load_shunt(); tr_load_active();} //DPAIR tr_review(); //CARD/nothing //void tr_accept(); //CARD/nothing void tr_unload() {untested(); tr_unload_shunt(); tr_unload_active();} //double tr_involts()const //DEV_VCCS //double tr_input()const //ELEMENT //double tr_involts_limited()const //DEV_VCCS //double tr_input_limited()const //ELEMENT //double tr_amps()const //ELEMENT //double tr_probe_num(CS&)const;//ELEMENT void ac_iwant_matrix() {ac_iwant_matrix_extended();} void ac_begin() {_ev = _y0.f0; _acg = _m0.c1;} void do_ac(); void ac_load() {ac_load_loss(); ac_load_active();} //COMPLEX ac_involts()const //DEV_VCCS //COMPLEX ac_amps()const; //ELEMENT //XPROBE ac_probe_ext(CS&)const;//ELEMENT }; /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/ void DEV_ADMITTANCE::precalc() { _y0.f0 = LINEAR; _y0.f1 = value(); _y1 = _y0; _m0.c1 = _y0.f1; _m0.c0 = 0.; _m1 = _m0; assert(_loss0 == 0.); assert(_loss1 == 0.); set_constant(!has_tr_eval()); set_converged(!has_tr_eval()); } /*--------------------------------------------------------------------------*/ bool DEV_ADMITTANCE::do_tr() { if (using_tr_eval()) { _y0.x = _m0.x = tr_involts_limited(); // _y0.x = tr_input_limited(); //assert(_y0.x == _m0.x); _y0.f0 = _m0.c1 * _m0.x + _m0.c0; //BUG// patch for diode tr_eval(); assert(_y0.f0 != LINEAR); store_values(); q_load(); _m0 = CPOLY1(_y0); }else{ assert(_y0.f0 == LINEAR); assert(_y0.f1 == value()); assert(_m0.c1 == _y0.f1); assert(_m0.c0 == 0.); assert(_y1 == _y0); assert(converged()); } return converged(); } /*--------------------------------------------------------------------------*/ void DEV_ADMITTANCE::do_ac() { if (using_ac_eval()) {untested(); ac_eval(); _acg = _ev; }else{ assert(_ev == _y0.f1); assert(has_tr_eval() || _ev == double(value())); } assert(_acg == _ev); ac_load(); } /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/ void DEV_VCG::precalc() { _loss1 = _loss0 = 0; _y0.f0 = 0; _y0.f1 = value(); _y1 = _y0; _m0.c1 = 0.; _m0.c0 = 0.; _m1 = _m0; assert(!is_constant()); set_not_converged(); } /*--------------------------------------------------------------------------*/ bool DEV_VCG::do_tr() { _y0.x = tr_input_limited(); tr_eval(); assert(_y0.f0 != LINEAR); store_values(); q_load(); _loss0 = _y0.f0; _m0.x = tr_outvolts(); // fake _m0.c1 = _y0.f1 * tr_outvolts(); _m0.c0 = -_y0.x * _m0.c1; return converged(); } /*--------------------------------------------------------------------------*/ void DEV_VCG::do_ac() { if (using_ac_eval()) { ac_eval(); _acg = _ev * _m0.x; _ev *= _y0.x; }else{ assert(_ev == _y0.f0); assert(_acg == _m0.c1); } ac_load(); } /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/ DEV_ADMITTANCE p1; DEV_VCCS p2; DEV_VCG p3; DISPATCHER::INSTALL d1(&device_dispatcher, "Y,admittance", &p1), d2(&device_dispatcher, "G,vccs", &p2), d3(&device_dispatcher, "vcg", &p3); } /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/