Come riconoscere i bug di overflow in Matlab?

Question

Ho trascorso parte di ieri e oggi rintracciando un bug in qualche codice Matlab. Avevo pensato che il mio problema fosse l'indicizzazione (con molte strutture che non ho definito e mi sto ancora abituando), ma si è rivelato un bug di overflow. Ho perso questo per un motivo ben preciso:

>> uint8(2) - uint8(1) 

ans = 

    1 

>> uint8(2) - uint8(2) 

ans = 

    0 

>> uint8(2) - uint8(3) 

ans = 

    0 

mi sarei aspettato l'ultimo ad essere qualcosa di simile a -1 (o 255). Nel mezzo di un grande vettore, gli errati 0 s erano difficili da rilevare, ma uno 255 si sarebbe distinto facilmente.

Eventuali suggerimenti su come rilevare facilmente questi problemi in futuro? (Idealmente, mi piacerebbe disattivare il controllo di overflow per farlo funzionare come C.) Cambiare in double funziona, ovviamente, ma se non mi rendo conto che è un uint8 per cominciare, che non aiuta.

Answer 1

Si può iniziare accendendo avvertimenti interi:

intwarning('on') 

Questo vi darà un avviso quando integer overflow aritmetici.

Attenzione però, come indicato in here, questo rallenta l'aritmetica dei numeri interi, quindi usalo solo durante il debug.

Answer 2

Vedere Funzione INTWARNING per controllare gli avvisi sulle operazioni con numero intero.

http://www.mathworks.com/access/helpdesk/help/techdoc/ref/intwarning.html

Answer 3

A partire dalla versione R2010b e più tardi, la funzione INTWARNING has been removed, insieme a these warning messages for integer math and conversion:

• MATLAB:intConvertNaN
• MATLAB:intConvertNonIntVal
• MATLAB:intConvertOverflow
• MATLAB:intMathOverflow

Quindi utilizzare INTWARNING non è più un'opzione valida per determinare quando si verificano overflow integer. Un'alternativa è utilizzare la funzione CLASS per verificare la classe dei dati e riformularli di conseguenza prima di eseguire l'operazione. Ecco un esempio:

if strcmp(class(data),'uint8') %# Check if data is a uint8 
    data = double(data);   %# Convert data to a double 
end 

È anche possibile utilizzare la funzione ISA così:

if ~isa(data,'single') %# Check if data is not a single 
    data = single(data); %# Convert data to a single 
end