def abscissae(quadn):
while len(fejer.abscissae) <= quadn:
fejer.extend()
for x in fejer.abscissae[quadn]:
yield centroid + width * x
def abscissae(quadn):
while len(fejer.abscissae) <= quadn:
fejer.extend()
for x in fejer.abscissae[quadn]:
t = (x + 1) / 2
yield domain.position - (1 - t) / t
def weights(quadn):
while len(fejer.weights) <= quadn:
fejer.extend()
for w in fejer.weights[quadn]:
yield width * w
def weights(quadn):
while len(fejer.weights) <= quadn:
fejer.extend()
for x, w in zip(fejer.abscissae[quadn], fejer.weights[quadn]):
t = (x + 1) / 2
yield 0.5 * w / t**2