LCOV - code coverage report
Current view: top level - src/components/verilog - vcresistor.core.cpp (source / functions) Hit Total Coverage
Test: qucs-core-0.0.19 Code Coverage Lines: 0 164 0.0 %
Date: 2015-01-05 16:01:02 Functions: 0 28 0.0 %
Legend: Lines: hit not hit | Branches: + taken - not taken # not executed Branches: 0 198 0.0 %

           Branch data     Line data    Source code
       1                 :            : /*
       2                 :            :  * vcresistor.core.cpp - device implementations for vcresistor module
       3                 :            :  *
       4                 :            :  * This is free software; you can redistribute it and/or modify
       5                 :            :  * it under the terms of the GNU General Public License as published by
       6                 :            :  * the Free Software Foundation; either version 2, or (at your option)
       7                 :            :  * any later version.
       8                 :            :  *
       9                 :            :  */
      10                 :            : 
      11                 :            : #if HAVE_CONFIG_H
      12                 :            : #include <config.h>
      13                 :            : #endif
      14                 :            : 
      15                 :            : #include "vcresistor.analogfunction.h"
      16                 :            : #include "component.h"
      17                 :            : #include "device.h"
      18                 :            : #include "vcresistor.core.h"
      19                 :            : 
      20                 :            : #ifndef CIR_vcresistor
      21                 :            : #define CIR_vcresistor -1
      22                 :            : #endif
      23                 :            : 
      24                 :            : // external nodes
      25                 :            : #define C1 0
      26                 :            : #define C2 1
      27                 :            : #define P1 2
      28                 :            : #define P2 3
      29                 :            : // internal nodes
      30                 :            : 
      31                 :            : // useful macro definitions
      32                 :            : #define NP(node) real (getV (node))
      33                 :            : #define BP(pnode,nnode) (NP(pnode) - NP(nnode))
      34                 :            : #define _load_static_residual2(pnode,nnode,current)\
      35                 :            :         _rhs[pnode] -= current;\
      36                 :            :         _rhs[nnode] += current;
      37                 :            : #define _load_static_augmented_residual2(pnode,nnode,current)\
      38                 :            :         _rhs[pnode] -= current;\
      39                 :            :         _rhs[nnode] += current;
      40                 :            : #define _load_static_residual1(node,current)\
      41                 :            :         _rhs[node] -= current;
      42                 :            : #define _load_static_augmented_residual1(node,current)\
      43                 :            :         _rhs[node] -= current;
      44                 :            : #define _load_static_jacobian4(pnode,nnode,vpnode,vnnode,conductance)\
      45                 :            :         _jstat[pnode][vpnode] += conductance;\
      46                 :            :         _jstat[nnode][vnnode] += conductance;\
      47                 :            :         _jstat[pnode][vnnode] -= conductance;\
      48                 :            :         _jstat[nnode][vpnode] -= conductance;\
      49                 :            :         if (doHB) {\
      50                 :            :         _ghs[pnode] += conductance * BP(vpnode,vnnode);\
      51                 :            :         _ghs[nnode] -= conductance * BP(vpnode,vnnode);\
      52                 :            :         } else {\
      53                 :            :         _rhs[pnode] += conductance * BP(vpnode,vnnode);\
      54                 :            :         _rhs[nnode] -= conductance * BP(vpnode,vnnode);\
      55                 :            :         }
      56                 :            : #define _load_static_jacobian2p(node,vpnode,vnnode,conductance)\
      57                 :            :         _jstat[node][vpnode] += conductance;\
      58                 :            :         _jstat[node][vnnode] -= conductance;\
      59                 :            :         if (doHB) {\
      60                 :            :         _ghs[node] += conductance * BP(vpnode,vnnode);\
      61                 :            :         } else {\
      62                 :            :         _rhs[node] += conductance * BP(vpnode,vnnode);\
      63                 :            :         }
      64                 :            : #define _load_static_jacobian2s(pnode,nnode,node,conductance)\
      65                 :            :         _jstat[pnode][node] += conductance;\
      66                 :            :         _jstat[nnode][node] -= conductance;\
      67                 :            :         if (doHB) {\
      68                 :            :         _ghs[pnode] += conductance * NP(node);\
      69                 :            :         _ghs[nnode] -= conductance * NP(node);\
      70                 :            :         } else {\
      71                 :            :         _rhs[pnode] += conductance * NP(node);\
      72                 :            :         _rhs[nnode] -= conductance * NP(node);\
      73                 :            :         }
      74                 :            : #define _load_static_jacobian1(node,vnode,conductance)\
      75                 :            :         _jstat[node][vnode] += conductance;\
      76                 :            :         if (doHB) {\
      77                 :            :         _ghs[node] += conductance * NP(vnode);\
      78                 :            :         } else {\
      79                 :            :         _rhs[node] += conductance * NP(vnode);\
      80                 :            :         }
      81                 :            : #define _load_dynamic_residual2(pnode,nnode,charge)\
      82                 :            :         if (doTR) _charges[pnode][nnode] += charge;\
      83                 :            :         if (doHB) {\
      84                 :            :         _qhs[pnode] -= charge;\
      85                 :            :         _qhs[nnode] += charge;\
      86                 :            :         }
      87                 :            : #define _load_dynamic_residual1(node,charge)\
      88                 :            :         if (doTR) _charges[node][node] += charge;\
      89                 :            :         if (doHB) {\
      90                 :            :         _qhs[node] -= charge;\
      91                 :            :         }
      92                 :            : #define _load_dynamic_jacobian4(pnode,nnode,vpnode,vnnode,capacitance)\
      93                 :            :         if (doAC) {\
      94                 :            :         _jdyna[pnode][vpnode] += capacitance;\
      95                 :            :         _jdyna[nnode][vnnode] += capacitance;\
      96                 :            :         _jdyna[pnode][vnnode] -= capacitance;\
      97                 :            :         _jdyna[nnode][vpnode] -= capacitance;\
      98                 :            :         }\
      99                 :            :         if (doTR) {\
     100                 :            :         _caps[pnode][nnode][vpnode][vnnode] += capacitance;\
     101                 :            :         }\
     102                 :            :         if (doHB) {\
     103                 :            :         _chs[pnode] += capacitance * BP(vpnode,vnnode);\
     104                 :            :         _chs[nnode] -= capacitance * BP(vpnode,vnnode);\
     105                 :            :         }
     106                 :            : #define _load_dynamic_jacobian2s(pnode,nnode,vnode,capacitance)\
     107                 :            :         if (doAC) {\
     108                 :            :         _jdyna[pnode][vnode] += capacitance;\
     109                 :            :         _jdyna[nnode][vnode] -= capacitance;\
     110                 :            :         }\
     111                 :            :         if (doTR) {\
     112                 :            :         _caps[pnode][nnode][vnode][vnode] += capacitance;\
     113                 :            :         }\
     114                 :            :         if (doHB) {\
     115                 :            :         _chs[pnode] += capacitance * NP(vnode);\
     116                 :            :         _chs[nnode] -= capacitance * NP(vnode);\
     117                 :            :         }
     118                 :            : #define _load_dynamic_jacobian2p(node,vpnode,vnnode,capacitance)\
     119                 :            :         if (doAC) {\
     120                 :            :         _jdyna[node][vpnode] += capacitance;\
     121                 :            :         _jdyna[node][vnnode] -= capacitance;\
     122                 :            :         }\
     123                 :            :         if (doTR) {\
     124                 :            :         _caps[node][node][vpnode][vnnode] += capacitance;\
     125                 :            :         }\
     126                 :            :         if (doHB) {\
     127                 :            :         _chs[node] += capacitance * BP(vpnode,vnnode);\
     128                 :            :         }
     129                 :            : #define _load_dynamic_jacobian1(node,vnode,capacitance)\
     130                 :            :         if (doAC) {\
     131                 :            :         _jdyna[node][vnode] += capacitance;\
     132                 :            :         }\
     133                 :            :         if (doTR) {\
     134                 :            :         _caps[node][node][vnode][vnode] += capacitance;\
     135                 :            :         }\
     136                 :            :         if (doHB) {\
     137                 :            :         _chs[node] += capacitance * NP(vnode);\
     138                 :            :         }
     139                 :            : 
     140                 :            : #define _save_whitenoise1(n1,pwr,type)\
     141                 :            :         _white_pwr[n1][n1] += pwr;
     142                 :            : #define _save_whitenoise2(n1,n2,pwr,type)\
     143                 :            :         _white_pwr[n1][n2] += pwr;
     144                 :            : #define _save_flickernoise1(n1,pwr,exp,type)\
     145                 :            :         _flicker_pwr[n1][n1] += pwr;\
     146                 :            :         _flicker_exp[n1][n1] += exp;
     147                 :            : #define _save_flickernoise2(n1,n2,pwr,exp,type)\
     148                 :            :         _flicker_pwr[n1][n2] += pwr;\
     149                 :            :         _flicker_exp[n1][n2] += exp;
     150                 :            : #define _load_whitenoise2(n1,n2,pwr)\
     151                 :            :         cy (n1,n2) -= pwr/kB/T0; cy (n2,n1) -= pwr/kB/T0;\
     152                 :            :         cy (n1,n1) += pwr/kB/T0; cy (n2,n2) += pwr/kB/T0;
     153                 :            : #define _load_whitenoise1(n1,pwr)\
     154                 :            :         cy (n1,n1) += pwr/kB/T0;
     155                 :            : #define _load_flickernoise2(n1,n2,pwr,exp)\
     156                 :            :         cy (n1,n2) -= pwr*pow(_freq,-exp)/kB/T0;\
     157                 :            :         cy (n2,n1) -= pwr*pow(_freq,-exp)/kB/T0;\
     158                 :            :         cy (n1,n1) += pwr*pow(_freq,-exp)/kB/T0;\
     159                 :            :         cy (n2,n2) += pwr*pow(_freq,-exp)/kB/T0;
     160                 :            : #define _load_flickernoise1(n1,pwr,exp)\
     161                 :            :         cy (n1,n1) += pwr*pow(_freq,-exp)/kB/T0;
     162                 :            : 
     163                 :            : // derivative helper macros
     164                 :            : // transcendental LRM p. 59
     165                 :            : #define m00_cos(v00,x)          v00 = cos(x);
     166                 :            : #define m10_cos(v10,v00,x)      v10 = (-sin(x));
     167                 :            : #define m00_sin(v00,x)          v00 = sin(x);
     168                 :            : #define m10_sin(v10,v00,x)      v10 = (cos(x));
     169                 :            : #define m00_tan(v00,x)          v00 = tan(x);
     170                 :            : #define m10_tan(v10,v00,x)      v10 = (1.0/cos(x)/cos(x));
     171                 :            : #define m00_cosh(v00,x)         v00 = cosh(x);
     172                 :            : #define m10_cosh(v10,v00,x)     v10 = (sinh(x));
     173                 :            : #define m00_sinh(v00,x)         v00 = sinh(x);
     174                 :            : #define m10_sinh(v10,v00,x)     v10 = (cosh(x));
     175                 :            : #define m00_tanh(v00,x)         v00 = tanh(x);
     176                 :            : #define m10_tanh(v10,v00,x)     v10 = (1.0/cosh(x)/cosh(x));
     177                 :            : #define m00_acos(v00,x)         v00 = acos(x);
     178                 :            : #define m10_acos(v10,v00,x)     v10 = (-1.0/sqrt(1-x*x));
     179                 :            : #define m00_asin(v00,x)         v00 = asin(x);
     180                 :            : #define m10_asin(v10,v00,x)     v10 = (+1.0/sqrt(1-x*x));
     181                 :            : #define m00_atan(v00,x)         v00 = atan(x);
     182                 :            : #define m10_atan(v10,v00,x)     v10 = (+1.0/(1+x*x));
     183                 :            : #define m00_hypot(v00,x,y)      v00 = sqrt((x)*(x)+(y)*(y));
     184                 :            : #define m10_hypot(v10,v00,x,y)  v10 = (x)/(v00);
     185                 :            : #define m11_hypot(v11,v00,x,y)  v11 = (y)/(v00);
     186                 :            : #define m00_atan2(v00,x,y)      v00 = atan2(x,y);
     187                 :            : // TODO atan2 derivatives ?
     188                 :            : #define m00_acosh(v00,x)        v00 = acosh(x);
     189                 :            : #define m10_acosh(v10,v00,x)    v10 = (1.0/(sqrt(x-1)*sqrt(x+1)));
     190                 :            : #define m00_asinh(v00,x)        v00 = asinh(x);
     191                 :            : #define m10_asinh(v10,v00,x)    v10 = (1.0/(sqrt(x*x+1)));
     192                 :            : #define m00_atanh(v00,x)        v00 = atanh(x);
     193                 :            : #define m10_atanh(v10,v00,x)    v10 = (1.0/(1-x*x));
     194                 :            : 
     195                 :            : 
     196                 :            : // standard functions LRM p.58
     197                 :            : #define m00_logE(v00,x)         v00 = log(x);
     198                 :            : #define m10_logE(v10,v00,x)     v10 = (1.0/x);
     199                 :            : #define m00_log10(v00,x)        v00 = log10(x);
     200                 :            : #define m10_log10(v10,v00,x)    v10 = (1.0/x/M_LN10);
     201                 :            : #define m00_exp(v00,x)          v00 = exp(x);
     202                 :            : #define m10_exp(v10,v00,x)      v10 = v00;
     203                 :            : #define m00_sqrt(v00,x)         v00 = sqrt(x);
     204                 :            : #define m10_sqrt(v10,v00,x)     v10 = (0.5/v00);
     205                 :            : #define m00_min(v00,x,y)        v00 = ((x)<(y))?(x):(y);
     206                 :            : #define m10_min(v10,v00,x,y)    v10 = ((x)<(y))?1.0:0.0;
     207                 :            : #define m11_min(v11,v00,x,y)    v11 = ((x)<(y))?0.0:1.0;
     208                 :            : #define m00_max(v00,x,y)        v00 = ((x)>(y))?(x):(y);
     209                 :            : #define m10_max(v10,v00,x,y)    v10 = ((x)>(y))?1.0:0.0;
     210                 :            : #define m11_max(v11,v00,x,y)    v11 = ((x)>(y))?0.0:1.0;
     211                 :            : #define m00_pow(v00,x,y)        v00 = pow(x,y);
     212                 :            : #define m10_pow(v10,v00,x,y)    v10 = (x==0.0)?0.0:(v00)*(y)/(x);
     213                 :            : #define m11_pow(v11,v00,x,y)    v11 = (x==0.0)?0.0:(log(x)*(v00));
     214                 :            : #define m00_abs(v00,x)          v00 = ((x)<(0)?(-(x)):(x));
     215                 :            : #define m10_abs(v10,v00,x)      v10 = (((x)>=0)?(+1.0):(-1.0));
     216                 :            : #define m00_floor(v00,x)        v00 = floor(x);
     217                 :            : #define m10_floor(v10,v00,x)    v10 = 1.0;
     218                 :            : 
     219                 :            : #define m00_ceil(v00,x)         v00 = ceil(x);
     220                 :            : // TODO ceil derivative, needed?
     221                 :            : 
     222                 :            : // analog operator, LRM p.61
     223                 :            : #define m00_limexp(v00,x)       v00 = ((x)<80.0?exp(x):exp(80.0)*(x-79.0));
     224                 :            : #define m10_limexp(v10,v00,x)   v10 = ((x)<80.0?(v00):exp(80.0));
     225                 :            : 
     226                 :            : // analog kernel parameter system functions, LRM p.215
     227                 :            : #define m00_vt(x)               (kBoverQ*(x))
     228                 :            : #define m10_vt(x)               (kBoverQ)
     229                 :            : 
     230                 :            : // extra functions (?)
     231                 :            : #define m00_div(v00,v10,x,y)    double v10=1/(y); double v00=(x)*v10;
     232                 :            : #define m10_div(v10,v00,vv,x,y)
     233                 :            : #define m11_div(v11,v00,vv,x,y) double v11 = -v00*vv;
     234                 :            : #define m00_mult(v00,v10,v11,x,y) double v10=(x); double v11=(y); double v00=v10*v11;
     235                 :            : #define m00_add(v00,x,y)        double v00=(x)+(y);
     236                 :            : 
     237                 :            : // second derivatives
     238                 :            : #define m20_logE(v00)           (-1.0/v00/v00)
     239                 :            : #define m20_exp(v00)            exp(v00)
     240                 :            : #define m20_limexp(v00)         ((v00)<80.0?exp(v00):0.0)
     241                 :            : #define m20_sqrt(v00)           (-0.25/(v00)/sqrt(v00))
     242                 :            : #define m20_abs(v00)           0.0
     243                 :            : #define m20_pow(x,y)            ((y)*((y)-1.0)*pow(x,y)/(x)/(x))
     244                 :            : 
     245                 :            : 
     246                 :            : // simulator specific definitions
     247                 :            : #define _modelname              "vcresistor"
     248                 :            : #define _instancename           getName()
     249                 :            : #define _circuit_temp           (getPropertyDouble("Temp")+273.15)
     250                 :            : #define _param_given(p)         (isPropertyGiven(p)?1:0)
     251                 :            : 
     252                 :            : 
     253                 :            : // $vt and $vt() functions
     254                 :            : #define _vt_nom                 (kBoverQ*_circuit_temp)
     255                 :            : 
     256                 :            : using namespace qucs::device;
     257                 :            : using qucs::matrix;
     258                 :            : 
     259                 :            : /* Device constructor. */
     260                 :          0 : vcresistor::vcresistor() : circuit (4)
     261                 :            : {
     262                 :          0 :   type = CIR_vcresistor;
     263                 :          0 : }
     264                 :            : 
     265                 :            : /* Initialization of model. */
     266                 :          0 : void vcresistor::initModel (void)
     267                 :            : {
     268                 :            :   // create internal nodes
     269                 :            : 
     270                 :            :   // get device model parameters
     271                 :          0 :   loadVariables ();
     272                 :            :   // evaluate global model equations
     273                 :          0 :   initializeModel ();
     274                 :            :   // evaluate initial step equations
     275                 :          0 :   initialStep ();
     276                 :            :   // evaluate global instance equations
     277                 :          0 :   initializeInstance ();
     278                 :          0 : }
     279                 :            : 
     280                 :            : /* Initialization of DC analysis. */
     281                 :          0 : void vcresistor::initDC (void)
     282                 :            : {
     283                 :          0 :   allocMatrixMNA ();
     284                 :          0 :   initModel ();
     285                 :          0 :   pol = 1;
     286                 :          0 :   restartDC ();
     287                 :          0 :   doAC = 1;
     288                 :          0 :   doTR = 0;
     289                 :          0 :   doHB = 0;
     290                 :          0 : }
     291                 :            : 
     292                 :            : /* Run when DC is restarted (fallback algorithms). */
     293                 :          0 : void vcresistor::restartDC (void)
     294                 :            : {
     295                 :          0 : }
     296                 :            : 
     297                 :            : /* Initialize Verilog-AMS code. */
     298                 :          0 : void vcresistor::initVerilog (void)
     299                 :            : {
     300                 :            :   // initialization of noise variables
     301                 :            : 
     302                 :            :   int i1, i2, i3, i4;
     303                 :            : 
     304                 :            :   // zero charges
     305         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     306         [ #  # ]:          0 :   for (i2 = 0; i2 < 4; i2++) {
     307                 :          0 :     _charges[i1][i2] = 0.0;
     308                 :            :   } }
     309                 :            : 
     310                 :            :   // zero capacitances
     311         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     312         [ #  # ]:          0 :   for (i2 = 0; i2 < 4; i2++) {
     313         [ #  # ]:          0 :   for (i3 = 0; i3 < 4; i3++) {
     314         [ #  # ]:          0 :   for (i4 = 0; i4 < 4; i4++) {
     315                 :          0 :     _caps[i1][i2][i3][i4] = 0.0;
     316                 :            :   } } } }
     317                 :            : 
     318                 :            :   // zero right hand side, static and dynamic jacobian
     319         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     320                 :          0 :     _rhs[i1] = 0.0;
     321                 :          0 :     _qhs[i1] = 0.0;
     322                 :          0 :     _chs[i1] = 0.0;
     323                 :          0 :     _ghs[i1] = 0.0;
     324         [ #  # ]:          0 :     for (i2 = 0; i2 < 4; i2++) {
     325                 :          0 :       _jstat[i1][i2] = 0.0;
     326                 :          0 :       _jdyna[i1][i2] = 0.0;
     327                 :            :     }
     328                 :            :   }
     329                 :          0 : }
     330                 :            : 
     331                 :            : /* Load device model input parameters. */
     332                 :          0 : void vcresistor::loadVariables (void)
     333                 :            : {
     334                 :          0 :   gain = getPropertyDouble ("gain");
     335                 :          0 : }
     336                 :            : 
     337                 :            : /* #define's for translated code */
     338                 :            : #undef  _DDT
     339                 :            : #define _DDT(q) q
     340                 :            : #define _DYNAMIC
     341                 :            : #define _DERIVATE
     342                 :            : #define _DDX
     343                 :            : #define _DERIVATEFORDDX
     344                 :            : 
     345                 :            : /* Evaluate Verilog-AMS equations in model initialization. */
     346                 :          0 : void vcresistor::initializeModel (void)
     347                 :            : {
     348                 :          0 : }
     349                 :            : 
     350                 :            : /* Evaluate Verilog-AMS equations in instance initialization. */
     351                 :          0 : void vcresistor::initializeInstance (void)
     352                 :            : {
     353                 :          0 : }
     354                 :            : 
     355                 :            : /* Evaluate Verilog-AMS equations in initial step. */
     356                 :          0 : void vcresistor::initialStep (void)
     357                 :            : {
     358                 :          0 : }
     359                 :            : 
     360                 :            : /* Evaluate Verilog-AMS equations in final step. */
     361                 :          0 : void vcresistor::finalStep (void)
     362                 :            : {
     363                 :          0 : }
     364                 :            : 
     365                 :            : /* Evaluate Verilog-AMS equations in analog block. */
     366                 :          0 : void vcresistor::calcVerilog (void)
     367                 :            : {
     368                 :            : 
     369                 :            : /* ----------------- evaluate verilog analog equations -------------------- */
     370                 :            : double resistance;
     371                 :            : #if defined(_DERIVATE)
     372                 :            : double resistance_VC1_C2;
     373                 :            : #endif
     374                 :            : #if defined(_DERIVATE)
     375                 :          0 : resistance_VC1_C2=(gain);
     376                 :            : #endif
     377         [ #  # ]:          0 : resistance=(BP(C1,C2)*gain);
     378         [ #  # ]:          0 : if
     379                 :            : (resistance>0.0)
     380                 :            : {
     381 [ #  # ][ #  # ]:          0 : _load_static_residual2(P1,P2,(BP(P1,P2)/resistance));
     382                 :            : #if defined(_DERIVATE)
     383 [ #  # ][ #  # ]:          0 : _load_static_jacobian4(P1,P2,C1,C2,(-BP(P1,P2)*resistance_VC1_C2/resistance/resistance));
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
                 [ #  # ]
     384 [ #  # ][ #  # ]:          0 : _load_static_jacobian4(P1,P2,P1,P2,(1/resistance));
         [ #  # ][ #  # ]
                 [ #  # ]
     385                 :            : #endif
     386                 :            : }
     387                 :            : else
     388                 :            : {
     389 [ #  # ][ #  # ]:          0 : _load_static_residual2(P1,P2,(BP(P1,P2)*1e6));
     390                 :            : #if defined(_DERIVATE)
     391 [ #  # ][ #  # ]:          0 : _load_static_jacobian4(P1,P2,P1,P2,(1e6));
         [ #  # ][ #  # ]
                 [ #  # ]
     392                 :            : #endif
     393                 :            : }
     394                 :            : 
     395                 :            : /* ------------------ end of verilog analog equations --------------------- */
     396                 :            : 
     397                 :            : /* ------------------ evaluate verilog noise equations -------------------- */
     398                 :            : 
     399                 :            : /* ------------------- end of verilog noise equations --------------------- */
     400                 :          0 : }
     401                 :            : 
     402                 :            : /* Perform DC iteration. */
     403                 :          0 : void vcresistor::calcDC (void)
     404                 :            : {
     405                 :            :   // evaluate Verilog code
     406                 :          0 :   initVerilog ();
     407                 :          0 :   calcVerilog ();
     408                 :            : 
     409                 :            :   // fill right hand side and static jacobian
     410         [ #  # ]:          0 :   for (int i1 = 0; i1 < 4; i1++) {
     411         [ #  # ]:          0 :     setI (i1, _rhs[i1]);
     412         [ #  # ]:          0 :     for (int i2 = 0; i2 < 4; i2++) {
     413         [ #  # ]:          0 :       setY (i1, i2, _jstat[i1][i2]);
     414                 :            :     }
     415                 :            :   }
     416                 :          0 : }
     417                 :            : 
     418                 :            : /* Save operating points. */
     419                 :          0 : void vcresistor::saveOperatingPoints (void)
     420                 :            : {
     421                 :            :   // save global instance operating points
     422                 :          0 : }
     423                 :            : 
     424                 :            : /* Load operating points. */
     425                 :          0 : void vcresistor::loadOperatingPoints (void)
     426                 :            : {
     427                 :          0 : }
     428                 :            : 
     429                 :            : /* Calculate operating points. */
     430                 :          0 : void vcresistor::calcOperatingPoints (void)
     431                 :            : {
     432                 :          0 : }
     433                 :            : 
     434                 :            : /* Initialization of AC analysis. */
     435                 :          0 : void vcresistor::initAC (void)
     436                 :            : {
     437                 :          0 :   allocMatrixMNA ();
     438                 :          0 : }
     439                 :            : 
     440                 :            : /* Perform AC calculations. */
     441                 :          0 : void vcresistor::calcAC (nr_double_t frequency)
     442                 :            : {
     443         [ #  # ]:          0 :   setMatrixY (calcMatrixY (frequency));
     444                 :          0 : }
     445                 :            : 
     446                 :            : /* Compute Y-matrix for AC analysis. */
     447                 :          0 : matrix vcresistor::calcMatrixY (nr_double_t frequency)
     448                 :            : {
     449                 :          0 :   _freq = frequency;
     450                 :          0 :   saveOperatingPoints ();
     451                 :          0 :   matrix y (4);
     452                 :            : 
     453         [ #  # ]:          0 :   for (int i1 = 0; i1 < 4; i1++) {
     454         [ #  # ]:          0 :     for (int i2 = 0; i2 < 4; i2++) {
     455                 :          0 :       y (i1,i2) = nr_complex_t (_jstat[i1][i2], _jdyna[i1][i2] * 2 * M_PI * _freq);
     456                 :            :     }
     457                 :            :   }
     458                 :            : 
     459                 :          0 :   return y;
     460                 :            : }
     461                 :            : 
     462                 :            : /* Initialization of S-parameter analysis. */
     463                 :          0 : void vcresistor::initSP (void)
     464                 :            : {
     465                 :          0 :   allocMatrixS ();
     466                 :          0 : }
     467                 :            : 
     468                 :            : /* Perform S-parameter calculations. */
     469                 :          0 : void vcresistor::calcSP (nr_double_t frequency)
     470                 :            : {
     471 [ #  # ][ #  # ]:          0 :   setMatrixS (ytos (calcMatrixY (frequency)));
         [ #  # ][ #  # ]
                 [ #  # ]
     472                 :          0 : }
     473                 :            : 
     474                 :            : /* Initialization of transient analysis. */
     475                 :          0 : void vcresistor::initTR (void)
     476                 :            : {
     477                 :          0 :   setStates (2 * 4 * 4);
     478                 :          0 :   initDC ();
     479                 :          0 : }
     480                 :            : 
     481                 :            : /* Perform transient analysis iteration step. */
     482                 :          0 : void vcresistor::calcTR (nr_double_t)
     483                 :            : {
     484                 :          0 :   doHB = 0;
     485                 :          0 :   doAC = 1;
     486                 :          0 :   doTR = 1;
     487                 :          0 :   calcDC ();
     488                 :            : 
     489                 :            :   int i1, i2, i3, i4, state;
     490                 :            : 
     491                 :            :   // 2-node charge integrations
     492         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     493         [ #  # ]:          0 :   for (i2 = 0; i2 < 4; i2++) {
     494                 :          0 :     state = 2 * (i2 + 4 * i1);
     495         [ #  # ]:          0 :     if (i1 != i2)
     496         [ #  # ]:          0 :     if (_charges[i1][i2] != 0.0)
     497                 :          0 :       transientCapacitanceQ (state, i1, i2, _charges[i1][i2]);
     498                 :            :   } }
     499                 :            : 
     500                 :            :   // 1-node charge integrations
     501         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     502                 :          0 :     state = 2 * (i1 + 4 * i1);
     503         [ #  # ]:          0 :     if (_charges[i1][i1] != 0.0)
     504                 :          0 :       transientCapacitanceQ (state, i1, _charges[i1][i1]);
     505                 :            :   }
     506                 :            : 
     507                 :            :   // charge: 2-node, voltage: 2-node
     508         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     509         [ #  # ]:          0 :   for (i2 = 0; i2 < 4; i2++) {
     510         [ #  # ]:          0 :   if (i1 != i2)
     511         [ #  # ]:          0 :   for (i3 = 0; i3 < 4; i3++) {
     512         [ #  # ]:          0 :   for (i4 = 0; i4 < 4; i4++) {
     513         [ #  # ]:          0 :     if (i3 != i4)
     514         [ #  # ]:          0 :     if (_caps[i1][i2][i3][i4] != 0.0)
     515 [ #  # ][ #  # ]:          0 :       transientCapacitanceC (i1, i2, i3, i4, _caps[i1][i2][i3][i4], BP(i3,i4));
     516                 :            :   } } } }
     517                 :            : 
     518                 :            :   // charge: 2-node, voltage: 1-node
     519         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     520         [ #  # ]:          0 :   for (i2 = 0; i2 < 4; i2++) {
     521         [ #  # ]:          0 :   if (i1 != i2)
     522         [ #  # ]:          0 :   for (i3 = 0; i3 < 4; i3++) {
     523         [ #  # ]:          0 :     if (_caps[i1][i2][i3][i3] != 0.0)
     524         [ #  # ]:          0 :       transientCapacitanceC2Q (i1, i2, i3, _caps[i1][i2][i3][i3], NP(i3));
     525                 :            :   } } }
     526                 :            : 
     527                 :            :   // charge: 1-node, voltage: 2-node
     528         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     529         [ #  # ]:          0 :   for (i3 = 0; i3 < 4; i3++) {
     530         [ #  # ]:          0 :   for (i4 = 0; i4 < 4; i4++) {
     531         [ #  # ]:          0 :     if (i3 != i4)
     532         [ #  # ]:          0 :     if (_caps[i1][i1][i3][i4] != 0.0)
     533 [ #  # ][ #  # ]:          0 :       transientCapacitanceC2V (i1, i3, i4, _caps[i1][i1][i3][i4], BP(i3,i4));
     534                 :            :   } } }
     535                 :            : 
     536                 :            :   // charge: 1-node, voltage: 1-node
     537         [ #  # ]:          0 :   for (i1 = 0; i1 < 4; i1++) {
     538         [ #  # ]:          0 :   for (i3 = 0; i3 < 4; i3++) {
     539         [ #  # ]:          0 :     if (_caps[i1][i1][i3][i3] != 0.0)
     540         [ #  # ]:          0 :       transientCapacitanceC (i1, i3, _caps[i1][i1][i3][i3], NP(i3));
     541                 :            :   } }
     542                 :          0 : }
     543                 :            : 
     544                 :            : /* Compute Cy-matrix for AC noise analysis. */
     545                 :          0 : matrix vcresistor::calcMatrixCy (nr_double_t frequency) 
     546                 :            : {
     547                 :          0 :   _freq = frequency;
     548                 :          0 :   matrix cy (4);
     549                 :            : 
     550                 :            : 
     551                 :          0 :   return cy;
     552                 :            : }
     553                 :            : 
     554                 :            : /* Perform AC noise computations. */
     555                 :          0 : void vcresistor::calcNoiseAC (nr_double_t frequency) 
     556                 :            : {
     557         [ #  # ]:          0 :   setMatrixN (calcMatrixCy (frequency));
     558                 :          0 : }
     559                 :            : 
     560                 :            : /* Perform S-parameter noise computations. */
     561                 :          0 : void vcresistor::calcNoiseSP (nr_double_t frequency) 
     562                 :            : {
     563 [ #  # ][ #  # ]:          0 :   setMatrixN (cytocs (calcMatrixCy (frequency) * z0, getMatrixS ()));
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
                 [ #  # ]
     564                 :          0 : }
     565                 :            : 
     566                 :            : /* Initialization of HB analysis. */
     567                 :          0 : void vcresistor::initHB (int)
     568                 :            : {
     569                 :          0 :   initDC ();
     570                 :          0 :   allocMatrixHB ();
     571                 :          0 : }
     572                 :            : 
     573                 :            : /* Perform HB analysis. */
     574                 :          0 : void vcresistor::calcHB (int)
     575                 :            : {
     576                 :          0 :   doHB = 1;
     577                 :          0 :   doAC = 1;
     578                 :          0 :   doTR = 0;
     579                 :            : 
     580                 :            :   // jacobian dI/dV and currents get filled
     581                 :          0 :   calcDC ();
     582                 :          0 :   saveOperatingPoints ();
     583                 :            : 
     584                 :            :   // fill in HB matrices
     585         [ #  # ]:          0 :   for (int i1 = 0; i1 < 4; i1++) {
     586         [ #  # ]:          0 :     setQ  (i1, _qhs[i1]); // charges
     587         [ #  # ]:          0 :     setCV (i1, _chs[i1]); // jacobian dQ/dV * V
     588         [ #  # ]:          0 :     setGV (i1, _ghs[i1]); // jacobian dI/dV * V
     589         [ #  # ]:          0 :     for (int i2 = 0; i2 < 4; i2++) {
     590         [ #  # ]:          0 :       setQV (i1, i2, _jdyna[i1][i2]); // jacobian dQ/dV
     591                 :            :     }
     592                 :            :   }
     593                 :          0 : }
     594                 :            : 
     595                 :            : #include "vcresistor.defs.h"
     596                 :            : 

Generated by: LCOV version 1.11