def test_bigmat(self): sci.eval("x=rand(1000,800)") y = sci.read("x") sci.write("xx", y) sci.eval("dist = sum((x - xx).^2)") dist = sci.read("dist") assert(dist[0] == 0)
def csim(P,I,D,n0,n1,n2,n3,d0,d1,d2,d3,u): code_string1 = "s = %s;" # If the plant simulation is not required, the block functions as a controller block. if n0 == 0 and n1 == 0 and n2 == 0 and n3 == 0 and d0 == 0 and d1 == 0 and d2 == 0 and d3 == 0: code_string2 = "Gc = syslin('c'," + str(P*I) + "*s + " + str(D) + "*s^2 + 1," + str(I) + "*s);" code_string3 = "r = tf2ss(Gc);" code_string4 = "u = " + str(u) + ";" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 # If the Controller simulation is not required, the block functions as a plant block elif P == 0 and I ==0 and D == 0: code_string2 = "G = syslin('c'," + str(n0) + "*s^3 + " + str(n1) + "*s^2 + " + str(n2) + "*s + "+ str(n3)+ "," + str(d0) + "*s^3 + " + str(d1) + "*s^2 + " + str(d2) + "*s + " + str(d3) + ");" code_string3 = "r = tf2ss(G);" code_string4 = "u = " + str(u) + ";" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 # If combined plant and controller dynamics are required else: code_string2 = "Gc = syslin('c'," + str(P*I) + "*s + " + str(D) + "*s^2 + 1," + str(I) + "*s);" code_string3 = "G = syslin('c'," + str(n0) + "*s^3 + " + str(n1) + "*s^2 + " + str(n2) + "*s + "+ str(n3)+ "," + str(d0) + "*s^3 + " + str(d1) + "*s^2 + " + str(d2) + "*s + " + str(d3) + ");" code_string4 = "r = tf2ss(G*Gc);" code_string5 = "u = " + str(u)+ ";" code_string6 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 + code_string6 sciscipy.eval(code_string) y = sciscipy.read("y") return y
def test_bigmat(self): sci.eval("x=rand(1000,800)") y = sci.read("x") sci.write("xx", y) sci.eval("dist = sum((x - xx).^2)") dist = sci.read("dist") assert (dist[0] == 0)
def nplot(self): a1,a2,a3,a4=self.a.split(',') b1,b2,b3,b4=self.b.split(',') c1,c2,c3,c4=self.c.split(',') d1,d2,d3,d4=self.d.split(',') inp = "s=poly(0,'s');h=syslin('c',("+str(a1)+"*s^3+"+str(a2)+"*s^2+"+str(a3)+"*s+"+str(a4)+")/("+str(b1)+"*s^3+"+str(b2)+"*s^2+"+str(b3)+"*s+"+str(b4)+"));h1=h*syslin('c',("+str(c1)+"*s^3+"+str(c2)+"*s^2+"+str(c3)+"*s+"+str(c4)+")/("+str(d1)+"*s^3+"+str(d2)+"*s^2+"+str(d3)+"*s+"+str(d4)+"));clf(); nyquist(h1)" sciscipy.eval(inp)
def csim(P,I,D,string1,u): code_string1 = "s = %s;" code_string_u = "u = "+str(u)+";" #If the controller simulation is not required, the block functions as a plant block. if P == 0 and I == 0 and D==0: code_string2 = "G = syslin('c'," + str(string1) + ");" code_string3 = "r = tf2ss(G);" code_string4 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string_u + code_string4 # If the Plant simulation is not required, the block functions as a controller block. elif string1 == "": code_string2 = "Gc = syslin('c'," + str(P*I+D)+"*s,("+str(I)+")*s);" code_string3 = "r = tf2ss(Gc);" code_string4 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string_u + code_string4 # Combining the plant and controller dynamics else: code_string2 = "Gc=syslin('c',("+str(P*I+D)+"*s)"+","+str(I)+"*s);" code_string3 = "G = syslin('c'," + str(string1) + ");" code_string4 = "r=tf2ss(G*Gc);" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string_u + code_string5 import sciscipy sciscipy.eval(code_string) y = sciscipy.read("y") return y
def csim(P, I, D, string1, u): code_string1 = "s = %s;" code_string_u = "u = " + str(u) + ";" #If the controller simulation is not required, the block functions as a plant block. if P == 0 and I == 0 and D == 0: code_string2 = "G = syslin('c'," + str(string1) + ");" code_string3 = "r = tf2ss(G);" code_string4 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string_u + code_string4 # If the Plant simulation is not required, the block functions as a controller block. elif string1 == "": code_string2 = "Gc = syslin('c'," + str(P * I + D) + "*s,(" + str(I) + ")*s);" code_string3 = "r = tf2ss(Gc);" code_string4 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string_u + code_string4 # Combining the plant and controller dynamics else: code_string2 = "Gc=syslin('c',(" + str(P * I + D) + "*s)" + "," + str( I) + "*s);" code_string3 = "G = syslin('c'," + str(string1) + ");" code_string4 = "r=tf2ss(G*Gc);" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string_u + code_string5 import sciscipy sciscipy.eval(code_string) y = sciscipy.read("y") return y
def plot(self,b,c): string1 = "s=%s; h=syslin('c'," string2 = str(b[4])+"*s^4+"+str(b[3])+"*s^3+"+str(b[2])+"*s^2+"+str(b[1])+"*s+"+str(b[0])+"," string3 = str(c[4])+"*s^4+"+str(c[3])+"*s^3+"+str(c[2])+"*s^2+"+str(c[1])+"*s+"+str(c[0])+");" string4 = "plzr(h);" string = string1 + string2 +string3 + string4 sciscipy.eval(string)
def test_mean(self): """ [test_call] Testing mean """ mean1 = self.sci.mean([[1, 2],[3, 4]]) sciscipy.eval("mean1 = mean([1,2;3,4])") mean2 = sciscipy.read("mean1") comp = mean1 == mean2 assert(comp)
def test_strcat(self): """ [test_call] Testing strcat """ strcat1 = self.sci.strcat(["1", "4"], "x") sciscipy.eval("strcat1 = strcat(['1', '4'], 'x')") strcat2 = sciscipy.read("strcat1") comp = strcat1 == strcat2 assert(comp)
def test_length(self): """ [test_call] Testing length """ strlength1 = self.sci.length(["3ch","5char","plenty of char"]) sciscipy.eval("strlength = length(['3ch','5char','plenty of char'])") strlength2 = sciscipy.read("strlength") for l1, l2 in zip(strlength1, strlength2): self.assertEquals(l1, l2)
def run_scilab_cmd(cmd_str): """ Defines the Scilab start command (with error handle) """ new_cmd = "_ier_ = execstr('%s', 'errcatch'); _er_msg_ = lasterror() ;" % cmd_str eval(new_cmd) ier = read("_ier_") if ier != 0 and ier != [0]: lasterror = read("_er_msg_") raise ScilabError, lasterror
def test_spec(self): """ [test_call] Testing spec """ spec1 = self.sci.spec([[1, 2],[3, 4]]) sciscipy.eval("spec1 = spec([1,2;3,4])") spec2 = sciscipy.read("spec1") for l1, l2 in zip(spec1, spec2): self.assertAlmostEqual(l1, l2)
def test_readwrite(self): sci.eval("x=[1,2,3 ; 4,5,6]") y = sci.read("x") sci.write("z", y) w = sci.read("z") if numpy_is_avail: assert (numpy.alltrue(numpy.equal(y, w))) else: assert (y == w)
def polyrange(pol, x_range, xlabel, ylabel): code_string1 = "x = " + x_range + ";" code_string2 = "y = " + pol + ";" code_string3 = "plot(x,y);" code_string4 = "xtitle('Function Plot','" + xlabel + "','" + ylabel + "');" code_string5 = "xgrid();" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 sciscipy.eval(code_string)
def find_roots(self, n): string1 = "s=%s;" string2 = "h2=" + str(n[3]) + "*s^3+" + str(n[2]) + "*s^2+" + str( n[1] ) + "*s+" + str( n[0] ) + ";[E1]=roots(h2);a11=real(E1(1));a12=imag(E1(1));b11 = real(E1(2));b12 = imag(E1(2));c11 = real(E1(3));c12 = imag(E1(3)); d11 = real(E1(4)); d12 = imag(E1(4)); e11 = real(E1(5)); e12 = imag(E1(5));" string = string1 + string2 print "This is string\n", string sciscipy.eval(string) try: self.a11 = sciscipy.read("a11") print "I am a11", self.a11 except TypeError: self.a11 = 0 try: self.a12 = sciscipy.read("a12") print "I am a12", self.a12 except TypeError: self.a12 = 0 try: self.b11 = sciscipy.read("b11") print "I am b11", self.b11 except TypeError: self.b11 = 0 try: self.b12 = sciscipy.read("b12") except TypeError: self.b12 = 0 try: self.c11 = sciscipy.read("c11") print "I am c11\n", self.c11 except TypeError: self.c11 = 0 try: self.c12 = sciscipy.read("c12") print "I am c12\n", self.c12 except TypeError: self.c12 = 0 try: self.d11 = sciscipy.read("d11") print "I am d11\n", self.d11 except TypeError: self.d11 = 0 try: self.d12 = sciscipy.read("d12") except TypeError: self.d12 = 0 try: self.e11 = sciscipy.read("e11") except TypeError: self.e11 = 0 try: self.e12 = sciscipy.read("e12") except TypeError: self.e12 = 0
def test_readwrite(self): sci.eval("x=[1,2,3 ; 4,5,6]") y = sci.read("x") sci.write("z", y) w = sci.read("z") if numpy_is_avail: assert(numpy.alltrue(numpy.equal(y, w))) else: assert(y == w)
def polyrange(pol,x_range,xlabel,ylabel): code_string1 = "x = " + x_range + ";" code_string2 = "y = " + pol + ";" code_string3 = "plot(x,y);" code_string4 = "xtitle('Function Plot','" + xlabel + "','" + ylabel + "');" code_string5 = "xgrid();" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 sciscipy.eval(code_string)
def test_read_tlist(self): eval("x=tlist(['test','a','b'],12,'item')") x=read('x') if numpy_is_avail: num = numpy.array(12) else: num = 12 py_x = {'__tlist_name': 'test', 'a': num, 'b': ['item']} assert x == py_x, str(py_x) + " != tlist(['test','a','b'],12,'item')"
def test_complex(self): sci.eval("x=[1+11*%i, 2+222*%i, 3+333*%i ; 4+444*%i , 5+55*%i, 6+66*%i]") y = sci.read("x") sci.write("z", y) w = sci.read("z") if numpy_is_avail: assert(numpy.alltrue(numpy.equal(y, w))) else: assert(y == w)
def test_readwrite1dT(self): sci.eval("x=rand(100, 1)") sci.eval("my_sum=sum(x)") x = sci.read("x") my_sum = 0 for i in range(len(x)): my_sum += x[i] my_other_sum = sci.read("my_sum") assert (my_other_sum[0] == my_sum)
def test_readwrite1dT(self): sci.eval("x=rand(100, 1)") sci.eval("my_sum=sum(x)") x = sci.read("x") my_sum = 0 for i in range(len(x)): my_sum += x[i] my_other_sum = sci.read("my_sum") assert(my_other_sum[0] == my_sum)
def test_complex(self): sci.eval( "x=[1+11*%i, 2+222*%i, 3+333*%i ; 4+444*%i , 5+55*%i, 6+66*%i]") y = sci.read("x") sci.write("z", y) w = sci.read("z") if numpy_is_avail: assert (numpy.alltrue(numpy.equal(y, w))) else: assert (y == w)
def find_roots(self,n): string1 = "s=%s;" string2 = "h2="+str(n[3])+"*s^3+"+str(n[2])+"*s^2+"+str(n[1])+"*s+"+str(n[0])+";[E1]=roots(h2);a11=real(E1(1));a12=imag(E1(1));b11 = real(E1(2));b12 = imag(E1(2));c11 = real(E1(3));c12 = imag(E1(3)); d11 = real(E1(4)); d12 = imag(E1(4)); e11 = real(E1(5)); e12 = imag(E1(5));" string = string1+string2 print "This is string\n", string sciscipy.eval(string) try: self.a11 = sciscipy.read("a11") print "I am a11", self.a11 except TypeError: self.a11 = 0 try: self.a12 = sciscipy.read("a12") print "I am a12", self.a12 except TypeError: self.a12 = 0 try: self.b11 = sciscipy.read("b11") print "I am b11", self.b11 except TypeError: self.b11 = 0 try: self.b12 = sciscipy.read("b12") except TypeError: self.b12 = 0 try: self.c11 = sciscipy.read("c11") print "I am c11\n", self.c11 except TypeError: self.c11 = 0 try: self.c12 = sciscipy.read("c12") print "I am c12\n", self.c12 except TypeError: self.c12 = 0 try: self.d11 = sciscipy.read("d11") print "I am d11\n", self.d11 except TypeError: self.d11 = 0 try: self.d12 = sciscipy.read("d12") except TypeError: self.d12 = 0 try: self.e11 = sciscipy.read("e11") except TypeError: self.e11 = 0 try: self.e12 = sciscipy.read("e12") except TypeError: self.e12 = 0
def file_import(self): import sciscipy f = open(self.path) x = f.read() x = x.split("\n") code_string = "" for i in range(0, len(x)): code_string += x[i] sciscipy.eval(code_string) self.ret_array = sciscipy.read(self.var_name)
def file_import(self): import sciscipy f = open(self.path) x = f.read() x = x.split("\n") code_string = "" for i in range(0,len(x)): code_string += x[i] sciscipy.eval(code_string) self.ret_array = sciscipy.read(self.var_name)
def nplot(self): a1, a2, a3, a4 = self.a.split(',') b1, b2, b3, b4 = self.b.split(',') c1, c2, c3, c4 = self.c.split(',') d1, d2, d3, d4 = self.d.split(',') inp = "s=poly(0,'s');h=syslin('c',(" + str(a1) + "*s^3+" + str( a2) + "*s^2+" + str(a3) + "*s+" + str(a4) + ")/(" + str( b1) + "*s^3+" + str(b2) + "*s^2+" + str(b3) + "*s+" + str( b4) + "));h1=h*syslin('c',(" + str(c1) + "*s^3+" + str( c2) + "*s^2+" + str(c3) + "*s+" + str( c4) + ")/(" + str(d1) + "*s^3+" + str( d2) + "*s^2+" + str(d3) + "*s+" + str( d4) + "));clf(); nyquist(h1)" sciscipy.eval(inp)
def csim(P,I,D,n0,n1,d0,d1,u): code_string1 = "s=%s;" code_string2 = "Gc=syslin('c',("+str(P*I+D)+"*s)"+","+str(I)+"*s);" code_string3 = "G=syslin(" code_string4 = "'c'"+","+str(n0)+"*s"+"+"+str(n1)+","+str(d0)+"*s"+"+"+str(d1)+");" code_string5 = "r=tf2ss(G*Gc);" code_string6 = "u="+str((u))+";" code_string7 = "y=csim(u,1:length(u),r)" code_string = code_string1+code_string2+code_string3+code_string4+code_string5+code_string6+code_string7 import sciscipy sciscipy.eval(code_string) y = sciscipy.read("y") return y
def find_scilab_type(var_name): """ Find the scilab type of var_name @param var_name: name of a scilab variable @type var_name: string @return: type(var_name) """ if type(var_name) != type(""): raise TypeError, "var_name must be a string" run_scilab_cmd("_tmp1_ = type(" + var_name + ")") res = read("_tmp1_") eval("clear _tmp1_") return res[0]
def find_resp(self,b,c): if (self.itype == 11): typo = "'imp'" elif (self.itype == 12): typo = "'step'" else: typo = "t" string1 = "s=%s; h=syslin('c'," string2 = str(b[4])+"*s^4+"+str(b[3])+"*s^3+"+str(b[2])+"*s^2+"+str(b[1])+"*s+"+str(b[0])+"," string3 = str(c[4])+"*s^4+"+str(c[3])+"*s^3+"+str(c[2])+"*s^2+"+str(c[1])+"*s+"+str(c[0])+");" string4 = "t="+str(self.tstart)+":"+str(self.tstep)+":"+str(self.tstop)+";" string5 = "deff('u=input(t)','u=50');" string6 = "r=tf2ss(h); a=csim("+typo+",t,r);" string = string1+string2+string3+string4+string5+string6 sciscipy.eval(string) self.a = sciscipy.read("a")
def csim(P, I, D, n0, n1, d0, d1, u): code_string1 = "s=%s;" code_string2 = "Gc=syslin('c',(" + str(P * I + D) + "*s)" + "," + str(I) + "*s);" code_string3 = "G=syslin(" code_string4 = "'c'" + "," + str(n0) + "*s" + "+" + str(n1) + "," + str( d0) + "*s" + "+" + str(d1) + ");" code_string5 = "r=tf2ss(G*Gc);" code_string6 = "u=" + str((u)) + ";" code_string7 = "y=csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 + code_string6 + code_string7 import sciscipy sciscipy.eval(code_string) y = sciscipy.read("y") return y
def discrete_sim(P,I,D,n0,n1,st,d0,d1,u): code_string1 = "s=%s;" code_string2 = "Gc=syslin("+str(st)+",("+str(P*I+D)+"*s)"+","+str(I)+"*s);" code_string3 = "G=syslin(" code_string4 = str(st)+","+str(n0)+"*s"+"+"+str(n1)+","+str(d0)+"*s"+"+"+str(d1)+");" code_string5 = "r=tf2ss(G*Gc);" code_string6 = "u="+str((u))+";" code_string7 = "y=dsimul(r,u)" code_string = code_string1+code_string2+code_string3+code_string4+code_string5+code_string6+code_string7 # Check complete_code_string #print code_string import sciscipy sciscipy.eval(code_string) y = sciscipy.read("y") return y
def discrete_sim(P, I, D, n0, n1, st, d0, d1, u): code_string1 = "s=%s;" code_string2 = "Gc=syslin(" + str(st) + ",(" + str( P * I + D) + "*s)" + "," + str(I) + "*s);" code_string3 = "G=syslin(" code_string4 = str(st) + "," + str(n0) + "*s" + "+" + str(n1) + "," + str( d0) + "*s" + "+" + str(d1) + ");" code_string5 = "r=tf2ss(G*Gc);" code_string6 = "u=" + str((u)) + ";" code_string7 = "y=dsimul(r,u)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 + code_string6 + code_string7 # Check complete_code_string #print code_string import sciscipy sciscipy.eval(code_string) y = sciscipy.read("y") return y
def csim(n0,n1,st,d0,d1,u): string1 = "s=%s;" string2 = "h=syslin(" string3 = str(st)+","+str(n0)+"*s"+"+"+str(n1)+","+str(d0)+"*s"+"+"+str(d1)+");" string4 = "r=tf2ss(h);" string5 = "u="+str(u)+";" string5 = u+";" string6 = "y=csim(u,1:length(u),r)" string = string1+string2+string3+string4+string5+string6 # Check complete_string #print string import sciscipy sciscipy.eval(string) y = sciscipy.read("y") #plt.plot(y) #plt.show() return y
def find_resp(self,b,c): if (self.itype == 11): typo = "'imp'" elif (self.itype == 12): typo = "'step'" else: typo = "t" # print "self b",self.b string1 = "s=%s; h=syslin('c'," string2 = str(b[4])+"*s^4+"+str(b[3])+"*s^3+"+str(b[2])+"*s^2+"+str(b[1])+"*s+"+str(b[0])+"," string3 = str(c[4])+"*s^4+"+str(c[3])+"*s^3+"+str(c[2])+"*s^2+"+str(c[1])+"*s+"+str(c[0])+");" string4 = "t="+str(self.tstart)+":"+str(self.tstep)+":"+str(self.tstop)+";" string5 = "deff('u=input(t)','u=50');" string6 = "r=tf2ss(h); a=csim("+typo+",t,r);" string = string1+string2+string3+string4+string5+string6 print "I am strin g",string sciscipy.eval(string) self.a = sciscipy.read("a") print "value of a\n",self.a
def calc_op(n0,n1,n2, d0,d1,d2): string1 = "s=%s; h=syslin('c'," string2 = str(n0)+"*s^2+"+str(n1)+"*s+"+str(n2)+"," string3 = str(d0)+"*s^2+"+str(d1)+"*s+"+str(d2)+");" string4 = "t=0:0.01:10;" string5 = "deff('u=input(t)','u=50');" string6 = "r=tf2ss(h); a=csim(input,t,r);" full_string = string1+string2+string3+string4+string5+string6 print "ex:" print full_string import sciscipy sciscipy.eval(full_string) b = sciscipy.read("a") t = sciscipy.read("t") print "output is ",b return b
def csim(n0, n1, st, d0, d1, u): string1 = "s=%s;" string2 = "h=syslin(" string3 = str(st) + "," + str(n0) + "*s" + "+" + str(n1) + "," + str( d0) + "*s" + "+" + str(d1) + ");" string4 = "r=tf2ss(h);" string5 = "u=" + str(u) + ";" string5 = u + ";" string6 = "y=csim(u,1:length(u),r)" string = string1 + string2 + string3 + string4 + string5 + string6 # Check complete_string #print string import sciscipy sciscipy.eval(string) y = sciscipy.read("y") #plt.plot(y) #plt.show() return y
def calc_op(n0, n1, n2, d0, d1, d2): string1 = "s=%s; h=syslin('c'," string2 = str(n0) + "*s^2+" + str(n1) + "*s+" + str(n2) + "," string3 = str(d0) + "*s^2+" + str(d1) + "*s+" + str(d2) + ");" string4 = "t=0:0.01:10;" string5 = "deff('u=input(t)','u=50');" string6 = "r=tf2ss(h); a=csim(input,t,r);" full_string = string1 + string2 + string3 + string4 + string5 + string6 print "ex:" print full_string import sciscipy sciscipy.eval(full_string) b = sciscipy.read("a") t = sciscipy.read("t") print "output is ", b return b
def csim(P, I, D, n0, n1, n2, n3, d0, d1, d2, d3, u): code_string1 = "s = %s;" # If the plant simulation is not required, the block functions as a controller block. if n0 == 0 and n1 == 0 and n2 == 0 and n3 == 0 and d0 == 0 and d1 == 0 and d2 == 0 and d3 == 0: code_string2 = "Gc = syslin('c'," + str( P * I) + "*s + " + str(D) + "*s^2 + 1," + str(I) + "*s);" code_string3 = "r = tf2ss(Gc);" code_string4 = "u = " + str(u) + ";" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 # If the Controller simulation is not required, the block functions as a plant block elif P == 0 and I == 0 and D == 0: code_string2 = "G = syslin('c'," + str(n0) + "*s^3 + " + str( n1) + "*s^2 + " + str(n2) + "*s + " + str(n3) + "," + str( d0) + "*s^3 + " + str(d1) + "*s^2 + " + str( d2) + "*s + " + str(d3) + ");" code_string3 = "r = tf2ss(G);" code_string4 = "u = " + str(u) + ";" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 # If combined plant and controller dynamics are required else: code_string2 = "Gc = syslin('c'," + str( P * I) + "*s + " + str(D) + "*s^2 + 1," + str(I) + "*s);" code_string3 = "G = syslin('c'," + str(n0) + "*s^3 + " + str( n1) + "*s^2 + " + str(n2) + "*s + " + str(n3) + "," + str( d0) + "*s^3 + " + str(d1) + "*s^2 + " + str( d2) + "*s + " + str(d3) + ");" code_string4 = "r = tf2ss(G*Gc);" code_string5 = "u = " + str(u) + ";" code_string6 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 + code_string6 sciscipy.eval(code_string) y = sciscipy.read("y") return y
def find_output_param(macro_name): """ Find out the number of output param of macro_name First we look in the __known_func dico to see if we have a special case for that macro. If not, we use macrovar for type 13 functions. Otherwise, we return 1. @param macro_name: the name of a scilab macro @type macro_name: string @return: number of ouput param of macro_name @rtype: integer """ if type(macro_name) != type(""): raise TypeError, "macro_name must be a string" if macro_name in __known_func.keys(): return __known_func[macro_name] if find_scilab_type(macro_name) == 13: eval("_tmp1_ = macrovar(" + macro_name + ");") eval("_tmp2_ = length(length(_tmp1_(2)))") res = read("_tmp2_") eval("clear _tmp1_, _tmp2_") return int(res[0]) return 1
import sciscipy import scilab code_string1 = "s = %s;" code_string2 = "Gc = syslin('c'," + str( 2 * 1) + "*s + " + str(2) + "*s^2 + 1," + str(3) + "*s);" code_string3 = "r = tf2ss(Gc);" code_string4 = "u = " + str(3) + ";" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 sciscipy.eval(code_string) y = sciscipy.read("y")
def nplot(self): inp = "s=poly(0,'s'); n="+str(self.a)+";"+"d="+str(self.b)+";h=syslin('c',n./d); plzr(h);" sciscipy.eval(inp)
def scilab_evaluate(request): try: user_id = request.session['user_id'] except: return HttpResponseRedirect("/login") all_code = request.POST.get('scilab_code') if not all_code: return HttpResponseRedirect("/scilab_view") all_code = all_code.replace(" ","") graphics_mode = request.POST.get('graphicsmode') all_code = all_code.replace("//Type Code Here","") filter_for_system = re.compile("unix_g|unix_x|unix_w|unix_s") if not (filter_for_system.findall(all_code)): split_code = all_code.split() plot_filter = re.compile("plot2d\(.*\)") function_filter = re.compile("deff\(.*\)") variable =[] expression = [] the_data_set = {} output = "" the_variable = "" graphs = [] links = [] for i in range(0,len(split_code)): plot_data = plot_filter.findall(split_code[i]) function_data = function_filter.findall(split_code[i]) if function_data: #a=sciscipy.eval(split_code[i]) #re_fnname = re.compile("\[.*\]=.*\(.*,.\)") #fn_data = re_fnname.findall(split_code[i]) return_variables = scilab_instances(split_code,all_code) #eval("function_name = fn_data[0].split('=')[1]") return render_to_response('default.html',{'input':all_code, 'output':return_variables , 'username':request.session['username']}) #return_variables = eval('sci.function_name') #print return_variables #return HttpResponse(return_variables) split_more = split_code[i].split("=") if (len(split_more)>1): expression.append(split_more[1]) evaluated_code = sciscipy.eval(split_code[i]) the_variable = split_code[i].split("=")[0] the_variable_value = sciscipy.read(the_variable) the_data_set[the_variable] = the_variable_value output = output +str(the_variable) + " =\n " + " " + str(the_variable_value) + "\n" if plot_data: coordinates = plot_data[0].split("(")[1].split(")")[0].split(",") x=coordinates[0] y=coordinates[1] pylab.plot(the_data_set[x],the_data_set[y]) graphs.append('simple_plot' + str(i)) cwd = str(os.getcwd()) + "/graphs/" + str(request.session['user_id']) cwdsf = cwd + str(graphs[-1]) if not os.path.exists(cwd): os.makedirs(cwd) pylab.savefig(cwdsf) for graph in graphs: p = canvas.Canvas(cwd+str(graph)+".pdf",pagesize=letter) links.append(str(request.META['SERVER_NAME'])+":"+str(request.META['SERVER_PORT']) + str(graph)) p.drawImage(cwd+str(graph)+".png", 1*inch,1*inch, width=5*inch,height=5*inch,mask=None) p.showPage() p.save() return render_to_response('default.html',{'input':all_code, 'output':output , "graphs":graphs,'username':request.session['username'],"links":links}) else: return render_to_response('default.html',{'input':all_code,'output':"error"})
def nplot(self): inp = "s=poly(0,'s');h=syslin('c',("+self.a+")/("+self.b+"));clf();bode(h,0.1,100)" sciscipy.eval(inp)
def nplot(self): inp = "s=poly(0,'s'); n=" + str(self.a) + ";" + "d=" + str( self.b) + ";h=syslin('c',n./d); plzr(h);" sciscipy.eval(inp)
def scipoll(): HOW_LONG = 0.1 # sec while 1: eval("") time.sleep(HOW_LONG)
def work(self, input_items, output_items): import sciscipy out_eval_string = "exec('" + self.path + "', -1)" sciscipy.eval(out_eval_string)
import sciscipy import scilab code_string1 = "s = %s;" code_string2 = "Gc = syslin('c'," + str(2*1) + "*s + " + str(2) + "*s^2 + 1," + str(3) + "*s);" code_string3 = "r = tf2ss(Gc);" code_string4 = "u = " + str(3) + ";" code_string5 = "y = csim(u,1:length(u),r)" code_string = code_string1 + code_string2 + code_string3 + code_string4 + code_string5 sciscipy.eval(code_string) y = sciscipy.read("y")