Python kdot_helper Examples

Example #1

File: dot.py Project: kostyfisik/accupy

def fdot(x, y):
    '''Algorithm 5.10. Dot product algorithm in K-fold working precision,
    K >= 3.
    '''
    xx = x.reshape(-1, x.shape[-1])
    yy = y.reshape(y.shape[0], -1)
    r = _accupy.kdot_helper(xx,
                            yy).reshape((-1, ) + x.shape[:-1] + y.shape[1:])
    return fsum(r)

Example #2

def fdot(x, y):
    """Algorithm 5.10. Dot product algorithm in K-fold working precision,
    K >= 3.
    """
    xx = x.reshape(-1, x.shape[-1])
    yy = y.reshape(y.shape[0], -1)

    xx = numpy.ascontiguousarray(xx)
    yy = numpy.ascontiguousarray(yy)

    r = _accupy.kdot_helper(xx,
                            yy).reshape((-1, ) + x.shape[:-1] + y.shape[1:])
    return fsum(r)