Problem modeling faults - Full block matrix for faults or drift #804
-
Hi all, I am having a hard time figuring out how to deal with a problem modeling three faults. This has been troubling me for a while, but as the faults have amost the same azimuth, I was putting azimuth & dip (180,85) for the three (in my rotated model - see #776). Nonetheless, now I wanted to use their real azimuth, but then I get an error. First I thought it was something related to my rotated model but even when I do not rotate the model I still get the problem. My faults have 2 points each (2 surfaces, 2 orientations), which worked fine if they are all 180 azimuth and 85 dip. The error I get when trying to put different azimuths for each is the huge text below. Any idea what is happening here and how can I solve it? Thanks! EDIT1: I checked both #733 and #583, but in my case the error only arrives when I change those three parameters (azimuth of each fault), so I do not know how I can step-wise increase the complexity of the model. If I change one azimuth it already give me this error. scan_perform.pyx in theano.scan_module.scan_perform.perform() ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/gof/op.py in rval(p, i, o, n) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/tensor/slinalg.py in perform(self, node, inputs, output_storage) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/scipy/linalg/basic.py in solve(a, b, sym_pos, lower, overwrite_a, overwrite_b, debug, check_finite, assume_a, transposed) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/scipy/linalg/basic.py in _solve_check(n, info, lamch, rcond) LinAlgError: Matrix is singular. Backtrace when the node is created(use Theano flag traceback.limit=N to make it longer): Debugprint of the apply node: HINT: Use the Theano flag 'exception_verbosity=high' for a debugprint and storage map footprint of this apply node. During handling of the above exception, another exception occurred: LinAlgError Traceback (most recent call last) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/gempy/gempy_api.py in compute_model(model, output, at, compute_mesh, reset_weights, reset_scalar, reset_block, sort_surfaces, debug, set_solutions, **kwargs) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/compile/function_module.py in call(self, *args, **kwargs) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/gof/link.py in raise_with_op(node, thunk, exc_info, storage_map) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/six.py in reraise(tp, value, tb) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/compile/function_module.py in call(self, *args, **kwargs) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/scan_module/scan_op.py in rval(p, i, o, n, allow_gc) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/scan_module/scan_op.py in p(node, args, outs) scan_perform.pyx in theano.scan_module.scan_perform.perform() ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/gof/link.py in raise_with_op(node, thunk, exc_info, storage_map) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/six.py in reraise(tp, value, tb) scan_perform.pyx in theano.scan_module.scan_perform.perform() ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/gof/op.py in rval(p, i, o, n) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/theano/tensor/slinalg.py in perform(self, node, inputs, output_storage) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/scipy/linalg/basic.py in solve(a, b, sym_pos, lower, overwrite_a, overwrite_b, debug, check_finite, assume_a, transposed) ~/mambaforge/envs/gempy/lib/python3.7/site-packages/scipy/linalg/basic.py in _solve_check(n, info, lamch, rcond) LinAlgError: Matrix is singular. Backtrace when the node is created(use Theano flag traceback.limit=N to make it longer): Debugprint of the apply node: HINT: Use the Theano flag 'exception_verbosity=high' for a debugprint and storage map footprint of this apply node. Backtrace when the node is created(use Theano flag traceback.limit=N to make it longer): Debugprint of the apply node: Inner graphs of the scan ops: for{cpu,Looping}.0 [id A] <TensorType(float64, 4D)> ''
for{cpu,Looping}.1 [id A] <TensorType(float64, matrix)> ''
for{cpu,Looping}.2 [id A] <TensorType(float64, 3D)> ''
for{cpu,Looping}.3 [id A] <TensorType(float64, 3D)> ''
for{cpu,Looping}.4 [id A] <TensorType(bool, 3D)> ''
for{cpu,Looping}.5 [id A] <TensorType(bool, 3D)> ''
for{cpu,Looping}.6 [id A] <TensorType(float64, 4D)> ''
for{cpu,Looping}.7 [id A] <TensorType(int64, vector)> ''
Storage map footprint:
|
Beta Was this translation helpful? Give feedback.
Replies: 1 comment
-
Ok, for some reason sharing my frustration here somehow solved my problem, haha! I have no idea what I did but not the error does not appear anymore. So, marked as solved. |
Beta Was this translation helpful? Give feedback.
Ok, for some reason sharing my frustration here somehow solved my problem, haha! I have no idea what I did but not the error does not appear anymore. So, marked as solved.