Mercurial > forge

comparison extra/NaN/inst/covm.m @ 12640:de98e4cb9248 octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
check for sparse matrices and and convert to full if needed
author schloegl
date Thu, 18 Jun 2015 15:09:49 +0000
parents e0a8637557ab
children
comparison
equal deleted inserted replaced
12639:061284acaabd 12640:de98e4cb9248
141 Y = Y - ones(r1,1)*(S2./N2); 141 Y = Y - ones(r1,1)*(S2./N2);
142 end; 142 end;
143 end; 143 end;
144 end; 144 end;
145 145
146 if issparse(X) || issparse(Y),
147 fprintf(2,'sumskipnan: sparse matrix converted to full matrix\n');
148 X=full(X);
149 Y=full(Y);
150 end;
151
146 [CC,NN] = covm_mex(real(X), real(Y), FLAG_NANS_OCCURED, W); 152 [CC,NN] = covm_mex(real(X), real(Y), FLAG_NANS_OCCURED, W);
147 %% complex matrices 153 %% complex matrices
148 if ~isreal(X) && ~isreal(Y) 154 if ~isreal(X) && ~isreal(Y)
149 [iCC,inn] = covm_mex(imag(X), imag(Y), FLAG_NANS_OCCURED, W); 155 [iCC,inn] = covm_mex(imag(X), imag(Y), FLAG_NANS_OCCURED, W);
150 CC = CC + iCC; 156 CC = CC + iCC;