ss_ion_bal_adf11.py eigenvector error
I was running ss_ion_bal_adf11.py and got the following error:
```
Traceback (most recent call last):
File ~/anaconda3/lib/python3.11/site-packages/spyder_kernels/py3compat.py:356 in compat_exec
exec(code, globals, locals)
File ~/Documents/au physics/colradpy/examples/ss_ion_bal_adf11.py:35
ion.solve_time_independent()#solve the ionization balance matrix in steady state
File ~/Documents/au physics/colradpy/colradpy/ionization_balance_class.py:607 in solve_time_independent
self.data['processed']['eigen_vec'] = solve_matrix_exponential_steady_state(
File ~/Documents/au physics/colradpy/colradpy/solve_matrix_exponential.py:117 in solve_matrix_exponential_steady_state
ev = eigenvectors.transpose(0,1,3,2)[index]#egienvectors sorted on eigenvals
ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (3,) + inhomogeneous part.
```
I peeked into the code and determined that index is initialized as a list of inhomogenous zero arrays. I'm not sure how this ever works, perhaps this is a deprecated feature? However, I didn't write the code or spend >5 min attempting to understand what goes on, as I suspect I would need to spend >5 weeks to debug this on my own.
I did pull up an old copy of colradpy on a different machine since I changed file reading parts on another machine. Without these mods, I get the same error. This reinforces the deprecated feature narrative.
关闭于 2024-06-03 2 条评论