pyeq 0.54.9
There is a shape error when lambda_stf is not 0. and constraints are added.
This constraint does seem however to be so important.
Traceback (most recent call last):
File "/usr/local/geodesy/anaconda38/bin/pyaks.py", line 1056, in
model = pyeq.regularization.laplace.add_laplace_cons(model)
File "/usr/local/geodesy/anaconda38/lib/python3.8/site-packages/pyeq/regularization/laplace/add_laplace_cons.py", line 98, in add_laplace_cons
model = pyeq.regularization.laplace.make_normal_dlo_space_time( model, ndlos, ndlot)
File "/usr/local/geodesy/anaconda38/lib/python3.8/site-packages/pyeq/regularization/laplace/make_normal_dlo_space_time.py", line 182, in make_normal_dlo_space_time
model.N[:nfault, :3 nfault] += np.hstack([2 one, -3 * one, one])
ValueError: operands could not be broadcast together with shapes (533,1066) (533,1599) (533,1066)
pyeq 0.54.9 There is a shape error when lambda_stf is not 0. and constraints are added. This constraint does seem however to be so important.
Traceback (most recent call last): File "/usr/local/geodesy/anaconda38/bin/pyaks.py", line 1056, in
model = pyeq.regularization.laplace.add_laplace_cons(model)
File "/usr/local/geodesy/anaconda38/lib/python3.8/site-packages/pyeq/regularization/laplace/add_laplace_cons.py", line 98, in add_laplace_cons
model = pyeq.regularization.laplace.make_normal_dlo_space_time( model, ndlos, ndlot)
File "/usr/local/geodesy/anaconda38/lib/python3.8/site-packages/pyeq/regularization/laplace/make_normal_dlo_space_time.py", line 182, in make_normal_dlo_space_time
model.N[:nfault, :3 nfault] += np.hstack([2 one, -3 * one, one])
ValueError: operands could not be broadcast together with shapes (533,1066) (533,1599) (533,1066)