Using sweep names and values with famMap()

I should also have mentioned that while ocnPrint() can print the data in this test case, I need a different print format, and  in the future I may also want to do something that ocnPrint() can't handle (such as assigning a sweep-based custom plot label to a plotted waveform).

Best Regards,

Stephen Greenwood

Stephen,

I've attached a function that you might be able to modify to do what you want. The function igfDescendAndSortFamMeasurements() takes a family of measurements, such as OP("da.M1" "vds") when evaluated on a parametric analysis (ADEXL), and returns a list of values sorted based on a string constructed from the sweep names and values.

I originally wrote this so that I could create a table of mixed measurements (dcOpInfo and tran) where each type of measurement had the same order relative to the sweep parameters.

The example below shows the dc operating point Vds extracted from a 3 parameter sweep. I believe the infomation you need for your purpose is available in data table (in x->swpList) and you just need to reformat it.

Hope this helps.

(I'll be out of email contact for the next two weeks, so any further replies will need to wait until early April).

Ivars

Example of output (function definition below):

x = igfDescendAndSortFamMeasurement(OP("da.M1" "vds"))

x->??
(sortList
(("VDD12 1.08 temperature -40.0 vincm 0.5" 0.5417283)
("VDD12 1.08 temperature -40.0 vincm 0.6" 0.4519738)
("VDD12 1.08 temperature -40.0 vincm 0.7" 0.36213)
("VDD12 1.08 temperature -40.0 vincm 0.8" 0.2721419)
("VDD12 1.08 temperature 125.0 vincm 0.5" 0.5636597)
("VDD12 1.08 temperature 125.0 vincm 0.6" 0.4758717)
("VDD12 1.08 temperature 125.0 vincm 0.7" 0.3882314)
("VDD12 1.08 temperature 125.0 vincm 0.8" 0.3005837)
("VDD12 1.32 temperature -40.0 vincm 0.5" 0.7779688)
("VDD12 1.32 temperature -40.0 vincm 0.6" 0.6880795)
("VDD12 1.32 temperature -40.0 vincm 0.7" 0.5980505)
("VDD12 1.32 temperature -40.0 vincm 0.8" 0.5078215)
("VDD12 1.32 temperature 125.0 vincm 0.5" 0.7994591)
("VDD12 1.32 temperature 125.0 vincm 0.6" 0.7115096)
("VDD12 1.32 temperature 125.0 vincm 0.7" 0.6236678)
("VDD12 1.32 temperature 125.0 vincm 0.8" 0.5357573)
) swpList
(((("VDD12" 1.08)
("temperature" -40.0)
("vincm" 0.5)
) 0.5417283
)
((("VDD12" 1.08)
("temperature" 125.0)
("vincm" 0.5)
) 0.5636597
)
((("VDD12" 1.32)
("temperature" -40.0)
("vincm" 0.5)
) 0.7779688
)
((("VDD12" 1.32)
("temperature" 125.0)
("vincm" 0.5)
) 0.7994591
)
((("VDD12" 1.08)
("temperature" -40.0)
("vincm" 0.6)
) 0.4519738
)
((("VDD12" 1.08)
("temperature" 125.0)
("vincm" 0.6)
) 0.4758717
)
((("VDD12" 1.32)
("temperature" -40.0)
("vincm" 0.6)
) 0.6880795
)
((("VDD12" 1.32)
("temperature" 125.0)
("vincm" 0.6)
) 0.7115096
)
((("VDD12" 1.08)
("temperature" -40.0)
("vincm" 0.7)
) 0.36213
)
((("VDD12" 1.08)
("temperature" 125.0)
("vincm" 0.7)
) 0.3882314
)
((("VDD12" 1.32)
("temperature" -40.0)
("vincm" 0.7)
) 0.5980505
)
((("VDD12" 1.32)
("temperature" 125.0)
("vincm" 0.7)
) 0.6236678
)
((("VDD12" 1.08)
("temperature" -40.0)
("vincm" 0.8)
) 0.2721419
)
((("VDD12" 1.08)
("temperature" 125.0)
("vincm" 0.8)
) 0.3005837
)
((("VDD12" 1.32)
("temperature" -40.0)
("vincm" 0.8)
) 0.5078215
)
((("VDD12" 1.32)
("temperature" 125.0)
("vincm" 0.8)
) 0.5357573
)
) meas
(0.5417283 0.4519738 0.36213 0.2721419 0.5636597
0.4758717 0.3882314 0.3005837 0.7779688 0.6880795
0.5980505 0.5078215 0.7994591 0.7115096 0.6236678
0.5357573
) swpTbl table:igfTable
)



;; When different types of measurements are made on a family of results, for
;; example a riseTime() vs. a pv(), the returned list/waveform of results
;; can be in different orders since the descent path through the family heirarchy
;; appears to be different for simulation results compared to simulation parameters.
;; To allow the different type of measurements to be correctly grouped, the sweep
;; parameters are recorded during the descent through the family heirarchy for each
;; measurement. Then the sweep parameters are sorted alphanumerical and then a sort
;; key consisting of concatenation of all the sweep parameter names and values is created.
;; This sweep key is then used to sort the measurement values. This ensures that all
;; measurements will be ordered in the same way (since the sweep parameters are constant
;; across each job).
;;
;; Usage example:
;;
;; x = igfDescendAndSortFamMeasurement(OP("x1.x2.M1" "vds"))
;;
;; A list of the sorted measurements can be obtained using:
;; x->meas
;;
;; The sorting list can be seen using:
;; x->sortList
;;
;; A list of the sweep parameters plus measurement can be obtained:
;; x->swpList
;;
;; use x->?? to see the full table
;;

procedure( igfDescendAndSortFamMeasurement(w @key (data ncons(nil)) )
let( (name val_list next_wave meas smeas)
when(!data->swpTbl
data->swpTbl = makeTable("igfTable")
);when
when(!data->meas ;; list of individual measurements
data->meas = nil
)
when(!data->swpList ;; list of measurements plus the sweep name/values
data->swpList = nil
)
when(!data->sortList ;; sorted list of measurements
data->sortList = nil
)
name = famGetSweepName(w)
if(name
val_list = famGetSweepValues(w)
val_list = nil
)
if(val_list
then
foreach(val val_list
data->swpTbl[name] = val
if(famIsFamily(w)
then
next_wave = famValue(w val)
else
next_wave = nil
meas = famValue(w val)
)
if(next_wave
then
data = igfDescendAndSortFamMeasurement(next_wave ?data data )
else
when(meas == nil
meas = "--"
)
if(drIsWaveform(meas)
m = igfVecToList(drGetWaveformYVec(meas))
m = meas
)
data->meas = append1(data->meas m)
data->swpList = append1(data->swpList list(tableToList(data->swpTbl) m))
strList = nil
keys = sort(setof(k data->swpTbl t) nil ) ;; extract and alpha sort all keys
foreach( key keys strList = append1(strList sprintf(nil "%s %L" key data->swpTbl[key])))
data->sortList = append1(data->sortList list(buildString(strList) m) )
);if
);foreach
smeas = sortcar(data->sortList nil)
data->meas = foreach(mapcar m smeas cadr(m))
else
data->meas = w
);if
`data
);let
);procedure

;; Convert a waveform vector into a list. Based on code by Andrew Beckett
procedure( igfVecToList( vec )
let( (len i vecList)
len = drVectorLength(vec)
for( i 0 sub1(len)
vecList = tconc(vecList drGetElem(vec i) )
);for
car(vecList)
);let
);procedure( igfVecToList

Thanks, Ivars! I will go through your code. Even if I can't use it directly I can try to adapt the concepts in it.

Much obliged,

Stephen

Stephen,

You might find this useful too. It takes a family, and produces a list of abFamilySignal objects, for each leaf waveform. As part of doing this, it constructs a name which contains the sweep name and value. I did this as part of some code which generates a "VCSV" (ViVA CSV) file from arbitrary results - I've not included the entire code, but the key part is here. May not be quite what you want, but it may be illustrative of an approach you could take.

Regards,

Andrew.

;------------------------------------------------------------------------
; Container structure to keep the cleaned signal name,
; length, whether it's complex etc
;------------------------------------------------------------------------
(defstruct abFamilySignal wave name len isComplex)

/*****************************************************************
*                                                                *
* (abUnwrapFamily fam @optional topName top (level 0) (name "")) *
*                                                                *
*    Convert a family or waveform into a flat list of single     *
*  waveforms, wrapped by an abFamilySignal container structure   *
*                                                                *
*****************************************************************/

(defun abUnwrapFamily (fam @optional topName top (level 0) (name ""))
  (let (dim sweepName)
    (unless top (setq top fam))
    (setq dim (if (or (famIsFamily top) (drIsWaveform top))
                (drWaveformGetDimension top)
                0))
    (cond
      ((and (famIsFamily top) (lessp level (sub1 dim)))
       (setq sweepName (famGetSweepName top level))
       (foreach mapcan sweepValue (famGetSweepValues fam)
                (abUnwrapFamily 
                  (famValue fam sweepValue) topName top (add1 level) 
                  (sprintf nil "%s (%s=%Y)" name sweepName sweepValue))
                )
       )
      (t (list (make_abFamilySignal 
                 ?wave fam 
                 ; the messing with strcat in format is to stop SCCS doing 
                 ; keyword substitution!
                 ?name (sprintf nil (strcat "%Y" "%s")
                                (or topName (famGetExpr top)) name))))
      )
    )
  )

Hi Andrew,

After playing with your code and reading some key sections of the sklangref, I am past the point that was hanging me up. I understand these forums are not your "day job" and I appreciate the help!

Best Regards,

Stephen