card_list.do_ac() card_list.ac_load() print "Loaded AC-matrix", gnucap.get_complex_array(gnucap.cvar.SIM_ac, n) print gnucap.bsmatrix_to_array_c(acx) ## Solve acx.lu_decomp() acx.fbsub(gnucap.cvar.SIM_ac) print "rhs after", gnucap.get_complex_array(gnucap.cvar.SIM_ac, n) def setup(self, cmd): pass def sweep(self): pass myac = MyAC() d0=gnucap.install_command("myac", myac) d1=gnucap.install_command("ac", myac) gnucap.command("store ac vm(2)") gnucap.command("myac") w= gnucap.SIM.find_wave("vm(2)") x,y = gnucap.wave_to_arrays(w) print x,y
class MyAC2(gnucap.SIM): def do_it(self, cmd, scope): print("HELLOWORLD2") def setup(self, cmd): pass def sweep(self): pass ac = MyAC() ac2 = MyAC2() print("install1") d1 = gnucap.install_command("myac1", ac) print("install2") d2 = gnucap.install_command("myac", ac2) d3 = gnucap.install_command("myac", ac) d4 = gnucap.install_command("ac", ac) gnucap.command("set trace") gnucap.command("ac 1 2 * 2") gnucap.command("myac1 1 2 * 2") gnucap.command("myac:0 1 2 * 2") gnucap.command("myac 1 2 * 2") print("side effects?") del (d2) try: gnucap.command("myac:0 1 2 * 2") # bad command? yes. this was d2
i[idx] = ir[1] - 1 raw[idx] *= -1 elif ir[0] == -1 and ir[1] == i[idx]: raw[idx] = 0. i[idx] = 0 elif ir[0] == -1 and i[idx] >= ir[1]: i[idx] -= 1 elif ir[0] != -1 and i[idx] > ir[0]: i[idx] -= 1 assert (i[idx] >= 0) set_io(j, ir, raw) set_io(i, jr, raw) f = coo_matrix((raw, (i, j))) M = f.todense() R = np.real(M) I = np.imag(M) self._out << "zeroes\n" E = eig(R, I, homogeneous_eigvals=True) E0 = E[0] self.sort_and_output(E0) sys.stdout.flush() sim = spice_pz() d1 = install_command("pz", sim)
class mysim(SIM): def do_it(self, cmd, scope): # print("running mysim...") out("running mysim...\n") def setup(self, cmd): pass #out = IO.mstdout; #out.reset(); # BUG// don't know why this is needed def sweep(self): pass sim = mysim() d1 = install_command("test", sim) #del(d1) #what does it do? d1 = install_command("test", sim) out("installing again\n") d2 = install_command("test", sim) out("installing again\n") d3 = install_command("mysim", sim) #BUG: that does not work as expected. # it just breaks "mysim". why? #d3 = install_command("mysim", sim) sys.stdout.flush()
c = COMPONENT() print("ERROR, component must be pure") assert(False) except AttributeError: print("ignoring AttributeError") pass except RuntimeError: pass class MyAC(gnucap.SIM): def do_it(self, cmd, scope): print("HELLOWORLD") stdout.flush() ac = MyAC() c1 = gnucap.install_command("myac", ac) stdout.flush() # test: "already installed" c2 = gnucap.install_command("myac", ac) stdout.flush() print("running ac") stdout.flush() gnucap.command("ac 1 2 * 2") stdout.flush() print("running myac") stdout.flush() gnucap.command("myac 1 2 * 2")
ac = s._ac for i in range(n + 1): ac[i] = 0 s.count_iterations(iTOTAL) cl = CARD_LIST().card_list_() cl.do_ac() cl.ac_load() if (0): print("M", self.sim_()._acx[0][0], self.sim_()._acx[0][1], self.sim_()._acx[0][2]) print("M", self.sim_()._acx[1][0], self.sim_()._acx[1][1], self.sim_()._acx[1][2]) print("M", self.sim_()._acx[2][0], self.sim_()._acx[2][1], self.sim_()._acx[2][2]) acx.lu_decomp() acx.fbsub_(self.sim_()._ac) sim = mysim() d1 = install_command("mysim", sim)