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