Hola Alfonso, las funciones para calcular el SSI y el LI las había probado solo con números, pero fsolve parece ser que no funciona con matrices y por eso te daba el error. Con vectores 1-D en teoría sí funciona, pero el cálculo terminaba killed por falta de memoria. Lo que he hecho ha sido aplicar un ravel() a las matrices de entrada para convertirlas en vectores 1-D y resolver dentro de un for fsolve tantas veces como valores hay en los vectores 1-D. He incluido una línea de código para que se muestre por pantalla el % del cálculo completado porque tarda unos cuantos minutos. He visto también que habías hecho cambios en las funciones y creo que los he respetado en los commits.
Hola Alfonso, las funciones para calcular el SSI y el LI las había probado solo con números, pero fsolve parece ser que no funciona con matrices y por eso te daba el error. Con vectores 1-D en teoría sí funciona, pero el cálculo terminaba killed por falta de memoria. Lo que he hecho ha sido aplicar un ravel() a las matrices de entrada para convertirlas en vectores 1-D y resolver dentro de un for fsolve tantas veces como valores hay en los vectores 1-D. He incluido una línea de código para que se muestre por pantalla el % del cálculo completado porque tarda unos cuantos minutos. He visto también que habías hecho cambios en las funciones y creo que los he respetado en los commits.