From eaeeff1629aa4f4ad3de5b397887a237d19167b7 Mon Sep 17 00:00:00 2001 From: maaijke Date: Mon, 30 Oct 2023 11:29:16 +0100 Subject: [PATCH] added check on y-axis plot --- losoto/operations/structure.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/losoto/operations/structure.py b/losoto/operations/structure.py index 6f939f3c..9bd11ce8 100644 --- a/losoto/operations/structure.py +++ b/losoto/operations/structure.py @@ -112,7 +112,7 @@ def run( soltab, doUnwrap=False, refAnt='', plotName='', ndiv=1 ): #par = np.dot(np.linalg.inv(np.dot(A,A.T)),np.dot(A,np.log10(variance[myselect]))) mask = variance[myselect].mask A = np.vstack([np.log10(D2[myselect][~mask]), np.ones(len(D2[myselect][~mask]))]) - par = np.linalg.lstsq( A.T, np.log10(variance[myselect][~mask]) )[0] + par = np.linalg.lstsq( A.T, np.log10(variance[myselect][~mask]) ,rcond=None)[0] S0 = 10**(-1*par[1]/par[0]) logging.info(r't%i: beta=%.2f - R_diff=%.2f km' % (i, par[0], S0/1.e3)) variances.append(variance) @@ -149,7 +149,7 @@ def run( soltab, doUnwrap=False, refAnt='', plotName='', ndiv=1 ): ax.set_xscale('log') ax.set_yscale('log') - ymin = np.min(variance[myselect]) + ymin = np.max(1e-9,np.min(variance[myselect])) ymax = np.max(variance[myselect]) ax.set_xlim(xmin=0.1,xmax=3) ax.set_ylim(ymin,ymax)