23 #ifndef INTEGRATOR6581_H
24 #define INTEGRATOR6581_H
32 #include "siddefs-fp.h"
162 const unsigned short* vcr_Vg;
163 const unsigned short* vcr_n_Ids_term;
164 const unsigned short* opamp_rev;
166 unsigned int Vddt_Vw_2;
178 const unsigned short Vddt;
179 const unsigned short nVt;
180 const unsigned short nVmin;
181 const unsigned short n_snake;
184 Integrator6581(
const unsigned short* vcr_Vg,
const unsigned short* vcr_n_Ids_term,
185 const unsigned short* opamp_rev,
unsigned short Vddt,
unsigned short nVt,
186 unsigned short nVmin,
unsigned short n_snake,
double N16) :
188 vcr_n_Ids_term(vcr_n_Ids_term),
189 opamp_rev(opamp_rev),
204 void setVw(
unsigned short Vw) { Vddt_Vw_2 = ((Vddt - Vw) * (Vddt - Vw)) >> 1; }
206 int solve(
int vi)
const;
211 #if RESID_INLINING || defined(INTEGRATOR_CPP)
217 int Integrator6581::solve(
int vi)
const
227 const unsigned int Vgst = Vddt - vx;
228 const unsigned int Vgdt = Vddt - vi;
230 const unsigned int Vgst_2 = Vgst * Vgst;
231 const unsigned int Vgdt_2 = Vgdt * Vgdt;
234 const int n_I_snake = n_snake * (
static_cast<int>(Vgst_2 - Vgdt_2) >> 15);
238 const int Vg =
static_cast<int>(vcr_Vg[(Vddt_Vw_2 + (Vgdt_2 >> 1)) >> 16]);
239 const int Vp = (Vg - nVt) / n;
240 const int kVg =
static_cast<int>(Vp) - nVmin;
243 const int Vgs = (vx < kVg) ? kVg - vx : 0;
244 assert(Vgs < (1 << 16));
245 const int Vgd = (vi < kVg) ? kVg - vi : 0;
246 assert(Vgd < (1 << 16));
249 const unsigned int If =
static_cast<unsigned int>(vcr_n_Ids_term[Vgs]) << 15;
250 const unsigned int Ir =
static_cast<unsigned int>(vcr_n_Ids_term[Vgd]) << 15;
251 const int n_I_vcr = (If - Ir) * n;
255 const double gamma = 1.0;
256 const double phi = 0.8;
257 const double Ut = 26.0e-3;
258 const double nVp = Vp / N16;
259 n = 1. + (gamma / (2 * sqrt(nVp + phi + 4*Ut)));
260 assert((n > 1.2) && (n < 1.8));
264 vc += n_I_snake + n_I_vcr;
267 const int tmp = (vc >> 15) + (1 << 15);
268 assert(tmp < (1 << 16));
272 return vx - (vc >> 14);
Definition: Integrator6581.h:160