def netlist_solver(cir_netlist):
net = ntl.Network(cir_netlist)
net_solve(net)
net.branch_voltage()
net.branch_current()
net.branch_power()
net.print()
def solve_for_v_i(log, bin_num):
with open("C:\\Users\\Nolan\\Desktop\\vanilla.txt", "r+",
encoding='utf8') as f:
data = f.readlines()
combined_data = []
for i in range(len(bin_num)):
ranges = [[0, 8], [8, 16], [16, 24], [24, 32]]
counter = 0
if int(bin_num[i]) == 0:
for things in range(ranges[i][0], ranges[i][1]):
combined_data.append(''.join([
data[things].strip(),
str(i + 4), ' ',
str(counter + 1), ' ',
str(log[things]), '\n'
]))
counter += 1
else:
for things in range(ranges[i][0], ranges[i][1]):
combined_data.append(''.join([
data[things].strip(),
str(1), ' ',
str(counter + 1), ' ',
str(log[things]), '\n'
]))
counter += 1
with open("C:\\Users\\Nolan\\Desktop\\input.net", 'w+',
encoding='utf8') as f:
f.writelines(combined_data)
f.writelines(data[32:])
net = ntl.Network('C:/Users/Nolan/Desktop/input.net')
try:
net_solve(net)
except:
return (0, 0, False)
resistor_num_list = list(range(1, 33))
resistor_list = []
currents = []
for items in resistor_num_list:
resistor_list.append('R' + str(items))
currents.append(net.get_current('R' + str(items))[0])
resistor_list_dividers = list(range(33, 41))
voltages = []
for items in resistor_list_dividers:
voltages.append(net.get_voltage('R' + str(items))[0])
return (voltages, currents, True)
# import SpicePy modules
import netlist as ntl
from netsolve import net_solve
import matplotlib.pyplot as plt
plt.ion()
# read netlist
net = ntl.Network('tran_network1.net')
# compute the circuit solution
net_solve(net)
# plot results
net.plot()
import numpy as np
# check rms
def check(rms, toll):
if rms < toll:
print('OK')
else:
print('*FAILED*')
print('Start benchmark:')
# DC network
print(' * Testing DC network ...', end=' ',flush=True)
net = ntl.Network('../demo/dc_network.net')
net_solve(net)
ref = np.loadtxt('./dc_network/solution.txt')
rms = np.sqrt(np.sum((ref - net.x) ** 2)) / ref.size
check(rms, 1e-10)
# OP network
print(' * Testing OP network ...', end=' ',flush=True)
net = ntl.Network('../demo/op_network.net')
net_solve(net)
ref = np.loadtxt('./op_network/solution.txt')
rms = np.sqrt(np.sum((ref - net.x) ** 2)) / ref.size
check(rms, 1e-5)
# AC network (single frequency)
print(' * Testing AC network (single frequency) ...', end=' ',flush=True)
# import SpicePy modules
import netlist as ntl
from netsolve import net_solve
# read netlist
net = ntl.Network('dc_network.net')
# compute the circuit solution
net_solve(net)
# compute branch quantities
net.branch_voltage()
net.branch_current()
net.branch_power()
# print results
net.print()
def get_solution(fname, update, context):
"""
'get_solution' computes the solution of a network using SpicePy
:param fname: filename with the netlist
:param update: bot update
:param context: CallbackContext
:return:
* mex: solution formatted in a string
"""
try:
# create network and solve it
net = ntl.Network(fname)
# check if number of nodes exceeds the limit
NMAX = 40
if net.node_num > NMAX:
mex = "Your netlist includes more than {:d} nodes.\n".format(
net.node_num)
mex += "*The maximum allowed number on this bot is {:d}.*\n".format(
NMAX)
mex += "Please reduce the number of nodes or take a look to the computational core of this bot "
mex += "that does not have this limitation:\n"
mex += "[SpicePy project](https://github.com/giaccone/SpicePy)"
else:
# limit sample for .tran to 2000 (max)
if net.analysis[0].lower() == '.tran':
Nsamples = float(net.convert_unit(net.analysis[2])) / float(
net.convert_unit(net.analysis[1]))
if Nsamples > 2000:
step = float(net.convert_unit(net.analysis[2])) / 1999
net.analysis[1] = '{:.3e}'.format(step)
mex = "Your netlits defines a '.tran' analysis with *{:d}* samples\n".format(
int(Nsamples))
mex += "Since this bot runs on a limited hardware shared by many users\n"
mex += "The analysis has been limited to *2000* samples:\n"
mex += "`.tran " + net.analysis[1] + " " + net.analysis[
-1] + "`"
context.bot.send_message(
chat_id=update.message.chat_id,
text=mex,
parse_mode=telegram.ParseMode.MARKDOWN)
# limit sample for .ac to 2000 (max)
elif net.analysis[0].lower() == '.ac':
# get frequencies
net.frequency_span()
if not np.isscalar(net.f):
# get Nsamples
Nsamples = len(net.f)
# limit di 2000 max
if Nsamples > 2000:
scale = 2000 / Nsamples
old_analysys = "`" + " ".join(net.analysis) + "`"
net.analysis[2] = str(
int(
np.ceil(
scale *
float(net.convert_unit(net.analysis[2])))))
new_analysys = "`" + " ".join(net.analysis) + "`"
mex = "Your netlits defines a '.tran' analysis with *{:d}* samples\n".format(
int(Nsamples))
mex += "Since this bot runs on a limited hardware shared by many users\n"
mex += "The analysis has been limited to *2000* samples:\n"
mex += "original analysis: " + old_analysys + "\n"
mex += "ner analysis: " + new_analysys + "\n"
context.bot.send_message(
chat_id=update.message.chat_id,
text=mex,
parse_mode=telegram.ParseMode.MARKDOWN)
# get configurations
fname = './users/' + str(update.message.chat_id) + '.cnf'
fid = open(fname, 'r')
flag = fid.readline()[:-1] # read nodal_pot conf
nodal_pot = flag == 'True'
flag = fid.readline()[:-1] # read polar conf
polar = flag == 'True'
flag = fid.readline() # read dB conf
dB = flag == 'True'
if net.analysis[0] == '.op':
# forcepolar to False for .op problems
polar = False
# solve the network
net_solve(net)
# .op and .ac (single-frequency): prepare mex to be printed
if (net.analysis[0].lower() == '.op') | (
(net.analysis[0].lower() == '.ac') & (np.isscalar(net.f))):
# get branch quantities
net.branch_voltage()
net.branch_current()
net.branch_power()
# prepare message
mex = net.print(polar=polar, message=True)
mex = mex.replace(
'==============================================\n'
' branch quantities'
'\n==============================================\n',
'*branch quantities*\n`')
mex = mex.replace(
'----------------------------------------------', '')
mex += '`'
# if the user wants node potentials add it to mex
if nodal_pot:
# create local dictionary node-number 2 node-label
num2node_label = {
num: name
for name, num in net.node_label2num.items()
if name != '0'
}
# compute the node potentials
mex0 = '*node potentials*\n`'
for num in num2node_label:
voltage = net.get_voltage('(' + num2node_label[num] +
')')[0]
if polar:
mex0 += 'v(' + num2node_label[
num] + ') = {:10.4f} V < {:10.4f}°\n'.format(
np.abs(voltage),
np.angle(voltage) * 180 / np.pi)
else:
mex0 += 'v(' + num2node_label[
num] + ') = {:10.4f} V\n'.format(voltage)
# add newline
mex0 += '`\n\n'
# add node potentials before branch quantities
mex = mex0 + mex
elif net.analysis[0].lower() == '.tran':
hf = net.plot(to_file=True,
filename='./users/tran_plot_' +
str(update.message.chat_id) + '.png',
dpi_value=150)
mex = None
plt.close(hf)
elif net.analysis[0].lower() == '.ac':
hf = net.bode(to_file=True,
decibel=dB,
filename='./users/bode_plot_' +
str(update.message.chat_id) + '.png',
dpi_value=150)
mex = None
if isinstance(hf, list):
for fig in hf:
plt.close(fig)
else:
plt.close(hf)
# Log every time a network is solved
# To make stat it is saved the type of network and the UserID
StatLog.info('Analysis: ' + net.analysis[0] + ' - UserID: ' +
str(update.effective_user.id))
return net, mex
except:
# set network to None
net = None
# read network with issues
wrong_net = ''
with open(fname) as f:
for line in f:
wrong_net += line
wrong_net = wrong_net.replace('\n', ' / ')
# log error
SolverLog.error('UserID: ' + str(update.effective_user.id) +
' - Netlist error: ' + wrong_net)
return net, "*Something went wrong with your netlist*.\nPlease check the netlist format."
# import SciPy modules
import netlist as ntl
from netsolve import net_solve
import matplotlib.pyplot as plt
plt.ion()
# read netlist
net = ntl.Network('ac_low_high_pass_filter.net')
# compute the circuit solution
net_solve(net)
# plot bode diagrams
net.bode()
# import SpicePy modules
import netlist as ntl
from netsolve import net_solve
# read netlist
net = ntl.Network('op_network.net')
# compute the circuit solution
net_solve(net)
# compute branch quantities
net.branch_voltage()
net.branch_current()
net.branch_power()
# print results
net.print()
# import SpicePy modules
import netlist as ntl
from netsolve import net_solve
# read netlist
net = ntl.Network('VCVS_and_CCCS.net')
# compute the circuit solution
net_solve(net)
# compute branch quantities
net.branch_voltage()
net.branch_current()
net.branch_power()
# print results
net.print()
# import SciPy modules
import netlist as ntl
from netsolve import net_solve
import matplotlib.pyplot as plt
plt.ion()
# read netlist
net = ntl.Network('ac_band_pass_filter.net')
# compute the circuit solution
net_solve(net)
# plot bode diagrams
net.bode()
# import SpicePy modules
import netlist as ntl
from netsolve import net_solve
# read netlist
net = ntl.Network('ac_single_frequency.net')
# compute the circuit solution
net_solve(net)
# compute branch quantities
net.branch_voltage()
net.branch_current()
net.branch_power()
# print results
net.print()
# print results in polar format
net.print(polar=True)
# -*- coding: utf-8 -*-
"""
@author: motanov
"""
import netlist as ntl
from netsolve import net_solve
"""
Prototype : just solves the circuit using SpicePy library
To do:
Image to netlist converter ( photo of circuit -> netlist of given circuit )
"""
net = ntl.Network('network.net')
net_solve(net)
net.branch_voltage()
net.branch_current()
net.branch_power()
net.print()
# import SpicePy modules
import netlist as ntl
from netsolve import net_solve
# read netlist
net = ntl.Network('VCCS_and_CCVS.net')
# compute the circuit solution
net_solve(net)
# compute branch quantities
net.branch_voltage()
net.branch_current()
net.branch_power()
# print results
net.print()
