Commit 39a0925b authored by Ryan Gutenkunst's avatar Ryan Gutenkunst
Browse files

Protect from_phi_linalg from xx != yy

parent 05e317e2
......@@ -1391,6 +1391,8 @@ class Spectrum(numpy.ma.masked_array):
raw: If True, return data as a numpy array, not a Spectrum object
"""
if not (len(xx) == len(yy) and np.allclose(xx,yy)):
raise ValueError('Must have xx==yy to use linear algebra calculation of FS from phi.')
dbeta1_xx, dbeta2_xx = cached_dbeta(nx, xx)
dbeta1_yy, dbeta2_yy = cached_dbeta(ny, yy)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment