/*$Id: d_cap.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. *------------------------------------------------------------------ * capacitance devices: * self-capacitance (C device) * trans-capacitance (non-spice charge transfer device) *------------------------------------------------------------------ * capacitor models * y.x = volts, y.f0 = coulombs, ev = y.f1 = farads * q = y history in time * i.x = volts, i.f0 = amps, i.f1 = mhos * m.x = volts, m.c0 = amps, acg = m.c1 = mhos */ //testing=script 2006.07.17 #include "l_dispatcher.h" #include "e_storag.h" extern DISPATCHER device_dispatcher; /*--------------------------------------------------------------------------*/ namespace { /*--------------------------------------------------------------------------*/ class DEV_CAPACITANCE : public STORAGE { protected: explicit DEV_CAPACITANCE(const DEV_CAPACITANCE& p) :STORAGE(p) {} public: explicit DEV_CAPACITANCE() :STORAGE() {} protected: // override virtual char id_letter()const {return 'C';} const char* dev_type()const {return "capacitor";} 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_CAPACITANCE(*this);} //void parse_spice(CS&); //ELEMENT //void print(OMSTREAM,LANGUAGE)const; //ELEMENT //void elabo1(); //COMPONENT //void map_nodes(); //ELEMENT //void precalc() //STORAGE void tr_iwant_matrix() {tr_iwant_matrix_passive();} //void dc_begin(); //STORAGE //void tr_begin(); //STORAGE //void tr_restore(); //STORAGE //void dc_advance(); //STORAGE //void tr_advance(); //STORAGE //bool tr_needs_eval(); //STORAGE //void tr_queue_eval(); //ELEMENT bool do_tr(); void tr_load() {tr_load_passive();} //DPAIR tr_review(); //STORAGE //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; void ac_iwant_matrix() {ac_iwant_matrix_passive();} void ac_begin() {_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_TRANSCAP : public DEV_CAPACITANCE { private: explicit DEV_TRANSCAP(const DEV_TRANSCAP& p) :DEV_CAPACITANCE(p){} public: explicit DEV_TRANSCAP() :DEV_CAPACITANCE() {} private: // override virtual char id_letter()const {untested();return '\0';} const char* dev_type()const {return "tcap";} 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 {untested();return true;} bool f_is_value()const {untested();return true;} CARD* clone()const {return new DEV_TRANSCAP(*this);} //void parse_spice(CS&); //ELEMENT //void print(OMSTREAM,LANGUAGE)const; //ELEMENT //void elabo1(); //COMPONENT //void map_nodes(); //ELEMENT //void precalc() //STORAGE void tr_iwant_matrix() {tr_iwant_matrix_active();} //void dc_begin(); //STORAGE //void tr_begin(); //STORAGE //void tr_restore(); //STORAGE //void dc_advance(); //STORAGE //void tr_advance(); //STORAGE //bool tr_needs_eval(); //STORAGE //void tr_queue_eval(); //ELEMENT //bool do_tr(); //DEV_CAPACITANCE void tr_load() {tr_load_active();} //DPAIR tr_review(); //STORAGE //void tr_accept(); //CARD/nothing //void tr_unload(); //DEV_CAPACITANCE 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;//DEV_CAPACITANCE void ac_iwant_matrix() {ac_iwant_matrix_active();} //void ac_begin(); //DEV_CAPACITANCE //void do_ac(); //DEV_CAPACITANCE void ac_load() {untested(); ac_load_active();} //COMPLEX ac_involts()const; //DEV_CAPACITANCE //COMPLEX ac_amps()const; //ELEMENT //XPROBE ac_probe_ext(CS&)const;//ELEMENT }; /*--------------------------------------------------------------------------*/ //BUG// doesn't model dynamic effects of control. class DEV_VCCAP : public DEV_CAPACITANCE { private: explicit DEV_VCCAP(const DEV_VCCAP& p) :DEV_CAPACITANCE(p) {} public: explicit DEV_VCCAP() :DEV_CAPACITANCE() {} private: // override virtual char id_letter()const {untested();return '\0';} const char* dev_type()const {return "vccap";} 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;} bool f_is_value()const {untested();return true;} CARD* clone()const {return new DEV_VCCAP(*this);} //void parse_spice(CS&); //ELEMENT //void print(OMSTREAM,LANGUAGE)const; //ELEMENT //void elabo1(); //COMPONENT //void map_nodes(); //ELEMENT //void precalc() //STORAGE void tr_iwant_matrix() {tr_iwant_matrix_extended();} //void dc_begin(); //STORAGE //void tr_begin(); //STORAGE //void tr_restore(); //STORAGE //void dc_advance(); //STORAGE //void tr_advance(); //STORAGE //bool tr_needs_eval(); //STORAGE //void tr_queue_eval(); //ELEMENT bool do_tr(); //void tr_load(); //DEV_CAPACITANCE //DPAIR tr_review(); //STORAGE //void tr_accept(); //CARD/nothing //void tr_unload(); //DEV_CAPACITANCE 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;//DEV_CAPACITANCE void ac_iwant_matrix() {ac_iwant_matrix_extended();} //void ac_begin(); //DEV_CAPACITANCE //void do_ac(); //DEV_CAPACITANCE //void ac_load(); //DEV_CAPACITANCE //COMPLEX ac_involts()const; //DEV_CAPACITANCE //COMPLEX ac_amps()const; //ELEMENT //XPROBE ac_probe_ext(CS&)const;//ELEMENT }; /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/ bool DEV_CAPACITANCE::do_tr() { if (using_tr_eval()) { _y0.x = tr_input_limited(); tr_eval(); }else{ _y0.x = tr_input(); // tr_involts(); assert(_y0.f1 == value()); _y0.f0 = _y0.x * _y0.f1; assert(converged()); } store_values(); q_load(); _q[0] = _y0; _i[0] = differentiate(); _m0 = CPOLY1(_i[0]); return converged(); } /*--------------------------------------------------------------------------*/ void DEV_CAPACITANCE::do_ac() { if (using_ac_eval()) { ac_eval(); }else{ assert(_ev == _y0.f1); assert(has_tr_eval() || _ev == double(value())); } _acg = _ev * SIM::jomega; ac_load(); } /*--------------------------------------------------------------------------*/ double DEV_CAPACITANCE::tr_probe_num(CS& cmd)const { if (cmd.pmatch("Qcap") || cmd.pmatch("CHarge")) { return _q[0].f0; }else if (cmd.pmatch("Capacitance")) { return _q[0].f1; }else if (cmd.pmatch("DCDT")) {untested(); return (_q[0].f1 - _q[1].f1) / _dt; }else if (cmd.pmatch("DC")) {untested(); return (_q[0].f1 - _q[1].f1); }else if (cmd.pmatch("DQDT")) { return (_q[0].f0 - _q[1].f0) / _dt; }else if (cmd.pmatch("DQ")) { return (_q[0].f0 - _q[1].f0); }else{ return STORAGE::tr_probe_num(cmd); } } /*--------------------------------------------------------------------------*/ bool DEV_VCCAP::do_tr() { _y0.x = tr_input_limited(); tr_eval(); store_values(); q_load(); _q[0].x = tr_outvolts(); _q[0].f1 = _y0.f0; // self capacitance _q[0].f0 = _q[0].x * _q[0].f1; // charge _i[0] = differentiate(); _m0.x = _i[0].x; _m0.c1 = _i[0].f1; _m0.c0 = _i[0].f0 - _i[0].x * _i[0].f1; return converged(); } /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/ DEV_CAPACITANCE p1; DEV_TRANSCAP p2; DEV_VCCAP p3; DISPATCHER::INSTALL d1(&device_dispatcher, "C,capacitor", &p1), d2(&device_dispatcher, "tcap,tcapacitor", &p2), d3(&device_dispatcher, "vccap", &p3); } /*--------------------------------------------------------------------------*/ /*--------------------------------------------------------------------------*/