Source Code for Module gmisclib.lreg_fill

 1  import g_pipe 
 2  import Num 
 3   
 4  LREG = "/home/gpk/misctts/lreg_fit2" 
 5   
 6 -def fill(f, wt, s): 
 7          assert len(f)==len(wt) 
 8          pi, po = g_pipe.popen2(LREG, ["lreg_fit2", "-smooth", "%g" % s]) 
 9   
10          for tfw in zip(f, wt): 
11                  pi.write("%g %g\n" % tfw) 
12          pi.close() 
13          o = map(float, po.readlines()) 
14          po.close() 
15          assert len(o) == len(f) 
16          tmp = Num.asarray(o, Num.Float) 
17          return tmp 
18