plotiterator.py 文件源码

def solint_numpy_indexing(dsref):
    start = time.time()

    dsref.as_numarray()
    tms = numpy.unique(dsref.x)

    # check if there is something to be averaged at all
    if len(tms)==len(dsref.x):
        return time.time() - start

    newds = dataset()
    for tm in tms:
        newds.append(tm, numpy.average(dsref.y[numpy.where(dsref.x==tm)]) )
    dsref.x = newds.x
    dsref.y = newds.y
    return time.time() - start