From 370e52736f3b75c8a96ad0d5ca7134623c6f97a6 Mon Sep 17 00:00:00 2001 From: h_vogt Date: Sat, 4 Feb 2017 10:41:35 +0100 Subject: [PATCH] MC_ring.sp, replace variables by vectors in the loop --- examples/Monte_Carlo/MC_ring.sp | 90 ++++++++++++++++++--------------- 1 file changed, 49 insertions(+), 41 deletions(-) diff --git a/examples/Monte_Carlo/MC_ring.sp b/examples/Monte_Carlo/MC_ring.sp index 58e5c141f..f63ae0467 100644 --- a/examples/Monte_Carlo/MC_ring.sp +++ b/examples/Monte_Carlo/MC_ring.sp @@ -31,7 +31,7 @@ cout buf ss 0.2pF .options noacct .control save buf $ we just need buf, save memory by more than 10x - let mc_runs = 10 $ number of runs for monte carlo + let mc_runs = 30 $ number of runs for monte carlo let run = 0 $ number of actual run set curplot = new $ create a new plot set curplottitle = "Transient outputs" @@ -60,39 +60,38 @@ cout buf ss 0.2pF * of the BSIM3 model for the NMOS and PMOS transistors. * We may obtain the nominal values (nom) by manually extracting them from * the parameter set. Here we get them automatically and store them into -* variables. This has the advantage that you may change the parameter set +* vectors. This has the advantage that you may change the parameter set * without having to look up the values again. - set n1vth0=@n1[vth0] - set n1u0=@n1[u0] - set n1tox=@n1[tox] - set n1lint=@n1[lint] - set n1wint=@n1[wint] - set p1vth0=@p1[vth0] - set p1u0=@p1[u0] - set p1tox=@p1[tox] - set p1lint=@p1[lint] - set p1wint=@p1[wint] + let n1vth0=@n1[vth0] + let n1u0=@n1[u0] + let n1tox=@n1[tox] + let n1lint=@n1[lint] + let n1wint=@n1[wint] + let p1vth0=@p1[vth0] + let p1u0=@p1[u0] + let p1tox=@p1[tox] + let p1lint=@p1[lint] + let p1wint=@p1[wint] + * * run the simulation loop dowhile run <= mc_runs - * without the reset switch there is some strange drift - * towards lower and lower frequencies - reset * run=0 simulates with nominal parameters if run > 0 - altermod @n1[vth0]=gauss($n1vth0, 0.1, 3) - altermod @n1[u0]=gauss($n1u0, 0.05, 3) - altermod @n1[tox]=gauss($n1tox, 0.1, 3) - altermod @n1[lint]=gauss($n1lint, 0.1, 3) - altermod @n1[wint]=gauss($n1wint, 0.1, 3) - altermod @p1[vth0]=gauss($p1vth0, 0.1, 3) - altermod @p1[u0]=gauss($p1u0, 0.1, 3) - altermod @p1[tox]=gauss($p1tox, 0.1, 3) - altermod @p1[lint]=gauss($p1lint, 0.1, 3) - altermod @p1[wint]=gauss($p1wint, 0.1, 3) + setplot $max_fft + altermod @n1[vth0] = gauss(n1vth0, 0.1, 3) + altermod @n1[u0] = gauss(n1u0, 0.05, 3) + altermod @n1[tox] = gauss(n1tox, 0.1, 3) + altermod @n1[lint] = gauss(n1lint, 0.1, 3) + altermod @n1[wint] = gauss(n1wint, 0.1, 3) + altermod @p1[vth0] = gauss(p1vth0, 0.1, 3) + altermod @p1[u0] = gauss(p1u0, 0.1, 3) + altermod @p1[tox] = gauss(p1tox, 0.1, 3 ) + altermod @p1[lint] = gauss(p1lint, 0.1, 3) + altermod @p1[wint] = gauss(p1wint, 0.1, 3) end - tran 15p 50n 0 -* select stop and step so that number of data points after linearization is not too + tran 15p 100n 0 +* select stop and step so that number of data points after linearization is not too * close to 8192, which would yield varying number of line length and thus scale for fft. * * We have to figure out what to do if a single simulation will not converge. @@ -104,8 +103,10 @@ cout buf ss 0.2pF set run ="$&run" $ create a variable from the vector set mc_runs ="$&mc_runs" $ create a variable from the vector echo simulation run no. $run of $mc_runs + set dt = $curplot * save the linearized data for having equal time scales for all runs linearize buf $ linearize only buf, no other vectors needed + destroy $dt $ delete the tran i plot set dt = $curplot $ store the current plot to dt (tran i+1) setplot $plot_out $ make 'plt_out' the active plot * firstly save the time scale once to become the default scale @@ -115,11 +116,14 @@ cout buf ss 0.2pF let vout{$run}={$dt}.buf $ store the output vector to plot 'plot_out' setplot $dt $ go back to the previous plot (tran i+1) fft buf $ run fft on vector buf + destroy $dt $ delete the tran i+1 plot let buf2=db(mag(buf)) * find the frequency where buf has its maximum of the fft signal meas sp fft_max MAX_AT buf2 from=0.1G to=0.7G * find the frequency where buf is -40dB at rising fft signal - meas sp fft_40 WHEN buf2=-40 RISE=1 from=0.1G to=0.7G + meas sp fft_40 WHEN buf2=-40 RISE=1 from=0.1G to=0.7G + echo + echo * store the fft vector set dt = $curplot $ store the current plot to dt (spec i) setplot $plot_fft $ make 'plot_fft' the active plot @@ -130,19 +134,15 @@ cout buf ss 0.2pF * store the measured value setplot $max_fft $ make 'max_fft' the active plot let maxffts[{$run}]={$dt}.fft_max - let halfffts[{$run}]={$dt}.fft_40 -* setplot $plot_out -* The following command does not work here. Why not? Probably not a real copy. -* destroy $dt $ save memory, we don't need this plot (spec) any more - setplot $dt $ go back to the previous plot + let halfffts[{$run}]={$dt}.fft_40 let run = run + 1 end ***** plotting ********************************************************** -* plot {$plot_out}.allv - plot {$plot_out}.vout0 $ just plot the tran output with nominal parameters -* setplot $plot_fft -* plot db(mag(ally)) xlimit .1G 1G ylimit -80 10 - plot db(mag({$plot_fft}.ally)) xlimit .1G 1G ylimit -80 10 + setplot $plot_out + plot vout0 ylabel 'RO output, original parameters' $ just plot the tran output with nominal parameters + setplot $plot_fft + settype decibel ally + plot db(mag(ally)) xlimit .1G 1G ylimit -80 10 ylabel 'fft output' * * create a histogram from vector maxffts setplot $max_fft $ make 'max_fft' the active plot @@ -171,9 +171,17 @@ cout buf ss 0.2pF end let run = run + 1 end + * plot the histogram set plotstyle=combplot - plot yvec-1 vs xvec $ subtract 1 because with started with unitvec containing ones + plot yvec-1 vs xvec xlabel 'oscillation frequency' ylabel 'bin count' $ subtract 1 because we started with unitvec containing ones + + * plot simulation series + set plotstyle=linplot + let xx = vector(mc_runsp) + settype frequency maxffts + plot maxffts vs xx xlabel 'iteration no.' ylabel 'RO frequency' + * calculate jitter let diff40 = (vecmax(halfffts) - vecmin(halfffts))*1e-6 echo @@ -193,11 +201,11 @@ cout buf ss 0.2pF +k3b=2.233 +vsat=86301.58 ua=6.47e-9 ub=4.23e-18 uc=-4.706281e-11 +rdsw=650 u0=388.3203 wr=1 -+a0=.3496967 ags=.1 b0=0.546 b1=1 ++a0=.3496967 ags=.1 b0=0.546 b1=1 +dwg=-6.0e-09 dwb=-3.56e-09 prwb=-.213 +keta=-3.605872e-02 a1=2.778747e-02 a2=.9 +voff=-6.735529e-02 nfactor=1.139926 cit=1.622527e-04 -+cdsc=-2.147181e-05 ++cdsc=-2.147181e-05 +cdscb=0 dvt0w=0 dvt1w=0 dvt2w=0 +cdscd=0 prwg=0 +eta0=1.0281729e-02 etab=-5.042203e-03