fiber_utils.py 文件源码

def get_norm_for_gauss_power(siglow=1.,sighigh=4.,powlow=1.5,powhigh=4.):
    S,P = np.meshgrid(np.linspace(siglow,sighigh,31), 
                      np.linspace(powlow,powhigh,26))
    x = np.linspace(-20,20,41)

    S = S.ravel()
    P = P.ravel()
    v = np.zeros(S.shape)
    for i,s in enumerate(S):
        v[i] = gauss_power(x, s, P[i]).sum()
    X = np.zeros((len(S),2))
    X[:,0] = S
    X[:,1] = P
    return LinearNDInterpolator(X,v, fill_value=1.)