разница между аналогичными разреженными структурами

Question

Моего кода работает с такого родом структуры

K>> mlf=sparse([],[],[],2^31+1,1); mlf(1)=1; mlf 

mlf = 

        (1,1)      1 

но он терпит неудачу с таким родом входы ниже, где я выбираю условие в МЗФЕ, которые больше нуля (я не могу понять, как это выбор делает ввод разным)

K>> mlf=sparse([],[],[],2^31+1,1); mlf(1)=1; mlf(mlf>0) 

ans = 

    (1,1)  1 

где единственное визуальное отличие - это некоторые вкладки/пробелы.

Пожалуйста, объясните, как они отличаются.

Answer 1

Я думаю, что ответ является размер результирующего массива, как Rody предложил:

>> mlf=sparse([],[],[],2^31+1,1); mlf(1)=1; size(mlf(mlf>0)) 
ans = 
    1  1 
>> mlf=sparse([],[],[],2^31+1,1); mlf(1)=1; size(mlf) 
ans = 
       2147483649       1 

* EDIT 1: Индексация работает правильно:

>> mlf(mlf>0) = 2 
mlf = 
        (1,1)      2 

Это функционально эквивалентно использованию find

>> mlf(find(mlf)) = 2 
mlf = 
        (1,1)      2 

Похоже, что хороший вывод о том, что display форматирует выход с достаточным пространством для элемента в (2147483649,1), , но только при индексировании для назначения этому элементу (думаю, lvalue vs rvalue).

* EDIT 2: Если вы собираетесь после этих элементов в полном (не разреженный) переменной, используйте full:

>> full(mlf(mlf>0)) 
ans = 
    1 

* EDIT 3: Для того, чтобы назначить до последнего элемента в соответствии с размеры mlf, а не до последнего ненулевого элемента,

>> mlf(numel(mlf))=77 
mlf = 
        (1,1)      1 
     (2147483649,1)      77 

* EDIT 4 : Для того, чтобы удалить отрицательные значения:

mlf(mlf<0)=0; % or mlf(find(mlf<0)) = 0; 

Если вы хотите, чтобы сделать копию и удалить негативы:

mlf2 = mlf; 
mlf2(mlf2<0) = 0; 
mlf3 = mlf; 
mlf3(mlf3>0) = 0; 

Тогда у вас есть mlf со всеми значениями, mlf2 только с позитивов и mlf3 только с негативов , Ключевым моментом в этом является то, что размер остается таким же, как у оригинала mlf, поэтому вы можете использовать такие вещи, как end, оригинальным способом, основанным на размере разреженного, ура!

Answer 2

Rody Oldenhuis рекомендовал whos

>> mlf=sparse([],[],[],2^31+1,1); mlf(1)=1; mlf=mlf(mlf>0) 

mlf = 

    (1,1)  1 


>> whos mlf 
    Name  Size   Bytes Class  Attributes 

    mlf  1x1    32 double sparse  

>> mlf=sparse([],[],[],2^31+1,1); mlf(1)=1; mlf 

mlf = 

        (1,1)      1 

>> whos mlf 
    Name    Size   Bytes Class  Attributes 

    mlf  2147483649x1    32 double sparse  

который показывает ключевую проблему: размер структур изменились. chappjc предоставил возможность решить эту проблему, введя новые переменные.