Update __LinearSolver__.py

Get pore data for rate() calculation instead of the result of last algorithm run.
Use rate in _calc_eff_prop() which works for non-cubic networks.