jameshensman/VFF

def make_Kuf_np(X, a, b, ms):

hhjontheway opened this issue · 0 comments

Why the code [omegas = omegas[omegas != 0] # don't compute zeros freq.] lies after 'Kuf_cos'? This may lead to a different length between Kuf_cos and Kuf_sin. Is it OK?