Without eval() I cannot read arr (the code throws an error) because then t inside arr[t] is still a tensor.
The idea of arr_theano = theano.shared(arr) sounds great! Now, I am reading arr_theano[t] inside my SIR-function but it will throw an error TypeError: Expected an integer.
Funny because in my DifferentialEquation-implementation I have
timeline =[0,1,2,3,4,5,6,7,8,9,10,11,12,13] (integer type). I also tried timeline which is numpy-int64 type.