Beispiel #1
0
 def test_legdiv(self) :
     for i in range(5) :
         for j in range(5) :
             msg = "At i=%d, j=%d" % (i, j)
             ci = [0]*i + [1]
             cj = [0]*j + [1]
             tgt = leg.legadd(ci, cj)
             quo, rem = leg.legdiv(tgt, ci)
             res = leg.legadd(leg.legmul(quo, ci), rem)
             assert_equal(trim(res), trim(tgt), err_msg=msg)
Beispiel #2
0
 def test_legdiv(self) :
     for i in range(5) :
         for j in range(5) :
             msg = "At i=%d, j=%d" % (i,j)
             ci = [0]*i + [1]
             cj = [0]*j + [1]
             tgt = leg.legadd(ci, cj)
             quo, rem = leg.legdiv(tgt, ci)
             res = leg.legadd(leg.legmul(quo, ci), rem)
             assert_equal(trim(res), trim(tgt), err_msg=msg)
Beispiel #3
0
def legdiv(c1, c2):
    from numpy.polynomial.legendre import legdiv
    return legdiv(c1, c2)
Beispiel #4
0
def legdiv(c1, c2):
    from numpy.polynomial.legendre import legdiv
    return legdiv(c1, c2)