Skip to content

Instantly share code, notes, and snippets.

Last active January 22, 2024 10:22
Show Gist options
  • Save snmishra/27dcc624b639c2626137 to your computer and use it in GitHub Desktop.
Save snmishra/27dcc624b639c2626137 to your computer and use it in GitHub Desktop.
Short python script to read ngspice raw binary files
# MIT license:
# See
# for a more complete library. Isotel's version is GPL licensed
from __future__ import division
import numpy as np
BSIZE_SP = 512 # Max size of a line of data; we don't want to read the
# whole file to find a line, in case file does not have
# expected structure.
MDATA_LIST = [b'title', b'date', b'plotname', b'flags', b'no. variables',
b'no. points', b'dimensions', b'command', b'option']
def rawread(fname: str):
"""Read ngspice binary raw files. Return tuple of the data, and the
plot metadata. The dtype of the data contains field names. This is
not very robust yet, and only supports ngspice.
>>> darr, mdata = rawread('')
>>> darr.dtype.names
>>> plot(np.real(darr['frequency']), np.abs(darr['v(out)']))
# Example header of raw file
# Title: rc band pass example circuit
# Date: Sun Feb 21 11:29:14 2016
# Plotname: AC Analysis
# Flags: complex
# No. Variables: 3
# No. Points: 41
# Variables:
# 0 frequency frequency grid=3
# 1 v(out) voltage
# 2 v(in) voltage
# Binary:
fp = open(fname, 'rb')
count = 0
arrs = []
plots = []
plot = {}
while (True):
mdata = fp.readline(BSIZE_SP).split(b':', maxsplit=1)
if len(mdata) == 2:
if mdata[0].lower() in MDATA_LIST:
plot[mdata[0].lower()] = mdata[1].strip()
if mdata[0].lower() == b'variables':
nvars = int(plot[b'no. variables'])
npoints = int(plot[b'no. points'])
plot['varnames'] = []
plot['varunits'] = []
for varn in range(nvars):
varspec = (fp.readline(BSIZE_SP).strip()
assert(varn == int(varspec[0]))
if mdata[0].lower() == b'binary':
rowdtype = np.dtype({'names': plot['varnames'],
'formats': [np.complex_ if b'complex'
in plot[b'flags']
else np.float_]*nvars})
# We should have all the metadata by now
arrs.append(np.fromfile(fp, dtype=rowdtype, count=npoints))
plot = {} # reset the plot dict
fp.readline() # Read to the end of line
return (arrs, plots)
if __name__ == '__main__':
arrs, plots = rawread('test.raw')
# Local Variables:
# mode: python
# End:
Copy link

That's not quite right, unfortunately.
Now the plot dict is reset between each line, instead of between separate plots.

Copy link

@jmspiewak How about now?

Copy link

That's exactly it, thanks.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment