This turns out to be a division by zero in this first part of your objective function:
sum{i in 0..&max_i., p in 0..&max_p., f in 0..&max_f.,r in 0..&max_r.}
(9*w_v[i,p,f,r]+1)*
(((e_v[i,p,f,r] - a_v[i,p,f,r])/(10 + a_v[i,p,f,r]) * 1000)**2) +
For data B, the value of a_v[8,10,3,7] is -10, yielding a zero denominator.
... View more